1 function W = fv_mass_matrix(U,grid,model)
2 %
function M = fv_mass_matrix([],grid,[])
4 %
function computing the l2 inner product matrix
for fv-functions
6 % By
this function, inner products between functions can be
7 % computed by their DOF vectors: K = U1
' * W * U2;
9 % most arguments are empty as they are only required for ldg functions
11 % Bernard Haasdonk 20.7.2006
13 % currently simple diagonal weighting of element values sufficient
14 W = sparse(1:grid.nelements,1:grid.nelements,grid.A);