Hello!
I'm having a weird problem with the method "triangulate_pointset()",
from the TetGenMeshInterface class. Essentially, it seems like it is
stalling after several calls. In my code, I have to calculate several
polyhedron intersections and mesh them. To do so, I built a
Mesh_Intersection class, which has the following method:
void Mesh_Intersection::triangulate_intersection(const
std::set<libMesh::Point> & input_points)
{
/*
* input_points : set containing the intersection
* polyhedron's vertices.
* m_libMesh_PolyhedronMesh : a libMesh::Mesh, member of the
Mesh_Intersection
* cleared when calling this method.
*/
m_libMesh_PolyhedronMesh.clear();
// Insert points into the mesh
for(std::set<libMesh::Point>::const_iterator it_set =
input_points.begin();
it_set != input_points.end();
++it_set)
{
m_libMesh_PolyhedronMesh.add_point(*it_set);
}
libMesh::TetGenMeshInterface
temp_tetgen_interface(m_libMesh_PolyhedronMesh);
temp_tetgen_interface.triangulate_pointset();
}
The code stalls at the last line. I already tested the input_points,
and they are valid. I'm using libMesh 1.0.0, compiled with Intel icpc v.
16.0.3 and Intel MPI v. 5.1.2, on Linux. The m_libMesh_PolyhedronMesh is
associated to a single processor. I don't think that it is a
communicator problem because it happens even with only one processor. I
think that the problem is associated to the compiler, since it works
well on a machine where I used gcc 4.9. Did anyone else have a similar
problem?
Thanks in advance,
Thiago Milanetto Schlittler
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
