void palAmpqk	(	double	ra,
		double	da,
		double	amprms[21],
		double *	rm,
		double *	dm
	)

Definition at line 85 of file palAmpqk.c.

References eraAnp(), eraC2s(), eraPdp(), eraPn(), eraS2c(), eraTrxp(), and i.

Referenced by PointingModel::CalcPointingPos(), palAmp(), and t_ampqk().

                            {

/* Local Variables: */
   double ab1;                 /* sqrt(1-v*v) where v=modulus of Earth vel */
   double abv[3];              /* Earth velocity wrt SSB (c, FK5) */
   double p1[3], p2[3], p3[3]; /* work vectors */
   double ab1p1, p1dv, p1dvp1, w;
   int i, j;

/* Unpack some of the parameters */
   ab1  = amprms[11];
   for( i = 0; i < 3; i++ ) {
      abv[i] = amprms[i + 8];
   }

/* Apparent RA,Dec to Cartesian */
   eraS2c( ra, da, p3 );

/* Precession and nutation */
   eraTrxp( (double(*)[3]) &amprms[12], p3, p2 );

/* Aberration */
   ab1p1 = ab1 + 1.0;
   for( i = 0; i < 3; i++ ) {
      p1[i] = p2[i];
   }
   for( j = 0; j < 2; j++ ) {
      p1dv = eraPdp( p1, abv );
      p1dvp1 = 1.0 + p1dv;
      w = 1.0 + p1dv / ab1p1;
      for( i = 0; i < 3; i++ ) {
         p1[i] = ( p1dvp1 * p2[i] - w * abv[i] ) / ab1;
      }
      eraPn( p1, &w, p3 );
      for( i = 0; i < 3; i++ ) {
         p1[i] = p3[i];
      }
   }

/* Mean RA,Dec */
   eraC2s( p1, rm, dm );
   *rm = eraAnp( *rm );
}

Here is the call graph for this function:

Here is the caller graph for this function: