Changeset 40456

Timestamp:

Jun 27, 2018, 1:05:36 PM (8 years ago)

Author:

eugene

Message:

use fp, not tp coordinates

File:

• branches/czw_branch/20170908/psastro/src/psastroLoadGlints.c (modified) (5 diffs)

branches/czw_branch/20170908/psastro/src/psastroLoadGlints.c

@@ -305 +305 @@
               double GLINT_RADIUS_OUTER = psMetadataLookupF32(&status, glintItem->data.md, "GLINT.RADIUS.OUTER");
               
-              //              double R_FPA = sqrt(pow(star->FP->y,2) + pow(star->FP->x,2));
-              double R_FPA = sqrt(pow(star->TP->y,2) + pow(star->TP->x,2));
+              // double R_FPA = sqrt(pow(star->FP->y,2) + pow(star->FP->x,2));
+              double R_FPA =  sqrt(pow(star->TP->y,2) + pow(star->TP->x,2));
               //              R_FPA *= 1.30;
               if (R_FPA < GLINT_RADIUS_INNER) { continue; }
@@ -420 +420 @@
               x_start->data.F32[0] = fp->x;
               y_start->data.F32[0] = fp->y;
-              x_start->data.F32[0] = tp->x;
-              y_start->data.F32[0] = tp->y;
+              //              x_start->data.F32[0] = tp->x;
+              //              y_start->data.F32[0] = tp->y;
 
               tp->x = 13.5 * (y1e / 0.015 + 12.78);
@@ -428 +428 @@
               x_end->data.F32[0]   = fp->x;
               y_end->data.F32[0]   = fp->y;
-              x_end->data.F32[0]   = tp->x;
-              y_end->data.F32[0]   = tp->y;
+              //              x_end->data.F32[0]   = tp->x;
+              //              y_end->data.F32[0]   = tp->y;
 
               x_start->data.F32[1] = x_end->data.F32[0];  
@@ -441 +441 @@
               x_start->data.F32[2] = fp->x;
               y_start->data.F32[2] = fp->y;
-              x_start->data.F32[2] = tp->x;
-              y_start->data.F32[2] = tp->y;
+              //              x_start->data.F32[2] = tp->x;
+              //              y_start->data.F32[2] = tp->y;
 
               tp->x = 13.5 * (y2e / 0.015 + 12.78);
@@ -449 +449 @@
               x_end->data.F32[2]   = fp->x;
               y_end->data.F32[2]   = fp->y;
-              x_end->data.F32[2]   = tp->x;
-              y_end->data.F32[2]   = tp->y;
+              //              x_end->data.F32[2]   = tp->x;
+              //              y_end->data.F32[2]   = tp->y;
               
               x_start->data.F32[3] = x_end->data.F32[2];