elp82b.c   elp82b.c 
skipping to change at line 162749 skipping to change at line 162749
const double x1 = rh * cos(r[0]); const double x1 = rh * cos(r[0]);
const double x2 = rh * sin(r[0]); const double x2 = rh * sin(r[0]);
double pw = t*(p1 + t*(p2 + t*(p3 + t*(p4 + t*p5)))); double pw = t*(p1 + t*(p2 + t*(p3 + t*(p4 + t*p5))));
double qw = t*(q1 + t*(q2 + t*(q3 + t*(q4 + t*q5)))); double qw = t*(q1 + t*(q2 + t*(q3 + t*(q4 + t*q5))));
const double pwq = pw * pw; const double pwq = pw * pw;
const double qwq = qw * qw; const double qwq = qw * qw;
const double pwqw = 2.0 * pw * qw; const double pwqw = 2.0 * pw * qw;
const double pw2 = 1.0 - 2.0 * pwq; const double pw2 = 1.0 - 2.0 * pwq;
const double qw2 = 1.0 - 2.0 * qwq; const double qw2 = 1.0 - 2.0 * qwq;
const double ra = 2.0 * sqrt(1.0 - pwq - qwq); const double h = 1.0 - pwq - qwq;
const double ra = (h > 0.0) ? (2.0 * sqrt(h)) : 0.0;
pw *= ra; pw *= ra;
qw *= ra; qw *= ra;
/* VSOP87 coordinates: */ /* VSOP87 coordinates: */
xyz[0] = pw2 *x1 + pwqw*x2 + pw*x3; xyz[0] = pw2 *x1 + pwqw*x2 + pw*x3;
xyz[1] = pwqw*x1 + qw2 *x2 - qw*x3; xyz[1] = pwqw*x1 + qw2 *x2 - qw*x3;
xyz[2] = -pw *x1 + qw *x2 + (pw2 + qw2 - 1.0)*x3; xyz[2] = -pw *x1 + qw *x2 + (pw2 + qw2 - 1.0)*x3;
/* /*
printf("Moon: %f %22.15f %22.15f %22.15f\n", printf("Moon: %f %22.15f %22.15f %22.15f\n",
 End of changes. 1 change blocks. 
1 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/