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
!> 
!> 

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.