From: Buğrahan Temür <bugrahan.temuer@tum.de>
Date: Wed, 15 Mar 2023 14:39:12 +0000 (+0100)
Subject: Fix bug in inverse JxW values
X-Git-Tag: v9.5.0-rc1~469^2~5
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=93729ce3ce192cd11b3d0790cff91e0900e0cc9d;p=dealii.git

Fix bug in inverse JxW values

wrong number of dofs for fe_degree=-1
---

diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h
index e2c5334900..320ee734b8 100644
--- a/include/deal.II/matrix_free/operators.h
+++ b/include/deal.II/matrix_free/operators.h
@@ -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,