]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Small changes.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 9 Apr 1998 09:24:32 +0000 (09:24 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 9 Apr 1998 09:24:32 +0000 (09:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@166 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/Attic/examples/convergence/convergence.cc
tests/big-tests/convergence/convergence.cc

index 5b17341eb4a3d21e06b7b9c0a385a8e8e73c7d30..cff66569ed6386f9424f61b3364b98297f221aaf 100644 (file)
@@ -292,28 +292,33 @@ void PoissonProblem<dim>::run (const unsigned int level) {
   cout << h1_error_per_cell.l2_norm() << endl;
   h1_error.push_back (h1_error_per_cell.l2_norm());
 
-  dVector l1_error_per_dof, l2_error_per_dof, linfty_error_per_dof;
-  dVector h1_seminorm_error_per_dof, h1_error_per_dof;
-  dof->distribute_cell_to_dof_vector (l1_error_per_cell, l1_error_per_dof);
-  dof->distribute_cell_to_dof_vector (l2_error_per_cell, l2_error_per_dof);
-  dof->distribute_cell_to_dof_vector (linfty_error_per_cell, linfty_error_per_dof);
-  dof->distribute_cell_to_dof_vector (h1_seminorm_error_per_cell, h1_seminorm_error_per_dof);
-  dof->distribute_cell_to_dof_vector (h1_error_per_cell, h1_error_per_dof);
-
-  string filename = "gnuplot.";
-  filename += ('0'+level);
-  cout << "    Writing error plots to <" << filename << ">..." << endl;
-
-  DataOut<dim> out;
-  ofstream gnuplot(filename.c_str());
-  fill_data (out);
-  out.add_data_vector (l1_error_per_dof, "L1-Error");
-  out.add_data_vector (l2_error_per_dof, "L2-Error");
-  out.add_data_vector (linfty_error_per_dof, "Linfty-Error");
-  out.add_data_vector (h1_seminorm_error_per_dof, "H1-seminorm-Error");
-  out.add_data_vector (h1_error_per_dof, "H1-Error");
-  out.write_gnuplot (gnuplot);
-  gnuplot.close ();
+  if (level<=5) 
+    {
+      dVector l1_error_per_dof, l2_error_per_dof, linfty_error_per_dof;
+      dVector h1_seminorm_error_per_dof, h1_error_per_dof;
+      dof->distribute_cell_to_dof_vector (l1_error_per_cell, l1_error_per_dof);
+      dof->distribute_cell_to_dof_vector (l2_error_per_cell, l2_error_per_dof);
+      dof->distribute_cell_to_dof_vector (linfty_error_per_cell, linfty_error_per_dof);
+      dof->distribute_cell_to_dof_vector (h1_seminorm_error_per_cell, h1_seminorm_error_per_dof);
+      dof->distribute_cell_to_dof_vector (h1_error_per_cell, h1_error_per_dof);
+
+      string filename = "gnuplot.";
+      filename += ('0'+level);
+      cout << "    Writing error plots to <" << filename << ">..." << endl;
+      
+      DataOut<dim> out;
+      ofstream gnuplot(filename.c_str());
+      fill_data (out);
+      out.add_data_vector (l1_error_per_dof, "L1-Error");
+      out.add_data_vector (l2_error_per_dof, "L2-Error");
+      out.add_data_vector (linfty_error_per_dof, "Linfty-Error");
+      out.add_data_vector (h1_seminorm_error_per_dof, "H1-seminorm-Error");
+      out.add_data_vector (h1_error_per_dof, "H1-Error");
+      out.write_gnuplot (gnuplot);
+      gnuplot.close ();
+    }
+  else
+    cout << "    Not writing error as grid." << endl;
   
   cout << endl;
 };
@@ -367,7 +372,7 @@ void PoissonProblem<dim>::print_history () const {
 int main () {
   PoissonProblem<2> problem;
 
-  for (unsigned int level=1; level<8; ++level)
+  for (unsigned int level=1; level<10; ++level)
     problem.run (level);
 
   cout << endl << "Printing convergence history to <gnuplot.history>..." << endl;
index 5b17341eb4a3d21e06b7b9c0a385a8e8e73c7d30..cff66569ed6386f9424f61b3364b98297f221aaf 100644 (file)
@@ -292,28 +292,33 @@ void PoissonProblem<dim>::run (const unsigned int level) {
   cout << h1_error_per_cell.l2_norm() << endl;
   h1_error.push_back (h1_error_per_cell.l2_norm());
 
-  dVector l1_error_per_dof, l2_error_per_dof, linfty_error_per_dof;
-  dVector h1_seminorm_error_per_dof, h1_error_per_dof;
-  dof->distribute_cell_to_dof_vector (l1_error_per_cell, l1_error_per_dof);
-  dof->distribute_cell_to_dof_vector (l2_error_per_cell, l2_error_per_dof);
-  dof->distribute_cell_to_dof_vector (linfty_error_per_cell, linfty_error_per_dof);
-  dof->distribute_cell_to_dof_vector (h1_seminorm_error_per_cell, h1_seminorm_error_per_dof);
-  dof->distribute_cell_to_dof_vector (h1_error_per_cell, h1_error_per_dof);
-
-  string filename = "gnuplot.";
-  filename += ('0'+level);
-  cout << "    Writing error plots to <" << filename << ">..." << endl;
-
-  DataOut<dim> out;
-  ofstream gnuplot(filename.c_str());
-  fill_data (out);
-  out.add_data_vector (l1_error_per_dof, "L1-Error");
-  out.add_data_vector (l2_error_per_dof, "L2-Error");
-  out.add_data_vector (linfty_error_per_dof, "Linfty-Error");
-  out.add_data_vector (h1_seminorm_error_per_dof, "H1-seminorm-Error");
-  out.add_data_vector (h1_error_per_dof, "H1-Error");
-  out.write_gnuplot (gnuplot);
-  gnuplot.close ();
+  if (level<=5) 
+    {
+      dVector l1_error_per_dof, l2_error_per_dof, linfty_error_per_dof;
+      dVector h1_seminorm_error_per_dof, h1_error_per_dof;
+      dof->distribute_cell_to_dof_vector (l1_error_per_cell, l1_error_per_dof);
+      dof->distribute_cell_to_dof_vector (l2_error_per_cell, l2_error_per_dof);
+      dof->distribute_cell_to_dof_vector (linfty_error_per_cell, linfty_error_per_dof);
+      dof->distribute_cell_to_dof_vector (h1_seminorm_error_per_cell, h1_seminorm_error_per_dof);
+      dof->distribute_cell_to_dof_vector (h1_error_per_cell, h1_error_per_dof);
+
+      string filename = "gnuplot.";
+      filename += ('0'+level);
+      cout << "    Writing error plots to <" << filename << ">..." << endl;
+      
+      DataOut<dim> out;
+      ofstream gnuplot(filename.c_str());
+      fill_data (out);
+      out.add_data_vector (l1_error_per_dof, "L1-Error");
+      out.add_data_vector (l2_error_per_dof, "L2-Error");
+      out.add_data_vector (linfty_error_per_dof, "Linfty-Error");
+      out.add_data_vector (h1_seminorm_error_per_dof, "H1-seminorm-Error");
+      out.add_data_vector (h1_error_per_dof, "H1-Error");
+      out.write_gnuplot (gnuplot);
+      gnuplot.close ();
+    }
+  else
+    cout << "    Not writing error as grid." << endl;
   
   cout << endl;
 };
@@ -367,7 +372,7 @@ void PoissonProblem<dim>::print_history () const {
 int main () {
   PoissonProblem<2> problem;
 
-  for (unsigned int level=1; level<8; ++level)
+  for (unsigned int level=1; level<10; ++level)
     problem.run (level);
 
   cout << endl << "Printing convergence history to <gnuplot.history>..." << endl;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.