Mercurial > repos > vipints > rdiff
diff rDiff/src/locfit/Source/lfev.c @ 0:0f80a5141704
version 0.3 uploaded
author | vipints |
---|---|
date | Thu, 14 Feb 2013 23:38:36 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rDiff/src/locfit/Source/lfev.c Thu Feb 14 23:38:36 2013 -0500 @@ -0,0 +1,44 @@ +#include "local.h" + +int lf_error; + + lfit lf; + design des; + +void calllf(x,y,n) +double *x, *y; +int n; +{ lfdata *lfd; + fitpt *fp; + int i; + + lfit_init(&lf); + + lfd = &(lf.lfd); + lfd->x[0] = x; + lfd->y = y; + lfd->n = n; + lfd->d = 1; + + startlf(&des,&lf,procv,1); +// startmodule(&lf,&des,"std",0); + + fp = &lf.fp; + for (i=0; i<fp->nv; i++) + printf("%8.5f %8.5f\n",evptx(fp,i,0),fp->coef[i]); +} + +int main() +{ double x[10], y[10]; + + x[0] = 0; x[1] = 1; x[2] = 2; x[3] = 3; x[4] = 4; + x[5] = 5; x[6] = 6; x[7] = 7; x[8] = 8; x[9] = 9; + + y[0] = 0.3692449; y[1] = 0.8194270; + y[2] = 1.6363139; y[3] =-0.9969944; + y[4] = 0.5359200; y[5] = 1.8642622; + y[6] = 0.3568127; y[7] = 0.4746753; + y[8] =-2.0038246; y[9] = 1.6636109; + + calllf(x,y,10); +}