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

Diff for /rpl/lapack/lapack/zpbtf2.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:53 version 1.8, 2011/07/22 07:38:18
Line 1 Line 1
       SUBROUTINE ZPBTF2( UPLO, N, KD, AB, LDAB, INFO )        SUBROUTINE ZPBTF2( UPLO, N, KD, AB, LDAB, INFO )
 *  *
 *  -- LAPACK routine (version 3.2) --  *  -- LAPACK routine (version 3.3.1) --
 *  -- 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..--
 *     November 2006  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 20 Line 20
 *  positive definite band matrix A.  *  positive definite band matrix A.
 *  *
 *  The factorization has the form  *  The factorization has the form
 *     A = U' * U ,  if UPLO = 'U', or  *     A = U**H * U ,  if UPLO = 'U', or
 *     A = L  * L',  if UPLO = 'L',  *     A = L  * L**H,  if UPLO = 'L',
 *  where U is an upper triangular matrix, U' is the conjugate transpose  *  where U is an upper triangular matrix, U**H is the conjugate transpose
 *  of U, and L is lower triangular.  *  of U, and L is lower triangular.
 *  *
 *  This is the unblocked version of the algorithm, calling Level 2 BLAS.  *  This is the unblocked version of the algorithm, calling Level 2 BLAS.
Line 52 Line 52
 *          if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for j<=i<=min(n,j+kd).  *          if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for j<=i<=min(n,j+kd).
 *  *
 *          On exit, if INFO = 0, the triangular factor U or L from the  *          On exit, if INFO = 0, the triangular factor U or L from the
 *          Cholesky factorization A = U'*U or A = L*L' of the band  *          Cholesky factorization A = U**H *U or A = L*L**H of the band
 *          matrix A, in the same storage format as A.  *          matrix A, in the same storage format as A.
 *  *
 *  LDAB    (input) INTEGER  *  LDAB    (input) INTEGER
Line 137 Line 137
 *  *
       IF( UPPER ) THEN        IF( UPPER ) THEN
 *  *
 *        Compute the Cholesky factorization A = U'*U.  *        Compute the Cholesky factorization A = U**H * U.
 *  *
          DO 10 J = 1, N           DO 10 J = 1, N
 *  *
Line 165 Line 165
    10    CONTINUE     10    CONTINUE
       ELSE        ELSE
 *  *
 *        Compute the Cholesky factorization A = L*L'.  *        Compute the Cholesky factorization A = L*L**H.
 *  *
          DO 20 J = 1, N           DO 20 J = 1, N
 *  *
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>