From 7408b4590010d11cb15d81e25aff66a58224948f Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 25 Aug 2022 09:26:02 +0200 Subject: [PATCH] MatrixFree: Fix bug in setup of new cell_index for MG --- include/deal.II/matrix_free/matrix_free.templates.h | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/include/deal.II/matrix_free/matrix_free.templates.h b/include/deal.II/matrix_free/matrix_free.templates.h index 4eef8291f0..e0272e9f8e 100644 --- a/include/deal.II/matrix_free/matrix_free.templates.h +++ b/include/deal.II/matrix_free/matrix_free.templates.h @@ -869,7 +869,9 @@ MatrixFree::internal_reinit( mf_cell_indices.resize((mg_level == numbers::invalid_unsigned_int) ? tria.n_active_cells() : - tria.n_cells(mg_level), + (mg_level < tria.n_levels() ? + tria.n_raw_cells(mg_level) : + 0), numbers::invalid_unsigned_int); for (unsigned int cell = 0; -- 2.39.5