LAPACK  3.6.1 LAPACK: Linear Algebra PACKage
Go to the documentation of this file.
1 *> \brief \b ZLADIV performs complex division in real arithmetic, avoiding unnecessary overflow.
2 *
3 * =========== DOCUMENTATION ===========
4 *
5 * Online html documentation available at
6 * http://www.netlib.org/lapack/explore-html/
7 *
8 *> \htmlonly
11 *> [TGZ]</a>
13 *> [ZIP]</a>
15 *> [TXT]</a>
16 *> \endhtmlonly
17 *
18 * Definition:
19 * ===========
20 *
21 * COMPLEX*16 FUNCTION ZLADIV( X, Y )
22 *
23 * .. Scalar Arguments ..
24 * COMPLEX*16 X, Y
25 * ..
26 *
27 *
28 *> \par Purpose:
29 * =============
30 *>
31 *> \verbatim
32 *>
33 *> ZLADIV := X / Y, where X and Y are complex. The computation of X / Y
34 *> will not overflow on an intermediary step unless the results
35 *> overflows.
36 *> \endverbatim
37 *
38 * Arguments:
39 * ==========
40 *
41 *> \param[in] X
42 *> \verbatim
43 *> X is COMPLEX*16
44 *> \endverbatim
45 *>
46 *> \param[in] Y
47 *> \verbatim
48 *> Y is COMPLEX*16
49 *> The complex scalars X and Y.
50 *> \endverbatim
51 *
52 * Authors:
53 * ========
54 *
55 *> \author Univ. of Tennessee
56 *> \author Univ. of California Berkeley
57 *> \author Univ. of Colorado Denver
58 *> \author NAG Ltd.
59 *
60 *> \date September 2012
61 *
62 *> \ingroup complex16OTHERauxiliary
63 *
64 * =====================================================================
65  COMPLEX*16 FUNCTION zladiv( X, Y )
66 *
67 * -- LAPACK auxiliary routine (version 3.4.2) --
68 * -- LAPACK is a software package provided by Univ. of Tennessee, --
69 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
70 * September 2012
71 *
72 * .. Scalar Arguments ..
73  COMPLEX*16 X, Y
74 * ..
75 *
76 * =====================================================================
77 *
78 * .. Local Scalars ..
79  DOUBLE PRECISION ZI, ZR
80 * ..
81 * .. External Subroutines ..
83 * ..
84 * .. Intrinsic Functions ..
85  INTRINSIC dble, dcmplx, dimag
86 * ..
87 * .. Executable Statements ..
88 *
89  CALL dladiv( dble( x ), dimag( x ), dble( y ), dimag( y ), zr,
90  \$ zi )
91  zladiv = dcmplx( zr, zi )
92 *
93  RETURN
94 *