You can probably fix this within gmsh. Gmsh will save all faces for which "Physical Line/Surface/Volume" has been assigned. Just dont assign a physical tag to the interior adjacent faces.
praveen On Mon, May 21, 2012 at 7:33 PM, Javier Muñoz <[email protected]>wrote: > 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 >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
