LAPACK  3.4.2 LAPACK: Linear Algebra PACKage
zlctes.f
Go to the documentation of this file.
1 *> \brief \b ZLCTES
2 *
3 * =========== DOCUMENTATION ===========
4 *
5 * Online html documentation available at
6 * http://www.netlib.org/lapack/explore-html/
7 *
8 * Definition:
9 * ===========
10 *
11 * LOGICAL FUNCTION ZLCTES( Z, D )
12 *
13 * .. Scalar Arguments ..
14 * COMPLEX*16 D, Z
15 * ..
16 *
17 *
18 *> \par Purpose:
19 * =============
20 *>
21 *> \verbatim
22 *>
23 *> ZLCTES returns .TRUE. if the eigenvalue Z/D is to be selected
24 *> (specifically, in this subroutine, if the real part of the
25 *> eigenvalue is negative), and otherwise it returns .FALSE..
26 *>
27 *> It is used by the test routine ZDRGES to test whether the driver
28 *> routine ZGGES succesfully sorts eigenvalues.
29 *> \endverbatim
30 *
31 * Arguments:
32 * ==========
33 *
34 *> \param[in] Z
35 *> \verbatim
36 *> Z is COMPLEX*16
37 *> The numerator part of a complex eigenvalue Z/D.
38 *> \endverbatim
39 *>
40 *> \param[in] D
41 *> \verbatim
42 *> D is COMPLEX*16
43 *> The denominator part of a complex eigenvalue Z/D.
44 *> \endverbatim
45 *
46 * Authors:
47 * ========
48 *
49 *> \author Univ. of Tennessee
50 *> \author Univ. of California Berkeley
51 *> \author Univ. of Colorado Denver
52 *> \author NAG Ltd.
53 *
54 *> \date November 2011
55 *
56 *> \ingroup complex16_eig
57 *
58 * =====================================================================
59  LOGICAL FUNCTION zlctes( Z, D )
60 *
61 * -- LAPACK test routine (version 3.4.0) --
62 * -- LAPACK is a software package provided by Univ. of Tennessee, --
63 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
64 * November 2011
65 *
66 * .. Scalar Arguments ..
67  COMPLEX*16 d, z
68 * ..
69 *
70 * =====================================================================
71 *
72 * .. Parameters ..
73 *
74  DOUBLE PRECISION zero, one
75  parameter( zero = 0.0d+0, one = 1.0d+0 )
76  COMPLEX*16 czero
77  parameter( czero = ( 0.0d+0, 0.0d+0 ) )
78 * ..
79 * .. Local Scalars ..
80  DOUBLE PRECISION zmax
81 * ..
82 * .. Intrinsic Functions ..
83  INTRINSIC abs, dble, dimag, max, sign
84 * ..
85 * .. Executable Statements ..
86 *
87  IF( d.EQ.czero ) THEN
88  zlctes = ( dble( z ).LT.zero )
89  ELSE
90  IF( dble( z ).EQ.zero .OR. dble( d ).EQ.zero ) THEN
91  zlctes = ( sign( one, dimag( z ) ).NE.
92  \$ sign( one, dimag( d ) ) )
93  ELSE IF( dimag( z ).EQ.zero .OR. dimag( d ).EQ.zero ) THEN
94  zlctes = ( sign( one, dble( z ) ).NE.
95  \$ sign( one, dble( d ) ) )
96  ELSE
97  zmax = max( abs( dble( z ) ), abs( dimag( z ) ) )
98  zlctes = ( ( dble( z ) / zmax )*dble( d )+
99  \$ ( dimag( z ) / zmax )*dimag( d ).LT.zero )
100  END IF
101  END IF
102 *
103  return
104 *
105 * End of ZLCTES
106 *
107  END