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
