test/poly.h

   1 #define MAIN(x) poly##x
   2 #include "test.h"
   3 #include "a/poly.h"
   4
   5 static void print_(a_float const *a, a_size n)
   6 {
   7     a_size i;
   8     (void)putchar('{');
   9     A_FORENUM(a_size, i, n)
  10     {
  11         printf("%" A_FLOAT_PRI "g", a[i]);
  12         if (i < n - 1) { (void)putchar(','); }
  13     }
  14     (void)putchar('}');
  15 }
  16
  17 int main(int argc, char *argv[]) /* NOLINT(misc-definitions-in-headers) */
  18 {
  19     a_float x = 2;
  20     a_float x4[] = {A_FLOAT_C(1.0), A_FLOAT_C(2.0), A_FLOAT_C(3.0), A_FLOAT_C(4.0)};
  21     print_(x4, 4);
  22     printf("=%" A_FLOAT_PRI "g\n", a_poly_eval(x4, 4, x));
  23     a_poly_swap(x4, 4);
  24     print_(x4, 4);
  25     printf("=%" A_FLOAT_PRI "g\n", a_poly_eval(x4, 4, x));
  26     print_(x4, 4);
  27     a_poly_swap(x4, 4);
  28     printf("=%" A_FLOAT_PRI "g\n", a_poly_evar(x4, 4, x));
  29     print_(x4, 3);
  30     printf("=%" A_FLOAT_PRI "g\n", a_poly_eval(x4, 3, x));
  31     a_poly_swap(x4, 3);
  32     print_(x4, 3);
  33     printf("=%" A_FLOAT_PRI "g\n", a_poly_eval(x4, 3, x));
  34     print_(x4, 3);
  35     printf("=%" A_FLOAT_PRI "g\n", a_poly_evar(x4, 3, x));
  36     a_poly_swap(x4, 3);
  37     (void)argc;
  38     (void)argv;
  39     return 0;
  40 }