Skip to content

Conversation

@alexbenedicto
Copy link
Collaborator

Currently when using generateFractures feature from mesh-doctor, the default precision for vtkPoints objects is float.
This can become an issue when the input mesh uses double precision and we want to create a new instance of this input mesh with the split points to output.

In certain cases, it will mess up the node mapping from the input to the output mesh and we will get cells with invalid nodes ordering, therefore leading to invalid/zero/negative volumes.

This PR solves this issue by checking the precision of points from the input mesh and applying this precision to the new vtkPoints objects of the output mesh.

@alexbenedicto alexbenedicto self-assigned this Dec 15, 2025
@alexbenedicto alexbenedicto added flag: ready for review type: bug test-geos-integration Triggers the testing of geosPythonPackages import and integration in GEOS CI labels Dec 15, 2025
Copy link

@jhuang2601 jhuang2601 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks, it fixed the reported negative volume issue.

@alexbenedicto alexbenedicto merged commit 1ac2be7 into main Dec 15, 2025
56 checks passed
@alexbenedicto alexbenedicto deleted the benedicto/fix/mesh-doctor-keep-input_mesh-precision branch December 15, 2025 21:32
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

flag: ready for review test-geos-integration Triggers the testing of geosPythonPackages import and integration in GEOS CI type: bug

Projects

None yet

Development

Successfully merging this pull request may close these issues.

5 participants