Re: [Getfem-users] Lumped mass matrix

[Roman Putanowicz]

Dear All,

Continuing my previous post:
> I see two other approaches :
>   a)  summing the rows to the assembled global matrix.
In the attachment you have the second version of the
assembly routine  asm_lumped_mass_matrix_v1() where
the lumping is done after assembly.
The only thing I wonder is, if it makes sense to check for
non-zero entries of the matrix and how to set sensibly the
threshold (which here is set arbitrarily to 1e-10).

Regards,

Roman

lumped_matrix.cxx
Description: Text Data