Adapt crashes for mixed-element 2D meshes with frozen quad layers
Description
PUMI adapt crashes when operating on 2D mixed-element meshes (tri + quad) if quad layers are frozen during adaptation.
This was observed while using PUMI in CFD workflow, where adaptation is intentionally restricted to triangle elements while quadrilateral boundary-layer elements are excluded.
Expected behavior
Adaptation should proceed safely on triangle elements in mixed-element 2D meshes even when quad layers are frozen.
Actual behavior
Adapt crashes in the adapt/maQuality path when encountering frozen quad elements in a mixed-element 2D mesh.
Versions
- PUMI v4.1.0 (release)
- Also reproducible prior to the fix on current develop
Notes
This use case is common in CFD workflows that employ quad boundary-layer meshes with triangular farfield regions, where only the triangular portion is adapted.
I have a proposed fix that adds a guard in maQuality to allow safe adaptation of triangles in this configuration without affecting existing workflows.