# To execute additional action apart from running the program, fill
# in this list:
-additional-run-action =
+additional-run-action = gnuplot make_ps
# To specify which files are to be deleted by "make clean" (apart from
# the usual ones: object files, executables, backups, etc), fill in the
template <int dim>
PoissonProblem<dim>::PoissonProblem () :
tria(0), dof(0), rhs(0),
- boundary_values(0), boundary(Point<dim>(),1) {};
+ boundary_values(0), boundary() {};
cout << " Making grid... ";
tria->set_boundary (&boundary);
- tria->create_hyper_ball (Point<dim>(),+1);
+ tria->create_hyper_ball ();
tria->refine_global (start_level);
cout << tria->n_active_cells() << " active cells." << endl;
--- /dev/null
+set xlabel "Number of degrees of freedom"
+set ylabel "Error"
+set data style linespoints
+set logscale xy
+
+set term postscript eps
+set output "error-estimation.eps"
+
+plot "gnuplot.history" using 1:2 title "L1 error","gnuplot.history" using 1:3 title "L2 error","gnuplot.history" using 1:4 title "Linfty error","gnuplot.history" using 1:5 title "H1 seminorm error","gnuplot.history" using 1:6 title "H1 error","gnuplot.history" using 1:7 title "Estimated H1 error"
+
# To execute additional action apart from running the program, fill
# in this list:
-additional-run-action =
+additional-run-action = gnuplot make_ps
# To specify which files are to be deleted by "make clean" (apart from
# the usual ones: object files, executables, backups, etc), fill in the
template <int dim>
PoissonProblem<dim>::PoissonProblem () :
tria(0), dof(0), rhs(0),
- boundary_values(0), boundary(Point<dim>(),1) {};
+ boundary_values(0), boundary() {};
cout << " Making grid... ";
tria->set_boundary (&boundary);
- tria->create_hyper_ball (Point<dim>(),+1);
+ tria->create_hyper_ball ();
tria->refine_global (start_level);
cout << tria->n_active_cells() << " active cells." << endl;
--- /dev/null
+set xlabel "Number of degrees of freedom"
+set ylabel "Error"
+set data style linespoints
+set logscale xy
+
+set term postscript eps
+set output "error-estimation.eps"
+
+plot "gnuplot.history" using 1:2 title "L1 error","gnuplot.history" using 1:3 title "L2 error","gnuplot.history" using 1:4 title "Linfty error","gnuplot.history" using 1:5 title "H1 seminorm error","gnuplot.history" using 1:6 title "H1 error","gnuplot.history" using 1:7 title "Estimated H1 error"
+