250
250
subroutine compute_mass(positions, mesh, mass, lumped_mass, density)
251
251
type(vector_field), intent(in) :: positions
253
253
type(csr_matrix), intent(inout) :: mass
254
254
type(scalar_field), intent(inout), optional :: lumped_mass
255
255
type(scalar_field), intent(inout), target, optional :: density