From 93729ce3ce192cd11b3d0790cff91e0900e0cc9d Mon Sep 17 00:00:00 2001 From: =?utf8?q?Bu=C4=9Frahan=20Tem=C3=BCr?= Date: Wed, 15 Mar 2023 15:39:12 +0100 Subject: [PATCH] Fix bug in inverse JxW values wrong number of dofs for fe_degree=-1 --- include/deal.II/matrix_free/operators.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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, -- 2.39.5