]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Avoid the use of GeometryInfo in FE_RaviartThomas. 18040/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Mon, 27 Jan 2025 22:10:00 +0000 (15:10 -0700)
committerWolfgang Bangerth <bangerth@colostate.edu>
Tue, 28 Jan 2025 03:56:43 +0000 (20:56 -0700)
source/fe/fe_raviart_thomas.cc

index 81e8c879adca122ae53cc93965598497d2676e67..3585e105a2f74ba7c9f5260387c1821e489ccb9b 100644 (file)
@@ -92,7 +92,9 @@ FE_RaviartThomas<dim>::FE_RaviartThomas(const unsigned int deg)
 
   // TODO[TL]: for anisotropic refinement we will probably need a table of
   // submatrices with an array for each refine case
-  FullMatrix<double> face_embeddings[GeometryInfo<dim>::max_children_per_face];
+  std::vector<FullMatrix<double>> face_embeddings(
+    this->reference_cell().face_reference_cell(0).template n_children<dim - 1>(
+      RefinementCase<dim - 1>::isotropic_refinement));
   for (auto &face_embedding : face_embeddings)
     face_embedding.reinit(this->n_dofs_per_face(face_no),
                           this->n_dofs_per_face(face_no));

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.