mirror/astronomy
1
0
Fork 0
mirror of https://github.com/cosinekitty/astronomy.git synced 2026-04-26 01:17:05 -04:00
Code Issues Packages Projects Releases Wiki Activity

Now running riseset test in 22 seconds: better _iau2000 loop.

Browse Source
This commit is contained in:
Don Cross
2019-07-04 19:08:21 -04:00
parent 0909aeeddd
commit ba244cdb0d
3 changed files with 18 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
source/python/astronomy.py
View File

@@ -469,6 +469,8 @@ _cls_t = [
[ 1290.0, 0.0, 0.0, -556.0, 0.0, 0.0]
]
_iaudata = list(zip(_nals_t, _cls_t))
class _iau2000b:
def __init__(self, time):
t = time.tt / 36525
@@ -479,12 +481,12 @@ class _iau2000b:
om = math.fmod((450160.398036 - t * 6962890.5431), _ASEC360) * _ASEC2RAD
dp = 0
de = 0
for i in range(77):
arg = math.fmod((_nals_t[i][0]*el + _nals_t[i][1]*elp + _nals_t[i][2]*f + _nals_t[i][3]*d + _nals_t[i][4]*om), _PI2)
for n, c in _iaudata:
arg = math.fmod((n[0]*el + n[1]*elp + n[2]*f + n[3]*d + n[4]*om), _PI2)
sarg = math.sin(arg)
carg = math.cos(arg)
dp += (_cls_t[i][0] + _cls_t[i][1] * t)*sarg + _cls_t[i][2]*carg
de += (_cls_t[i][3] + _cls_t[i][4] * t)*carg + _cls_t[i][5]*sarg
dp += (c[0] + c[1] * t)*sarg + c[2]*carg
de += (c[3] + c[4] * t)*carg + c[5]*sarg
self.dpsi = -0.000135 + (dp * 1.0e-7)
self.deps = +0.000388 + (de * 1.0e-7)