event perfs (i += 1) {
  static FILE * fp = fopen ("perfs", "w");
  if (i == 0)
    fprintf (fp,
	     "t dt mgp.i mgp.nrelax grid->tn perf.t perf.speed npe\n");
  fprintf (fp, "%g %g %d %d %ld %g %g %d\n", 
	   t, dt,
#if NH
	   mgp.i, mgp.nrelax,
#elif IMPLICIT_H
	   mgH.i, mgH.nrelax,
#else
	   0, 0,
#endif
	   grid->tn*nl, perf.t, perf.speed*nl, npe());
  fflush (fp);
}

event perf_plot (i = 10) {
  if (getenv ("DISPLAY"))
    popen ("gnuplot -e 'set term x11 noraise title perfs' "
	   "$BASILISK/layered/perfs.plot 2> /dev/null "
	   "& read dummy; kill $!", "w");
}