Skip to content
Snippets Groups Projects
Commit 7cd702a8 authored by Tom Fischer's avatar Tom Fischer Committed by Dmitri Naumov
Browse files

[PL/HT] Constify K_over_mu in staggered HT.

parent 7ed38ee9
No related branches found
No related tags found
No related merge requests found
......@@ -132,7 +132,7 @@ void StaggeredHTFEM<ShapeFunction, IntegrationMethod, GlobalDim>::
solid_phase
.property(MaterialPropertyLib::PropertyType::permeability)
.value(vars));
GlobalDimMatrixType K_over_mu = intrinsic_permeability / viscosity;
GlobalDimMatrixType const K_over_mu = intrinsic_permeability / viscosity;
// matrix assembly
local_M.noalias() +=
......@@ -269,7 +269,8 @@ void StaggeredHTFEM<ShapeFunction, IntegrationMethod, GlobalDim>::
.property(MaterialPropertyLib::PropertyType::permeability)
.value(vars));
GlobalDimMatrixType K_over_mu = intrinsic_permeability / viscosity;
GlobalDimMatrixType const K_over_mu =
intrinsic_permeability / viscosity;
GlobalDimVectorType const velocity =
material_properties.has_gravity
? GlobalDimVectorType(-K_over_mu * (dNdx * local_p_Eigen_type -
......
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