There are occasions in which duplicated vertices are actually intended.... Maybe you'd like to add a switch to the function on order to maintain the same behavior if this is what is intended...
Luca On Jun 1, 2012, at 1:58 PM, Wolfgang Bangerth <[email protected]> wrote: > > Javier, > >> //empty vector to consider all vertices >> std::vector< unsigned int> considered_vertices; >> GridTools::delete_duplicated_vertices (vertices, cells, subcelldata, >> considered_vertices, 1e-10); >> >> and recompile the library. >> >> This, apparently, eliminates the problem of Gmsh grids of duplicated >> surfaces that arises when the domain is formed by two (or more I guess) >> adjacent volumes. >> >> Would it be good to include GridTools::delete_duplicated_vertices per >> default in GridIn::read_mesh and perhaps others? > > Yes, that would be excellent. I don't know if this can happen with other > formats as well but if GMSH doesn't do this step itself it makes sense to > make this change in deal.II. Can you send me a diff file with the changes > between the old version and yours? > > Best & thanks > W. > > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.tamu.edu/~bangerth/ > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
