pclasizeheevx.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	pclasizeheevx (wknown, range, n, desca, vl, vu, il, iu, iseed, win, maxsize, vecsize, valsize)

Function/Subroutine Documentation

pclasizeheevx()

subroutine pclasizeheevx	(	logical	wknown,
		character	range,
		integer	n,
		integer, dimension( * )	desca,
		real	vl,
		real	vu,
		integer	il,
		integer	iu,
		integer, dimension( 4 )	iseed,
		real, dimension( * )	win,
		integer	maxsize,
		integer	vecsize,
		integer	valsize )

Definition at line 3 of file pclasizeheevx.f.

*
*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      LOGICAL            WKNOWN
      CHARACTER          RANGE
      INTEGER            IL, IU, MAXSIZE, N, VALSIZE, VECSIZE
      REAL               VL, VU
*     ..
*     .. Array Arguments ..
      INTEGER            DESCA( * ), ISEED( 4 )
      REAL               WIN( * )
*     ..
*
*  Purpose
*  =======
*
*  PCLASIZEHEEVX computes the amount of memory needed by PCHEEVX
*  to ensure:
*    1)  Orthogonal Eigenvectors
*    2)  Eigenvectors
*    3)  Eigenvalues
*
*  Arguments
*  =========
*
*  WKNOWN  (global input) INTEGER
*          .FALSE.:  WIN does not contain the eigenvalues
*          .TRUE.:   WIN does contain the eigenvalues
*
*  RANGE   (global input) CHARACTER*1
*          = 'A': all eigenvalues will be found.
*          = 'V': all eigenvalues in the interval [VL,VU]
*                 will be found.
*          = 'I': the IL-th through IU-th eigenvalues will be found.
*
*  N       (global input) INTEGER
*          Size of the matrix to be tested.  (global size)
*
*  DESCA   (global input) INTEGER array dimension ( DLEN_ )
*
*  VL      (global input/output ) REAL
*          If RANGE='V', the lower bound of the interval to be searched
*          for eigenvalues.  Not referenced if RANGE = 'A' or 'I'.
*          If VL > VU, RANGE='V' and WKNOWN = .TRUE., VL is set
*          to a random value near an entry in WIN
*
*  VU      (global input/output ) REAL
*          If RANGE='V', the upper bound of the interval to be searched
*          for eigenvalues.  Not referenced if RANGE = 'A' or 'I'.
*          If VL > VU, RANGE='V' and WKNOWN = .TRUE., VU is set
*          to a random value near an entry in WIN
*
*  IL      (global input/output ) INTEGER
*          If RANGE='I', the index (from smallest to largest) of the
*          smallest eigenvalue to be returned.  IL >= 1.
*          Not referenced if RANGE = 'A' or 'V'.
*          If IL < 0, RANGE='I' and WKNOWN = .TRUE., IL is set
*          to a random value from 1 to N
*
*  IU      (global input/output ) INTEGER
*          If RANGE='I', the index (from smallest to largest) of the
*          largest eigenvalue to be returned.  min(IL,N) <= IU <= N.
*          Not referenced if RANGE = 'A' or 'V'.
*          If IU < 0, RANGE='I' and WKNOWN = .TRUE., IU is set
*          to a random value from IL to N
*
*  ISEED   (global input/output) INTEGER array, dimension (4)
*          On entry, the seed of the random number generator; the array
*          elements must be between 0 and 4095, and ISEED(4) must be
*          odd.
*          On exit, the seed is updated.
*          ISEED is not touched unless IL, IU, VL or VU are modified.
*
*  WIN     (global input) REAL array, dimension (N)
*          If WKNOWN=1, WIN contains the eigenvalues of the matrix.
*
*  MAXSIZE (global output) INTEGER
*          Workspace required to guarantee that PCHEEVX will return
*          orthogonal eigenvectors.  IF WKNOWN=0, MAXSIZE is set to a
*          a value which guarantees orthogonality no matter what the
*          spectrum is.  If WKNOWN=1, MAXSIZE is set to a value which
*          guarantees orthogonality on a matrix with eigenvalues given
*          by WIN.
*
*  VECSIZE (global output) INTEGER
*          Workspace required to guarantee that PCHEEVX
*          will compute eigenvectors.
*
*  VALSIZE (global output) INTEGER
*          Workspace required to guarantee that PCHEEVX
*          will compute eigenvalues.
*
*
*     .. Parameters ..
      INTEGER            BLOCK_CYCLIC_2D, DLEN_, DTYPE_, CTXT_, M_, N_,
     $                   MB_, NB_, RSRC_, CSRC_, LLD_
      parameter( block_cyclic_2d = 1, dlen_ = 9, dtype_ = 1,
     $                   ctxt_ = 2, m_ = 3, n_ = 4, mb_ = 5, nb_ = 6,
     $                   rsrc_ = 7, csrc_ = 8, lld_ = 9 )
      REAL               TWENTY
      parameter( twenty = 20.0e0 )
*     ..
*     .. Local Scalars ..
*
      INTEGER            CLUSTERSIZE, I, ILMIN, IUMAX, MAXCLUSTERSIZE,
     $                   MQ0, MYCOL, MYIL, MYIU, MYROW, NB, NEIG, NN,
     $                   NP0, NPCOL, NPROW
      REAL               ANORM, EPS, ORFAC, SAFMIN, VLMIN, VUMAX
*     ..
*     .. External Functions ..
*
*
      LOGICAL            LSAME
      INTEGER            ICEIL, NUMROC
      REAL               PSLAMCH, SLARAN
      EXTERNAL           lsame, iceil, numroc, pslamch, slaran
*     ..
*     .. External Subroutines ..
      EXTERNAL           blacs_gridinfo
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          abs, int, max, real
*     ..
*     .. Executable Statements ..
*       This is just to keep ftnchek happy
      IF( block_cyclic_2d*csrc_*ctxt_*dlen_*dtype_*lld_*mb_*m_*nb_*n_*
     $    rsrc_.LT.0 )RETURN
*
      orfac = 1.0e-3
*
*
      CALL blacs_gridinfo( desca( ctxt_ ), nprow, npcol, myrow, mycol )
      eps = pslamch( desca( ctxt_ ), 'Precision' )
      safmin = pslamch( desca( ctxt_ ), 'Safe Minimum' )
      nb = desca( mb_ )
      nn = max( n, nb, 2 )
      np0 = numroc( nn, nb, 0, 0, nprow )
*
      valsize = 5*nn + 4*n
*
      IF( wknown ) THEN
         anorm = safmin / eps
         IF( n.GE.1 )
     $      anorm = max( abs( win( 1 ) ), abs( win( n ) ), anorm )
*
         IF( lsame( range, 'i' ) ) THEN
.LT.            IF( IL0 )
     $         IL = INT( SLARAN( ISEED )*REAL( N ) ) + 1
.LT.            IF( IU0 )
     $         IU = INT( SLARAN( ISEED )*REAL( N-IL ) ) + IL
.EQ.            IF( N0 )
     $         IU = 0
         ELSE IF( LSAME( RANGE, 'v' ) ) THEN
.GT.            IF( VLVU ) THEN
               MYIL = INT( SLARAN( ISEED )*REAL( N ) ) + 1
               MYIU = INT( SLARAN( ISEED )*REAL( N-MYIL ) ) + MYIL
               VL = WIN( MYIL ) + TWENTY*EPS*ABS( WIN( MYIL ) )
               VU = WIN( MYIU ) + TWENTY*EPS*ABS( WIN( MYIU ) )
               VU = MAX( VU, VL+EPS*TWENTY*ABS( VL )+SAFMIN )
            END IF
         END IF
*
      END IF
      IF( LSAME( RANGE, 'v' ) ) THEN
*
*     Compute ILMIN, IUMAX (based on VL, VU and WIN)
*
         IF( WKNOWN ) THEN
            VLMIN = VL - TWENTY*EPS*ANORM
            VUMAX = VU + TWENTY*EPS*ANORM
            ILMIN = 1
            IUMAX = 0
            DO 10 I = 1, N
.LT.               IF( WIN( I )VLMIN )
     $            ILMIN = ILMIN + 1
.LT.               IF( WIN( I )VUMAX )
     $            IUMAX = IUMAX + 1
   10       CONTINUE
         ELSE
            ILMIN = 1
            IUMAX = N
         END IF
      ELSE IF( LSAME( RANGE, 'i' ) ) THEN
         ILMIN = IL
         IUMAX = IU
      ELSE IF( LSAME( RANGE, 'a' ) ) THEN
         ILMIN = 1
         IUMAX = N
      END IF
*
      NEIG = IUMAX - ILMIN + 1
*
      MQ0 = NUMROC( MAX( NEIG, NB, 2 ), NB, 0, 0, NPCOL )
      VECSIZE = 4*N + MAX( 5*NN, NP0*MQ0 ) +
     $          ICEIL( NEIG, NPROW*NPCOL )*NN
*
      IF( WKNOWN ) THEN
         CLUSTERSIZE = 1
         MAXCLUSTERSIZE = 1
         DO 20 I = ILMIN + 1, IUMAX
.LT.            IF( ( WIN( I )-WIN( I-1 ) )ORFAC*2*ANORM ) THEN
               CLUSTERSIZE = CLUSTERSIZE + 1
.GT.               IF( CLUSTERSIZEMAXCLUSTERSIZE )
     $            MAXCLUSTERSIZE = CLUSTERSIZE
            ELSE
               CLUSTERSIZE = 1
            END IF
   20    CONTINUE
.GT.         IF( CLUSTERSIZEMAXCLUSTERSIZE )
     $      MAXCLUSTERSIZE = CLUSTERSIZE
      ELSE
         MAXCLUSTERSIZE = N
      END IF
*
      MAXSIZE = VECSIZE + MAX( ( MAXCLUSTERSIZE-1 ), 0 )*N
*
*
      RETURN
*
*     End of PCLASIZEHEEVX
*