1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
| #include "layered/check_eta.h"
event init (i = 0)
{
if (res_eta.i < 0)
res_eta = new scalar;
}
event gnuplots (i += 10)
{
static FILE * fp = popen ("gnuplot 2> /dev/null", "w");
if (i == 0)
fprintf (fp, "set term x11\n");
FILE * fp1 = fopen ("gnuplot", "w");
foreach()
fprintf (fp1, "%g %g %g %g\n", x, eta[], deta[], res_eta[]);
fclose (fp1);
fprintf (fp,
"p 'gnuplot' u 1:3 w l t 'etap - eta',"
" '' u 1:4 w p t 'res'\n");
fflush (fp);
}
|