cgd_read_parallel_test.c

.cgd read test file.
- Adaptivity

.cgd read test file.

This works in serial currently but not in parallel parallel compiling that doesn’t work: qcc -g -O2 -fopenmp -Wall cgd_read_parallel_test.c -o cgd_read_parallel_test.exe -lm serial compiling that does qcc -g -O2 -Wall cgd_read_parallel_test.c -o cgd_read_parallel_test.exe -lm

#include "saint-venant.h"
#include "cgd_read.h"
#include "input.h"

We then define a few useful macros and constants.

#define MAXLEVEL 8
#define MINLEVEL 4
#define ETAE     1e-2 // error on free surface elevation (1 cm)

int main()
{
 size (1.);
  origin (0.,0.);
#if QUADTREE
  // 32^2 grid points to start with
  init_grid (1 << MINLEVEL);
#else // Cartesian
  // 1024^2 grid points
  init_grid (1 << MAXLEVEL);
#endif

  run();
}


int adapt() {
#if QUADTREE
 
  astats s = adapt_wavelet ({eta}, (double[]){ETAE},
			    MAXLEVEL, MINLEVEL);
  fprintf (stderr, "# refined %d cells, coarsened %d cells\n", s.nf, s.nc);
  return s.nf;
#else // Cartesian
  return 0;
#endif
}

event init (i = 0)
{

  foreach(){
    zb[]=-1.;
    h[] = max(0., - zb[]);}
  boundary ({h,zb});

/* Open the file. */
  char *fname = "/data/basilisk_test/test.cgd";
  FILE *fidin=fopen(fname,"r");
  fprintf(stderr,"deformation\n");
  deformation_cgd_read (x=0., y=0., fid = fidin, iterate = adapt );
  fclose ( fidin );

}

event logfile (i=0) {
  stats s = statsf (h);
  norm n = normf (u.x);
  if (i == 0)
    fprintf (stderr, "t i h.min h.max h.sum u.x.rms u.x.max dt\n");
  fprintf (stderr, "%g %d %g %g %g %g %g %g\n", t, i, s.min, s.max, s.sum, 
	   n.rms, n.max, dt);
}

Adaptivity

And finally we apply our adapt() function at every timestep.

event do_adapt (i++) adapt();