pub fn compinv<T: Number>(c: T, d: T) -> (T, T)
Expand description

Inversion of complex numbers that avoids overflows.

Arguments

  • c - Complex number real part
  • d - Complex number imaginary part

It is implemented using the division algorithm, replacing the numerator with 1+0i and simplify terms. Michael Baudin, Robert L. Smith, A Robust Complex Division in Scilab, 2012, arXiv:1210.4539v2 [cs.MS]