From 961ceddbe28a78c301292d926a31c422afc9e6dc Mon Sep 17 00:00:00 2001
From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
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<dim> &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<dim>::vertices_per_face+n_additional_points;
 
   typename Triangulation<dim>::active_cell_iterator        cell=tria.begin_active();
   const typename Triangulation<dim>::active_cell_iterator  endc=tria.end();
@@ -628,8 +629,10 @@ void GridOut::write_gnuplot (const Triangulation<dim> &tria,
   if (mapping!=0)
     {
       std::vector<Point<dim-1> > boundary_points(n_points);
-      for (unsigned int i=0; i<n_points; ++i)
-	boundary_points[i](0)= 1.*(i+1)/(n_points+1);
+      boundary_points[0][0]=0;
+      boundary_points[n_points-1][0]=1;      
+      for (unsigned int i=1; i<n_points-1; ++i)
+        boundary_points[i](0)= 1.*i/(n_points-1);
 
       Quadrature<dim-1> quadrature(boundary_points);
       q_projector = new Quadrature<dim> (QProjector<dim>::project_to_all_faces(quadrature));
@@ -695,11 +698,6 @@ void GridOut::write_gnuplot (const Triangulation<dim> &tria,
 		      face = cell->face(face_no);
 		    if (face->at_boundary())
 		      {
-			out << face->vertex(0)
-			    << ' ' << cell->level()
-			    << ' ' << static_cast<unsigned int>(cell->material_id())
-			    << std::endl;
-
                                                          // compute
                                                          // offset of
                                                          // quadrature
@@ -735,12 +733,8 @@ void GridOut::write_gnuplot (const Triangulation<dim> &tria,
 			      << ' ' << cell->level()
 			      << ' ' << static_cast<unsigned int>(cell->material_id())
 			      << std::endl;
-			
-			out << face->vertex(1)
-			    << ' ' << cell->level()
-			    << ' ' << static_cast<unsigned int>(cell->material_id())
-			    << std::endl
-			    << std::endl
+
+			out << std::endl
 			    << std::endl;
 		      }
 		    else
-- 
2.39.5