83 palDs2tp (
double ra,
double dec,
double raz,
double decz,
84 double *xi,
double *eta,
int *j ) {
86 const double TINY = 1.0e-6;
107 denom = sdec * sdecz + cdec * cdecz * cradif;
112 }
else if (denom >= 0.) {
115 }
else if (denom > -TINY) {
123 *xi = cdec * sradif / denom;
124 *eta = (sdec * cdecz - cdec * sdecz * cradif) / denom;
void palDs2tp(double ra, double dec, double raz, double decz, double *xi, double *eta, int *j)