void eraAb	(	double	pnat[3],
		double	v[3],
		double	s,
		double	bm1,
		double	ppr[3]
	)

Definition at line 3 of file ab.c.

References eraPdp(), ERFA_SRS, and i.

Referenced by eraAtciq(), eraAtciqn(), eraAtciqz(), eraAticq(), eraAticqn(), and t_ab().

{
   int i;
   double pdv, w1, w2, r2, w, p[3], r;

   pdv = eraPdp(pnat, v);
   w1 = 1.0 + pdv/(1.0 + bm1);
   w2 = ERFA_SRS/s;
   r2 = 0.0;
   for (i = 0; i < 3; i++) {
      w = pnat[i]*bm1 + w1*v[i] + w2*(v[i] - pdv*pnat[i]);
      p[i] = w;
      r2 = r2 + w*w;
   }
   r = sqrt(r2);
   for (i = 0; i < 3; i++) {
      ppr[i] = p[i]/r;
   }

/* Finished. */

}

Here is the call graph for this function:

Here is the caller graph for this function: