--- rpl/lapack/lapack/dlaqp2.f	2010/08/06 15:28:41	1.3
+++ rpl/lapack/lapack/dlaqp2.f	2011/07/22 07:38:07	1.9
@@ -1,10 +1,10 @@
       SUBROUTINE DLAQP2( M, N, OFFSET, A, LDA, JPVT, TAU, VN1, VN2,
      $                   WORK )
 *
-*  -- LAPACK auxiliary routine (version 3.2) --
+*  -- LAPACK auxiliary routine (version 3.3.1) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     November 2006
+*  -- April 2011                                                      --
 *
 *     .. Scalar Arguments ..
       INTEGER            LDA, M, N, OFFSET
@@ -75,7 +75,7 @@
 *  Partial column norm updating strategy modified by
 *    Z. Drmac and Z. Bujanovic, Dept. of Mathematics,
 *    University of Zagreb, Croatia.
-*    June 2006.
+*  -- April 2011                                                      --
 *  For more details see LAPACK Working Note 176.
 *  =====================================================================
 *
@@ -88,7 +88,7 @@
       DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, DSWAP
+      EXTERNAL           DLARF, DLARFG, DSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -125,15 +125,15 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL DLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL DLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL DLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL DLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LE.N ) THEN
 *
-*           Apply H(i)' to A(offset+i:m,i+1:n) from the left.
+*           Apply H(i)**T to A(offset+i:m,i+1:n) from the left.
 *
             AII = A( OFFPI, I )
             A( OFFPI, I ) = ONE