Hi Wolfgang, Thinking about it a bit more, if I generate a grid with :
subdivided_hyper_rectangle<http://www.dealii.org/6.3.1/doxygen/deal.II/classGridGenerator.html#ac5d882db5e1fa09d42501b038ef4dc03> (Triangulation<http://www.dealii.org/6.3.1/doxygen/deal.II/classTriangulation.html>< dim > &tria, const std::vector< unsigned<http://www.dealii.org/6.3.1/doxygen/deal.II/classunsigned.html> int <http://www.dealii.org/6.3.1/doxygen/deal.II/classint.html> > &repetitions, const Point<http://www.dealii.org/6.3.1/doxygen/deal.II/classPoint.html>< dim > &p1, const Point<http://www.dealii.org/6.3.1/doxygen/deal.II/classPoint.html>< dim > &p2, const bool<http://www.dealii.org/6.3.1/doxygen/deal.II/classbool.html> colorize=false) I could then go through every possible voxel. If the data indicates that it is a cell, I could then use the list of vertices in the equivalent cell in the generated grid above and store this list in the appropriate place in: int cell_vertices[][GeometryInfo<dim>::vertices_per_cell <http://www.dealii.org/6.3.1/doxygen/deal.II/structGeometryInfo.html>] Is there a way to obtain the list/numbering of vertices of a cell? Is this a reasonable way to do it or is there a cleaner and better way? Many thanks Ted
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
