[BACK]Return to zheevr.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/zheevr.f between versions 1.14 and 1.22

version 1.14, 2012/12/14 14:22:47 version 1.22, 2023/08/07 08:39:23
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *> \htmlonly  *> \htmlonly
 *> Download ZHEEVR + dependencies   *> Download ZHEEVR + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zheevr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zheevr.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zheevr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zheevr.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zheevr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zheevr.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 21 Line 21
 *       SUBROUTINE ZHEEVR( JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU,  *       SUBROUTINE ZHEEVR( JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU,
 *                          ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK,  *                          ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK,
 *                          RWORK, LRWORK, IWORK, LIWORK, INFO )  *                          RWORK, LRWORK, IWORK, LIWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          JOBZ, RANGE, UPLO  *       CHARACTER          JOBZ, RANGE, UPLO
 *       INTEGER            IL, INFO, IU, LDA, LDZ, LIWORK, LRWORK, LWORK,  *       INTEGER            IL, INFO, IU, LDA, LDZ, LIWORK, LRWORK, LWORK,
Line 33 Line 33
 *       DOUBLE PRECISION   RWORK( * ), W( * )  *       DOUBLE PRECISION   RWORK( * ), W( * )
 *       COMPLEX*16         A( LDA, * ), WORK( * ), Z( LDZ, * )  *       COMPLEX*16         A( LDA, * ), WORK( * ), Z( LDZ, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 75 Line 75
 *> 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 155 Line 155
 *> \param[in] VL  *> \param[in] VL
 *> \verbatim  *> \verbatim
 *>          VL is DOUBLE PRECISION  *>          VL is DOUBLE PRECISION
   *>          If RANGE='V', the lower bound of the interval to
   *>          be searched for eigenvalues. VL < VU.
   *>          Not referenced if RANGE = 'A' or 'I'.
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in] VU  *> \param[in] VU
 *> \verbatim  *> \verbatim
 *>          VU is DOUBLE PRECISION  *>          VU is DOUBLE PRECISION
 *>          If RANGE='V', the lower and upper bounds of the interval to  *>          If RANGE='V', the upper bound of the interval to
 *>          be searched for eigenvalues. VL < VU.  *>          be searched for eigenvalues. VL < VU.
 *>          Not referenced if RANGE = 'A' or 'I'.  *>          Not referenced if RANGE = 'A' or 'I'.
 *> \endverbatim  *> \endverbatim
Line 168 Line 171
 *> \param[in] IL  *> \param[in] IL
 *> \verbatim  *> \verbatim
 *>          IL is INTEGER  *>          IL is INTEGER
   *>          If RANGE='I', the index of the
   *>          smallest eigenvalue to be returned.
   *>          1 <= IL <= IU <= N, if N > 0; IL = 1 and IU = 0 if N = 0.
   *>          Not referenced if RANGE = 'A' or 'V'.
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in] IU  *> \param[in] IU
 *> \verbatim  *> \verbatim
 *>          IU is INTEGER  *>          IU is INTEGER
 *>          If RANGE='I', the indices (in ascending order) of the  *>          If RANGE='I', the index of the
 *>          smallest and largest eigenvalues to be returned.  *>          largest eigenvalue to be returned.
 *>          1 <= IL <= IU <= N, if N > 0; IL = 1 and IU = 0 if N = 0.  *>          1 <= IL <= IU <= N, if N > 0; IL = 1 and IU = 0 if N = 0.
 *>          Not referenced if RANGE = 'A' or 'V'.  *>          Not referenced if RANGE = 'A' or 'V'.
 *> \endverbatim  *> \endverbatim
Line 203 Line 210
 *>          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 250 Line 257
 *>          The support of the eigenvectors in Z, i.e., the indices  *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th eigenvector  *>          indicating the nonzero elements in Z. The i-th eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through  *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
 *>          ISUPPZ( 2*i ).  *>          ISUPPZ( 2*i ). This is an output of ZSTEMR (tridiagonal
   *>          matrix). The support of the eigenvectors of A is typically
   *>          1:N because of the unitary transformations applied by ZUNMTR.
 *>          Implemented only for RANGE = 'A' or 'I' and IU - IL = N - 1  *>          Implemented only for RANGE = 'A' or 'I' and IU - IL = N - 1
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 324 Line 333
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  
 *> \date September 2012  
 *  *
 *> \ingroup complex16HEeigen  *> \ingroup complex16HEeigen
 *  *
Line 348 Line 355
      $                   ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK,       $                   ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK,
      $                   RWORK, LRWORK, IWORK, LIWORK, INFO )       $                   RWORK, LRWORK, IWORK, LIWORK, INFO )
 *  *
 *  -- LAPACK driver routine (version 3.4.2) --  *  -- 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..--
 *     September 2012  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOBZ, RANGE, UPLO        CHARACTER          JOBZ, RANGE, UPLO
Line 617 Line 623
      $                   IWORK, LIWORK, INFO )       $                   IWORK, LIWORK, INFO )
 *  *
 *           Apply unitary matrix used in reduction to tridiagonal  *           Apply unitary matrix used in reduction to tridiagonal
 *           form to eigenvectors returned by ZSTEIN.  *           form to eigenvectors returned by ZSTEMR.
 *  *
             IF( WANTZ .AND. INFO.EQ.0 ) THEN              IF( WANTZ .AND. INFO.EQ.0 ) THEN
                INDWKN = INDWK                 INDWKN = INDWK
Removed from v.1.14  
changed lines
  Added in v.1.22

CVSweb interface <joel.bertrand@systella.fr>