]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove calls to said functions on the library.
authorDavid Wells <drwells@email.unc.edu>
Sat, 30 Mar 2024 12:24:29 +0000 (08:24 -0400)
committerDavid Wells <drwells@email.unc.edu>
Mon, 1 Apr 2024 13:09:58 +0000 (09:09 -0400)
14 files changed:
source/base/qprojector.cc
source/fe/fe_abf.cc
source/fe/fe_bdm.cc
source/fe/fe_hermite.cc
source/fe/fe_nedelec.cc
source/fe/fe_nedelec_sz.cc
source/fe/fe_poly_tensor.cc
source/fe/fe_raviart_thomas.cc
source/fe/fe_rt_bubbles.cc
source/fe/mapping_cartesian.cc
source/fe/mapping_fe.cc
source/fe/mapping_fe_field.cc
source/fe/mapping_manifold.cc
source/fe/mapping_q.cc

index 0f6c8817e534cd0d15f799e98e45b05ca6e6d912..a9dda2cef17a5547c78fb663951177889e65ac3d 100644 (file)
@@ -1618,9 +1618,7 @@ QProjector<1>::DataSetDescriptor::subface(
   const ReferenceCell &reference_cell,
   const unsigned int   face_no,
   const unsigned int   subface_no,
-  const bool,
-  const bool,
-  const bool,
+  const unsigned char /*combined_orientation*/,
   const unsigned int n_quadrature_points,
   const internal::SubfaceCase<1>)
 {
@@ -1643,9 +1641,7 @@ QProjector<2>::DataSetDescriptor::subface(
   const ReferenceCell &reference_cell,
   const unsigned int   face_no,
   const unsigned int   subface_no,
-  const bool,
-  const bool,
-  const bool,
+  const unsigned char /*combined_orientation*/,
   const unsigned int n_quadrature_points,
   const internal::SubfaceCase<2>)
 {
@@ -1668,14 +1664,15 @@ QProjector<3>::DataSetDescriptor::subface(
   const ReferenceCell           &reference_cell,
   const unsigned int             face_no,
   const unsigned int             subface_no,
-  const bool                     face_orientation,
-  const bool                     face_flip,
-  const bool                     face_rotation,
+  const unsigned char            combined_orientation,
   const unsigned int             n_quadrature_points,
   const internal::SubfaceCase<3> ref_case)
 {
   const unsigned int dim = 3;
 
+  const auto [face_orientation, face_rotation, face_flip] =
+    internal::split_face_orientation(combined_orientation);
+
   Assert(reference_cell == ReferenceCells::Hexahedron, ExcNotImplemented());
   (void)reference_cell;
 
index 08519959e4f4760e3e62be4d435337a18600daa1..44b8fc65677e6add93b7bfae278e966ed17124f4 100644 (file)
@@ -623,7 +623,10 @@ FE_ABF<dim>::convert_generalized_support_point_values_to_dof_values(
           // TODO: Check what the face_orientation, face_flip and face_rotation
           // have to be in 3d
           unsigned int k = QProjector<dim>::DataSetDescriptor::face(
-            this->reference_cell(), face, false, false, false, n_face_points);
+            this->reference_cell(),
+            face,
+            ReferenceCell::default_combined_face_orientation(),
+            n_face_points);
           for (unsigned int i = 0; i < boundary_weights_abf.size(1); ++i)
             nodal_values[start_abf_dofs + i] +=
               n_orient * boundary_weights_abf(k + fp, i) *
index 2a8b2ff821dab19f62ea29ec2d0189a2465a0c6a..00d60cba06769a2bb59f512ca8156b7682af0d51 100644 (file)
@@ -388,14 +388,12 @@ FE_BDM<dim>::initialize_support_points(const unsigned int deg)
   for (unsigned int k = 0;
        k < face_points.size() * GeometryInfo<dim>::faces_per_cell;
        ++k)
-    this->generalized_support_points[k] = faces.point(
-      k +
-      QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
-                                               0,
-                                               true,
-                                               false,
-                                               false,
-                                               this->n_dofs_per_face(face_no)));
+    this->generalized_support_points[k] =
+      faces.point(k + QProjector<dim>::DataSetDescriptor::face(
+                        this->reference_cell(),
+                        0,
+                        ReferenceCell::default_combined_face_orientation(),
+                        this->n_dofs_per_face(face_no)));
 
   // Currently, for backward compatibility, we do not use moments, but
   // point values on faces in 2d. In 3d, this is impossible, since the
index e98507f646c05c9567e58123aee778dd18b959d5..6cd8ce03c15f7c9687214070a105e44eadc14cec 100644 (file)
@@ -887,9 +887,7 @@ FE_Hermite<dim, spacedim>::fill_fe_face_values(
     QProjector<dim>::DataSetDescriptor::face(
       ReferenceCells::get_hypercube<dim>(),
       face_no,
-      cell->face_orientation(face_no),
-      cell->face_flip(face_no),
-      cell->face_rotation(face_no),
+      cell->combined_face_orientation(face_no),
       quadrature[0].size());
 
   const UpdateFlags flags(fe_data.update_each);
index 6e2c95166e40489f08833c738c4622dfb48b7318..207d3a846b7aa48b7fa8bd75801b5ca8c95c93f8 100644 (file)
@@ -325,12 +325,11 @@ FE_Nedelec<2>::initialize_support_points(const unsigned int order)
                ++line)
             this->generalized_support_points[line * n_edge_points + q_point] =
               edge_quadrature.point(
-                QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
-                                                         line,
-                                                         true,
-                                                         false,
-                                                         false,
-                                                         n_edge_points) +
+                QProjector<dim>::DataSetDescriptor::face(
+                  this->reference_cell(),
+                  line,
+                  ReferenceCell::default_combined_face_orientation(),
+                  n_edge_points) +
                 q_point);
 
           for (unsigned int i = 0; i < order; ++i)
@@ -356,12 +355,11 @@ FE_Nedelec<2>::initialize_support_points(const unsigned int order)
         for (unsigned int q_point = 0; q_point < n_edge_points; ++q_point)
           this->generalized_support_points[line * n_edge_points + q_point] =
             edge_quadrature.point(
-              QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
-                                                       line,
-                                                       true,
-                                                       false,
-                                                       false,
-                                                       n_edge_points) +
+              QProjector<dim>::DataSetDescriptor::face(
+                this->reference_cell(),
+                line,
+                ReferenceCell::default_combined_face_orientation(),
+                n_edge_points) +
               q_point);
     }
 }
@@ -429,9 +427,7 @@ FE_Nedelec<3>::initialize_support_points(const unsigned int order)
                 QProjector<dim - 1>::DataSetDescriptor::face(
                   ReferenceCells::get_hypercube<dim - 1>(),
                   line,
-                  true,
-                  false,
-                  false,
+                  ReferenceCell::default_combined_face_orientation(),
                   n_edge_points) +
                 q_point);
 
@@ -497,12 +493,11 @@ FE_Nedelec<3>::initialize_support_points(const unsigned int order)
                                              GeometryInfo<dim>::lines_per_cell *
                                                n_edge_points] =
               face_quadrature.point(
-                QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
-                                                         face,
-                                                         true,
-                                                         false,
-                                                         false,
-                                                         n_face_points) +
+                QProjector<dim>::DataSetDescriptor::face(
+                  this->reference_cell(),
+                  face,
+                  ReferenceCell::default_combined_face_orientation(),
+                  n_face_points) +
                 q_point);
           }
 
@@ -530,9 +525,7 @@ FE_Nedelec<3>::initialize_support_points(const unsigned int order)
                 QProjector<dim - 1>::DataSetDescriptor::face(
                   ReferenceCells::get_hypercube<dim - 1>(),
                   line,
-                  true,
-                  false,
-                  false,
+                  ReferenceCell::default_combined_face_orientation(),
                   n_edge_points) +
                 q_point);
 
index a870d6acef1d5cebd09755fa99a927d1ef8a1aa5..0a4e04b2e2656a21fe7eeb10867242bfd1c6144a 100644 (file)
@@ -3142,9 +3142,8 @@ FE_NedelecSZ<dim, spacedim>::fill_fe_face_values(
   const auto         offset =
     QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
                                              face_no,
-                                             cell->face_orientation(face_no),
-                                             cell->face_flip(face_no),
-                                             cell->face_rotation(face_no),
+                                             cell->combined_face_orientation(
+                                               face_no),
                                              n_q_points);
 
   if (flags & update_values)
index 8f054174afad7dec0b8cafa653492042617ecfc1..bea29e6db01981669f9b008a97382fb78de860e6 100644 (file)
@@ -1080,9 +1080,8 @@ FE_PolyTensor<dim, spacedim>::fill_fe_face_values(
   const auto offset =
     QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
                                              face_no,
-                                             cell->face_orientation(face_no),
-                                             cell->face_flip(face_no),
-                                             cell->face_rotation(face_no),
+                                             cell->combined_face_orientation(
+                                               face_no),
                                              n_q_points);
 
   // TODO: Size assertions
@@ -1764,9 +1763,8 @@ FE_PolyTensor<dim, spacedim>::fill_fe_subface_values(
     QProjector<dim>::DataSetDescriptor::subface(this->reference_cell(),
                                                 face_no,
                                                 sub_no,
-                                                cell->face_orientation(face_no),
-                                                cell->face_flip(face_no),
-                                                cell->face_rotation(face_no),
+                                                cell->combined_face_orientation(
+                                                  face_no),
                                                 n_q_points,
                                                 cell->subface_case(face_no));
 
index 2b7f559397fb4762fba2f814e53c1b28830a4d34..f5c194b6bcfef70ec2c2eed8442532b476a75e7b 100644 (file)
@@ -211,9 +211,11 @@ FE_RaviartThomas<dim>::initialize_support_points(const unsigned int deg)
           // Enter the support point
           // into the vector
           this->generalized_support_points[current] = faces.point(
-            current +
-            QProjector<dim>::DataSetDescriptor::face(
-              this->reference_cell(), 0, true, false, false, n_face_points));
+            current + QProjector<dim>::DataSetDescriptor::face(
+                        this->reference_cell(),
+                        0,
+                        ReferenceCell::default_combined_face_orientation(),
+                        n_face_points));
         }
     }
 
index cbda9a1e9de294e124e5838d0a0f766a437df829..a9339f68b24e2117a7e5bc5dabd5a91e3121167a 100644 (file)
@@ -187,14 +187,12 @@ FE_RT_Bubbles<dim>::initialize_support_points(const unsigned int deg)
       for (unsigned int k = 0; k < this->n_dofs_per_face(face_no) *
                                      GeometryInfo<dim>::faces_per_cell;
            ++k)
-        this->generalized_support_points[k] = faces.point(
-          k + QProjector<dim>::DataSetDescriptor::face(this->reference_cell(),
-                                                       0,
-                                                       true,
-                                                       false,
-                                                       false,
-                                                       this->n_dofs_per_face(
-                                                         face_no)));
+        this->generalized_support_points[k] =
+          faces.point(k + QProjector<dim>::DataSetDescriptor::face(
+                            this->reference_cell(),
+                            0,
+                            ReferenceCell::default_combined_face_orientation(),
+                            this->n_dofs_per_face(face_no)));
 
       current =
         this->n_dofs_per_face(face_no) * GeometryInfo<dim>::faces_per_cell;
index f8842f8a1aaa5c2ab332cad55edfe61c7cfdec5c..85be4007eeac9f555ffb8375f881bd9640749652 100644 (file)
@@ -293,9 +293,7 @@ MappingCartesian<dim, spacedim>::maybe_update_face_quadrature_points(
       const auto offset = QProjector<dim>::DataSetDescriptor::face(
         ReferenceCells::get_hypercube<dim>(),
         face_no,
-        cell->face_orientation(face_no),
-        cell->face_flip(face_no),
-        cell->face_rotation(face_no),
+        cell->combined_face_orientation(face_no),
         quadrature_points.size());
 
 
@@ -327,9 +325,7 @@ MappingCartesian<dim, spacedim>::maybe_update_subface_quadrature_points(
         ReferenceCells::get_hypercube<dim>(),
         face_no,
         sub_no,
-        cell->face_orientation(face_no),
-        cell->face_flip(face_no),
-        cell->face_rotation(face_no),
+        cell->combined_face_orientation(face_no),
         quadrature_points.size(),
         cell->subface_case(face_no));
 
index 7a1322aa20e4b33b15da247729c623d67d72b594..9707bb7a811eb14082a4c46f465c6a171f6e1584 100644 (file)
@@ -1616,9 +1616,8 @@ MappingFE<dim, spacedim>::fill_fe_face_values(
     numbers::invalid_unsigned_int,
     QProjector<dim>::DataSetDescriptor::face(this->fe->reference_cell(),
                                              face_no,
-                                             cell->face_orientation(face_no),
-                                             cell->face_flip(face_no),
-                                             cell->face_rotation(face_no),
+                                             cell->combined_face_orientation(
+                                               face_no),
                                              quadrature),
     quadrature[quadrature.size() == 1 ? 0 : face_no],
     data,
@@ -1664,9 +1663,8 @@ MappingFE<dim, spacedim>::fill_fe_subface_values(
     QProjector<dim>::DataSetDescriptor::subface(this->fe->reference_cell(),
                                                 face_no,
                                                 subface_no,
-                                                cell->face_orientation(face_no),
-                                                cell->face_flip(face_no),
-                                                cell->face_rotation(face_no),
+                                                cell->combined_face_orientation(
+                                                  face_no),
                                                 quadrature.size(),
                                                 cell->subface_case(face_no)),
     quadrature,
index 9210ecbef4ede7b1d19ce16b4d99d855abc19ecb..5fa44c8548618544b82e67001fad6a609384711c 100644 (file)
@@ -1718,9 +1718,8 @@ MappingFEField<dim, spacedim, VectorType>::fill_fe_face_values(
       numbers::invalid_unsigned_int,
       QProjector<dim>::DataSetDescriptor::face(reference_cell,
                                                face_no,
-                                               cell->face_orientation(face_no),
-                                               cell->face_flip(face_no),
-                                               cell->face_rotation(face_no),
+                                               cell->combined_face_orientation(
+                                                 face_no),
                                                quadrature[0].size()),
       data,
       euler_dof_handler->get_fe(),
@@ -1759,9 +1758,8 @@ MappingFEField<dim, spacedim, VectorType>::fill_fe_subface_values(
     QProjector<dim>::DataSetDescriptor::subface(reference_cell,
                                                 face_no,
                                                 subface_no,
-                                                cell->face_orientation(face_no),
-                                                cell->face_flip(face_no),
-                                                cell->face_rotation(face_no),
+                                                cell->combined_face_orientation(
+                                                  face_no),
                                                 quadrature.size(),
                                                 cell->subface_case(face_no)),
     data,
index a9e72275d43db08b2f33dd727cdd295c472291c4..5ae5ffb3e633219a42e1bfb4aec0699d67e6bdc1 100644 (file)
@@ -1228,9 +1228,7 @@ MappingManifold<dim, spacedim>::fill_fe_face_values(
     QProjector<dim>::DataSetDescriptor::face(
       ReferenceCells::get_hypercube<dim>(),
       face_no,
-      cell->face_orientation(face_no),
-      cell->face_flip(face_no),
-      cell->face_rotation(face_no),
+      cell->combined_face_orientation(face_no),
       quadrature[0].size()),
     quadrature[0],
     data,
@@ -1264,9 +1262,7 @@ MappingManifold<dim, spacedim>::fill_fe_subface_values(
       ReferenceCells::get_hypercube<dim>(),
       face_no,
       subface_no,
-      cell->face_orientation(face_no),
-      cell->face_flip(face_no),
-      cell->face_rotation(face_no),
+      cell->combined_face_orientation(face_no),
       quadrature.size(),
       cell->subface_case(face_no)),
     quadrature,
index e5931b091f8b002356d590fcc7b8fde92121923c..f1631738add1e694b3a99e27ccc91fd1fd18fffe 100644 (file)
@@ -1090,9 +1090,7 @@ MappingQ<dim, spacedim>::fill_fe_subface_values(
       ReferenceCells::get_hypercube<dim>(),
       face_no,
       subface_no,
-      cell->face_orientation(face_no),
-      cell->face_flip(face_no),
-      cell->face_rotation(face_no),
+      cell->combined_face_orientation(face_no),
       quadrature.size(),
       cell->subface_case(face_no)),
     quadrature,

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.