From dbca423978e4411ef1a47979db1b73668d8b41ea Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Wed, 6 Dec 2017 14:34:15 +0100
Subject: [PATCH] [PL] LIE/HM; Don't recompute perm. for output.
---
.../LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h b/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h
index 413970d5952..6a009bdb979 100644
--- a/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h
+++ b/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h
@@ -384,8 +384,7 @@ void HydroMechanicsLocalAssemblerFracture<ShapeFunctionDisplacement,
effective_stress_prev, effective_stress, C, state);
// permeability
- double const local_k = b_m * b_m / 12;
- ip_data.permeability = local_k;
+ double const local_k = ip_data.permeability;
local_k_tensor.diagonal().head(GlobalDim - 1).setConstant(local_k);
GlobalDimMatrix const k = R.transpose() * local_k_tensor * R;
--
GitLab