Changeset 39238 for trunk/Ohana/src/relastro/src/FitAstromOps.c

Timestamp:

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

Author:

eugene

Message:

adding IRLS to relastro

File:

• trunk/Ohana/src/relastro/src/FitAstromOps.c (modified) (3 diffs)

trunk/Ohana/src/relastro/src/FitAstromOps.c

@@ -93 +93 @@
 
   /* do I need to do this as 2 2x2 matrix equations? */
-  fit->A = array_init (Nterms, Nterms);
-  fit->B = array_init (Nterms, 1);
+  fit->B   = array_init (Nterms, 1);
+  fit->A   = array_init (Nterms, Nterms);
+  fit->Cov = array_init (Nterms, Nterms);
+
+  ALLOCATE (fit->Beta, double, Nterms);
+  ALLOCATE (fit->Beta_prev, double, Nterms);
   fit->Nterms = Nterms;
 
@@ -106 +110 @@
   array_free (fit->A, fit->Nterms);
   array_free (fit->B, fit->Nterms);
+  array_free (fit->Cov, fit->Nterms);
+
+  free (fit->Beta);
+  free (fit->Beta_prev);
+
   free (fit);
   return;
@@ -126 +135 @@
   object->C_red  = 0.0;
   object->measure= -1;
+
+  object->Wx     = 1.0;
+  object->Wy     = 1.0;
+
+  object->rx     = 0.0;
+  object->ry     = 0.0;
+  object->u      = 0.0;
+
+  object->mask   = 0;
   return;
 }