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

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

version 1.7, 2010/12/21 13:53:49 version 1.8, 2011/07/22 07:38:17
Line 1 Line 1
       SUBROUTINE ZLAHEF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, INFO )        SUBROUTINE ZLAHEF( 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
 *  matrix A using the Bunch-Kaufman diagonal pivoting method. The  *  matrix A using the Bunch-Kaufman diagonal pivoting method. The
 *  partial factorization has the form:  *  partial 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**H U22**H )
 *  *
 *  A  =  ( L11  0 ) (  D   0  ) ( L11' L21' )  if UPLO = 'L'  *  A  =  ( L11  0 ) (  D   0  ) ( L11**H L21**H )  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 conjugate transpose of U.  *  Note that U**H denotes the conjugate transpose of U.
 *  *
 *  ZLAHEF is an auxiliary routine called by ZHETRF. It uses blocked code  *  ZLAHEF is an auxiliary routine called by ZHETRF. 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 344 Line 344
 *  *
 *        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**H = A11 - U12*W**H
 *  *
 *        computing blocks of NB columns at a time (note that conjg(W) is  *        computing blocks of NB columns at a time (note that conjg(W) is
 *        actually stored)  *        actually stored)
Line 593 Line 593
 *  *
 *        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**H = A22 - L21*W**H
 *  *
 *        computing blocks of NB columns at a time (note that conjg(W) is  *        computing blocks of NB columns at a time (note that conjg(W) is
 *        actually stored)  *        actually stored)
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>