sandbox/bugs/my_report_bug.c

#include "grid/multigrid3D.h"
#include "navier-stokes/centered.h"
#include "lambda2.h"
#include "navier-stokes/perfs.h"
#include "maxruntime.h"
const double MU = 0.01;
face vector av[];
double RUNID;            // For restarting purpose                                 
int maxlevel = 7;

int main (int argc, char * argv[])
{
  maxruntime (&argc, argv);
  if (argc > 1)
    RUNID = atoi(argv[1]);

  L0 = 2.*pi [1];
  dimensions(nx=4,ny=1,nz=1);
  foreach_dimension()
    periodic (right);

  mu[] = {MU,MU,MU};
  a = av;
  N = 2048;
  run();
}

event init (i = 0) {
  double u0 = 1., k = 1.;
  if (!restore (file = "restart"))
    foreach() {
      u.x[] = u0*(cos(k*y) + sin(k*z));
      u.y[] = u0*(sin(k*x) + cos(k*z));
      u.z[] = u0*(cos(k*x) + sin(k*y));
    }
}

event end(i = 5 + (RUNID * 10)){
  char name[80];
  sprintf (name, "restart");
  dump (file = name);
  return 1;
}