pub fn complex_div(a: Complex64, b: Complex64) -> Complex64
NumPy’s complex division (Smith’s scaled formula, including a division by zero), so quotients match NumPy to the last bit.