Changeset 1463 for trunk/psLib/src/astronomy/psAstrometry.c

Timestamp:

Aug 10, 2004, 1:59:41 PM (22 years ago)

Author:

gusciora

Message:

This is the first compiling version of astometry.c

File:

• trunk/psLib/src/astronomy/psAstrometry.c (modified) (13 diffs)

trunk/psLib/src/astronomy/psAstrometry.c

@@ -1 @@
-
 /** @file  psAstrometry.c
 *
@@ -9 +8 @@
 *  @author George Gusciora, MHPCC
 *
-*  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-08-09 23:34:57 $
+*  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-08-10 23:59:41 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -21 +20 @@
 #include "psMemory.h"
 
-static void grommitFree(psGrommit* grommit);
 static int checkValidChipCoords(double x, double y, psChip* tmpChip);
 static int checkValidImageCoords(double x, double y, psImage* tmpImage);
@@ -49 +47 @@
 }
 
+/*
+    Several members of the psGrommit data structure have no direct counterpart
+    in the psExposure structure.  How are we to determine them?
+*/
 psGrommit* psGrommitAlloc(const psExposure* exp)
 {
@@ -79 +81 @@
 }
 
-void p_psGrommitFree(psGrommit* grommit)
+void psGrommitFree(psGrommit* grommit)
 {
     psFree(grommit);
@@ -94 +96 @@
 }
 
+/*
+    XXX: do this
+*/
 int checkValidChipCoords(double x, double y, psChip* tmpChip)
 {
@@ -120 +125 @@
 }
 
+/*
+    XXX: do this
+*/
 int checkValidImageCoords(double x, double y, psImage* tmpImage)
 {
@@ -142 +150 @@
 psCell* psCellinChip(psCell* out, const psPlane* coord, const psChip* chip)
 {
+    int i = 0;
     psPlane* tmpCoord = NULL;
     psArray* cells;
@@ -155 +164 @@
     }
 
-    for (int i = 0; i < cells->n; i++) {
+    for (i = 0; i < cells->n; i++) {
         psCell* tmpCell = (psCell* ) cells->data[i];
         psArray* readouts = tmpCell->readouts;
@@ -210 +219 @@
 // transformation, as well as a few psPlane structs.  Can this be implemented
 // better?
-psSphere* psCoordCelltoSky(psSphere* out, const psPlane* in, const psCell* cell)
+psSphere* psCoordCelltoSky(psSphere* out,
+                           const psPlane* in,
+                           const psCell* cell)
 {
     psPlane* tmp1 = NULL;
@@ -216 +227 @@
     psFPA* parFPA = (cell->parent)->parent;
     psGrommit* tmpGrommit = NULL;
+    float XXX1 = 0.0;
+    float XXX2 = 0.0;
+
 
     tmp1 = psPlaneTransformApply(tmp1, cell->toFPA, in);
-    tmp2 = psPlaneTransformApply(tmp2, parFPA->toTangentPlane, tmp1);
+    // XXX:
+    //
+    tmp2 = psPlaneDistortApply(tmp2, parFPA->toTangentPlane, tmp1, XXX1, XXX2);
     tmpGrommit = psGrommitAlloc(parFPA->exposure);
-    tmp3 = psCoordTPtoSky(out, tmp2, psGrommit);
-
     psFree(tmp1);
     psFree(tmp2);
     psFree(tmpGrommit);
 
-    return (psCoordTPtoSky(out, tmp2, psGrommit));
-
+    return(psCoordTPtoSky(out, tmp2, tmpGrommit));
 }
 
@@ -247 +260 @@
 psPlane* psCoordSkytoTP(psPlane* out, const psSphere* in, const psGrommit* grommit)
 {
-    extern void sla_AOPQK(RAP, DAP, AOPRMS, AOB, ZOB, HOB, DOB, ROB);
-    double AOB;
-    double ZOB;
-    double HOB;
-    double DOB;
-    double ROB;
-
-    if (out == NULL) {
-        out = (psPlane* ) psAlloc(sizeof(psPlane));
-    }
-
-    sla_AOPQK(psSphere->r, psSphere->d, *grommit, &AOB, &ZOB, &HOB, &DOB, &ROB);
-    out->x = XXX;
-    out->y = XXX;
-    return (out);
-}
-
+    /*
+        extern void sla_AOPQK(RAP, DAP, AOPRMS, AOB, ZOB, HOB, DOB, ROB);
+        double AOB;
+        double ZOB;
+        double HOB;
+        double DOB;
+        double ROB;
+     
+        if (out == NULL) {
+            out = (psPlane* ) psAlloc(sizeof(psPlane));
+        }
+     
+        sla_AOPQK(psSphere->r, psSphere->d, *grommit, &AOB, &ZOB, &HOB, &DOB, &ROB);
+        out->x = XXX;
+        out->y = XXX;
+        return (out);
+    */
+    return(NULL);
+}
+
+
+/*
+    XXX: What are the XXX1 and XXX2 args supposed to be?
+*/
 psPlane* psCoordTPtoFPA(psPlane* out, const psPlane* in, const psFPA* fpa)
 {
-    return (psPlaneTransformApply(out, fpa->fromTangentPlane, in));
+    float XXX1 = 0.0;
+    float XXX2 = 0.0;
+
+    return (psPlaneDistortApply(out, fpa->fromTangentPlane, in, XXX1, XXX2));
 }
 
@@ -281 +304 @@
 psPlane* psCoordSkytoCell(psPlane* out, const psSphere* in, const psCell* cell)
 {
-    out = psCoordSkytoTP(out, in, tmpGrommit);
+    psGrommit* XXXGrommit = NULL;
+    psFPA *whichFPA = NULL;
+    psChip *whichChip = NULL;
+    psCell *whichCell = NULL;
+
+    out = psCoordSkytoTP(out, in, XXXGrommit);
     out = psCoordTPtoFPA(out, out, whichFPA);
     out = psCoordFPAtoChip(out, out, whichChip);