From d1a86a654ab4d7c15e0f16c05b79f932ddcb0fdc Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Mon, 1 Mar 2021 20:58:10 +0000
Subject: [PATCH] [md][1d3d] Avoid square root call in kernel evaluation
---
dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh b/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh
index bbc33fde6d..2fc693cf36 100644
--- a/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh
+++ b/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh
@@ -540,9 +540,9 @@ private:
// compute distance
auto proj = a; proj.axpy(t, ab);
- const auto r = (proj - point).two_norm();
+ const auto radiusSquared = (proj - point).two_norm2();
- if (r > rho)
+ if (radiusSquared > rho*rho)
return 0.0;
return 1.0/(M_PI*rho*rho);
--
GitLab