src/poly.c

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