7034
|
1 COMPLEX*16 FUNCTION ZLADIV( X, Y ) |
2329
|
2 * |
7034
|
3 * -- LAPACK auxiliary routine (version 3.1) -- |
|
4 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. |
|
5 * November 2006 |
2329
|
6 * |
|
7 * .. Scalar Arguments .. |
|
8 COMPLEX*16 X, Y |
|
9 * .. |
|
10 * |
|
11 * Purpose |
|
12 * ======= |
|
13 * |
|
14 * ZLADIV := X / Y, where X and Y are complex. The computation of X / Y |
|
15 * will not overflow on an intermediary step unless the results |
|
16 * overflows. |
|
17 * |
|
18 * Arguments |
|
19 * ========= |
|
20 * |
|
21 * X (input) COMPLEX*16 |
|
22 * Y (input) COMPLEX*16 |
|
23 * The complex scalars X and Y. |
|
24 * |
|
25 * ===================================================================== |
|
26 * |
|
27 * .. Local Scalars .. |
|
28 DOUBLE PRECISION ZI, ZR |
|
29 * .. |
|
30 * .. External Subroutines .. |
|
31 EXTERNAL DLADIV |
|
32 * .. |
|
33 * .. Intrinsic Functions .. |
|
34 INTRINSIC DBLE, DCMPLX, DIMAG |
|
35 * .. |
|
36 * .. Executable Statements .. |
|
37 * |
|
38 CALL DLADIV( DBLE( X ), DIMAG( X ), DBLE( Y ), DIMAG( Y ), ZR, |
|
39 $ ZI ) |
|
40 ZLADIV = DCMPLX( ZR, ZI ) |
|
41 * |
|
42 RETURN |
|
43 * |
|
44 * End of ZLADIV |
|
45 * |
|
46 END |