tests/kidneycontour.c

   1 #include <qpms/beyn.h>
   2 #include <stdio.h>
   3
   4 #define CPAIR(x) creal(x), cimag(x)
   5
   6 int main(int argc, char **argv)
   7 {
   8   double rRe = 2e3, rIm = 1.5e3, rounding = 0.2;
   9   complex double centre = 1e3 * I;
  10   size_t n = 100;
  11
  12   beyn_contour_t *c = beyn_contour_kidney(centre, rRe, rIm, rounding, n, BEYN_CONTOUR_HALFELLIPSE_IM_PLUS);
  13
  14   for(size_t i = 0; i < n; ++i)
  15     printf("%g\t%g\t%g\t%g\n", CPAIR(c->z_dz[i][0]), CPAIR(c->z_dz[i][1]));
  16
  17   free(c);
  18   return 0;
  19 }