%
% HW2
%

clear all;

gllpts = 8;

% Initialize operators
[Operators.D] = InitDerivativeGaussLobattoSEM(gllpts);

% 2D operators
[BL,ML]       = InitGaussLobattoLocalMassMatrix(gllpts);


% Weak derivative
Operators.Dwk = (Operators.D)'*diag(ML);
Operators.ML    = ML;

% Stiffness matrix
Operators.K = Operators.Dwk*diag(ML)*Operators.D;

% Construct M and K as in homework 1
% and solve the same problem.


