[Libmesh-users] Mass lumping and trapezoidal quadrature rule

[marco <in...@gm...>]

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