sandbox/acastillo/output_fields/tests_outputs/test_output_vtu_box.c

    Testing output_vtu_box()

    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_box.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, but only inside a region defined by box 
  #if dimension == 2
    coord box[2] = {{-0.25, -0.25}, {0.25, 0.25}};
  #else
    coord box[2] = {{-0.25, -0.25, -0.25}, {0.25, 0.25, 0.25}};
  #endif
  output_vtu_box({f,p}, {u}, "domain", box);

    Verify the written file against the analytical fields, and additionally check that no cell outside the requested box was emitted. The report goes to stderr, i.e. to the log diffed against test_output_vtu_box.ref.

      if (pid() == 0)
    system ("python3 ../test_output_vtu.py --check "
            "--box=-0.25,-0.25,0.25,0.25 domain.vtu 1>&2");
}