celestial-core 0.1.1-alpha.1

Core types and utilities for the Celestial astronomy library
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

# celestial-core

Low-level astronomical calculations for coordinate transformations.

[![Crates.io](https://img.shields.io/crates/v/celestial-core)](https://crates.io/crates/celestial-core)
[![Documentation](https://docs.rs/celestial-core/badge.svg)](https://docs.rs/celestial-core)
[![License: MIT OR Apache-2.0](https://img.shields.io/crates/l/celestial-core)](https://github.com/gaker/celestial)

Pure Rust implementation of IAU 2000/2006 standards for celestial mechanics: rotation matrices, nutation/precession models, angle handling, and geodetic conversions. No runtime FFI.

## Installation

```toml
[dependencies]
celestial-core = "0.1"
```

## Modules

| Module        | Purpose                                                   |
|---------------|-----------------------------------------------------------|
| `angle`       | Angle types, parsing (HMS/DMS), normalization, validation |
| `matrix`      | 3×3 rotation matrices and 3D vectors                      |
| `nutation`    | IAU 2000A/2000B/2006A nutation models                     |
| `precession`  | IAU 2000/2006 precession (Fukushima-Williams angles)      |
| `cio`         | CIO-based GCRS↔CIRS transformations                       |
| `obliquity`   | Mean obliquity of the ecliptic (IAU 1980, 2006)           |
| `location`    | Observer geodetic coordinates, geocentric conversion      |
| `constants`   | Astronomical constants (J2000, WGS84, unit conversions)   |

## Example

```rust
use celestial_core::nutation::NutationIAU2006A;
use celestial_core::constants::J2000_JD;

// Compute nutation at J2000.0
let nutation = NutationIAU2006A::new().compute(J2000_JD, 0.0).unwrap();
println!("Δψ = {:.6}″", nutation.delta_psi * 206264.806); // radians to arcsec
println!("Δε = {:.6}″", nutation.delta_eps * 206264.806);
```

## Features

- **`serde`** — Enables serialization for `Angle` and other types

## Design Notes

- **Two-part Julian Dates**: Functions accept `(jd1, jd2)` to preserve precision. Typically `jd1 = 2451545.0` (J2000.0) and `jd2` is days from epoch.
- **Radians internally**: All angular computations use radians. The `Angle` type provides conversion methods for degrees/HMS/DMS display.
- **Stateless models**: Nutation and precession calculators have no internal state. Call `compute(jd1, jd2)` with any epoch.

## License

Licensed under either of:

- Apache License, Version 2.0
- MIT License

## Contributing

See the [repository](https://github.com/gaker/celestial) for contribution guidelines.