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

version 1.7, 2010/12/21 13:53:52 version 1.8, 2011/07/22 07:38:18
Line 1 Line 1
       SUBROUTINE ZLASYF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, INFO )        SUBROUTINE ZLASYF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, 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 21 Line 21
 *  A using the Bunch-Kaufman diagonal pivoting method. The partial  *  A using the Bunch-Kaufman diagonal pivoting method. The partial
 *  factorization has the form:  *  factorization has the form:
 *  *
 *  A  =  ( I  U12 ) ( A11  0  ) (  I    0   )  if UPLO = 'U', or:  *  A  =  ( I  U12 ) ( A11  0  ) (  I       0    )  if UPLO = 'U', or:
 *        ( 0  U22 ) (  0   D  ) ( U12' U22' )  *        ( 0  U22 ) (  0   D  ) ( U12**T U22**T )
 *  *
 *  A  =  ( L11  0 ) ( D    0  ) ( L11' L21' )  if UPLO = 'L'  *  A  =  ( L11  0 ) ( D    0  ) ( L11**T L21**T )  if UPLO = 'L'
 *        ( L21  I ) ( 0   A22 ) (  0    I   )  *        ( L21  I ) ( 0   A22 ) (  0       I    )
 *  *
 *  where the order of D is at most NB. The actual order is returned in  *  where the order of D is at most NB. The actual order is returned in
 *  the argument KB, and is either NB or NB-1, or N if N <= NB.  *  the argument KB, and is either NB or NB-1, or N if N <= NB.
 *  Note that U' denotes the transpose of U.  *  Note that U**T denotes the transpose of U.
 *  *
 *  ZLASYF is an auxiliary routine called by ZSYTRF. It uses blocked code  *  ZLASYF is an auxiliary routine called by ZSYTRF. It uses blocked code
 *  (calling Level 3 BLAS) to update the submatrix A11 (if UPLO = 'U') or  *  (calling Level 3 BLAS) to update the submatrix A11 (if UPLO = 'U') or
Line 322 Line 322
 *  *
 *        Update the upper triangle of A11 (= A(1:k,1:k)) as  *        Update the upper triangle of A11 (= A(1:k,1:k)) as
 *  *
 *        A11 := A11 - U12*D*U12' = A11 - U12*W'  *        A11 := A11 - U12*D*U12**T = A11 - U12*W**T
 *  *
 *        computing blocks of NB columns at a time  *        computing blocks of NB columns at a time
 *  *
Line 546 Line 546
 *  *
 *        Update the lower triangle of A22 (= A(k:n,k:n)) as  *        Update the lower triangle of A22 (= A(k:n,k:n)) as
 *  *
 *        A22 := A22 - L21*D*L21' = A22 - L21*W'  *        A22 := A22 - L21*D*L21**T = A22 - L21*W**T
 *  *
 *        computing blocks of NB columns at a time  *        computing blocks of NB columns at a time
 *  *

Removed from v.1.7  
changed lines
  Added in v.1.8


CVSweb interface <joel.bertrand@systella.fr>