src/test/tag.c

    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    
    #include "utils.h"
    #include "tag.h"
    
    int main()
    {
      size (1.[0]);
      init_grid (64);
      X0 = Y0 = -0.5;
      scalar t[];
      foreach()
        t[] = sin(4*pi*x)*cos(4.*pi*y) > 0.5;
      int n = tag (t);
      foreach() {
        if (t[] == 0)
          t[] = nodata;
        else
          fprintf (qerr, "%g %g %g\n", x, y, t[]);
      }
      output_ppm (t, file = "t.ppm", min = 1, max = n, n = 256);
    }