sandbox/acastillo/output_fields/tests_outputs/test_output_vtu.c

    Testing output_vtu() and output_vtu_slice()

    Tests VTU output functions

    import sys, os
if os.path.exists('test_output_vtu.py'):
  sys.path.append('.')
else:
  sys.path.append('../')
import test_output_vtu
data = test_output_vtu.print_vtu_info('domain.vtu')
    (script) 
    #include "acastillo/output_fields/vtu/output_vtu.h"
#define MAXLEVEL 4
#define ASPECTRATIO 8
#define r2 (sq(x) + sq(y))

int main(){

  #if !TREE
    #if dimension == 2
      dimensions (nx = 1, ny = ASPECTRATIO);
    #else
      dimensions (nx = 1, ny = ASPECTRATIO, nz = 1);
    #endif
  #endif

  L0 = 1.0;
  X0 = Y0 = Z0 = -L0 / 2;
  N = 1 << (MAXLEVEL-1);
  init_grid(N);

  #if TREE
    double outer_radius = 0.25;
    double inner_radius = 0.1 ;
    refine(((r2 < sq(outer_radius)) && (r2 > sq(inner_radius))) && level < MAXLEVEL);
  #endif

  // Create test fields 
  scalar f[], p[];
  vector u[];
  foreach(){    
    f[] = point.level;
    p[] = pid();
    u.x[] = x;
    u.y[] = y;
    #if dimension == 3
      u.z[] = z;
    #endif
  }

  // And write a vtu file
  output_vtu({f,p}, {u}, "domain");

#if dimension > 2
  output_slice_vtu({f,p}, {u}, "slice_x", (coord){1,0,0}, 0);
  output_slice_vtu({f,p}, {u}, "slice_y", (coord){0,1,0}, 0);
  output_slice_vtu({f,p}, {u}, "slice_z", (coord){0,0,1}, 0);
#endif

}