Changeset 26764 for branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/burnutils.c

Timestamp:

Feb 2, 2010, 7:06:18 PM (17 years ago)

Author:

watersc1

Message:

new burntool code

File:

• branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/burnutils.c (modified) (1 diff)

branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/burnutils.c

@@ -290 +290 @@
    return(0);
 }
+
+/****************************************************************/
+/* Return of RMS relative to Fit y = ax + b */
+STATIC int linearrms(int npt, double *x, double *y, double a, double b, 
+                     double *rms)
+{
+   *rms = 0.0;
+   while(--npt >= 0) *rms += (y[npt]-a*x[npt]-b)*(y[npt]-a*x[npt]-b);
+   if(*rms > 0) *rms = sqrt(*rms);
+   return(0);
+}