From 9d07ff8002da4d5a598640d714c52967d26fefa3 Mon Sep 17 00:00:00 2001 From: David Wells Date: Sun, 9 Mar 2025 13:01:34 -0400 Subject: [PATCH] FE_PolyTensor: don't assume the default line orientation is 1. --- source/fe/fe_poly_tensor.cc | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/source/fe/fe_poly_tensor.cc b/source/fe/fe_poly_tensor.cc index 6b5c3f739f..f17a13ee80 100644 --- a/source/fe/fe_poly_tensor.cc +++ b/source/fe/fe_poly_tensor.cc @@ -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) -- 2.39.5