SRC/cladiv.f

   1 *> \brief \b CLADIV performs complex division in real arithmetic, avoiding unnecessary overflow.
   2 *
   3 *  =========== DOCUMENTATION ===========
   4 *
   5 * Online html documentation available at
   6 *            http://www.netlib.org/lapack/explore-html/
   7 *
   8 *> \htmlonly
   9 *> Download CLADIV + dependencies
  10 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/cladiv.f">
  11 *> [TGZ]</a>
  12 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/cladiv.f">
  13 *> [ZIP]</a>
  14 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/cladiv.f">
  15 *> [TXT]</a>
  16 *> \endhtmlonly
  17 *
  18 *  Definition:
  19 *  ===========
  20 *
  21 *       COMPLEX FUNCTION CLADIV( X, Y )
  22 *
  23 *       .. Scalar Arguments ..
  24 *       COMPLEX            X, Y
  25 *       ..
  26 *
  27 *
  28 *> \par Purpose:
  29 *  =============
  30 *>
  31 *> \verbatim
  32 *>
  33 *> CLADIV := X / Y, where X and Y are complex.  The computation of X / Y
  34 *> will not overflow on an intermediary step unless the results
  35 *> overflows.
  36 *> \endverbatim
  37 *
  38 *  Arguments:
  39 *  ==========
  40 *
  41 *> \param[in] X
  42 *> \verbatim
  43 *>          X is COMPLEX
  44 *> \endverbatim
  45 *>
  46 *> \param[in] Y
  47 *> \verbatim
  48 *>          Y is COMPLEX
  49 *>          The complex scalars X and Y.
  50 *> \endverbatim
  51 *
  52 *  Authors:
  53 *  ========
  54 *
  55 *> \author Univ. of Tennessee
  56 *> \author Univ. of California Berkeley
  57 *> \author Univ. of Colorado Denver
  58 *> \author NAG Ltd.
  59 *
  60 *> \date September 2012
  61 *
  62 *> \ingroup complexOTHERauxiliary
  63 *
  64 *  =====================================================================
  65       COMPLEX FUNCTION CLADIV( X, Y )
  66 *
  67 *  -- LAPACK auxiliary routine (version 3.4.2) --
  68 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
  69 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  70 *     September 2012
  71 *
  72 *     .. Scalar Arguments ..
  73       COMPLEX            X, Y
  74 *     ..
  75 *
  76 *  =====================================================================
  77 *
  78 *     .. Local Scalars ..
  79       REAL               ZI, ZR
  80 *     ..
  81 *     .. External Subroutines ..
  82       EXTERNAL           SLADIV
  83 *     ..
  84 *     .. Intrinsic Functions ..
  85       INTRINSIC          AIMAG, CMPLX, REAL
  86 *     ..
  87 *     .. Executable Statements ..
  88 *
  89       CALL SLADIV( REAL( X ), AIMAG( X ), REAL( Y ), AIMAG( Y ), ZR,
  90      $             ZI )
  91       CLADIV = CMPLX( ZR, ZI )
  92 *
  93       RETURN
  94 *
  95 *     End of CLADIV
  96 *
  97       END