--- rpl/lapack/lapack/zlaqp2.f	2010/08/06 15:28:57	1.3
+++ rpl/lapack/lapack/zlaqp2.f	2011/07/22 07:38:17	1.9
@@ -1,10 +1,10 @@
       SUBROUTINE ZLAQP2( 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.
 *  =====================================================================
 *
@@ -91,7 +91,7 @@
       COMPLEX*16         AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZLARF, ZLARFP, ZSWAP
+      EXTERNAL           ZLARF, ZLARFG, ZSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DCONJG, MAX, MIN, SQRT
@@ -128,15 +128,15 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL ZLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL ZLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL ZLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL ZLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN
 *
-*           Apply H(i)' to A(offset+i:m,i+1:n) from the left.
+*           Apply H(i)**H to A(offset+i:m,i+1:n) from the left.
 *
             AII = A( OFFPI, I )
             A( OFFPI, I ) = CONE