1 2 3 4 5 6 7 8
// complex_ops.rs use num_complex::Complex; use num_traits::Float; pub fn complex_exp<T: Float>(z: Complex<T>) -> Complex<T> { z.exp() } // TODO: complex matrix exponential, conjugate-pair folding (S4D-style)