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

Diff for /rpl/lapack/lapack/zsytrs_aa_2stage.f between versions 1.1 and 1.2

version 1.1, 2018/05/29 14:55:45 version 1.2, 2020/05/21 21:46:11
Line 36 Line 36
 *> \verbatim  *> \verbatim
 *>  *>
 *> ZSYTRS_AA_2STAGE solves a system of linear equations A*X = B with a complex  *> ZSYTRS_AA_2STAGE solves a system of linear equations A*X = B with a complex
 *> symmetric matrix A using the factorization A = U*T*U**T or  *> symmetric matrix A using the factorization A = U**T*T*U or
 *> A = L*T*L**T computed by ZSYTRF_AA_2STAGE.  *> A = L*T*L**T computed by ZSYTRF_AA_2STAGE.
 *> \endverbatim  *> \endverbatim
 *  *
Line 48 Line 48
 *>          UPLO is CHARACTER*1  *>          UPLO is CHARACTER*1
 *>          Specifies whether the details of the factorization are stored  *>          Specifies whether the details of the factorization are stored
 *>          as an upper or lower triangular matrix.  *>          as an upper or lower triangular matrix.
 *>          = 'U':  Upper triangular, form is A = U*T*U**T;  *>          = 'U':  Upper triangular, form is A = U**T*T*U;
 *>          = 'L':  Lower triangular, form is A = L*T*L**T.  *>          = 'L':  Lower triangular, form is A = L*T*L**T.
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 85 Line 85
 *>  *>
 *> \param[in] LTB  *> \param[in] LTB
 *> \verbatim  *> \verbatim
   *>          LTB is INTEGER
 *>          The size of the array TB. LTB >= 4*N.  *>          The size of the array TB. LTB >= 4*N.
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 207 Line 208
 *  *
       IF( UPPER ) THEN        IF( UPPER ) THEN
 *  *
 *        Solve A*X = B, where A = U*T*U**T.  *        Solve A*X = B, where A = U**T*T*U.
 *  *
          IF( N.GT.NB ) THEN           IF( N.GT.NB ) THEN
 *  *
 *           Pivot, P**T * B  *           Pivot, P**T * B -> B
 *  *
             CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, 1 )              CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, 1 )
 *  *
 *           Compute (U**T \P**T * B) -> B    [ (U**T \P**T * B) ]  *           Compute (U**T \ B) -> B    [ (U**T \P**T * B) ]
 *  *
             CALL ZTRSM( 'L', 'U', 'T', 'U', N-NB, NRHS, ONE, A(1, NB+1),              CALL ZTRSM( 'L', 'U', 'T', 'U', N-NB, NRHS, ONE, A(1, NB+1),
      $                 LDA, B(NB+1, 1), LDB)       $                 LDA, B(NB+1, 1), LDB)
Line 233 Line 234
             CALL ZTRSM( 'L', 'U', 'N', 'U', N-NB, NRHS, ONE, A(1, NB+1),              CALL ZTRSM( 'L', 'U', 'N', 'U', N-NB, NRHS, ONE, A(1, NB+1),
      $                  LDA, B(NB+1, 1), LDB)       $                  LDA, B(NB+1, 1), LDB)
 *  *
 *           Pivot, P * B  [ P * (U \ (T \ (U**T \P**T * B) )) ]  *           Pivot, P * B -> B  [ P * (U \ (T \ (U**T \P**T * B) )) ]
 *  *
             CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, -1 )              CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, -1 )
 *  *
Line 245 Line 246
 *  *
          IF( N.GT.NB ) THEN           IF( N.GT.NB ) THEN
 *  *
 *           Pivot, P**T * B  *           Pivot, P**T * B -> B
 *  *
             CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, 1 )              CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, 1 )
 *  *
 *           Compute (L \P**T * B) -> B    [ (L \P**T * B) ]  *           Compute (L \ B) -> B    [ (L \P**T * B) ]
 *  *
             CALL ZTRSM( 'L', 'L', 'N', 'U', N-NB, NRHS, ONE, A(NB+1, 1),              CALL ZTRSM( 'L', 'L', 'N', 'U', N-NB, NRHS, ONE, A(NB+1, 1),
      $                 LDA, B(NB+1, 1), LDB)       $                 LDA, B(NB+1, 1), LDB)
Line 267 Line 268
             CALL ZTRSM( 'L', 'L', 'T', 'U', N-NB, NRHS, ONE, A(NB+1, 1),              CALL ZTRSM( 'L', 'L', 'T', 'U', N-NB, NRHS, ONE, A(NB+1, 1),
      $                  LDA, B(NB+1, 1), LDB)       $                  LDA, B(NB+1, 1), LDB)
 *  *
 *           Pivot, P * B  [ P * (L**T \ (T \ (L \P**T * B) )) ]  *           Pivot, P * B -> B  [ P * (L**T \ (T \ (L \P**T * B) )) ]
 *  *
             CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, -1 )              CALL ZLASWP( NRHS, B, LDB, NB+1, N, IPIV, -1 )
 *  *
Removed from v.1.1  
changed lines
  Added in v.1.2

CVSweb interface <joel.bertrand@systella.fr>