Skip to content

Commit 99bb23d

Browse files
committed
Allow nodal quadrature on pyramids in Elem::has_invertible_map() tests
1 parent 7c7f467 commit 99bb23d

File tree

1 file changed

+3
-0
lines changed

1 file changed

+3
-0
lines changed

src/geom/elem.C

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2424,6 +2424,9 @@ bool Elem::has_invertible_map(Real /*tol*/) const
24242424
FEMap fe_map;
24252425
auto & jac = fe_map.get_jacobian();
24262426

2427+
// We have a separate check for is_singular_node() below, so in this
2428+
// case its "OK" to do nodal quadrature on pyramids.
2429+
qnodal.allow_nodal_pyramid_quadrature = true;
24272430
qnodal.init(this->type());
24282431
auto & qp = qnodal.get_points();
24292432
libmesh_assert_equal_to(qp.size(), this->n_nodes());

0 commit comments

Comments
 (0)