g_math 0.4.2

Multi-domain fixed-point arithmetic with geometric extension: Lie groups, manifolds, ODE solvers, tensors, fiber bundles — zero-float, 0 ULP transcendentals
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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111

//! Transcendental Functions — Fixed-Point Implementations
//!
//! 18 transcendental functions across 6 engines, computed at tier N+1 for full storage-tier precision.
//!
//! Dedicated engines: `exp`, `ln`, `sqrt`, `sin_cos`, `atan`, `pow`.
//! FASC-composed: `tan`, `asin`, `acos`, `sinh`, `cosh`, `tanh`, `asinh`, `acosh`, `atanh`, `atan2`.
//!
//! | Profile    | Format    | Decimals | ULP |
//! |------------|-----------|----------|-----|
//! | Embedded   | Q64.64    | 19       | 0   |
//! | Balanced   | Q128.128  | 38       | 0   |
//! | Scientific | Q256.256  | 77       | 0   |
//!
//! Strategy: compute at tier N+1 (double precision), downscale with rounding.
//! Profile-aware compile-time table selection via `GMATH_PROFILE` env var.

pub mod exp_tier_n_plus_1;
pub mod ln_tier_n_plus_1;
pub mod pow_tier_n_plus_1;
pub mod sqrt_tier_n_plus_1;
pub mod sin_cos_tier_n_plus_1;
pub mod atan_tier_n_plus_1;
// Public API - Export profile-aware wrappers
// i128 wrappers exist for all profiles that define them (q64_64, q128_128, q256_256, q32_32)
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use exp_tier_n_plus_1::exp_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use exp_tier_n_plus_1::{exp_binary_i256, exp_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use ln_tier_n_plus_1::ln_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use ln_tier_n_plus_1::{ln_binary_i256, ln_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256"))]
pub use pow_tier_n_plus_1::{pow_binary_i128, pow_binary_i256, pow_binary_i512};
pub use pow_tier_n_plus_1::pow_integer_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sqrt_tier_n_plus_1::sqrt_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sqrt_tier_n_plus_1::{sqrt_binary_i256, sqrt_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::sin_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::{sin_binary_i256, sin_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::cos_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::{cos_binary_i256, cos_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::pi_half_i128;
#[cfg(table_format = "q128_128")]
pub use sin_cos_tier_n_plus_1::pi_half_i256;
#[cfg(table_format = "q256_256")]
pub use sin_cos_tier_n_plus_1::pi_half_i512;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::atan_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::{atan_binary_i256, atan_binary_i512};
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::atan2_binary_i128;
#[cfg(any(table_format = "q64_64", table_format = "q128_128", table_format = "q256_256", table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::{atan2_binary_i256, atan2_binary_i512};

// Export I1024 exp/ln/sqrt/sin/cos/atan/atan2 for tier N+1 computation (scientific profile)
#[cfg(table_format = "q256_256")]
pub use exp_tier_n_plus_1::exp_binary_i1024;
#[cfg(table_format = "q256_256")]
pub use ln_tier_n_plus_1::ln_binary_i1024;
#[cfg(table_format = "q256_256")]
pub use sqrt_tier_n_plus_1::sqrt_binary_i1024;
#[cfg(table_format = "q256_256")]
pub use sin_cos_tier_n_plus_1::{sin_binary_i1024, cos_binary_i1024, pi_half_i1024};
#[cfg(table_format = "q256_256")]
pub use atan_tier_n_plus_1::{atan_binary_i1024, atan2_binary_i1024};

// Signed Q512.512 multiply helper - used by all Q512.512 transcendental implementations
#[cfg(table_format = "q256_256")]
pub(crate) use ln_tier_n_plus_1::multiply_i1024_q512_512;

// Compute-tier dispatch functions for StackEvaluator tier N+1 integration
#[cfg(table_format = "q256_256")]
pub use sin_cos_tier_n_plus_1::{sin_compute_tier_i1024, cos_compute_tier_i1024};
#[cfg(table_format = "q128_128")]
pub use sin_cos_tier_n_plus_1::{sin_compute_tier_i512, cos_compute_tier_i512};
#[cfg(any(table_format = "q64_64", table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::{sin_compute_tier_i256, cos_compute_tier_i256};

#[cfg(table_format = "q256_256")]
pub use atan_tier_n_plus_1::{atan_compute_tier_i1024, atan2_compute_tier_i1024};
#[cfg(table_format = "q128_128")]
pub use atan_tier_n_plus_1::{atan_compute_tier_i512, atan2_compute_tier_i512};
#[cfg(any(table_format = "q64_64", table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::{atan_compute_tier_i256, atan2_compute_tier_i256};

// Export i64 wrappers for Q32.32 / Q16.16 profiles
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use exp_tier_n_plus_1::exp_binary_i64;
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use ln_tier_n_plus_1::ln_binary_i64;
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use sqrt_tier_n_plus_1::sqrt_binary_i64;
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::{sin_binary_i64, cos_binary_i64};
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use sin_cos_tier_n_plus_1::{sin_compute_tier_i64, cos_compute_tier_i64, sincos_compute_tier_i64};
#[cfg(any(table_format = "q32_32", table_format = "q16_16"))]
pub use atan_tier_n_plus_1::{atan_binary_i64, atan_compute_tier_i64, atan2_compute_tier_i64};

// ln decomposition functions are internal to ln_tier_n_plus_1.rs (not public API)

// traits.rs and exp.rs archived — FASC uses exp_tier_n_plus_1 directly