]> https://gitweb.dealii.org/ - dealii.git/commitdiff
MappingP1: Small improvement to avoid some divisions 18636/head
authorMartin Kronbichler <martin.kronbichler@rub.de>
Mon, 7 Jul 2025 07:50:48 +0000 (09:50 +0200)
committerMartin Kronbichler <martin.kronbichler@rub.de>
Mon, 7 Jul 2025 07:50:48 +0000 (09:50 +0200)
source/fe/mapping_p1.cc

index f224a58b95ba3dc9a1e36b4b6bad720e8a10334e..96801aac15d73aaa392f2c69f9321ede35f4d239 100644 (file)
@@ -578,7 +578,7 @@ MappingP1<dim, spacedim>::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<dim, spacedim>::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 =

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.