diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h
index 91d919096c7e262a14fdad0f0fa4297860acf657..7cec15a3911e5e74310ee16d36128ff3e5bcb803 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 =