l1.c   l1.c 
skipping to change at line 929 skipping to change at line 929
for (j=2;j<4;j++) { for (j=2;j<4;j++) {
const struct L1Term *const begin = bp->lists[j].terms; const struct L1Term *const begin = bp->lists[j].terms;
const struct L1Term *p = begin + bp->lists[j].size; const struct L1Term *p = begin + bp->lists[j].size;
while (--p >= begin) { while (--p >= begin) {
const double d = p->phase + t*p->frequency; const double d = p->phase + t*p->frequency;
elem[2*j-2] += p->amplitude * cos(d); elem[2*j-2] += p->amplitude * cos(d);
elem[2*j-1] += p->amplitude * sin(d); elem[2*j-1] += p->amplitude * sin(d);
} }
} }
elem++;
for (j=0;j<5;j++) { for (j=0;j<5;j++) {
elem[j] += (*cheb++); elem[j] += (*cheb++);
} }
/* This is an artificial hack. Valery Lainey recommends /* This is an artificial hack. Valery Lainey recommends
not using polynoials outside of [-700,700]. not using polynoials outside of [-700,700].
I do the following: I do the following:
inside [-695,695] use polynomials inside [-695,695] use polynomials
outside [-705,705] do not use polynomials outside [-705,705] do not use polynomials
inbetween avoid jumping moons by partly using polynomials */ inbetween avoid jumping moons by partly using polynomials */
skipping to change at line 962 skipping to change at line 961
for (i=2;i<9;i++) { for (i=2;i<9;i++) {
const double ti = 2.0 * x * ti1 - ti2; const double ti = 2.0 * x * ti1 - ti2;
ti2 = ti1; ti2 = ti1;
ti1 = ti; ti1 = ti;
for (j=0;j<5;j++) { for (j=0;j<5;j++) {
elem[j] += use_polynomials * ((*cheb++) * ti); elem[j] += use_polynomials * ((*cheb++) * ti);
} }
} }
} }
elem[0] += t*bp->l; elem[1] += t*bp->l;
} }
/* /*
const double L1toJ2000[9] = { const double L1toJ2000[9] = {
9.994327653386544723e-01, 3.039594289062820484e-02,-1.449945596633516053 e-02, 9.994327653386544723e-01, 3.039594289062820484e-02,-1.449945596633516053 e-02,
-3.367710746976384242e-02, 9.020579123528089974e-01,-4.302991694091006926 e-01, -3.367710746976384242e-02, 9.020579123528089974e-01,-4.302991694091006926 e-01,
0.000000000000000000e+00, 4.305433885422950373e-01, 9.025698812737539884 e-01 0.000000000000000000e+00, 4.305433885422950373e-01, 9.025698812737539884 e-01
}; };
*/ */
 End of changes. 2 change blocks. 
2 lines changed or deleted 1 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/