/* lsfit.h  ultimately for 2, 3, and 4 dimensions          */
/*          Compute a multiple dimensional polynomial      */
/*          that is a least square fit of some data.       */
/*          Evaluate the polynomial at desired points.     */
/* use:   fit_init(n,m); n is power of polynomial          */
/*                       m is number of dimensions         */
/*        fit2_bpoly(nx, ny, xg, yg, u) build coefficients */
/*                   nx number of x coordinates i=0,nx-1   */
/*                   ny number of y coordinates j=0,ny-1   */
/*                   xg array of x coordinate values       */
/*                   yg array of y coordinate values       */
/*                   u[i*ny+j] value at point i,j          */
void fit_init(int n, int m); /* only one at a time */
void fit_free();             /* free allocated ram */
void fit2_data(double x, double y, double u, int last);
void fit3_data(double x, double y, double z, double u, int last);
void fit4_data(double x, double y, double z, double t, double u, int last);
void fit5_data(double x, double y, double z, double t, double u, double w,
               int last);
void fit6_data(double x, double y, double z, double t, double u, double v,
               double w, int last);
void fit2_bpoly(int nx, int ny, double x[], double y[], double u[]);
double eval2_poly(double x, double y);
double eval3_poly(double x, double y, double z);
double eval4_poly(double x, double y, double z, double t);
double eval5_poly(double x, double y, double z, double t, double u);
double eval6_poly(double x, double y, double z, double t, double u, double v);