Diff for /rpl/lapack/lapack/dgeqpf.f between versions 1.3 and 1.8

version 1.3, 2010/08/06 15:28:36 version 1.8, 2010/12/21 13:53:25
Line 1 Line 1
       SUBROUTINE DGEQPF( M, N, A, LDA, JPVT, TAU, WORK, INFO )        SUBROUTINE DGEQPF( M, N, A, LDA, JPVT, TAU, WORK, INFO )
 *  *
 *  -- LAPACK deprecated driver routine (version 3.2) --  *  -- LAPACK deprecated computational 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            INFO, LDA, M, N        INTEGER            INFO, LDA, M, N
Line 78 Line 78
 *  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 92 Line 92
       DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z        DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DGEQR2, DLARF, DLARFP, DORM2R, DSWAP, XERBLA        EXTERNAL           DGEQR2, DLARF, DLARFG, DORM2R, DSWAP, XERBLA
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT        INTRINSIC          ABS, MAX, MIN, SQRT
Line 182 Line 182
 *           Generate elementary reflector H(i)  *           Generate elementary reflector H(i)
 *  *
             IF( I.LT.M ) THEN              IF( I.LT.M ) THEN
                CALL DLARFP( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )                 CALL DLARFG( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
             ELSE              ELSE
                CALL DLARFP( 1, A( M, M ), A( M, M ), 1, TAU( M ) )                 CALL DLARFG( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
             END IF              END IF
 *  *
             IF( I.LT.N ) THEN              IF( I.LT.N ) THEN

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


CVSweb interface <joel.bertrand@systella.fr>