cube and cylinder "touching" about ftetwild HOT 2 CLOSED

Hi,

The simplest solution to try first is to combine the cylinder and the cube into a single mesh file and input it to fTetWild (using --input only would work). The interface of the cylinder and the cube would be preserved by the algorithm. To separate the output tets for the two objects, you can check if the centers of the tets inside certain object.

from ftetwild.

Thank you very much Yixin
this is exactly the solution I adopted some days ago: indeed I was writing for sharing it: Some details:

1. the two - initially separated, distinct - solids (step file) are loaded into a viewer using an OpenCascade step reader
2. a non-manifold solid is created using the BOPAlgo_Builder class. In the resulting non-maniforld the base of the cylinder is a shared face
3. the triangulation of the composite solid is created using an OpenCascade class (BRepMesh_IncrementalMesh)
4. fTwtWild creates the tetrahedralization of the triangle soup (in the example I posted the soup is a watertight triangulation)
5. the barycenter of each tetrahedron is computed,: a tetrahedron "belongs" to a solid if the winding number of its center of mass, with respect to the triangulation of that solid, is greater than, say 0.5 (for problematic, non watertight, triangulations)
   I attach a picture of a results, with a mesh refinement also applied to one face and a segment.
   Grazie mille
   Giovanni

from ftetwild.