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

version 1.7, 2010/12/21 13:53:47 version 1.8, 2011/07/22 07:38:15
Line 1 Line 1
       SUBROUTINE ZHPGST( ITYPE, UPLO, N, AP, BP, INFO )        SUBROUTINE ZHPGST( ITYPE, UPLO, N, AP, BP, 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 108 Line 108
       IF( ITYPE.EQ.1 ) THEN        IF( ITYPE.EQ.1 ) THEN
          IF( UPPER ) THEN           IF( UPPER ) THEN
 *  *
 *           Compute inv(U')*A*inv(U)  *           Compute inv(U**H)*A*inv(U)
 *  *
 *           J1 and JJ are the indices of A(1,j) and A(j,j)  *           J1 and JJ are the indices of A(1,j) and A(j,j)
 *  *
Line 131 Line 131
    10       CONTINUE     10       CONTINUE
          ELSE           ELSE
 *  *
 *           Compute inv(L)*A*inv(L')  *           Compute inv(L)*A*inv(L**H)
 *  *
 *           KK and K1K1 are the indices of A(k,k) and A(k+1,k+1)  *           KK and K1K1 are the indices of A(k,k) and A(k+1,k+1)
 *  *
Line 161 Line 161
       ELSE        ELSE
          IF( UPPER ) THEN           IF( UPPER ) THEN
 *  *
 *           Compute U*A*U'  *           Compute U*A*U**H
 *  *
 *           K1 and KK are the indices of A(1,k) and A(k,k)  *           K1 and KK are the indices of A(1,k) and A(k,k)
 *  *
Line 186 Line 186
    30       CONTINUE     30       CONTINUE
          ELSE           ELSE
 *  *
 *           Compute L'*A*L  *           Compute L**H *A*L
 *  *
 *           JJ and J1J1 are the indices of A(j,j) and A(j+1,j+1)  *           JJ and J1J1 are the indices of A(j,j) and A(j+1,j+1)
 *  *

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


CVSweb interface <joel.bertrand@systella.fr>