dladiv()

subroutine dladiv	(	double precision	A,
		double precision	B,
		double precision	C,
		double precision	D,
		double precision	P,
		double precision	Q
	)

DLADIV performs complex division in real arithmetic, avoiding unnecessary overflow.

Download DLADIV + dependencies [TGZ] [ZIP] [TXT]

Purpose:

 DLADIV performs complex division in  real arithmetic

                       a + i*b
            p + i*q = ---------
                       c + i*d

 The algorithm is due to Michael Baudin and Robert L. Smith
 and can be found in the paper
 "A Robust Complex Division in Scilab"

Parameters

[in]	A	A is DOUBLE PRECISION
[in]	B	B is DOUBLE PRECISION
[in]	C	C is DOUBLE PRECISION
[in]	D	D is DOUBLE PRECISION The scalars a, b, c, and d in the above expression.
[out]	P	P is DOUBLE PRECISION
[out]	Q	Q is DOUBLE PRECISION The scalars p and q in the above expression.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.