Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

RichardsMechanics: mass lumping #2590

merged 3 commits into from Aug 6, 2019
Changes from 1 commit
File filter...
Filter file types
Jump to…
Jump to file or symbol
Failed to load files and symbols.


Just for now


[PL] RM; Implement mass lumping.

Both, Picard and Newton implementations are done.
  • Loading branch information...
endJunction committed Aug 5, 2019
commit dc69d0ed4b40b43a2e2b144b711caf2c33cbf2c2
@@ -281,6 +281,13 @@ void RichardsMechanicsLocalAssembler<
.noalias() += N_p.transpose() * S_L * rho_LR * alpha *
identity2.transpose() * B * w;

if (_process_data.apply_mass_lumping)
auto Mpp = M.template block<pressure_size, pressure_size>(
pressure_index, pressure_index);
Mpp = Mpp.colwise().sum().eval().asDiagonal();

template <typename ShapeFunctionDisplacement, typename ShapeFunctionPressure,
@@ -553,6 +560,11 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
dNdx_p.transpose() * rho_LR * k_rel * rho_Ki_over_mu * b * w;

if (_process_data.apply_mass_lumping)
storage_p = storage_p.colwise().sum().eval().asDiagonal();

This comment has been minimized.

Copy link

wenqing Aug 6, 2019


In future, the name of this term is better to be unified both here and where in assemble().


// pressure equation, pressure part.
.template block<pressure_size, pressure_size>(pressure_index,
ProTip! Use n and p to navigate between commits in a pull request.
You can’t perform that action at this time.