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

version 1.4, 2010/08/06 15:32:28 version 1.5, 2010/08/07 13:18:07
Line 1 Line 1
       SUBROUTINE DLAQPS( M, N, OFFSET, NB, KB, A, LDA, JPVT, TAU, VN1,        SUBROUTINE DLAQPS( 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 91 Line 91
 *  Partial column norm updating strategy modified by  *  Partial column norm updating strategy modified by
 *    Z. Drmac and Z. Bujanovic, Dept. of Mathematics,  *    Z. Drmac and Z. Bujanovic, Dept. of Mathematics,
 *    University of Zagreb, Croatia.  *    University of Zagreb, Croatia.
 *    June 2006.  *     June 2010
 *  For more details see LAPACK Working Note 176.  *  For more details see LAPACK Working Note 176.
 *  =====================================================================  *  =====================================================================
 *  *
Line 104 Line 104
       DOUBLE PRECISION   AKK, TEMP, TEMP2, TOL3Z        DOUBLE PRECISION   AKK, TEMP, TEMP2, TOL3Z
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DGEMM, DGEMV, DLARFP, DSWAP        EXTERNAL           DGEMM, DGEMV, DLARFG, DSWAP
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, MAX, MIN, NINT, SQRT        INTRINSIC          ABS, DBLE, MAX, MIN, NINT, SQRT
Line 152 Line 152
 *        Generate elementary reflector H(k).  *        Generate elementary reflector H(k).
 *  *
          IF( RK.LT.M ) THEN           IF( RK.LT.M ) THEN
             CALL DLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )              CALL DLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE           ELSE
             CALL DLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )              CALL DLARFG( 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>