From 64d9e8b793db096ce58b2aabbac4a0538372cbe2 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Mon, 7 Jul 2025 09:50:48 +0200 Subject: [PATCH] MappingP1: Small improvement to avoid some divisions --- source/fe/mapping_p1.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/fe/mapping_p1.cc b/source/fe/mapping_p1.cc index f224a58b95..96801aac15 100644 --- a/source/fe/mapping_p1.cc +++ b/source/fe/mapping_p1.cc @@ -578,7 +578,7 @@ MappingP1::transform( auto transformation = data.linear_component; Assert(data.determinant > 0.0, ExcDivideByZero()); for (unsigned int d = 0; d < spacedim; ++d) - transformation[d] /= data.determinant; + transformation[d] *= 1.0 / data.determinant; for (unsigned int i = 0; i < output.size(); ++i) output[i] = apply_transformation(transformation, input[i]); return; @@ -696,7 +696,7 @@ MappingP1::transform( auto scaled_contravariant = data.linear_component; for (unsigned int d = 0; d < spacedim; ++d) - scaled_contravariant[d] /= data.determinant; + scaled_contravariant[d] *= 1.0 / data.determinant; for (unsigned int i = 0; i < output.size(); ++i) { const DerivativeForm<1, spacedim, dim> A = -- 2.39.5