Here is the output of $ diff grid_in.cc grid_in_new.cc > grid_in_diff.txt
859c859,862 < GridTools::delete_unused_vertices (vertices, cells, subcelldata); --- > > std::vector< unsigned int > considered_vertices; > GridTools::delete_duplicated_vertices (vertices, cells, subcelldata, considered_vertices, 1e-10); > // GridTools::delete_unused_vertices (vertices, cells, subcelldata); and I'm including the txt file as well. Regards Javier On Fri, 2012-06-01 at 14:02 +0200, Luca Heltai wrote: > 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
859c859,862 < GridTools::delete_unused_vertices (vertices, cells, subcelldata); --- > > std::vector< unsigned int > considered_vertices; > GridTools::delete_duplicated_vertices (vertices, cells, subcelldata, > considered_vertices, 1e-10); > // GridTools::delete_unused_vertices (vertices, cells, subcelldata);
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
