Changeset 39238 for trunk/Ohana/src/relastro/include/relastro.h

Timestamp:

Dec 8, 2015, 1:56:09 PM (11 years ago)

Author:

eugene

Message:

adding IRLS to relastro

File:

• trunk/Ohana/src/relastro/include/relastro.h (modified) (4 diffs)

trunk/Ohana/src/relastro/include/relastro.h

@@ -134 +134 @@
   double **A;
   double **B;
+  double **Cov;
+  double *Beta;
+  double *Beta_prev;
   int Nterms;
 } FitAstromData;
 
+// XXX do we need doubles for all of these?  I actually only have of order 100 of these
+// allocated at a time, so size is not an issue.
 typedef struct {
   double X, dX;
@@ -143 +148 @@
   double D, dD;
   double T, dT;
+  double Wx, Wy;
+  double rx, ry;
+  double u;
   double pR;
   double pD;
@@ -148 +156 @@
   double C_red;
   int measure;
+  int mask;
 } FitAstromPoint;
 
@@ -739 +748 @@
 int FitAstromResultSetPM (FitAstromResult *fit, int Nfit, Average *average);
 void AstromErrorSetLoop (int Nloop, int isImageMode);
+
+int FitPM_IRLS (FitAstromResult *fit, FitAstromData *data, FitAstromPoint *points, int Npoints, int VERBOSE);
+
+double MedianAbsDeviation(FitAstromPoint *points, int Npoints);
+
+int weighted_LS_PM (FitAstromResult *fit, FitAstromData *data, FitAstromPoint *points, int Npoints, int VERBOSE);
+
+double weight_cauchy (double x);
+double dpsi_cauchy (double x);