]> https://gitweb.dealii.org/ - dealii.git/commitdiff
FE_PolyTensor: don't assume the default line orientation is 1. 18212/head
authorDavid Wells <drwells@email.unc.edu>
Sun, 9 Mar 2025 17:01:34 +0000 (13:01 -0400)
committerDavid Wells <drwells@email.unc.edu>
Sun, 9 Mar 2025 17:02:20 +0000 (13:02 -0400)
source/fe/fe_poly_tensor.cc

index 6b5c3f739f2aa9874eacfddaefe50dc9ecc4a305..f17a13ee80e7c043e32cd6ac7b7aeab1395d8887 100644 (file)
@@ -166,7 +166,8 @@ namespace internal
         const unsigned int k = fe.tensor_degree() - 1;
 
         for (unsigned int l = 0; l < GeometryInfo<2>::lines_per_cell; ++l)
-          if (!(cell->line_orientation(l)) &&
+          if (cell->line_orientation(l) !=
+                numbers::default_geometric_orientation &&
               mapping_kind[0] == mapping_nedelec)
             {
               if (k == 0)
@@ -236,7 +237,8 @@ namespace internal
         // Here we adjust only the line (edge) dofs. The line dofs need only
         // sign adjustment. That is, no permutation of the line dofs is needed.
         for (unsigned int l = 0; l < GeometryInfo<3>::lines_per_cell; ++l)
-          if (!(cell->line_orientation(l)) &&
+          if (cell->line_orientation(l) !=
+                numbers::default_geometric_orientation &&
               mapping_kind[0] == mapping_nedelec)
             {
               if (k == 0)

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.