From 961ceddbe28a78c301292d926a31c422afc9e6dc Mon Sep 17 00:00:00 2001 From: hartmann Date: Mon, 20 Sep 2004 09:18:58 +0000 Subject: [PATCH] Simplification of code in GridOut<2>::write_gnuplot. No change of output. git-svn-id: https://svn.dealii.org/trunk@9645 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/grid/grid_out.cc | 22 ++++++++-------------- 1 file changed, 8 insertions(+), 14 deletions(-) diff --git a/deal.II/deal.II/source/grid/grid_out.cc b/deal.II/deal.II/source/grid/grid_out.cc index 29edc54814..0ff5a2e586 100644 --- a/deal.II/deal.II/source/grid/grid_out.cc +++ b/deal.II/deal.II/source/grid/grid_out.cc @@ -613,8 +613,9 @@ void GridOut::write_gnuplot (const Triangulation &tria, { AssertThrow (out, ExcIO()); - const unsigned int n_points= + const unsigned int n_additional_points= gnuplot_flags.n_boundary_face_points; + const unsigned int n_points=GeometryInfo::vertices_per_face+n_additional_points; typename Triangulation::active_cell_iterator cell=tria.begin_active(); const typename Triangulation::active_cell_iterator endc=tria.end(); @@ -628,8 +629,10 @@ void GridOut::write_gnuplot (const Triangulation &tria, if (mapping!=0) { std::vector > boundary_points(n_points); - for (unsigned int i=0; i quadrature(boundary_points); q_projector = new Quadrature (QProjector::project_to_all_faces(quadrature)); @@ -695,11 +698,6 @@ void GridOut::write_gnuplot (const Triangulation &tria, face = cell->face(face_no); if (face->at_boundary()) { - out << face->vertex(0) - << ' ' << cell->level() - << ' ' << static_cast(cell->material_id()) - << std::endl; - // compute // offset of // quadrature @@ -735,12 +733,8 @@ void GridOut::write_gnuplot (const Triangulation &tria, << ' ' << cell->level() << ' ' << static_cast(cell->material_id()) << std::endl; - - out << face->vertex(1) - << ' ' << cell->level() - << ' ' << static_cast(cell->material_id()) - << std::endl - << std::endl + + out << std::endl << std::endl; } else -- 2.39.5