[Nektar-users] EulerCFD in 1D with variable coefficients in mass and stiffness matrix

Dear Spencer, How can I implement or modify the EulerCFE code (only 1D is needed) to obtain a variant with mass and stiffness matrices with variable coefficients (i.e. M_A = \int_{-1}^{1} L_q L_p A(\xi,t) d\xi), as well as setting up similar matrices with variable coefficients for setting up source terms (e.g. M_s S_p). Best regards, Helmut