#include "embed.h"
#include "run.h"

#define GAUSS (exp(-(sq(x))) - y)

scalar s[];
int main() {
  L0 = 100.;
  X0 = -L0/2.;
  init_grid (1 << 8);
  run();
}

event init (t = 0) {
  s.refine = refine_embed_linear;
  refine (fabs(GAUSS) < 0.1 && level < 9); //This goes well
  foreach()
    printf ("%g %g\n", x, y);
  fflush (stdout);

  FILE * fp2 = fopen ("iterated", "w");
  refine (fabs(GAUSS) < 0.1 && level < 10 &&
	  fprintf(fp2, "%g %g\n", x, y) && !fflush(fp2));//This goes wrong
}