src/regress.c

   1 #include "a/regress.h"
   2 #include "a/math.h"
   3
   4 void a_regress_odm_(a_size n, a_float const *p, a_float *odm, a_float p_mean)
   5 {
   6     for (; n; --n) { *odm++ = *p++ - p_mean; }
   7 }
   8
   9 a_float a_regress_odm(a_size n, a_float const *p, a_float *odm)
  10 {
  11     a_float m = a_float_mean(p, n);
  12     a_regress_odm_(n, p, odm, m);
  13     return m;
  14 }