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 e2c5334900e6c531cd05f30e3caa612a1fa26cc9..320ee734b8e618310fc7b213c6e328a849c55beb 100644 (file)
--- 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,