version 1.6, 2010/08/13 21:03:53
|
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 |
* |
* |