VSOP was developed and is maintained (updated with the latest data) by the scientists at the Bureau des Longitudes in Paris.
VSOP2013,another version of VSOP, computed the positions of the planets directly at any moment, as well as their orbital elements with improved accuracy.
Original VSOP2013 Solution was write by FORTRAN 77 . It's too old to use.
This repo is not just programming language translation, it's refactoring of VSOP2013.
This thing is totally useless for myself. But I think someone might need this algorithm in the future.
VSOP2013 is much much slower than VSOP87. Maybe more accurate than 87?
Using multi-thread and precalculation technique to accelerate iteration speed.
It is much faster than origin algorithm.
GetPlanetPosition(EMB) on i9-12950HX, .NET 10:
| Calculator | Scalar | Vector256 |
|---|---|---|
| Calculator (double, full series) | 1,241 μs | 243 μs |
| CalculatorF (float, truncated) | 12.6 μs | 2.4 μs |
CalculatorF is for visualization and animation only; it trades accuracy for speed.
See the Gitbook for details.
- Use
VSOPResultclass to manage calculate results. - Use
VSOPTimeclass to manage time.
Easy to convert time by callingVSOPTime.UTC,VSOPTime.TAI,VSOPTime.TDB - Very high performance per solution using SIMD (
Vector128/256/512, auto detected) and multithread. - Useful Utility class. Convert Elliptic coordinates to cartesian and spherical
- Async Api included.
- precalculation on φ in terms, which gives 20%+ speed up.
- Optional term truncation:
new Calculator(SimdMode.Auto, truncation). CalculatorF: single precision calculator for visualization and animation, 100 to 500 times faster (low precision).- Validated against the official
VSOP2013.ctlreference values (max error 1e-10). - Use MemoryPack for binary serialize.
Initialization time becomes less than 10% of previous version. - Brotli compression on source data. ~300Mb -> ~43MB.
.NET 10 Runtime