--- rpl/lapack/lapack/dgeqpf.f	2010/08/06 15:32:23	1.4
+++ rpl/lapack/lapack/dgeqpf.f	2010/08/07 13:18:06	1.5
@@ -1,9 +1,9 @@
       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,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     November 2006
+*     June 2010
 *
 *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, M, N
@@ -78,7 +78,7 @@
 *  Partial column norm updating strategy modified by
 *    Z. Drmac and Z. Bujanovic, Dept. of Mathematics,
 *    University of Zagreb, Croatia.
-*    June 2006.
+*     June 2010
 *  For more details see LAPACK Working Note 176.
 *
 *  =====================================================================
@@ -92,7 +92,7 @@
       DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEQR2, DLARF, DLARFP, DORM2R, DSWAP, XERBLA
+      EXTERNAL           DGEQR2, DLARF, DLARFG, DORM2R, DSWAP, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -182,9 +182,9 @@
 *           Generate elementary reflector H(i)
 *
             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
-               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
 *
             IF( I.LT.N ) THEN