From a06e8608043ff3194fd07b11f3bc71c06d3c8aa6 Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Fri, 27 Jan 2017 10:32:21 +0100
Subject: [PATCH] [PL] HT: Use Eigen::norm().
---
ProcessLib/HT/HTFEM.h | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h
index 91d919096c7..7cec15a3911 100644
--- a/ProcessLib/HT/HTFEM.h
+++ b/ProcessLib/HT/HTFEM.h
@@ -183,8 +183,7 @@ public:
-perm_over_visc *
(sm.dNdx * p_nodal_values - density_water_T * b);
- double const velocity_magnitude =
- std::sqrt(velocity.transpose() * velocity);
+ double const velocity_magnitude = velocity.norm();
GlobalDimMatrixType const& I(
GlobalDimMatrixType::Identity(GlobalDim, GlobalDim));
GlobalDimMatrixType thermal_dispersivity =
--
GitLab