src/regress.c

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