Just using pos_equal_eps for stitching slices together should work fine and avoids loss in floating point accuracy in cases where intersects are very close to one another (fixing debug assert case). The slice_join_eps was originally added in the C++ code to mask over problems of floating point accuracy in the intersects code.
Also when constructing the aabb_index from the slice start points no fuzzy epsilon is applied (they are added as points), epsilon is only applied during the querying. If it is applied both in creation and querying then the epsilon is effectively doubled, which is not desired.
Just using pos_equal_eps for stitching slices together should work fine and avoids loss in floating point accuracy in cases where intersects are very close to one another (fixing debug assert case). The slice_join_eps was originally added in the C++ code to mask over problems of floating point accuracy in the intersects code.
Also when constructing the aabb_index from the slice start points no fuzzy epsilon is applied (they are added as points), epsilon is only applied during the querying. If it is applied both in creation and querying then the epsilon is effectively doubled, which is not desired.