Diff for /rpl/lapack/lapack/zlaqps.f between versions 1.4 and 1.5

version 1.4, 2010/08/06 15:32:44 version 1.5, 2010/08/07 13:18:09
Line 1 Line 1
       SUBROUTINE ZLAQPS( M, N, OFFSET, NB, KB, A, LDA, JPVT, TAU, VN1,        SUBROUTINE ZLAQPS( M, N, OFFSET, NB, KB, A, LDA, JPVT, TAU, VN1,
      $                   VN2, AUXV, F, LDF )       $                   VN2, AUXV, F, LDF )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary routine (version 3.2.2) --
 *  -- 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  *     June 2010
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            KB, LDA, LDF, M, N, NB, OFFSET        INTEGER            KB, LDA, LDF, M, N, NB, OFFSET
Line 103 Line 103
       COMPLEX*16         AKK        COMPLEX*16         AKK
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           ZGEMM, ZGEMV, ZLARFP, ZSWAP        EXTERNAL           ZGEMM, ZGEMV, ZLARFG, ZSWAP
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, DCONJG, MAX, MIN, NINT, SQRT        INTRINSIC          ABS, DBLE, DCONJG, MAX, MIN, NINT, SQRT
Line 157 Line 157
 *        Generate elementary reflector H(k).  *        Generate elementary reflector H(k).
 *  *
          IF( RK.LT.M ) THEN           IF( RK.LT.M ) THEN
             CALL ZLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )              CALL ZLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE           ELSE
             CALL ZLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )              CALL ZLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF           END IF
 *  *
          AKK = A( RK, K )           AKK = A( RK, K )

Removed from v.1.4  
changed lines
  Added in v.1.5


CVSweb interface <joel.bertrand@systella.fr>