src/LAPACK/zladiv.f

   1       COMPLEX*16     FUNCTION ZLADIV( X, Y )
   2 *
   3 *  -- LAPACK auxiliary routine (version 3.1) --
   4 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
   5 *     November 2006
   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