Skip to content
Snippets Groups Projects
Commit a06e8608 authored by Tom Fischer's avatar Tom Fischer
Browse files

[PL] HT: Use Eigen::norm().

parent c8c88eb1
No related branches found
No related tags found
No related merge requests found
...@@ -183,8 +183,7 @@ public: ...@@ -183,8 +183,7 @@ public:
-perm_over_visc * -perm_over_visc *
(sm.dNdx * p_nodal_values - density_water_T * b); (sm.dNdx * p_nodal_values - density_water_T * b);
double const velocity_magnitude = double const velocity_magnitude = velocity.norm();
std::sqrt(velocity.transpose() * velocity);
GlobalDimMatrixType const& I( GlobalDimMatrixType const& I(
GlobalDimMatrixType::Identity(GlobalDim, GlobalDim)); GlobalDimMatrixType::Identity(GlobalDim, GlobalDim));
GlobalDimMatrixType thermal_dispersivity = GlobalDimMatrixType thermal_dispersivity =
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment