pub const NUMERICAL_STABILITY_EPSILON: f64 = 1e-10;
Epsilon for numerical stability in division operations