Hello,
In your code, for each input cell, you add a copy of each of its vertices - so the mesh will consist of topologically isolated tet cells that have vertices that are not identical, but share one position.
You should only call add_vertex once for each vertex and re-use the returned VertexHandle for all cells that use this vertex.
Good luck!
Best Regards
Martin Heistermann
On 8/29/19 5:10 PM, S Vijai Kumar wrote:
Hello there,
I am not sure if this is the place to ask questions about Open Volume Mesh. But i couldn't find any other mailing list. Hence i am posting here.
I am trying to create a mesh of type *OpenVolumeMesh::GeometricTetrahedralMeshV3d HalfFaceMesh;*
To this, i am adding tetrahedrons as follows
*O**penVolumeMesh::CellPropertyT< unsigned int > TetSphereAssociativeProperty = HalfFaceMesh.request_cell_property< unsigned int >("TET_SPHERE_ASSOCIATION_ID");
*int num_vol_elements = nglib::Ng_GetNE(CurrentMesh);**
*unsigned int tetra_cnt = 1;**
*for(unsigned int tetra = 0; tetra < num_vol_elements; ++tetra)**
*{**
* int CoordinateArray[4];**
* nglib::Ng_GetVolumeElement(CurrentMesh, tetra + 1, CoordinateArray);**
* long v1 = CoordinateArray[0];**
* long v2 = CoordinateArray[1];**
* long v3 = CoordinateArray[2];**
* long v4 = CoordinateArray[3];**
* VertexHandle _v1 = HalfFaceMesh.add_vertex(meshPointToVector3D(v1));**
* VertexHandle _v2 = HalfFaceMesh.add_vertex(meshPointToVector3D(v2));**
* VertexHandle _v3 = HalfFaceMesh.add_vertex(meshPointToVector3D(v3));**
* VertexHandle _v4 = HalfFaceMesh.add_vertex(meshPointToVector3D(v4));**
* CellHandle current_cell = HalfFaceMesh.add_cell(_v1, _v2, _v3, _v4);**
* TetSphereAssociativeProperty[current_cell] = tetra_cnt;**
* ++tetra_cnt;**
*}*
The issue is when i try to get neighbors of a particular tetrahedron, i run into issues.
I have tried *CellCellIterator* and it gives me nothing. So i wrote the code to do it in the most naive way possible also (i know its quadratic).
*std::set< CellHandle > getUniqueNeighbouringTetsForGivenTets(CellHandle input_tetrahedron)**
*{**
* std::set< CellHandle > matching_cells;**
* std::set< FaceHandle > matching_faces = getFacesForGivenCell(input_tetrahedron);**
* HalfFaceMesh.enable_bottom_up_incidences(true);**
* for(auto mf_iter = matching_faces.begin(); mf_iter != matching_faces.end(); ++mf_iter)**
* {**
* for(CellIterator c_it = HalfFaceMesh.cells_begin(); c_it != HalfFaceMesh.cells_end(); ++c_it)**
* {**
* if(*c_it != input_tetrahedron)**
* {**
* std::set< FaceHandle > intermediate_matching_faces = getFacesForGivenCell(*c_it);**
* for(auto mff_iter = intermediate_matching_faces.begin(); mff_iter != intermediate_matching_faces.end(); ++mff_iter)**
* {**
* if(*mff_iter == *mf_iter)**
* {**
* matching_cells.emplace(*c_it);**
* }**
* }**
* }**
* }**
* }**
* return matching_cells;**
*}*
I still do not get any neighbours at all. It almost seems like the neighbors are not associated automatically and basically one of the half faces of every face is just some garbage with an id -1. Should i manually establish the half face connectivity between adjacent cells while building a mesh in open volume mesh ? This was not the case in OpenMesh where we simply add polygons and edges are automatically split into half edges.
I am stuck and any help would be appreciated. Thanks in advance.
Kind regards, S Vijai Kumar
Sent with ProtonMail https://protonmail.com Secure Email.
OpenMesh mailing list -- openmesh@lists.rwth-aachen.de To unsubscribe send an email to openmesh-leave@lists.rwth-aachen.de https://lists.rwth-aachen.de/postorius/lists/openmesh.lists.rwth-aachen.de