Diff for /rpl/lapack/lapack/zheevr_2stage.f between versions 1.2 and 1.6

version 1.2, 2017/06/17 11:06:47 version 1.6, 2023/08/07 08:39:23
Line 81 Line 81
 *> The desired accuracy of the output can be specified by the input  *> The desired accuracy of the output can be specified by the input
 *> parameter ABSTOL.  *> parameter ABSTOL.
 *>  *>
 *> For more details, see DSTEMR's documentation and:  *> For more details, see ZSTEMR's documentation and:
 *> - Inderjit S. Dhillon and Beresford N. Parlett: "Multiple representations  *> - Inderjit S. Dhillon and Beresford N. Parlett: "Multiple representations
 *>   to compute orthogonal eigenvectors of symmetric tridiagonal matrices,"  *>   to compute orthogonal eigenvectors of symmetric tridiagonal matrices,"
 *>   Linear Algebra and its Applications, 387(1), pp. 1-28, August 2004.  *>   Linear Algebra and its Applications, 387(1), pp. 1-28, August 2004.
Line 217 Line 217
 *>          eigenvalues are computed to high relative accuracy when  *>          eigenvalues are computed to high relative accuracy when
 *>          possible in future releases.  The current code does not  *>          possible in future releases.  The current code does not
 *>          make any guarantees about high relative accuracy, but  *>          make any guarantees about high relative accuracy, but
 *>          furutre releases will. See J. Barlow and J. Demmel,  *>          future releases will. See J. Barlow and J. Demmel,
 *>          "Computing Accurate Eigensystems of Scaled Diagonally  *>          "Computing Accurate Eigensystems of Scaled Diagonally
 *>          Dominant Matrices", LAPACK Working Note #7, for a discussion  *>          Dominant Matrices", LAPACK Working Note #7, for a discussion
 *>          of which matrices define their eigenvalues to high relative  *>          of which matrices define their eigenvalues to high relative
Line 354 Line 354
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date June 2016  
 *  
 *> \ingroup complex16HEeigen  *> \ingroup complex16HEeigen
 *  *
 *> \par Contributors:  *> \par Contributors:
Line 408 Line 406
 *  *
       IMPLICIT NONE        IMPLICIT NONE
 *  *
 *  -- LAPACK driver routine (version 3.7.0) --  *  -- LAPACK driver routine --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     June 2016  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOBZ, RANGE, UPLO        CHARACTER          JOBZ, RANGE, UPLO
Line 445 Line 442
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            LSAME        LOGICAL            LSAME
       INTEGER            ILAENV        INTEGER            ILAENV, ILAENV2STAGE
       DOUBLE PRECISION   DLAMCH, ZLANSY        DOUBLE PRECISION   DLAMCH, ZLANSY
       EXTERNAL           LSAME, ILAENV, DLAMCH, ZLANSY        EXTERNAL           LSAME, DLAMCH, ZLANSY, ILAENV, ILAENV2STAGE
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DCOPY, DSCAL, DSTEBZ, DSTERF, XERBLA, ZDSCAL,        EXTERNAL           DCOPY, DSCAL, DSTEBZ, DSTERF, XERBLA, ZDSCAL,
Line 471 Line 468
       LQUERY = ( ( LWORK.EQ.-1 ) .OR. ( LRWORK.EQ.-1 ) .OR.        LQUERY = ( ( LWORK.EQ.-1 ) .OR. ( LRWORK.EQ.-1 ) .OR.
      $         ( LIWORK.EQ.-1 ) )       $         ( LIWORK.EQ.-1 ) )
 *  *
       KD     = ILAENV( 17, 'DSYTRD_2STAGE', JOBZ, N, -1, -1, -1 )        KD     = ILAENV2STAGE( 1, 'ZHETRD_2STAGE', JOBZ, N, -1, -1, -1 )
       IB     = ILAENV( 18, 'DSYTRD_2STAGE', JOBZ, N, KD, -1, -1 )        IB     = ILAENV2STAGE( 2, 'ZHETRD_2STAGE', JOBZ, N, KD, -1, -1 )
       LHTRD  = ILAENV( 19, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 )        LHTRD  = ILAENV2STAGE( 3, 'ZHETRD_2STAGE', JOBZ, N, KD, IB, -1 )
       LWTRD  = ILAENV( 20, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 )        LWTRD  = ILAENV2STAGE( 4, 'ZHETRD_2STAGE', JOBZ, N, KD, IB, -1 )
       LWMIN  = N + LHTRD + LWTRD        LWMIN  = N + LHTRD + LWTRD
       LRWMIN = MAX( 1, 24*N )        LRWMIN = MAX( 1, 24*N )
       LIWMIN = MAX( 1, 10*N )        LIWMIN = MAX( 1, 10*N )

Removed from v.1.2  
changed lines
  Added in v.1.6


CVSweb interface <joel.bertrand@systella.fr>