Meshfix().repair() crashes Python code with no error message

[plebbyd]

I am iterating through a set of .STL files and repairing the meshes if they are broken using MeshFix().repair(), and the program is terminating without any error message, warning message, exception or any other indication that there is an issue.

One thing that I have noticed is that it is terminating on larger files.

[SterlingButters]

Same issue here but not just with the repair() method

[SerdarHelli]

me to

[rewanaa]

Same issue

[JfosterBGC]

I am running into this same issue in certain cases. I am trying to create a closed volume from the clipping a mesh and an ellips.
I have attached an offending example below. In this case, I am wondering if it is because the surfaces are so small (1.95 m max thickness. the other surfaces that worked well had thicknesses if 5 m or greater

c19_SEG_Flare Addition_Remesh.zip

import numpy as np
import pandas as pd
import pyvista as pv
import pymeshfix as pmf

row = {'MPA ID': '01_SEGully 4-0',
 'direction': 'DS',
 'Axis Point X': 377453.56,
 'Axis Point Y': 2089177.45,
 'Axis Point Z': 288.66,
 'Radius X,Z': 16.98,
 'Aspect Ratio': 0.94,
 'Max Depth': 1.950000048,
 'Rotation °': 63.9,
 'idx': 4964}`

top_surf_fn = 'c19_SEG_Flare Addition_Remesh.obj'
top_surf = top_surface_remesh = pv.read(top_surf_fn)
top_surf = top_surf.triangulate()

aspect_ratio = row["Aspect Ratio"]
radius_x = row["Radius X,Z"]
radius_z = radius_x
radius_y = radius_x * aspect_ratio

axis_x = row["Axis Point X"]
axis_y = row["Axis Point Y"]
axis_z = row["Axis Point Z"]

rotation = row["Rotation °"]

#Create Ellips based on the above Geometry 
ellipse = pv.ParametricEllipsoid(radius_x, radius_y, radius_z)
ellipse = ellipse.translate((axis_x, axis_y, axis_z))
ellipse = ellipse.rotate_z(angle=rotation, point=(ellipse.center_of_mass()))

# Clip Ellipse to Top Surface and Keep Largest
ellipse_btm = ellipse.clip_surface(top_surf)
ellipse_btm = ellipse_btm.extract_largest()
ellipse_top = top_surf.clip_surface(ellipse, invert = False)
ellipse_top = ellipse_top.extract_largest()

slide_mass = pv.merg
e([ellipse_btm,ellipse_top])
slide_mass = pmf.MeshFix(slide_mass)
slide_mass.repair()
slide_mass = slide_mass.mesh
is_manifold = slide_mass.is_manifold # Is the slide mass volume closed i.e. no holes
print (is_manifold)