From 10e698bf9d3e5b06cc24f1ebba07d298301c206b Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Fri, 19 Feb 2021 16:11:06 +0100
Subject: [PATCH] [PL/RM] Comment the whole block for empty if-cnd.
---
.../RichardsMechanicsFEM-impl.h | 22 +++++++++----------
1 file changed, 11 insertions(+), 11 deletions(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 98f57b37cdc..ca36115eb74 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -995,17 +995,17 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
.noalias() +=
N_u_op.transpose() * phi * rho_LR * dS_L_dp_cap * b * N_p * w;
- if (medium->hasProperty(MPL::PropertyType::saturation_micro))
- {
- // For the swelling stress with double structure model the
- // corresponding Jacobian u-p entry would be:
- // -B.transpose() *
- // dsigma_sw_dS_L_m* dS_L_m_dp_cap_m*(p_L_m - p_L_m_prev) /
- // p_cap_dot_ip / dt* N_p* w;
- // but it does not improve convergence and sometimes worsen it.
- }
- else if (solid_phase.hasProperty(
- MPL::PropertyType::swelling_stress_rate))
+ // For the swelling stress with double structure model the corresponding
+ // Jacobian u-p entry would be required, but it does not improve
+ // convergence and sometimes worsens it:
+ // if (medium->hasProperty(MPL::PropertyType::saturation_micro))
+ // {
+ // -B.transpose() *
+ // dsigma_sw_dS_L_m* dS_L_m_dp_cap_m*(p_L_m - p_L_m_prev) /
+ // p_cap_dot_ip / dt* N_p* w;
+ // }
+ if (!medium->hasProperty(MPL::PropertyType::saturation_micro) &&
+ solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate))
{
using DimMatrix = Eigen::Matrix<double, 3, 3>;
auto const dsigma_sw_dS_L =
--
GitLab