Changeset 12050

Timestamp:

Feb 25, 2007, 4:33:39 PM (19 years ago)

Author:

eugene

Message:

more pm work

Location:

branches/dvo-mods-2007-02/Ohana/src/relastro

Files:

• include/relastro.h (modified) (1 diff)
• src/FitPM.c (modified) (4 diffs)
• src/FitPMandPar.c (added)
• src/UpdateObjects.c (modified) (1 diff)

branches/dvo-mods-2007-02/Ohana/src/relastro/include/relastro.h

@@ -4 +4 @@
 # include <signal.h>
 
-/* # define GRID_V1 */
-# define GRID_V2
-
-# if (0)
 typedef struct {
-  double xmin, xmax, ymin, ymax;
-  int style, ptype, ltype, etype, color;
-  double lweight, size;
-} Graphdata;
-# endif
+  double Ro, dRo;
+  double Do, dDo;
+  
+  double uR, duR;
+  double uD, duD;
+
+  double p, dp;
+} PMfit;
 
 typedef struct {

branches/dvo-mods-2007-02/Ohana/src/relastro/src/FitPM.c

@@ -2 +2 @@
 
 /* do we want an init function which does the alloc and a clear function to free? */
-
 FitPM (double *X, double *dX, double *Y, double *dY, double *T, double *pR, double *pD, int Npts) {
 
+  PFfit *fit;
   double **A, **B;
 
-  ALLOCATE (A, double *, 5);
-  ALLOCATE (B, double *, 5);
-  for (i = 0; i < 5; i++) {
-    ALLOCATE (A[i], double, 5);
+  /* do I need to do this as 2 2x2 matrix equations? */
+  ALLOCATE (A, double *, 4);
+  ALLOCATE (B, double *, 4);
+  for (i = 0; i < 4; i++) {
+    ALLOCATE (A[i], double, 4);
     ALLOCATE (B[i], double, 1);
-    bzero (A[i], 5*sizeof(double));
+    bzero (A[i], 4*sizeof(double));
     bzero (B[i], 1*sizeof(double));
   }
@@ -30 +31 @@
     Ty2 += SQ(T[i])*wy;
     
-    PR += pR[i]*wx;
-    PD += pD[i]*wy;
-    
-    PRT += pR[i]*T[i]*wx;
-    PDT += pD[i]*T[i]*wy;
-    
-    PRX += pR[i]*X[i]*wx;
-    PDY += pD[i]*Y[i]*wy;
-    
-    PR2 += SQ(pR[i])*wx;
-    PD2 += SQ(pD[i])*wy;
-
     Xs += X[i]*wx;
     Ys += Y[i]*wy;
@@ -51 +40 @@
   A[0][0] = Wx;
   A[0][1] = Tx;
-  A[0][4] = PR;
 
   A[1][0] = Tx;
   A[1][1] = Tx2;
-  A[1][4] = PRT;
 
   A[2][2] = Wy;
   A[2][3] = Ty;
-  A[2][4] = PD;
 
   A[3][2] = Ty;
   A[3][3] = Ty2;
-  A[3][4] = PDT;
-
-  A[4][0] = PR;
-  A[4][1] = PRT;
-  A[4][2] = PD;
-  A[4][3] = PDT;
-  A[4][4] = PR2 + PD2;
 
   B[0][0] = Xs;
@@ -75 +54 @@
   B[2][0] = Ys;
   B[3][0] = YT;
-  B[4][0] = PRX + PDY;
 
-  dgaussj (A, 5, B, 1);
+  dgaussj (A, 4, B, 1);
 
-  Ro = B[0];
-  uR = B[1];
-  Do = B[2];
-  uD = B[3];
-  p  = B[4];
+  ALLOCATE (fit, PMfit, 1);
+
+  fit[0].Ro = B[0];
+  fit[0].uR = B[1];
+  fit[0].Do = B[2];
+  fit[0].uD = B[3];
+  fit[0].p  = 0.0;
   
-  dRo = sqrt(A[0][0]);
-  duR = sqrt(A[1][1]);
-  dDo = sqrt(A[2][2]);
-  duD = sqrt(A[3][3]);
-  dp  = sqrt(A[4][4]);
+  fit[0].dRo = sqrt(A[0][0]);
+  fit[0].duR = sqrt(A[1][1]);
+  fit[0].dDo = sqrt(A[2][2]);
+  fit[0].duD = sqrt(A[3][3]);
+  fit[0].dp  = 0.0;
   
-  for (i = 0; i < 5; i++) {
+  for (i = 0; i < 4; i++) {
     free (A[i]);
     free (B[i]);

branches/dvo-mods-2007-02/Ohana/src/relastro/src/UpdateObjects.c

@@ -86 +86 @@
       }
 
+      if (FitProperMotion) {
+          fit = FitPMandPar (listX, dlistX, listY, dlistY, listT, pR, pD, N);
+      }
+
       if (AverageCoords) {
           liststats (listR, dlistR, N, &statsR);