Implement positivity preserving with lumped mass. See #353
Implement positivity preserving with lumped mass. See #353