Skip to content

Commit 11dd1c9

Browse files
committed
ref - Ensure at-points operators skip cells with no points
1 parent 9b9f611 commit 11dd1c9

File tree

1 file changed

+1
-0
lines changed

1 file changed

+1
-0
lines changed

backends/ref/ceed-ref-operator.c

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1054,6 +1054,7 @@ static int CeedOperatorApplyAddAtPoints_Ref(CeedOperator op, CeedVector in_vec,
10541054
// Setup points for element
10551055
CeedCallBackend(CeedElemRestrictionApplyAtPointsInElement(rstr_points, e, CEED_NOTRANSPOSE, point_coords, impl->point_coords_elem, request));
10561056
CeedCallBackend(CeedElemRestrictionGetNumPointsInElement(rstr_points, e, &num_points));
1057+
if (num_points <= 0) continue;
10571058

10581059
// Input basis apply
10591060
CeedCallBackend(CeedOperatorInputBasisAtPoints_Ref(e, num_points_offset, num_points, qf_input_fields, op_input_fields, num_input_fields, in_vec,

0 commit comments

Comments
 (0)