void eraPnm00b(double date1, double date2, double rbpn[3][3])
/*
** - - - - - - - - - -
** e r a P n m 0 0 b
** - - - - - - - - - -
**
** Form the matrix of precession-nutation for a given date (including
** frame bias), equinox-based, IAU 2000B model.
**
** Given:
** date1,date2 double TT as a 2-part Julian Date (Note 1)
**
** Returned:
** rbpn double[3][3] bias-precession-nutation matrix (Note 2)
**
** Notes:
**
** 1) The TT date date1+date2 is a Julian Date, apportioned in any
** convenient way between the two arguments. For example,
** JD(TT)=2450123.7 could be expressed in any of these ways,
** among others:
**
** date1 date2
**
** 2450123.7 0.0 (JD method)
** 2451545.0 -1421.3 (J2000 method)
** 2400000.5 50123.2 (MJD method)
** 2450123.5 0.2 (date & time method)
**
** The JD method is the most natural and convenient to use in
** cases where the loss of several decimal digits of resolution
** is acceptable. The J2000 method is best matched to the way
** the argument is handled internally and will deliver the
** optimum resolution. The MJD method and the date & time methods
** are both good compromises between resolution and convenience.
**
** 2) The matrix operates in the sense V(date) = rbpn * V(GCRS), where
** the p-vector V(date) is with respect to the true equatorial triad
** of date date1+date2 and the p-vector V(GCRS) is with respect to
** the Geocentric Celestial Reference System (IAU, 2000).
**
** 3) The present function is faster, but slightly less accurate (about
** 1 mas), than the eraPnm00a function.
**
** Called:
** eraPn00b bias/precession/nutation, IAU 2000B
**
** Reference:
**
** IAU: Trans. International Astronomical Union, Vol. XXIVB; Proc.
** 24th General Assembly, Manchester, UK. Resolutions B1.3, B1.6.
** (2000)
**
** Copyright (C) 2013-2016, NumFOCUS Foundation.
** Derived, with permission, from the SOFA library. See notes at end of file.
*/
{
double dpsi, deps, epsa, rb[3][3], rp[3][3], rbp[3][3], rn[3][3];
/* Obtain the required matrix (discarding other results). */
eraPn00b(date1, date2, &dpsi, &deps, &epsa, rb, rp, rbp, rn, rbpn);
return;
}
void eraNum00b(double date1, double date2, double rmatn[3][3])
/*
** - - - - - - - - - -
** e r a N u m 0 0 b
** - - - - - - - - - -
**
** Form the matrix of nutation for a given date, IAU 2000B model.
**
** Given:
** date1,date2 double TT as a 2-part Julian Date (Note 1)
**
** Returned:
** rmatn double[3][3] nutation matrix
**
** Notes:
**
** 1) The TT date date1+date2 is a Julian Date, apportioned in any
** convenient way between the two arguments. For example,
** JD(TT)=2450123.7 could be expressed in any of these ways,
** among others:
**
** date1 date2
**
** 2450123.7 0.0 (JD method)
** 2451545.0 -1421.3 (J2000 method)
** 2400000.5 50123.2 (MJD method)
** 2450123.5 0.2 (date & time method)
**
** The JD method is the most natural and convenient to use in
** cases where the loss of several decimal digits of resolution
** is acceptable. The J2000 method is best matched to the way
** the argument is handled internally and will deliver the
** optimum resolution. The MJD method and the date & time methods
** are both good compromises between resolution and convenience.
**
** 2) The matrix operates in the sense V(true) = rmatn * V(mean), where
** the p-vector V(true) is with respect to the true equatorial triad
** of date and the p-vector V(mean) is with respect to the mean
** equatorial triad of date.
**
** 3) The present function is faster, but slightly less accurate (about
** 1 mas), than the eraNum00a function.
**
** Called:
** eraPn00b bias/precession/nutation, IAU 2000B
**
** Reference:
**
** Explanatory Supplement to the Astronomical Almanac,
** P. Kenneth Seidelmann (ed), University Science Books (1992),
** Section 3.222-3 (p114).
**
** Copyright (C) 2013, NumFOCUS Foundation.
** Derived, with permission, from the SOFA library. See notes at end of file.
*/
{
double dpsi, deps, epsa, rb[3][3], rp[3][3], rbp[3][3], rbpn[3][3];
/* Obtain the required matrix (discarding other results). */
eraPn00b(date1, date2,
&dpsi, &deps, &epsa, rb, rp, rbp, rmatn, rbpn);
return;
}
