src/poly.c

   1 #define LIBA_POLY_C
   2 #include "a/poly.h"
   3
   4 a_float *a_poly_swap(a_float *a, a_size n)
   5 {
   6     a_float *b = a, *c = a + n;
   7     for (a_float x; b < --c; ++b)
   8     {
   9         x = *b;
  10         *b = *c;
  11         *c = x;
  12     }
  13     return a;
  14 }
  15
  16 a_float a_poly_eval_(a_float const *a, a_float const *b, a_float x)
  17 {
  18     a_float y;
  19     for (y = *--b; b > a;)
  20     {
  21         y = y * x + *--b;
  22     }
  23     return y;
  24 }
  25
  26 a_float a_poly_evar_(a_float const *a, a_float const *b, a_float x)
  27 {
  28     a_float y;
  29     for (y = *a; ++a < b;)
  30     {
  31         y = y * x + *a;
  32     }
  33     return y;
  34 }