Changeset 39610 for trunk/Ohana/src/opihi/include/astro.h

Timestamp:

Jun 24, 2016, 5:07:24 PM (10 years ago)

Author:

eugene

Message:

modify irls fitting to apply the multiply the modified weight by the error-based weight

File:

• trunk/Ohana/src/opihi/include/astro.h (modified) (2 diffs)

trunk/Ohana/src/opihi/include/astro.h

@@ -20 +20 @@
   double *dX;
   double *dY;
-  double *Wx;
-  double *Wy;
+  double *Wx; // modified weight based on distance from fit
+  double *Wy; // modified weight based on distance from fit
+  double *Qx; // raw error-based weight (1/dX^2)
+  double *Qy;
+  double *qx; // modified error-based weight (1/dX^2) * Wx
+  double *qy;
   int *index;
   int Npts;
@@ -59 +63 @@
 
 int FitPMonly_IRLS (PlxFit *fit, double *X, double *dX, double *Y, double *dY, double *T, int Npts, double binning_step, int VERBOSE);
-int FitPMandPar_IRLS (PlxFit *fit, double *X, double *dX, double *Y, double *dY, double *T, double *pR, double *pD, double *Wx, double *Wy, int Npts, int max_iterations, double outlier_limit, double binning_step, int VERBOSE);
+int FitPMandPar_IRLS (PlxFit *fit, double *X, double *dX, double *Y, double *dY, double *T, double *pR, double *pD, double *Wx, double *Wy, double *Qx, double *Qy, double *qx, double *qy, int Npts, int max_iterations, double outlier_limit, double binning_step, int VERBOSE);
 
 int weighted_LS_PLX (double *T, double *pR, double *pD, double *X, double *WX, double *Y, double *WY, int Npts, double **A, double **B, int VERBOSE);