s06.c | s06.c | |||
---|---|---|---|---|
skipping to change at line 230 | skipping to change at line 230 | |||
/* Number of terms in the series */ | /* Number of terms in the series */ | |||
static const int NS0 = (int) (sizeof s0 / sizeof (TERM)); | static const int NS0 = (int) (sizeof s0 / sizeof (TERM)); | |||
static const int NS1 = (int) (sizeof s1 / sizeof (TERM)); | static const int NS1 = (int) (sizeof s1 / sizeof (TERM)); | |||
static const int NS2 = (int) (sizeof s2 / sizeof (TERM)); | static const int NS2 = (int) (sizeof s2 / sizeof (TERM)); | |||
static const int NS3 = (int) (sizeof s3 / sizeof (TERM)); | static const int NS3 = (int) (sizeof s3 / sizeof (TERM)); | |||
static const int NS4 = (int) (sizeof s4 / sizeof (TERM)); | static const int NS4 = (int) (sizeof s4 / sizeof (TERM)); | |||
/*--------------------------------------------------------------------*/ | /*--------------------------------------------------------------------*/ | |||
/* Interval between fundamental epoch J2000.0 and current date (JC). */ | /* Interval between fundamental epoch J2000.0 and current date (JC). */ | |||
t = ((date1 - DJ00) + date2) / DJC; | t = ((date1 - ERFA_DJ00) + date2) / ERFA_DJC; | |||
/* Fundamental Arguments (from IERS Conventions 2003) */ | /* Fundamental Arguments (from IERS Conventions 2003) */ | |||
/* Mean anomaly of the Moon. */ | /* Mean anomaly of the Moon. */ | |||
fa[0] = eraFal03(t); | fa[0] = eraFal03(t); | |||
/* Mean anomaly of the Sun. */ | /* Mean anomaly of the Sun. */ | |||
fa[1] = eraFalp03(t); | fa[1] = eraFalp03(t); | |||
/* Mean longitude of the Moon minus that of the ascending node. */ | /* Mean longitude of the Moon minus that of the ascending node. */ | |||
skipping to change at line 311 | skipping to change at line 311 | |||
a += (double)s4[i].nfa[j] * fa[j]; | a += (double)s4[i].nfa[j] * fa[j]; | |||
} | } | |||
w4 += s4[i].s * sin(a) + s4[i].c * cos(a); | w4 += s4[i].s * sin(a) + s4[i].c * cos(a); | |||
} | } | |||
s = (w0 + | s = (w0 + | |||
(w1 + | (w1 + | |||
(w2 + | (w2 + | |||
(w3 + | (w3 + | |||
(w4 + | (w4 + | |||
w5 * t) * t) * t) * t) * t) * DAS2R - x*y/2.0; | w5 * t) * t) * t) * t) * t) * ERFA_DAS2R - x*y/2.0; | |||
return s; | return s; | |||
} | } | |||
/*---------------------------------------------------------------------- | /*---------------------------------------------------------------------- | |||
** | ** | |||
** | ** | |||
** Copyright (C) 2013, NumFOCUS Foundation. | ** Copyright (C) 2013, NumFOCUS Foundation. | |||
** All rights reserved. | ** All rights reserved. | |||
** | ** | |||
End of changes. 2 change blocks. | ||||
2 lines changed or deleted | 2 lines changed or added | |||
This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/ |