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
	)

Definition at line 3 of file pmsafe.c.

References eraSeps(), and eraStarpm().

Referenced by t_pmsafe().

{

/* Minimum allowed parallax (arcsec) */
   const double PXMIN = 5e-7;

/* Factor giving maximum allowed transverse speed of about 1% c */
   const double F = 326.0;

   int jpx, j;
   double pm, px1a;

/* Proper motion in one year (radians). */
   pm = eraSeps(ra1, dec1, ra1+pmr1, dec1+pmd1);

/* Override the parallax to reduce the chances of a warning status. */
   jpx = 0;
   px1a = px1;
   pm *= F;
   if (px1a < pm) {jpx = 1; px1a = pm;}
   if (px1a < PXMIN) {jpx = 1; px1a = PXMIN;}

/* Carry out the transformation using the modified parallax. */
   j = eraStarpm(ra1, dec1, pmr1, pmd1, px1a, rv1,
                 ep1a, ep1b, ep2a, ep2b,
                 ra2, dec2, pmr2, pmd2, px2, rv2);

/* Revise and return the status. */
   if ( !(j%2) ) j += jpx;
   return j;

/* Finished. */

}

Here is the call graph for this function:

Here is the caller graph for this function: