Mass lumping is defined by using trapezoidal quadrature rule. In this way, the weights of a trapezoidal rule should be the integral over the element of each basis function.
I noticed that this is not always true for hexahedral elements. Am I missing something or could there be a bug? Thank you _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
