]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix bug in inverse JxW values
authorBuğrahan Temür <bugrahan.temuer@tum.de>
Wed, 15 Mar 2023 14:39:12 +0000 (15:39 +0100)
committerBuğrahan Temür <bugrahan.temuer@tum.de>
Wed, 15 Mar 2023 14:39:12 +0000 (15:39 +0100)
wrong number of dofs for fe_degree=-1

include/deal.II/matrix_free/operators.h

index e2c5334900e6c531cd05f30e3caa612a1fa26cc9..320ee734b8e618310fc7b213c6e328a849c55beb 100644 (file)
@@ -1057,7 +1057,7 @@ namespace MatrixFreeOperators
       (fe_degree > -1) ?
         Utilities::pow(fe_degree + 1, dim) :
         Utilities::pow(fe_eval.get_shape_info().data.front().fe_degree + 1,
-                       dim - 1);
+                       dim);
 
     Assert(inverse_jxw.size() > 0 &&
              inverse_jxw.size() % dofs_per_component_on_cell == 0,

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.