deallog << std::endl;
}
- for (const auto face : tria.active_face_iterators())
+ for (const auto &face : tria.active_face_iterators())
{
if (face->at_boundary())
{
deallog << "Number of vertices: " << tria.get_vertices().size() << std::endl;
deallog << "Number of cells: " << tria.n_cells() << std::endl;
+ for (const auto &cell : tria.active_cell_iterators())
+ deallog << "cell " << cell->index()
+ << " type = " << cell->reference_cell().to_string()
+ << " volume = " << cell->measure() << std::endl;
+
GridOut go;
go.write_vtk(tria, deallog.get_file_stream());
}
DEAL::face center = 3.50000 7.00000 face boundary id = 0
DEAL::Number of vertices: 8
DEAL::Number of cells: 5
+DEAL::cell 0 type = Quad volume = 10.0000
+DEAL::cell 1 type = Quad volume = 10.0000
+DEAL::cell 2 type = Quad volume = 10.0000
+DEAL::cell 3 type = Quad volume = 10.0000
+DEAL::cell 4 type = Quad volume = 9.00000
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII
DEAL::face center = 3.50000 7.00000 3.50000 face boundary id = 0
DEAL::Number of vertices: 8
DEAL::Number of cells: 5
+DEAL::cell 0 type = Quad volume = 14.1421
+DEAL::cell 1 type = Quad volume = 14.1421
+DEAL::cell 2 type = Quad volume = 14.1421
+DEAL::cell 3 type = Quad volume = 14.1421
+DEAL::cell 4 type = Quad volume = 12.7279
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII
DEAL::face center = 8.00000 1.33333 3.66667 face boundary id = 0
DEAL::Number of vertices: 13
DEAL::Number of cells: 4
+DEAL::cell 0 type = Hex volume = 27.0000
+DEAL::cell 1 type = Pyramid volume = 6.00000
+DEAL::cell 2 type = Tet volume = 2.00000
+DEAL::cell 3 type = Wedge volume = -15.0000
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII
DEAL::face center = 8.00000 1.33333 3.66667 face boundary id = 0
DEAL::Number of vertices: 29
DEAL::Number of cells: 4
+DEAL::cell 0 type = Hex volume = 27.0000
+DEAL::cell 1 type = Pyramid volume = 6.00000
+DEAL::cell 2 type = Tet volume = 2.00000
+DEAL::cell 3 type = Wedge volume = -15.0000
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII
DEAL::face center = 8.00000 1.33333 3.66667 face boundary id = 6
DEAL::Number of vertices: 13
DEAL::Number of cells: 4
+DEAL::cell 0 type = Hex volume = 27.0000
+DEAL::cell 1 type = Pyramid volume = 6.00000
+DEAL::cell 2 type = Tet volume = 2.00000
+DEAL::cell 3 type = Wedge volume = -15.0000
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII
DEAL::face center = 8.00000 1.33333 3.66667 face manifold id = 6
DEAL::Number of vertices: 13
DEAL::Number of cells: 4
+DEAL::cell 0 type = Hex volume = 27.0000
+DEAL::cell 1 type = Pyramid volume = 6.00000
+DEAL::cell 2 type = Tet volume = 2.00000
+DEAL::cell 3 type = Wedge volume = -15.0000
# vtk DataFile Version 3.0
Triangulation generated with deal.II
ASCII