From 3569f2744c31c9103293806310dfaa4a03fb3b69 Mon Sep 17 00:00:00 2001
From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 27 Nov 2001 17:38:09 +0000
Subject: [PATCH] Include all data of the Shapes*.output files into the
 mapping.output file.

git-svn-id: https://svn.dealii.org/trunk@5281 0785d39b-7218-0410-832d-ea1e28bc413d
---
 tests/fe/shapes.cc | 44 +++++++++++++++++++++-----------------------
 1 file changed, 21 insertions(+), 23 deletions(-)

diff --git a/tests/fe/shapes.cc b/tests/fe/shapes.cc
index b7da61371f..3956cfbf7e 100644
--- a/tests/fe/shapes.cc
+++ b/tests/fe/shapes.cc
@@ -42,10 +42,8 @@ plot_shape_functions(Mapping<dim>& mapping,
   FEValues<dim> fe(mapping, finel, q, UpdateFlags(update_values));
   fe.reinit(c);
   
-  sprintf(fname, "Shapes%dd-%s.output", dim, name);
-  std::ofstream gnuplot(fname);
-  gnuplot.setf(std::ios::fixed);
-  gnuplot.precision (PRECISION);
+  sprintf(fname, "Shapes%dd-%s", dim, name);
+  deallog.push(fname);
 
   unsigned int k=0;
   for (unsigned int mz=0;mz<=((dim>2) ? div : 0) ;++mz)
@@ -54,19 +52,20 @@ plot_shape_functions(Mapping<dim>& mapping,
 	{
 	  for (unsigned int mx=0;mx<=div;++mx)
 	    {
-	      gnuplot << q.point(k);
+	      deallog << q.point(k);
 	  
 	      for (unsigned int i=0;i<finel.dofs_per_cell;++i)
 		{
-		  gnuplot << " " << fe.shape_value(i,k) + 1.;
+		  deallog << " " << fe.shape_value(i,k) + 1.;
 		}
-	      gnuplot << std::endl;
+	      deallog << std::endl;
 	      k++;
 	    }
-	  gnuplot << std::endl;
+	  deallog << std::endl;
 	}
-      gnuplot << std::endl;
+      deallog << std::endl;
     }
+  deallog.pop();
 }
 
 
@@ -98,10 +97,8 @@ plot_face_shape_functions(Mapping<dim>& mapping,
   FESubfaceValues<dim> sub(mapping, finel, q, UpdateFlags(update_values
     | update_q_points));
 
-  sprintf(fname, "ShapesFace%dd-%s.output", dim, name);
-  std::ofstream gnuplot(fname);
-  gnuplot.setf(std::ios::fixed);
-  gnuplot.precision (PRECISION);
+  sprintf(fname, "ShapesFace%dd-%s", dim, name);
+  deallog.push(fname);
   
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     {
@@ -114,18 +111,18 @@ plot_face_shape_functions(Mapping<dim>& mapping,
 	    {
 	      for (unsigned int mx=0;mx<=div;++mx)
 		{
-		  gnuplot << fe.quadrature_point(k);
+		  deallog << fe.quadrature_point(k);
 		  
 		  for (unsigned int i=0;i<finel.dofs_per_cell;++i)
 		    {
-		      gnuplot << " " << fe.shape_value(i,k) + 1.;
+		      deallog << " " << fe.shape_value(i,k) + 1.;
 		    }
-		  gnuplot << std::endl;
+		  deallog << std::endl;
 		  k++;
 		}
-	      gnuplot << std::endl;
+	      deallog << std::endl;
 	    }
-	  gnuplot << std::endl;
+	  deallog << std::endl;
 	} else {
 	  for (unsigned int s=0;s<GeometryInfo<dim>::subfaces_per_face; ++s)
 	    {
@@ -136,21 +133,22 @@ plot_face_shape_functions(Mapping<dim>& mapping,
 		{
 		  for (unsigned int mx=0;mx<=div;++mx)
 		    {
-		      gnuplot << sub.quadrature_point(k);
+		      deallog << sub.quadrature_point(k);
 		      
 		      for (unsigned int i=0;i<finel.dofs_per_cell;++i)
 			{
-			  gnuplot << " " << sub.shape_value(i,k) + 1.;
+			  deallog << " " << sub.shape_value(i,k) + 1.;
 			}
-		      gnuplot << std::endl;
+		      deallog << std::endl;
 		      k++;
 		    }
-		  gnuplot << std::endl;
+		  deallog << std::endl;
 		}
-	      gnuplot << std::endl;	      
+	      deallog << std::endl;	      
 	    }
 	}
     }
+  deallog.pop();
 }
 
 
-- 
2.39.5