src/test/hf1.c

    1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

    #include "fractions.h"
#include "curvature.h"
#include "utils.h"

int main()
{
  origin (-0.5, -0.5);
  init_grid (16);
  scalar c[];
  fraction (c, min(sq(x + 0.1) + sq(y + 0.22) - sq(0.18),
		   sq(x - 0.1) + sq(y - 0.22) - sq(0.18)));
  output_cells (stdout);
  output_facets (c, stdout);
  vector h[];
  heights (c, h);
  scalar kappa[];
  curvature (c, kappa);
  foreach() {
    if (h.x[] != nodata)
      fprintf (qerr, "%g %g %g %g hx\n", x, y, height(h.x[]), Delta);
    if (h.y[] != nodata)
      fprintf (qerr, "%g %g %g %g hy\n", x, y, height(h.y[]), Delta);
    if (kappa[] != nodata)
      fprintf (qerr, "%g %g %g kappa\n", x, y, kappa[]);
  }
}