Dear all I'm trying to import a mesh generated with Gmsh. It is composed by two volumes with adjacent faces. Deal.ii imports the mesh without problem but then I realize that it has assigned a boundary indicator of 0 to the cell faces between the volumes. So, deal.ii is not recognizing them as interior faces. I have search for a solution using Gmsh, but apparently it is still an issue. In deal.ii I found this command for deleting duplicated nodes
void GridTools::delete_duplicated_vertices ( std::vector< Point< spacedim > > & all_vertices,std::vector< CellData< dim > > & cells,SubCellData &subcelldata,std::vector< unsigned int > & considered_vertices, const double tol = 1e-12) that I think could help me solve the problem (as I have read that duplicated interfaces in Gmsh might be the cause). Unfortunately I haven't found any example of usage in the tutorials or in the mailing list. It says that the function can be called by some GridIn::read_* functions, how is it called? and the parameters, I'm not sure how to define them, I found something similar in step 14: std::vector<Point<3> > vertices (triangulation_3d.get_vertices()); std::vector<CellData<3> > cells (triangulation_3d.n_cells(), CellData<3>()); std::vector< unsigned int > considered_vertices(); and then maybe I can pass them as: GridTools::delete_duplicated_vertices(vertices,cells,SubCellData(),considered_vertices,1e-12); As you can see, I'm a bit lost. Any help would be highly appreciated. Regards Javier _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
