Floating point complex division

Complex floating point division isn't quite trivial. The naive method may encounter overflow, underflow, or loss of precision due to the intermediate floating point results. A simple remedy[1] is:

Calculates (a + ib)/(c + id) and returns a two-element JS array with the real and imaginary components.

result is an optional Array or typed array into which the output is inserted.

var cdiv = require('complex-division');

var result = cdiv( 1, 2, 3, 4 );

// result = [ 0.44, 0.08 ]

[1] Press, William H. Numerical Recipes 3rd Edition: The Art of Scientific Computing. Cambridge University Press, 2007.

(c) 2015 Ricky Reusser. MIT License