IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 5306 for trunk/psLib/src/astro


Ignore:
Timestamp:
Oct 13, 2005, 10:23:57 AM (21 years ago)
Author:
drobbin
Message:

Fixed psSphereRotAlloc and test. Correctly verified calculations.

Location:
trunk/psLib/src/astro
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/psLib/src/astro/psEarthOrientation.h

    r5012 r5306  
    88*  @author Robert Daniel DeSonia, MHPCC
    99*
    10 *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
    11 *  @date $Date: 2005-09-13 00:50:04 $
     10*  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
     11*  @date $Date: 2005-10-13 20:23:57 $
    1212*
    1313*  Copyright 2005 Maui High Performance Computing Center, University of Hawaii
    … …  
    2121#include "psSphereOps.h"
    2222
     23/**
     24 *
     25 *
     26 */
    2327typedef struct
    2428{
    25     double x;
    26     double y;
    27     double s;
     29    double x;                          ///<
     30    double y;                          ///<
     31    double s;                          ///<
    2832}
    2933psEarthPole;
    3034
     35/**
     36 *
     37 *
     38 */
    3139typedef enum {
    32     PS_PRECESS_ROUGH,
    33     PS_PRECESS_COMPLETE,
    34     PS_PRECESS_IAU2000A,
     40    PS_PRECESS_ROUGH,                  ///<
     41    PS_PRECESS_COMPLETE,               ///<
     42    PS_PRECESS_IAU2000A,               ///<
    3543} psPrecessMethod;
    3644
     45/**
     46 *
     47 *
     48 *
     49 */
    3750psSphere *psAberration(
    38     psSphere *apparent,
    39     const psSphere *actual,
    40     const psSphere *direction,
    41     double speed
     51    psSphere *apparent,                ///<
     52    const psSphere *actual,            ///<
     53    const psSphere *direction,         ///<
     54    double speed                       ///<
    4255);
    4356
     57/**
     58 *
     59 *
     60 *
     61 */
    4462psSphere *psGravityDeflection(
    45     psSphere *apparent,
    46     psSphere *actual,
    47     psSphere *sun
     63    psSphere *apparent,                ///<
     64    psSphere *actual,                  ///<
     65    psSphere *sun                      ///<
    4866);
    4967
     68/**
     69 *
     70 *
     71 */
    5072double psEOC_ParallaxFactor(
    51     const psSphere *coords,
    52     const psTime *time
     73    const psSphere *coords,            ///<
     74    const psTime *time                 ///<
    5375);
    5476
     77/**
     78 *
     79 *
     80 */
    5581psEarthPole *psEOC_PrecessionModel(
    56     const psTime *time
     82    const psTime *time                 ///<
    5783);
    5884
     85/**
     86 *
     87 *
     88 */
    5989psEarthPole *psEOC_PrecessionCorr(
    60     const psTime *time,
    61     psTimeBulletin bulletin
     90    const psTime *time,                ///<
     91    psTimeBulletin bulletin            ///<
    6292);
    6393
     94/**
     95 *
     96 *
     97 */
    6498psSphereRot *psSphereRot_CEOtoGCRS(
    65     const psEarthPole *pole
     99    const psEarthPole *pole            ///<
    66100);
    67101
     102/**
     103 *
     104 *
     105 */
    68106psSphereRot *psSphereRot_TEOtoCEO(
    69     const psTime *time
     107    const psTime *time                 ///<
    70108);
    71109
     110/**
     111 *
     112 *
     113 */
    72114psEarthPole *psEOC_GetPolarMotion(
    73     const psTime *time,
    74     psTimeBulletin bulletin
     115    const psTime *time,                ///<
     116    psTimeBulletin bulletin            ///<
    75117);
    76118
     119/**
     120 *
     121 *
     122 */
    77123psEarthPole *psEOC_PolarTideCorr(
    78     const psTime *time
     124    const psTime *time                 ///<
    79125);
    80126
     127/**
     128 *
     129 *
     130 *
     131 */
    81132psEarthPole *psEOC_NutationCorr(
    82     psTime *time
     133    psTime *time                       ///<
    83134);
    84135
     136/**
     137 *
     138 *
     139 *
     140 */
    85141psSphereRot *psSphereRot_ITRStoTEO(
    86     const psEarthPole *motion
     142    const psEarthPole *motion          ///<
    87143);
    88144

  • trunk/psLib/src/astro/psSphereOps.c

    r4898 r5306  
    77 *  @author Robert DeSonia, MHPCC
    88 *
    9  *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
    10  *  @date $Date: 2005-08-30 01:14:10 $
     9 *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
     10 *  @date $Date: 2005-10-13 20:23:57 $
    1111 *
    1212 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
    … …  
    153153
    154154    // following came from ADD
    155     out->q0 = b3*a0 + b2*a1 - b1*a2 + b0*a3;
    156     out->q1 = b3*a1 - b2*a0 + b1*a3 + b0*a2;
    157     out->q2 = b3*a2 + b2*a3 + b1*a0 - b0*a1;
    158     out->q3 = b3*a3 - b3*a2 - b1*a1 - b0*a0;
     155    //    out->q0 = b3*a0 + b2*a1 - b1*a2 + b0*a3;
     156    //    out->q1 = b3*a1 - b2*a0 + b1*a3 + b0*a2;
     157    //    out->q2 = b3*a2 + b2*a3 + b1*a0 - b0*a1;
     158    out->q0 = a3*b0 + a0*b3 + a1*b2 - a2*b1;
     159    out->q1 = a3*b1 - a0*b2 + a1*b3 + a2*b0;
     160    out->q2 = a3*b2 + a0*b1 - a1*b0 + a2*b3;
     161
     162    out->q3 = b3*a3 - b2*a2 - b1*a1 - b0*a0;
    159163
    160164    return out;

  • trunk/psLib/src/astro/psSphereOps.h

    r4937 r5306  
    77 *  @author Robert DeSonia, MHPCC
    88 *
    9  *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
    10  *  @date $Date: 2005-08-31 22:28:35 $
     9 *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
     10 *  @date $Date: 2005-10-13 20:23:57 $
    1111 *
    1212 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
    … …  
    7373);
    7474
     75psSphereRot* psSphereRotAlloc2(
     76    double alphaP,                      ///< north pole latitude
     77    double deltaP,                      ///< north pole longitude
     78    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     79);
     80
     81psSphereRot* psSphereRotAlloc3(
     82    double alphaP,                      ///< north pole latitude
     83    double deltaP,                      ///< north pole longitude
     84    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     85);
     86
     87psSphereRot* psSphereRotAlloc4(
     88    double alphaP,                      ///< north pole latitude
     89    double deltaP,                      ///< north pole longitude
     90    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     91);
     92
     93psSphereRot* psSphereRotAlloc5(
     94    double alphaP,                      ///< north pole latitude
     95    double deltaP,                      ///< north pole longitude
     96    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     97);
     98
     99psSphereRot* psSphereRotAlloc6(
     100    double alphaP,                      ///< north pole latitude
     101    double deltaP,                      ///< north pole longitude
     102    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     103);
     104
     105psSphereRot* psSphereRotAlloc7(
     106    double alphaP,                      ///< north pole latitude
     107    double deltaP,                      ///< north pole longitude
     108    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     109);
     110
     111psSphereRot* psSphereRotAlloc8(
     112    double alphaP,                      ///< north pole latitude
     113    double deltaP,                      ///< north pole longitude
     114    double phiP                         ///< defines the longitude in the input system of the equatorial intersection between the two systems (e.g, the first point of Ares).
     115);
     116
    75117
    76118/** Checks the type of a particular pointer.
Note: See TracChangeset for help on using the changeset viewer.