From e9bf477d5f4047299b15c62a90e9fd6ed07632cb Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Tue, 14 Jan 2020 13:45:22 +0100
Subject: [PATCH] [PL/LF] getFlux(): Init flux vec with zeros.
In 2d case the z component would be arbitrary when uninitialized.
---
ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
index 8e09f59a246..63a60f77371 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
@@ -87,7 +87,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::getFlux(
const double mu =
_material_properties.getViscosity(pressure, _reference_temperature);
- Eigen::Vector3d flux;
+ Eigen::Vector3d flux(0.0, 0.0, 0.0);
flux.head<GlobalDim>() =
-permeability / mu * shape_matrices.dNdx *
Eigen::Map<const NodalVectorType>(local_x.data(), local_x.size());
--
GitLab