From: Peter Munch Date: Sun, 22 Nov 2020 06:00:57 +0000 (+0100) Subject: Work on some asserts in MatrixFree X-Git-Tag: v9.3.0-rc1~869^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11227%2Fhead;p=dealii.git Work on some asserts in MatrixFree --- diff --git a/include/deal.II/matrix_free/evaluation_kernels.h b/include/deal.II/matrix_free/evaluation_kernels.h index 036b628f59..dff4ba1e00 100644 --- a/include/deal.II/matrix_free/evaluation_kernels.h +++ b/include/deal.II/matrix_free/evaluation_kernels.h @@ -2504,6 +2504,9 @@ namespace internal const unsigned int face_orientation, const Table<2, unsigned int> &orientation_map) { + Assert(data.element_type != MatrixFreeFunctions::tensor_none, + ExcNotImplemented()); + constexpr unsigned int static_dofs_per_face = fe_degree > -1 ? Utilities::pow(fe_degree + 1, dim - 1) : numbers::invalid_unsigned_int; @@ -2594,6 +2597,9 @@ namespace internal const unsigned int face_orientation, const Table<2, unsigned int> &orientation_map) { + Assert(data.element_type != MatrixFreeFunctions::tensor_none, + ExcNotImplemented()); + if (face_orientation) adjust_for_face_orientation(dim, n_components, diff --git a/include/deal.II/matrix_free/mapping_info.templates.h b/include/deal.II/matrix_free/mapping_info.templates.h index 02ef441e58..47e0d5c252 100644 --- a/include/deal.II/matrix_free/mapping_info.templates.h +++ b/include/deal.II/matrix_free/mapping_info.templates.h @@ -376,12 +376,17 @@ namespace internal flag &= quad[my_q][q].get_tensor_basis()[0] == quad[my_q][q].get_tensor_basis()[i]; } -#ifdef DEAL_II_WITH_SIMPLEX_SUPPORT + if (flag == false) - cell_data[my_q].descriptor[q].initialize(quad[my_q][q], - update_default); - else + { +#ifdef DEAL_II_WITH_SIMPLEX_SUPPORT + cell_data[my_q].descriptor[q].initialize(quad[my_q][q], + update_default); +#else + Assert(false, ExcNotImplemented()); #endif + } + else cell_data[my_q].descriptor[q].initialize( quad[my_q][q].get_tensor_basis()[0], update_default); }