1 function M = fv_local_mass_matrix_rect(qdeg,model)
2 %
function M = fv_local_mass_matrix_rect(qdeg,model)
4 %
function computing the local mass matrix of a scalar fv basis by
5 % numerical integration of degree qdeg.
7 % Result is a numlocal_base_function x num_local_base_functions
9 % @verbatim@f$\int_{\mbox{ref}_{\mbox{tria}}} \hat phi_i \hat phi_j@f$ \endverbatim.
11 % as only first order fv elements are supported currently,
12 % the
"matrix" is a simply 1/2
for triangles and 1
for rectangles
14 % Bernard Haasdonk 31.8.2009