if (y<x)
if (y<1-x)
- middle(1) = 0.03*sin(6*3.141592*middle(0));
+ middle(1) = 0.04*sin(6*3.141592*middle(0));
else
- middle(0) = 1+0.03*sin(6*3.141592*middle(1));
+ middle(0) = 1+0.04*sin(6*3.141592*middle(1));
else
if (y<1-x)
- middle(0) = 0.03*sin(6*3.141592*middle(1));
+ middle(0) = 0.04*sin(6*3.141592*middle(1));
else
- middle(1) = 1+0.03*sin(6*3.141592*middle(0));
+ middle(1) = 1+0.04*sin(6*3.141592*middle(0));
return middle;
};
template <int dim>
void PoissonProblem<dim>::run (ParameterHandler &prm) {
- cout << "Test case = " << prm.get ("Test run") << endl
+ cout << "Test case = " << prm.get ("Test run")
<< endl;
cout << " Making grid..." << endl;
cout << " Solving..." << endl;
solve ();
- cout << " Writing to file <" << prm.get("Output file") << ">..." << endl;
+ cout << " Writing to file <" << prm.get("Output file") << ">..."
+ << endl
+ << endl;
+
DataOut<dim> out;
-
ofstream gnuplot(prm.get("Output file"));
fill_data (out);
out.write_gnuplot (gnuplot);
--- /dev/null
+set data style lines
+set noxtics
+set noytics
+set noztics
+set nokey
+set para
+set term postscript eps
+
+!echo " Making <fig.zoom_in.eps>"
+set output "fig.zoom_in.eps"
+splot "gnuplot.zoom_in"
+
+!echo " Making <fig.ball.eps>"
+set output "fig.ball.eps"
+splot "gnuplot.ball"
+
+!echo " Making <fig.curved_line.eps>"
+set output "fig.curved_line.eps"
+splot "gnuplot.curved_line"
+
+!echo " Making <fig.random.eps>"
+set output "fig.random.eps"
+splot "gnuplot.random"
+
+
if (y<x)
if (y<1-x)
- middle(1) = 0.03*sin(6*3.141592*middle(0));
+ middle(1) = 0.04*sin(6*3.141592*middle(0));
else
- middle(0) = 1+0.03*sin(6*3.141592*middle(1));
+ middle(0) = 1+0.04*sin(6*3.141592*middle(1));
else
if (y<1-x)
- middle(0) = 0.03*sin(6*3.141592*middle(1));
+ middle(0) = 0.04*sin(6*3.141592*middle(1));
else
- middle(1) = 1+0.03*sin(6*3.141592*middle(0));
+ middle(1) = 1+0.04*sin(6*3.141592*middle(0));
return middle;
};
template <int dim>
void PoissonProblem<dim>::run (ParameterHandler &prm) {
- cout << "Test case = " << prm.get ("Test run") << endl
+ cout << "Test case = " << prm.get ("Test run")
<< endl;
cout << " Making grid..." << endl;
cout << " Solving..." << endl;
solve ();
- cout << " Writing to file <" << prm.get("Output file") << ">..." << endl;
+ cout << " Writing to file <" << prm.get("Output file") << ">..."
+ << endl
+ << endl;
+
DataOut<dim> out;
-
ofstream gnuplot(prm.get("Output file"));
fill_data (out);
out.write_gnuplot (gnuplot);