diff --git a/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh b/dumux/multidomain/embedded/couplingmanager1d3d_kernel.hh index bbc33fde6de9b612c9ec6aa1ec32635900f33aa2..2fc693cf36cb5d0ada8716b0718851c26038adab 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);