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

Diff for /rpl/lapack/lapack/zsysv_aa.f between versions 1.4 and 1.5

version 1.4, 2018/05/29 07:18:36 version 1.5, 2020/05/21 21:46:10
Line 42 Line 42
 *> matrices.  *> matrices.
 *>  *>
 *> Aasen's algorithm is used to factor A as  *> Aasen's algorithm is used to factor A as
 *>    A = U * T * U**T,  if UPLO = 'U', or  *>    A = U**T * T * U,  if UPLO = 'U', or
 *>    A = L * T * L**T,  if UPLO = 'L',  *>    A = L * T * L**T,  if UPLO = 'L',
 *> where U (or L) is a product of permutation and unit upper (lower)  *> where U (or L) is a product of permutation and unit upper (lower)
 *> triangular matrices, and T is symmetric tridiagonal. The factored  *> triangular matrices, and T is symmetric tridiagonal. The factored
Line 86 Line 86
 *>  *>
 *>          On exit, if INFO = 0, the tridiagonal matrix T and the  *>          On exit, if INFO = 0, the tridiagonal matrix T and the
 *>          multipliers used to obtain the factor U or L from the  *>          multipliers used to obtain the factor U or L from the
 *>          factorization A = U*T*U**T or A = L*T*L**T as computed by  *>          factorization A = U**T*T*U or A = L*T*L**T as computed by
 *>          ZSYTRF.  *>          ZSYTRF.
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 221 Line 221
          LWKOPT_SYTRS = INT( WORK(1) )           LWKOPT_SYTRS = INT( WORK(1) )
          LWKOPT = MAX( LWKOPT_SYTRF, LWKOPT_SYTRS )           LWKOPT = MAX( LWKOPT_SYTRF, LWKOPT_SYTRS )
          WORK( 1 ) = LWKOPT           WORK( 1 ) = LWKOPT
          IF( LWORK.LT.LWKOPT .AND. .NOT.LQUERY ) THEN  
             INFO = -10  
          END IF  
       END IF        END IF
 *  *
       IF( INFO.NE.0 ) THEN        IF( INFO.NE.0 ) THEN
Line 233 Line 230
          RETURN           RETURN
       END IF        END IF
 *  *
 *     Compute the factorization A = U*T*U**T or A = L*T*L**T.  *     Compute the factorization A = U**T*T*U or A = L*T*L**T.
 *  *
       CALL ZSYTRF_AA( UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO )        CALL ZSYTRF_AA( UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO )
       IF( INFO.EQ.0 ) THEN        IF( INFO.EQ.0 ) THEN
Removed from v.1.4  
changed lines
  Added in v.1.5

CVSweb interface <joel.bertrand@systella.fr>