3 int eraPmsafe(
double ra1,
double dec1,
double pmr1,
double pmd1,
4 double px1,
double rv1,
5 double ep1a,
double ep1b,
double ep2a,
double ep2b,
6 double *ra2,
double *dec2,
double *pmr2,
double *pmd2,
7 double *px2,
double *rv2)
114 const double PXMIN = 5e-7;
117 const double F = 326.0;
123 pm =
eraSeps(ra1, dec1, ra1+pmr1, dec1+pmd1);
129 if (px1a < pm) {jpx = 1; px1a = pm;}
130 if (px1a < PXMIN) {jpx = 1; px1a = PXMIN;}
133 j =
eraStarpm(ra1, dec1, pmr1, pmd1, px1a, rv1,
134 ep1a, ep1b, ep2a, ep2b,
135 ra2, dec2, pmr2, pmd2, px2, rv2);
138 if ( !(j%2) ) j += jpx;
double eraSeps(double al, double ap, double bl, double bp)
int eraStarpm(double ra1, double dec1, double pmr1, double pmd1, double px1, double rv1, double ep1a, double ep1b, double ep2a, double ep2b, double *ra2, double *dec2, double *pmr2, double *pmd2, double *px2, double *rv2)
int eraPmsafe(double ra1, double dec1, double pmr1, double pmd1, double px1, double rv1, double ep1a, double ep1b, double ep2a, double ep2b, double *ra2, double *dec2, double *pmr2, double *pmd2, double *px2, double *rv2)