FACT++  1.0
void eraPn06 ( double  date1,
double  date2,
double  dpsi,
double  deps,
double *  epsa,
double  rb[3][3],
double  rp[3][3],
double  rbp[3][3],
double  rn[3][3],
double  rbpn[3][3] 
)

Definition at line 3 of file pn06.c.

References eraCr(), eraFw2m(), eraPfw06(), eraRxr(), eraTr(), ERFA_DJM0, and ERFA_DJM00.

Referenced by eraPn06a(), and t_pn06().

99 {
100  double gamb, phib, psib, eps, r1[3][3], r2[3][3], rt[3][3];
101 
102 /* Bias-precession Fukushima-Williams angles of J2000.0 = frame bias. */
103  eraPfw06(ERFA_DJM0, ERFA_DJM00, &gamb, &phib, &psib, &eps);
104 
105 /* B matrix. */
106  eraFw2m(gamb, phib, psib, eps, r1);
107  eraCr(r1, rb);
108 
109 /* Bias-precession Fukushima-Williams angles of date. */
110  eraPfw06(date1, date2, &gamb, &phib, &psib, &eps);
111 
112 /* Bias-precession matrix. */
113  eraFw2m(gamb, phib, psib, eps, r2);
114  eraCr(r2, rbp);
115 
116 /* Solve for precession matrix. */
117  eraTr(r1, rt);
118  eraRxr(r2, rt, rp);
119 
120 /* Equinox-based bias-precession-nutation matrix. */
121  eraFw2m(gamb, phib, psib + dpsi, eps + deps, r1);
122  eraCr(r1, rbpn);
123 
124 /* Solve for nutation matrix. */
125  eraTr(r2, rt);
126  eraRxr(r1, rt, rn);
127 
128 /* Obliquity, mean of date. */
129  *epsa = eps;
130 
131  return;
132 
133 }
void eraPfw06(double date1, double date2, double *gamb, double *phib, double *psib, double *epsa)
Definition: pfw06.c:3
#define ERFA_DJM00
Definition: erfam.h:93
void eraFw2m(double gamb, double phib, double psi, double eps, double r[3][3])
Definition: fw2m.c:3
#define ERFA_DJM0
Definition: erfam.h:90
void eraCr(double r[3][3], double c[3][3])
Definition: cr.c:3
void eraRxr(double a[3][3], double b[3][3], double atb[3][3])
Definition: rxr.c:3
void eraTr(double r[3][3], double rt[3][3])
Definition: tr.c:3

+ Here is the call graph for this function:

+ Here is the caller graph for this function: