1
*> \brief \b DLADIV performs complex division in real arithmetic, avoiding unnecessary overflow.
3
* =========== DOCUMENTATION ===========
5
* Online html documentation available at
6
* http://www.netlib.org/lapack/explore-html/
9
*> Download DLADIV + dependencies
10
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dladiv.f">
12
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dladiv.f">
14
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dladiv.f">
21
* SUBROUTINE DLADIV( A, B, C, D, P, Q )
23
* .. Scalar Arguments ..
24
* DOUBLE PRECISION A, B, C, D, P, Q
33
*> DLADIV performs complex division in real arithmetic
36
*> p + i*q = ---------
39
*> The algorithm is due to Robert L. Smith and can be found
40
*> in D. Knuth, The art of Computer Programming, Vol.2, p.195
48
*> A is DOUBLE PRECISION
53
*> B is DOUBLE PRECISION
58
*> C is DOUBLE PRECISION
63
*> D is DOUBLE PRECISION
64
*> The scalars a, b, c, and d in the above expression.
69
*> P is DOUBLE PRECISION
74
*> Q is DOUBLE PRECISION
75
*> The scalars p and q in the above expression.
81
*> \author Univ. of Tennessee
82
*> \author Univ. of California Berkeley
83
*> \author Univ. of Colorado Denver
86
*> \date September 2012
88
*> \ingroup auxOTHERauxiliary
90
* =====================================================================
1
91
SUBROUTINE DLADIV( A, B, C, D, P, Q )
3
* -- LAPACK auxiliary routine (version 2.0) --
4
* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
5
* Courant Institute, Argonne National Lab, and Rice University
93
* -- LAPACK auxiliary routine (version 3.4.2) --
94
* -- LAPACK is a software package provided by Univ. of Tennessee, --
95
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
8
98
* .. Scalar Arguments ..
9
99
DOUBLE PRECISION A, B, C, D, P, Q
13
* $Id: dladiv.f 19697 2010-10-29 16:57:34Z d3y133 $
18
* DLADIV performs complex division in real arithmetic
24
* The algorithm is due to Robert L. Smith and can be found
25
* in D. Knuth, The art of Computer Programming, Vol.2, p.195
30
* A (input) DOUBLE PRECISION
31
* B (input) DOUBLE PRECISION
32
* C (input) DOUBLE PRECISION
33
* D (input) DOUBLE PRECISION
34
* The scalars a, b, c, and d in the above expression.
36
* P (output) DOUBLE PRECISION
37
* Q (output) DOUBLE PRECISION
38
* The scalars p and q in the above expression.
40
102
* =====================================================================
42
104
* .. Local Scalars ..