db/d43/zladiv_8f_source.html

00001       COMPLEX*16     FUNCTION ZLADIV( X, Y )
00002 *
00003 *  -- LAPACK auxiliary routine (version 3.2) --
00004 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
00005 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
00006 *     November 2006
00007 *
00008 *     .. Scalar Arguments ..
00009       COMPLEX*16         X, Y
00010 *     ..
00011 *
00012 *  Purpose
00013 *  =======
00014 *
00015 *  ZLADIV := X / Y, where X and Y are complex.  The computation of X / Y
00016 *  will not overflow on an intermediary step unless the results
00017 *  overflows.
00018 *
00019 *  Arguments
00020 *  =========
00021 *
00022 *  X       (input) COMPLEX*16
00023 *  Y       (input) COMPLEX*16
00024 *          The complex scalars X and Y.
00025 *
00026 *  =====================================================================
00027 *
00028 *     .. Local Scalars ..
00029       DOUBLE PRECISION   ZI, ZR
00030 *     ..
00031 *     .. External Subroutines ..
00032       EXTERNAL           DLADIV
00033 *     ..
00034 *     .. Intrinsic Functions ..
00035       INTRINSIC          DBLE, DCMPLX, DIMAG
00036 *     ..
00037 *     .. Executable Statements ..
00038 *
00039       CALL DLADIV( DBLE( X ), DIMAG( X ), DBLE( Y ), DIMAG( Y ), ZR,
00040      $             ZI )
00041       ZLADIV = DCMPLX( ZR, ZI )
00042 *
00043       RETURN
00044 *
00045 *     End of ZLADIV
00046 *
00047       END