mirror/astronomy
1
0
Fork 0
mirror of https://github.com/cosinekitty/astronomy.git synced 2026-05-19 14:27:52 -04:00
Code Issues Packages Projects Releases Wiki Activity

riseset 17.8 seconds: eliminate more 0 terms.

Browse Source
This commit is contained in:
Don Cross
2019-07-04 21:59:47 -04:00
parent eb7b1c18fc
commit 2a6e4fa8c9
3 changed files with 207 additions and 207 deletions
Download Patch File Download Diff File Expand all files Collapse all files

198
source/python/astronomy.py
View File

@@ -353,7 +353,7 @@ class _iau2000b:
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (711159.0 + 73.0*t)*sarg + -872.0*carg
de += (-6750.0 + 0.0*t)*carg + 358.0*sarg
de += (-6750.0)*carg + 358.0*sarg
arg = math.fmod(2.0*f + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -383,17 +383,17 @@ class _iau2000b:
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (156994.0 + 10.0*t)*sarg + -168.0*carg
de += (-1235.0 + 0.0*t)*carg + 82.0*sarg
de += (-1235.0)*carg + 82.0*sarg
arg = math.fmod(el + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (63110.0 + 63.0*t)*sarg + 27.0*carg
de += (-33228.0 + 0.0*t)*carg + -9.0*sarg
de += (-33228.0)*carg + -9.0*sarg
arg = math.fmod(-el + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-57976.0 + -63.0*t)*sarg + -189.0*carg
de += (31429.0 + 0.0*t)*carg + -75.0*sarg
de += (31429.0)*carg + -75.0*sarg
arg = math.fmod(-el + 2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -403,7 +403,7 @@ class _iau2000b:
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-51613.0 + -42.0*t)*sarg + 129.0*carg
de += (26366.0 + 0.0*t)*carg + 78.0*sarg
de += (26366.0)*carg + 78.0*sarg
arg = math.fmod(-2.0*el + 2.0*f + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -413,7 +413,7 @@ class _iau2000b:
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (63384.0 + 11.0*t)*sarg + -150.0*carg
de += (-1220.0 + 0.0*t)*carg + 29.0*sarg
de += (-1220.0)*carg + 29.0*sarg
arg = math.fmod(2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -422,13 +422,13 @@ class _iau2000b:
arg = math.fmod(-2.0*elp + 2.0*f - 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (32481.0 + 0.0*t)*sarg + 0.0*carg
de += (-13870.0 + 0.0*t)*carg + 0.0*sarg
dp += (32481.0)*sarg + 0.0*carg
de += (-13870.0)*carg + 0.0*sarg
arg = math.fmod(-2.0*el + 2.0*d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-47722.0 + 0.0*t)*sarg + -18.0*carg
de += (477.0 + 0.0*t)*carg + -25.0*sarg
dp += (-47722.0)*sarg + -18.0*carg
de += (477.0)*carg + -25.0*sarg
arg = math.fmod(2.0*el + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -437,23 +437,23 @@ class _iau2000b:
arg = math.fmod(el + 2.0*f - 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (28593.0 + 0.0*t)*sarg + -1.0*carg
dp += (28593.0)*sarg + -1.0*carg
de += (-12338.0 + 10.0*t)*carg + -3.0*sarg
arg = math.fmod(-el + 2.0*f + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (20441.0 + 21.0*t)*sarg + 10.0*carg
de += (-10758.0 + 0.0*t)*carg + -3.0*sarg
de += (-10758.0)*carg + -3.0*sarg
arg = math.fmod(2.0*el, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (29243.0 + 0.0*t)*sarg + -74.0*carg
de += (-609.0 + 0.0*t)*carg + 13.0*sarg
dp += (29243.0)*sarg + -74.0*carg
de += (-609.0)*carg + 13.0*sarg
arg = math.fmod(2.0*f, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (25887.0 + 0.0*t)*sarg + -66.0*carg
de += (-550.0 + 0.0*t)*carg + 11.0*sarg
dp += (25887.0)*sarg + -66.0*carg
de += (-550.0)*carg + 11.0*sarg
arg = math.fmod(elp + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -463,7 +463,7 @@ class _iau2000b:
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (15164.0 + 10.0*t)*sarg + 11.0*carg
de += (-8001.0 + 0.0*t)*carg + -1.0*sarg
de += (-8001.0)*carg + -1.0*sarg
arg = math.fmod(2.0*elp + 2.0*f - 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -472,23 +472,23 @@ class _iau2000b:
arg = math.fmod(-2.0*f + 2.0*d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (21783.0 + 0.0*t)*sarg + 13.0*carg
de += (-167.0 + 0.0*t)*carg + 13.0*sarg
dp += (21783.0)*sarg + 13.0*carg
de += (-167.0)*carg + 13.0*sarg
arg = math.fmod(el - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-12873.0 + -10.0*t)*sarg + -37.0*carg
de += (6953.0 + 0.0*t)*carg + -14.0*sarg
de += (6953.0)*carg + -14.0*sarg
arg = math.fmod(-elp + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-12654.0 + 11.0*t)*sarg + 63.0*carg
de += (6415.0 + 0.0*t)*carg + 26.0*sarg
de += (6415.0)*carg + 26.0*sarg
arg = math.fmod(-el + 2.0*f + 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-10204.0 + 0.0*t)*sarg + 25.0*carg
de += (5222.0 + 0.0*t)*carg + 15.0*sarg
dp += (-10204.0)*sarg + 25.0*carg
de += (5222.0)*carg + 15.0*sarg
arg = math.fmod(2.0*elp, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
@@ -497,213 +497,213 @@ class _iau2000b:
arg = math.fmod(el + 2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-7691.0 + 0.0*t)*sarg + 44.0*carg
de += (3268.0 + 0.0*t)*carg + 19.0*sarg
dp += (-7691.0)*sarg + 44.0*carg
de += (3268.0)*carg + 19.0*sarg
arg = math.fmod(-2.0*el + 2.0*f, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-11024.0 + 0.0*t)*sarg + -14.0*carg
de += (104.0 + 0.0*t)*carg + 2.0*sarg
dp += (-11024.0)*sarg + -14.0*carg
de += (104.0)*carg + 2.0*sarg
arg = math.fmod(elp + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (7566.0 + -21.0*t)*sarg + -11.0*carg
de += (-3250.0 + 0.0*t)*carg + -5.0*sarg
de += (-3250.0)*carg + -5.0*sarg
arg = math.fmod(2.0*f + 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-6637.0 + -11.0*t)*sarg + 25.0*carg
de += (3353.0 + 0.0*t)*carg + 14.0*sarg
de += (3353.0)*carg + 14.0*sarg
arg = math.fmod(-elp + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-7141.0 + 21.0*t)*sarg + 8.0*carg
de += (3070.0 + 0.0*t)*carg + 4.0*sarg
de += (3070.0)*carg + 4.0*sarg
arg = math.fmod(2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-6302.0 + -11.0*t)*sarg + 2.0*carg
de += (3272.0 + 0.0*t)*carg + 4.0*sarg
de += (3272.0)*carg + 4.0*sarg
arg = math.fmod(el + 2.0*f - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (5800.0 + 10.0*t)*sarg + 2.0*carg
de += (-3045.0 + 0.0*t)*carg + -1.0*sarg
de += (-3045.0)*carg + -1.0*sarg
arg = math.fmod(2.0*el + 2.0*f - 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (6443.0 + 0.0*t)*sarg + -7.0*carg
de += (-2768.0 + 0.0*t)*carg + -4.0*sarg
dp += (6443.0)*sarg + -7.0*carg
de += (-2768.0)*carg + -4.0*sarg
arg = math.fmod(-2.0*el + 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-5774.0 + -11.0*t)*sarg + -15.0*carg
de += (3041.0 + 0.0*t)*carg + -5.0*sarg
de += (3041.0)*carg + -5.0*sarg
arg = math.fmod(2.0*el + 2.0*f + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-5350.0 + 0.0*t)*sarg + 21.0*carg
de += (2695.0 + 0.0*t)*carg + 12.0*sarg
dp += (-5350.0)*sarg + 21.0*carg
de += (2695.0)*carg + 12.0*sarg
arg = math.fmod(-elp + 2.0*f - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-4752.0 + -11.0*t)*sarg + -3.0*carg
de += (2719.0 + 0.0*t)*carg + -3.0*sarg
de += (2719.0)*carg + -3.0*sarg
arg = math.fmod(-2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-4940.0 + -11.0*t)*sarg + -21.0*carg
de += (2720.0 + 0.0*t)*carg + -9.0*sarg
de += (2720.0)*carg + -9.0*sarg
arg = math.fmod(-el - elp + 2.0*d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (7350.0 + 0.0*t)*sarg + -8.0*carg
de += (-51.0 + 0.0*t)*carg + 4.0*sarg
dp += (7350.0)*sarg + -8.0*carg
de += (-51.0)*carg + 4.0*sarg
arg = math.fmod(2.0*el - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (4065.0 + 0.0*t)*sarg + 6.0*carg
de += (-2206.0 + 0.0*t)*carg + 1.0*sarg
dp += (4065.0)*sarg + 6.0*carg
de += (-2206.0)*carg + 1.0*sarg
arg = math.fmod(el + 2.0*d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (6579.0 + 0.0*t)*sarg + -24.0*carg
de += (-199.0 + 0.0*t)*carg + 2.0*sarg
dp += (6579.0)*sarg + -24.0*carg
de += (-199.0)*carg + 2.0*sarg
arg = math.fmod(elp + 2.0*f - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (3579.0 + 0.0*t)*sarg + 5.0*carg
de += (-1900.0 + 0.0*t)*carg + 1.0*sarg
dp += (3579.0)*sarg + 5.0*carg
de += (-1900.0)*carg + 1.0*sarg
arg = math.fmod(el - elp, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (4725.0 + 0.0*t)*sarg + -6.0*carg
de += (-41.0 + 0.0*t)*carg + 3.0*sarg
dp += (4725.0)*sarg + -6.0*carg
de += (-41.0)*carg + 3.0*sarg
arg = math.fmod(-2.0*el + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-3075.0 + 0.0*t)*sarg + -2.0*carg
de += (1313.0 + 0.0*t)*carg + -1.0*sarg
dp += (-3075.0)*sarg + -2.0*carg
de += (1313.0)*carg + -1.0*sarg
arg = math.fmod(3.0*el + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-2904.0 + 0.0*t)*sarg + 15.0*carg
de += (1233.0 + 0.0*t)*carg + 7.0*sarg
dp += (-2904.0)*sarg + 15.0*carg
de += (1233.0)*carg + 7.0*sarg
arg = math.fmod(-elp + 2.0*d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (4348.0 + 0.0*t)*sarg + -10.0*carg
de += (-81.0 + 0.0*t)*carg + 2.0*sarg
dp += (4348.0)*sarg + -10.0*carg
de += (-81.0)*carg + 2.0*sarg
arg = math.fmod(el - elp + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-2878.0 + 0.0*t)*sarg + 8.0*carg
de += (1232.0 + 0.0*t)*carg + 4.0*sarg
dp += (-2878.0)*sarg + 8.0*carg
de += (1232.0)*carg + 4.0*sarg
arg = math.fmod(d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-4230.0 + 0.0*t)*sarg + 5.0*carg
de += (-20.0 + 0.0*t)*carg + -2.0*sarg
dp += (-4230.0)*sarg + 5.0*carg
de += (-20.0)*carg + -2.0*sarg
arg = math.fmod(-el - elp + 2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-2819.0 + 0.0*t)*sarg + 7.0*carg
de += (1207.0 + 0.0*t)*carg + 3.0*sarg
dp += (-2819.0)*sarg + 7.0*carg
de += (1207.0)*carg + 3.0*sarg
arg = math.fmod(-el + 2.0*f, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-4056.0 + 0.0*t)*sarg + 5.0*carg
de += (40.0 + 0.0*t)*carg + -2.0*sarg
dp += (-4056.0)*sarg + 5.0*carg
de += (40.0)*carg + -2.0*sarg
arg = math.fmod(-elp + 2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-2647.0 + 0.0*t)*sarg + 11.0*carg
de += (1129.0 + 0.0*t)*carg + 5.0*sarg
dp += (-2647.0)*sarg + 11.0*carg
de += (1129.0)*carg + 5.0*sarg
arg = math.fmod(-2.0*el + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-2294.0 + 0.0*t)*sarg + -10.0*carg
de += (1266.0 + 0.0*t)*carg + -4.0*sarg
dp += (-2294.0)*sarg + -10.0*carg
de += (1266.0)*carg + -4.0*sarg
arg = math.fmod(el + elp + 2.0*f + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (2481.0 + 0.0*t)*sarg + -7.0*carg
de += (-1062.0 + 0.0*t)*carg + -3.0*sarg
dp += (2481.0)*sarg + -7.0*carg
de += (-1062.0)*carg + -3.0*sarg
arg = math.fmod(2.0*el + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (2179.0 + 0.0*t)*sarg + -2.0*carg
de += (-1129.0 + 0.0*t)*carg + -2.0*sarg
dp += (2179.0)*sarg + -2.0*carg
de += (-1129.0)*carg + -2.0*sarg
arg = math.fmod(-el + elp + d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (3276.0 + 0.0*t)*sarg + 1.0*carg
de += (-9.0 + 0.0*t)*carg + 0.0*sarg
dp += (3276.0)*sarg + 1.0*carg
de += (-9.0)*carg + 0.0*sarg
arg = math.fmod(el + elp, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-3389.0 + 0.0*t)*sarg + 5.0*carg
de += (35.0 + 0.0*t)*carg + -2.0*sarg
dp += (-3389.0)*sarg + 5.0*carg
de += (35.0)*carg + -2.0*sarg
arg = math.fmod(el + 2.0*f, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (3339.0 + 0.0*t)*sarg + -13.0*carg
de += (-107.0 + 0.0*t)*carg + 1.0*sarg
dp += (3339.0)*sarg + -13.0*carg
de += (-107.0)*carg + 1.0*sarg
arg = math.fmod(-el + 2.0*f - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-1987.0 + 0.0*t)*sarg + -6.0*carg
de += (1073.0 + 0.0*t)*carg + -2.0*sarg
dp += (-1987.0)*sarg + -6.0*carg
de += (1073.0)*carg + -2.0*sarg
arg = math.fmod(el + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-1981.0 + 0.0*t)*sarg + 0.0*carg
de += (854.0 + 0.0*t)*carg + 0.0*sarg
dp += (-1981.0)*sarg + 0.0*carg
de += (854.0)*carg + 0.0*sarg
arg = math.fmod(-el + d, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (4026.0 + 0.0*t)*sarg + -353.0*carg
de += (-553.0 + 0.0*t)*carg + -139.0*sarg
dp += (4026.0)*sarg + -353.0*carg
de += (-553.0)*carg + -139.0*sarg
arg = math.fmod(2.0*f + d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (1660.0 + 0.0*t)*sarg + -5.0*carg
de += (-710.0 + 0.0*t)*carg + -2.0*sarg
dp += (1660.0)*sarg + -5.0*carg
de += (-710.0)*carg + -2.0*sarg
arg = math.fmod(-el + 2.0*f + 4.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-1521.0 + 0.0*t)*sarg + 9.0*carg
de += (647.0 + 0.0*t)*carg + 4.0*sarg
dp += (-1521.0)*sarg + 9.0*carg
de += (647.0)*carg + 4.0*sarg
arg = math.fmod(-el + elp + d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (1314.0 + 0.0*t)*sarg + 0.0*carg
de += (-700.0 + 0.0*t)*carg + 0.0*sarg
dp += (1314.0)*sarg + 0.0*carg
de += (-700.0)*carg + 0.0*sarg
arg = math.fmod(-2.0*elp + 2.0*f - 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-1283.0 + 0.0*t)*sarg + 0.0*carg
de += (672.0 + 0.0*t)*carg + 0.0*sarg
dp += (-1283.0)*sarg + 0.0*carg
de += (672.0)*carg + 0.0*sarg
arg = math.fmod(el + 2.0*f + 2.0*d + om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (-1331.0 + 0.0*t)*sarg + 8.0*carg
de += (663.0 + 0.0*t)*carg + 4.0*sarg
dp += (-1331.0)*sarg + 8.0*carg
de += (663.0)*carg + 4.0*sarg
arg = math.fmod(-2.0*el + 2.0*f + 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (1383.0 + 0.0*t)*sarg + -2.0*carg
de += (-594.0 + 0.0*t)*carg + -2.0*sarg
dp += (1383.0)*sarg + -2.0*carg
de += (-594.0)*carg + -2.0*sarg
arg = math.fmod(-el + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (1405.0 + 0.0*t)*sarg + 4.0*carg
de += (-610.0 + 0.0*t)*carg + 2.0*sarg
dp += (1405.0)*sarg + 4.0*carg
de += (-610.0)*carg + 2.0*sarg
arg = math.fmod(el + elp + 2.0*f - 2.0*d + 2.0*om, _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (1290.0 + 0.0*t)*sarg + 0.0*carg
de += (-556.0 + 0.0*t)*carg + 0.0*sarg
dp += (1290.0)*sarg + 0.0*carg
de += (-556.0)*carg + 0.0*sarg
self.dpsi = -0.000135 + (dp * 1.0e-7)
self.deps = +0.000388 + (de * 1.0e-7)