From: David Wells Date: Sun, 9 Mar 2025 17:01:34 +0000 (-0400) Subject: FE_PolyTensor: don't assume the default line orientation is 1. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F18212%2Fhead;p=dealii.git FE_PolyTensor: don't assume the default line orientation is 1. --- 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)