logical function zslect

(

complex*16

Z

)

ZSLECT

Purpose:

 ZSLECT returns .TRUE. if the eigenvalue Z is to be selected,
 otherwise it returns .FALSE.
 It is used by ZCHK41 to test if ZGEES successfully sorts eigenvalues,
 and by ZCHK43 to test if ZGEESX successfully sorts eigenvalues.

 The common block /SSLCT/ controls how eigenvalues are selected.
 If SELOPT = 0, then ZSLECT return .TRUE. when real(Z) is less than
 zero, and .FALSE. otherwise.
 If SELOPT is at least 1, ZSLECT returns SELVAL(SELOPT) and adds 1
 to SELOPT, cycling back to 1 at SELMAX.

Parameters

[in]

Z

Z is COMPLEX*16 The eigenvalue Z.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Date

June 2016

Definition at line 58 of file zslect.f.

*
*  -- LAPACK test routine (version 3.6.1) --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     June 2016
*
*     .. Scalar Arguments ..
      COMPLEX*16         z
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      DOUBLE PRECISION   zero
      parameter                ( zero = 0.0d0 )
*     ..
*     .. Local Scalars ..
      INTEGER            i
      DOUBLE PRECISION   rmin, x
*     ..
*     .. Scalars in Common ..
      INTEGER            seldim, selopt
*     ..
*     .. Arrays in Common ..
      LOGICAL            selval( 20 )
      DOUBLE PRECISION   selwi( 20 ), selwr( 20 )
*     ..
*     .. Common blocks ..
      COMMON             / sslct / selopt, seldim, selval, selwr, selwi
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          abs, dble, dcmplx
*     ..
*     .. Executable Statements ..
*
      IF( selopt.EQ.0 ) THEN
         zslect = ( dble( z ).LT.zero )
      ELSE
         rmin = abs( z-dcmplx( selwr( 1 ), selwi( 1 ) ) )
         zslect = selval( 1 )
         DO 10 i = 2, seldim
            x = abs( z-dcmplx( selwr( i ), selwi( i ) ) )
            IF( x.LE.rmin ) THEN
               rmin = x
               zslect = selval( i )
            END IF
   10    CONTINUE
      END IF
      RETURN
*
*     End of ZSLECT
*