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