Changeset 39387

Timestamp:

Feb 24, 2016, 10:43:23 AM (10 years ago)

Author:

eugene

Message:

remove warning at 0,360 boundary

Location:

trunk/Ohana/src/relastro/src

Files:

• UpdateObjects.c (modified) (4 diffs)
• fitstk.c (modified) (2 diffs)

trunk/Ohana/src/relastro/src/UpdateObjects.c

@@ -205 +205 @@
     // project Ro, Do back to RA,DEC
     XY_to_RD (&fitPar.Ro, &fitPar.Do, fitPar.Ro, fitPar.Do, &fitStats->coords);
-    if (fabs(fitPar.Ro) < 0.01) fprintf (stderr, "watch out for 0,360 boundary\n");
 
     average[0].flags |= ID_OBJ_FIT_PAR;
@@ -275 +274 @@
     // project Ro, Do back to RA,DEC
     XY_to_RD (&fitPM.Ro, &fitPM.Do, fitPM.Ro, fitPM.Do, &fitStats->coords);
-    if (fabs(fitPM.Ro) < 0.01) fprintf (stderr, "watch out for 0,360 boundary\n");
 
     average[0].flags |= ID_OBJ_FIT_PM;
@@ -363 +361 @@
     // project Ro, Do back to RA,DEC
     XY_to_RD (&fitPos.Ro, &fitPos.Do, fitPos.Ro, fitPos.Do, &fitStats->coords);
-    if (fabs(fitPos.Ro) < 0.01) fprintf (stderr, "watch out for 0,360 boundary\n");
 
     average[0].flags |= ID_OBJ_FIT_AVE;
@@ -378 +375 @@
     // project Ro, Do back to RA,DEC
     XY_to_RD (&fitPos.Ro, &fitPos.Do, fitPos.Ro, fitPos.Do, &fitStats->coords);
-    if (fabs(fitPos.Ro) < 0.01) fprintf (stderr, "watch out for 0,360 boundary\n");
 
     average[0].flags |= ID_OBJ_USE_AVE;

trunk/Ohana/src/relastro/src/fitstk.c

@@ -120 +120 @@
   for (i = 0; i < Nstars; i++) {
 
-    double Ro  = 360.0*(drand48_cnt() - 0.5);
+    // double Ro  = 360.0*drand48_cnt();
+    double Ro  = 0.1*drand48_cnt();
+
     double Phi = 2.0*(drand48_cnt() - 0.5);
     double Do  = DEG_RAD * asin(Phi);
@@ -193 +195 @@
     catalog->measureT[i].D = Do + (dD + plx*pD + uD*(Tyears - MJD_REF_YRS))/3600;
 
+    // force some points to loop over the boundary:
+    if (drand48() > 0.5) {
+      catalog->measureT[i].R += 360.0;
+    }
+    if (drand48() > 0.8) {
+      catalog->measureT[i].R += 360.0;
+    }
+
     catalog->measureT[i].dXccd = ToShortPixels(POS_ERROR);
     catalog->measureT[i].dYccd = ToShortPixels(POS_ERROR);