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

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

version 1.7, 2010/12/21 13:53:33 version 1.8, 2011/07/22 07:38:08
Line 1 Line 1
       SUBROUTINE DLASYF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, INFO )        SUBROUTINE DLASYF( 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
 *  using the Bunch-Kaufman diagonal pivoting method. The partial  *  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.
Line 313 Line 313
 *  *
 *        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 536 Line 536
 *  *
 *        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>