FACT++  1.0
double eraEors ( double  rnpb[3][3],
double  s 
)

Definition at line 3 of file eors.c.

Referenced by eraApci13(), eraApco13(), eraEo06a(), eraGst06(), and t_eors().

39 {
40  double x, ax, xs, ys, zs, p, q, eo;
41 
42 /* Evaluate Wallace & Capitaine (2006) expression (16). */
43  x = rnpb[2][0];
44  ax = x / (1.0 + rnpb[2][2]);
45  xs = 1.0 - ax * x;
46  ys = -ax * rnpb[2][1];
47  zs = -x;
48  p = rnpb[0][0] * xs + rnpb[0][1] * ys + rnpb[0][2] * zs;
49  q = rnpb[1][0] * xs + rnpb[1][1] * ys + rnpb[1][2] * zs;
50  eo = ((p != 0) || (q != 0)) ? s - atan2(q, p) : s;
51 
52  return eo;
53 
54 }

+ Here is the caller graph for this function: