sidereon-core 1.3.3

Numerical astrodynamics propagation core plus the GNSS domain layer (SP3, broadcast ephemeris, multi-GNSS positioning, RTK/PPP, ionosphere/troposphere, DOP) behind a default-on gnss feature
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86

"""Binary64 operations rounded from a high-precision mpmath calculation.

Every helper returns a Python ``float`` immediately.  This is deliberate: the
broadcast reference recipe is a sequence of binary64 operations, not one
unrounded arbitrary-precision expression.
"""

from __future__ import annotations

import mpmath as mp


MPMATH_PRECISION_DIGITS = 200
mp.mp.dps = MPMATH_PRECISION_DIGITS


def _mp(value: float) -> mp.mpf:
    return mp.mpf(value)


def _f64(value: mp.mpf) -> float:
    return float(value)


def add(a: float, b: float) -> float:
    return _f64(_mp(a) + _mp(b))


def sub(a: float, b: float) -> float:
    return _f64(_mp(a) - _mp(b))


def mul(a: float, b: float) -> float:
    return _f64(_mp(a) * _mp(b))


def div(a: float, b: float) -> float:
    return _f64(_mp(a) / _mp(b))


def sqrt(a: float) -> float:
    return _f64(mp.sqrt(_mp(a)))


def sin(a: float) -> float:
    return _f64(mp.sin(_mp(a)))


def cos(a: float) -> float:
    return _f64(mp.cos(_mp(a)))


def tan(a: float) -> float:
    return _f64(mp.tan(_mp(a)))


def asin(a: float) -> float:
    return _f64(mp.asin(_mp(a)))


def acos(a: float) -> float:
    return _f64(mp.acos(_mp(a)))


def atan(a: float) -> float:
    return _f64(mp.atan(_mp(a)))


def atan2(y: float, x: float) -> float:
    return _f64(mp.atan2(_mp(y), _mp(x)))


def exp(a: float) -> float:
    return _f64(mp.exp(_mp(a)))


def ln(a: float) -> float:
    return _f64(mp.log(_mp(a)))


def log10(a: float) -> float:
    return _f64(mp.log10(_mp(a)))


def pow(a: float, b: float) -> float:
    return _f64(mp.power(_mp(a), _mp(b)))