]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Allow to include FEFaceEvaluation in std::vector
authorMartin Kronbichler <martin.kronbichler@rub.de>
Mon, 19 Aug 2024 19:28:34 +0000 (21:28 +0200)
committerMartin Kronbichler <martin.kronbichler@rub.de>
Tue, 20 Aug 2024 07:57:03 +0000 (09:57 +0200)
include/deal.II/matrix_free/fe_evaluation.h

index fda4197cce2d30e6f8f25d4876eac8bc5ce2a255..5246d5ac29c0ce289cc917d9deb904d2ddcef100 100644 (file)
@@ -2517,8 +2517,16 @@ inline FEEvaluationBase<dim,
           other.mapped_geometry->get_fe_values().get_mapping(),
           other.mapped_geometry->get_quadrature(),
           other.mapped_geometry->get_fe_values().get_update_flags());
-      this->mapping_data = &this->mapped_geometry->get_data_storage();
-      this->cell         = 0;
+
+      if constexpr (is_face == false)
+        this->mapping_data = &this->mapped_geometry->get_data_storage();
+      else
+        Assert(false,
+               ExcNotImplemented("On-the-fly geometry evaluation with "
+                                 "face evaluators is not currently "
+                                 "implemented!"));
+
+      this->cell = 0;
 
       this->jacobian =
         this->mapped_geometry->get_data_storage().jacobians[0].begin();
@@ -2585,8 +2593,16 @@ operator=(const FEEvaluationBase<dim,
           other.mapped_geometry->get_fe_values().get_mapping(),
           other.mapped_geometry->get_quadrature(),
           other.mapped_geometry->get_fe_values().get_update_flags());
-      this->cell         = 0;
-      this->mapping_data = &this->mapped_geometry->get_data_storage();
+
+      if constexpr (is_face == false)
+        this->mapping_data = &this->mapped_geometry->get_data_storage();
+      else
+        Assert(false,
+               ExcNotImplemented("On-the-fly geometry evaluation with "
+                                 "face evaluators is not currently "
+                                 "implemented!"));
+      this->cell = 0;
+
       this->jacobian =
         this->mapped_geometry->get_data_storage().jacobians[0].begin();
       this->J_value =

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.