Skip to content

Conversation

@pbehne
Copy link
Contributor

@pbehne pbehne commented Nov 14, 2025

Eliminated repeat accumulations, switched c-style for loops to range-based.

DO NOT MERGE until #31917 is merged.

Refs. #31795.

Given a cut line, the XYMeshCutter deletes elements on one side of the
line.
When the line cuts through quad elements, those elements are converted
to tris and refined to form a new mesh boundary coincident with the cut
line.
To determine whether the cut line cuts through a given element, each
vertex of the element is plugged into the line equation f(x, y) = ax + by + c.
An (x, y) point is on the line if f(x, y) = 0.
Vertices residing on one side of the line have positive f(x, y) and
vertices residing on the other side of the line have negative f(x, y).
If two vertices of an element have differently signed f values, the
element is deemed to be cut by the line and is marked for conversion to
tri elements and further refinement.

This commit patches a bug in the above logic.
Previously, the mesh generator did not check if a given vertex lies on
the cut line.
These vertices were assigned to one side of the line.
In the case of a quad that is not cut by the line, but has an edge
coincident with the line, this can cause those elements to be
incorrectly marked for tri conversion and refinement.
While the final cut mesh was correct, this results in an irrelevant, confusing
warning that quad elements are being converted to tri elements.
This commit adds checks to determine and ignore vertices that lie on the
cut line.
These vertices are no longer used when determining whether an element is cut
by the line.
Ref. idaholab#31795
@moosebuild
Copy link
Contributor

Job Test, step Results summary on 714a7c9 wanted to post the following:

Framework test summary

Compared against b79a48d in job civet.inl.gov/job/3380566.

No change

Modules test summary

Compared against b79a48d in job civet.inl.gov/job/3380566.

No change

@moosebuild
Copy link
Contributor

Job Documentation, step Docs: sync website on 714a7c9 wanted to post the following:

View the site here

This comment will be updated on new commits.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants