94 void palMappa(
double eq,
double date,
double amprms[21] ){
99 const double GR2 = 2.0 * 9.87063e-9;
103 double ebd[ 3 ], ehd[ 3 ], eh[ 3 ], e, vn[ 3 ], vm;
106 memset( amprms, 0, 21*
sizeof( *amprms ) );
112 palEvp( date, eq, ebd, &rms[ 1 ], ehd, eh );
115 eraPn( eh, &e, &rms[ 4 ] );
121 for( i = 0; i < 3; i++ ) {
124 eraPn( &rms[8], &vm, vn );
125 amprms[ 11 ] = sqrt( 1.0 - vm*vm );
128 palPrenut( eq, date, (
double(*)[ 3 ]) &rms[ 12 ] );
void palEvp(double date, double deqx, double dvb[3], double dpb[3], double dvh[3], double dph[3])
static const double PAL__CR
void palPrenut(double epoch, double date, double rmatpn[3][3])
void palMappa(double eq, double date, double amprms[21])
static const double PAL__MJD0
void eraPn(double p[3], double *r, double u[3])
double eraEpj(double dj1, double dj2)