Function wright_omega 

pub fn wright_omega(x: f64) -> f64

Compute the Wright Omega function ω(x), the real branch of the solution to ω + ln(ω) = x.

Uses a piecewise initial approximation followed by Halley refinement steps, achieving near-double-precision accuracy across the full real domain.

Reference: D’Angelo & Fontana, “Efficient Computer Music Clipping Audio Effects Using SPICE-Like Lumped Models”, ICMC 2010.

Examples

// ω(0) ≈ 0.5671
assert!((wright_omega(0.0) - 0.5671).abs() < 1e-4);