01:       COMPLEX*16     FUNCTION ZLADIV( X, Y )
02: *
03: *  -- LAPACK auxiliary routine (version 3.2) --
04: *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
05: *     November 2006
06: *
07: *     .. Scalar Arguments ..
08:       COMPLEX*16         X, Y
09: *     ..
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
47: