@kirma, OK, thanks. The only workaround that I found for this issue is that it seems to be scale dependent (a precision issue). If you use `m = OpenCascadeShapeSurfaceMeshToBoundaryMesh@
OpenCascadeShapeBooleanRegion@
RegionDifference[Ball[{0, 0, 0}, 100], Ball[{0, 100, 100}, 100]]` things work. You'd then need to rescale the result.
OpenCascadeShapeBooleanRegion@
RegionDifference[Ball[{0, 0, 0}, 100], Ball[{0, 100, 100}, 100]]` things work. You'd then need to rescale the result.