sladiv()

subroutine sladiv	(	real	a,
		real	b,
		real	c,
		real	d,
		real	p,
		real	q )

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

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

Purpose:

!>
!> SLADIV 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 REAL !>
[in]	B	!> B is REAL !>
[in]	C	!> C is REAL !>
[in]	D	!> D is REAL !> The scalars a, b, c, and d in the above expression. !>
[out]	P	!> P is REAL !>
[out]	Q	!> Q is REAL !> The scalars p and q in the above expression. !>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.