85 void palAmpqk (
double ra,
double da,
double amprms[21],
double *rm,
91 double p1[3], p2[3], p3[3];
92 double ab1p1, p1dv, p1dvp1, w;
97 for( i = 0; i < 3; i++ ) {
98 abv[
i] = amprms[i + 8];
105 eraTrxp( (
double(*)[3]) &rms[12], p3, p2 );
109 for( i = 0; i < 3; i++ ) {
112 for( j = 0; j < 2; j++ ) {
115 w = 1.0 + p1dv / ab1p1;
116 for( i = 0; i < 3; i++ ) {
117 p1[
i] = ( p1dvp1 * p2[
i] - w * abv[
i] ) / ab1;
120 for( i = 0; i < 3; i++ ) {
double eraPdp(double a[3], double b[3])
void eraS2c(double theta, double phi, double c[3])
void eraC2s(double p[3], double *theta, double *phi)
void eraTrxp(double r[3][3], double p[3], double trp[3])
void eraPn(double p[3], double *r, double u[3])
void palAmpqk(double ra, double da, double amprms[21], double *rm, double *dm)