mirror/astronomy
1
0
Fork 0
mirror of https://github.com/cosinekitty/astronomy.git synced 2026-05-24 08:45:55 -04:00
Code Issues Packages Projects Releases Wiki Activity

Python problem might be in trig functions

Browse Source
This commit is contained in:
Don Cross
2024-05-28 17:25:02 -04:00
parent 5dbede1a4b
commit e75bf335c4
4 changed files with 21 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
source/python/README.md
View File

@@ -3567,12 +3567,3 @@ latitude, longitude, and elevation for that vector.
The geographic latitude, longitude, and elevation above sea level
that corresponds to the given equatorial vector.
---
<a name="getenv"></a>
### getenv(key, default=None)
Get an environment variable, return None if it doesn't exist.
The optional second argument can specify an alternate default.
key, default and the result are str.

10
source/python/astronomy/astronomy.py
View File

@@ -36,9 +36,11 @@ import datetime
import enum
import re
import abc
from os import getenv
import os
from typing import Any, List, Tuple, Optional, Union, Callable, Dict
_ProblemDebug:bool = ('1' == os.getenv('ASTRONOMY_ENGINE_PYTHON_PROBLEM'))
def _cbrt(x: float) -> float:
'''Returns the cube root of x.'''
y = (x ** (1.0 / 3.0)) if (x >= 0.0) else -((-x) ** (1.0 / 3.0))
@@ -3181,8 +3183,6 @@ def _VsopSphereToRect(lon: float, lat: float, rad: float) -> _TerseVector:
rad * math.sin(lat)
)
_ProblemDebug:bool = ('1' == getenv('ASTRONOMY_ENGINE_PYTHON_PROBLEM'))
def _CalcVsop(model: _vsop_model_t, time: Time) -> Vector:
t = time.tt / _DAYS_PER_MILLENNIUM
lon = _VsopFormula(model.lon, t, True)
@@ -3191,6 +3191,10 @@ def _CalcVsop(model: _vsop_model_t, time: Time) -> Vector:
eclip = _VsopSphereToRect(lon, lat, rad)
if _ProblemDebug:
print('_CalcVsop: lon={:0.16g}, lat={:0.16g}, rad={:0.16g}, eclip=({:0.16g}, {:0.16g}, {:0.16g})'.format(lon, lat, rad, eclip.x, eclip.y, eclip.z))
print('_CalcVsop: cos(lon) = {:0.16g}'.format(math.cos(lon)))
print('_CalcVsop: sin(lon) = {:0.16g}'.format(math.sin(lon)))
print('_CalcVsop: cos(lat) = {:0.16g}'.format(math.cos(lat)))
print('_CalcVsop: sin(lat) = {:0.16g}'.format(math.sin(lat)))
return _VsopRotate(eclip).ToAstroVector(time)
class _body_state_t: