| #ifndef __WDONG_LOCAL__ | |
| #define __WDONG_LOCAL__ | |
| /* y = a + b x */ | |
| static inline void least_squares (double *a, double *b, int n, double sxx, double sxy, double sx, double sy) | |
| { | |
| *b = (n * sxy - sx*sy) / (n * sxx - sx*sx); | |
| *a = (sy - *b * sx) / n; | |
| } | |
| int localdim (double a, double b, double *alpha, double *beta, double W, int M, int L, int K); | |
| #endif |