From: Martin Kronbichler Date: Mon, 7 Jul 2025 07:50:48 +0000 (+0200) Subject: MappingP1: Small improvement to avoid some divisions X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=64d9e8b793db096ce58b2aabbac4a0538372cbe2;p=dealii.git MappingP1: Small improvement to avoid some divisions --- 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 =