Skip to content
Snippets Groups Projects
Commit 426472f3 authored by Christoph Lehmann's avatar Christoph Lehmann Committed by Dmitry Yu. Naumov
Browse files

[PL/RM] Reduced variable scope (sigma_eff, rho)

parent 0ca4756b
No related branches found
No related tags found
No related merge requests found
......@@ -1045,12 +1045,16 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
.noalias() += B.transpose() * C * B * w;
}
auto const& sigma_eff = _ip_data[ip].sigma_eff;
double const rho = *std::get<Density>(CD);
auto const& b = _process_data.specific_body_force;
local_rhs.template segment<displacement_size>(displacement_index)
.noalias() -=
(B.transpose() * sigma_eff - N_u_op(N_u).transpose() * rho * b) * w;
{
auto const& sigma_eff = _ip_data[ip].sigma_eff;
double const rho = *std::get<Density>(CD);
local_rhs.template segment<displacement_size>(displacement_index)
.noalias() -= (B.transpose() * sigma_eff -
N_u_op(N_u).transpose() * rho * b) *
w;
}
//
// displacement equation, pressure part
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment