lumped_mass(3rheolef) rheolef-6.1 lumped_mass(3rheolef)
NAME
lumped_mass -- lumped L2 scalar product
SYNOPSIS
form(const space& V, const space& V, "lumped_mass");
form(const space& M, const space& V, "lumped_mass");
form (const space& V, const space& V, "lumped_mass", const domain& gamma);
form_diag(const space& V, "mass");
EXAMPLE
The use of lumped mass form write:
form m(V, "lumped_mass");
or (see also mass(3)):
form_diag md(V, "mass");
The lumped procedure sums all extra-diagonal terms on the diagonal: let us denote M the original mass matrix, then, the lumped mass matrix
is the diagonal matrix defined by:
n
----
M1(i,i) = / M(i,j)
----
j=1
SEE ALSO
mass(3)
rheolef-6.1 rheolef-6.1 lumped_mass(3rheolef)
