--- rpl/lapack/lapack/zhetrs2.f 2010/12/21 13:53:47 1.2 +++ rpl/lapack/lapack/zhetrs2.f 2011/07/22 07:38:15 1.3 @@ -1,10 +1,10 @@ SUBROUTINE ZHETRS2( UPLO, N, NRHS, A, LDA, IPIV, B, LDB, $ WORK, INFO ) * -* -- LAPACK PROTOTYPE routine (version 3.3.0) -- +* -- LAPACK PROTOTYPE routine (version 3.3.1) -- * * -- Written by Julie Langou of the Univ. of TN -- -* November 2010 +* -- April 2011 -- * * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- @@ -21,9 +21,9 @@ * Purpose * ======= * -* ZHETRS2 solves a system of linear equations A*X = B with a real -* Hermitian matrix A using the factorization A = U*D*U**T or -* A = L*D*L**T computed by ZSYTRF and converted by ZSYCONV. +* ZHETRS2 solves a system of linear equations A*X = B with a complex +* Hermitian matrix A using the factorization A = U*D*U**H or +* A = L*D*L**H computed by ZHETRF and converted by ZSYCONV. * * Arguments * ========= @@ -118,9 +118,9 @@ * IF( UPPER ) THEN * -* Solve A*X = B, where A = U*D*U'. +* Solve A*X = B, where A = U*D*U**H. * -* P' * B +* P**T * B K=N DO WHILE ( K .GE. 1 ) IF( IPIV( K ).GT.0 ) THEN @@ -140,11 +140,11 @@ END IF END DO * -* Compute (U \P' * B) -> B [ (U \P' * B) ] +* Compute (U \P**T * B) -> B [ (U \P**T * B) ] * - CALL ZTRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','U','N','U',N,NRHS,ONE,A,LDA,B,LDB) * -* Compute D \ B -> B [ D \ (U \P' * B) ] +* Compute D \ B -> B [ D \ (U \P**T * B) ] * I=N DO WHILE ( I .GE. 1 ) @@ -169,11 +169,11 @@ I = I - 1 END DO * -* Compute (U' \ B) -> B [ U' \ (D \ (U \P' * B) ) ] +* Compute (U**H \ B) -> B [ U**H \ (D \ (U \P**T * B) ) ] * - CALL ZTRSM('L','U','C','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','U','C','U',N,NRHS,ONE,A,LDA,B,LDB) * -* P * B [ P * (U' \ (D \ (U \P' * B) )) ] +* P * B [ P * (U**H \ (D \ (U \P**T * B) )) ] * K=1 DO WHILE ( K .LE. N ) @@ -196,9 +196,9 @@ * ELSE * -* Solve A*X = B, where A = L*D*L'. +* Solve A*X = B, where A = L*D*L**H. * -* P' * B +* P**T * B K=1 DO WHILE ( K .LE. N ) IF( IPIV( K ).GT.0 ) THEN @@ -218,11 +218,11 @@ ENDIF END DO * -* Compute (L \P' * B) -> B [ (L \P' * B) ] +* Compute (L \P**T * B) -> B [ (L \P**T * B) ] * - CALL ZTRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','L','N','U',N,NRHS,ONE,A,LDA,B,LDB) * -* Compute D \ B -> B [ D \ (L \P' * B) ] +* Compute D \ B -> B [ D \ (L \P**T * B) ] * I=1 DO WHILE ( I .LE. N ) @@ -245,11 +245,11 @@ I = I + 1 END DO * -* Compute (L' \ B) -> B [ L' \ (D \ (L \P' * B) ) ] +* Compute (L**H \ B) -> B [ L**H \ (D \ (L \P**T * B) ) ] * - CALL ZTRSM('L','L','C','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','L','C','U',N,NRHS,ONE,A,LDA,B,LDB) * -* P * B [ P * (L' \ (D \ (L \P' * B) )) ] +* P * B [ P * (L**H \ (D \ (L \P**T * B) )) ] * K=N DO WHILE ( K .GE. 1 )