pub fn divide<R>(x: Complex<R>, y: Complex<R>) -> Complex<R>where
R: Native,Expand description
x / y for complex numbers, inline in the caller’s loop, by num_complex’s formula. Two real arguments divide as reals, so a real zero divisor gives ±∞ or NaN with an imaginary part of 0. Any other zero divisor gives NaN parts. |y|² overflows beyond about 1e154 at f64 and 1e19 at f32, and underflows below their inverses, where the quotient loses its digits. The decompositions’ Scalar::quotient uses Smith’s algorithm, which avoids that at the cost of speed in a vector loop.