Step 1 of xlarfp : back out the old changes (need to run testing before going on)

36 files changed, 82 insertions, 82 deletions

diff --git a/SRC/cgelq2.f b/SRC/cgelq2.f
index 529c6a3f..583f5eee 100644
--- a/SRC/cgelq2.f
+++ b/SRC/cgelq2.f
@@ -74,7 +74,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARF, CLARFP, XERBLA
+      EXTERNAL           CLACGV, CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -104,7 +104,7 @@
 *
          CALL CLACGV( N-I+1, A( I, I ), LDA )
          ALPHA = A( I, I )
-         CALL CLARFP( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
+         CALL CLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *
diff --git a/SRC/cgeql2.f b/SRC/cgeql2.f
index c0aea5b9..bf0bf8f0 100644
--- a/SRC/cgeql2.f
+++ b/SRC/cgeql2.f
@@ -76,7 +76,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, XERBLA
+      EXTERNAL           CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG, MAX, MIN
@@ -106,7 +106,7 @@
 *        A(1:m-k+i-1,n-k+i)
 *
          ALPHA = A( M-K+I, N-K+I )
-         CALL CLARFP( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
+         CALL CLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
 *
 *        Apply H(i)' to A(1:m-k+i,1:n-k+i-1) from the left
 *
diff --git a/SRC/cgeqpf.f b/SRC/cgeqpf.f
index d63127f6..8e85f319 100644
--- a/SRC/cgeqpf.f
+++ b/SRC/cgeqpf.f
@@ -96,7 +96,7 @@
       COMPLEX            AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEQR2, CLARF, CLARFP, CSWAP, CUNM2R, XERBLA
+      EXTERNAL           CGEQR2, CLARF, CLARFG, CSWAP, CUNM2R, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CMPLX, CONJG, MAX, MIN, SQRT
@@ -186,7 +186,7 @@
 *           Generate elementary reflector H(i)
 *
             AII = A( I, I )
-            CALL CLARFP( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
+            CALL CLARFG( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
      $                   TAU( I ) )
             A( I, I ) = AII
 *
diff --git a/SRC/cgeqr2.f b/SRC/cgeqr2.f
index 135a167c..7697baf2 100644
--- a/SRC/cgeqr2.f
+++ b/SRC/cgeqr2.f
@@ -74,7 +74,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, XERBLA
+      EXTERNAL           CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG, MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL CLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL CLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *
diff --git a/SRC/cgerq2.f b/SRC/cgerq2.f
index 50b61927..14f1ecbc 100644
--- a/SRC/cgerq2.f
+++ b/SRC/cgerq2.f
@@ -76,7 +76,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARF, CLARFP, XERBLA
+      EXTERNAL           CLACGV, CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -107,7 +107,7 @@
 *
          CALL CLACGV( N-K+I, A( M-K+I, 1 ), LDA )
          ALPHA = A( M-K+I, N-K+I )
-         CALL CLARFP( N-K+I, ALPHA, A( M-K+I, 1 ), LDA,
+         CALL CLARFG( N-K+I, ALPHA, A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right
diff --git a/SRC/claqp2.f b/SRC/claqp2.f
index b65f6157..9c7d99bf 100644
--- a/SRC/claqp2.f
+++ b/SRC/claqp2.f
@@ -91,7 +91,7 @@
       COMPLEX            AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, CSWAP
+      EXTERNAL           CLARF, CLARFG, CSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CONJG, MAX, MIN, SQRT
@@ -128,10 +128,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL CLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL CLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL CLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL CLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN
diff --git a/SRC/claqps.f b/SRC/claqps.f
index cefdf044..30daffbc 100644
--- a/SRC/claqps.f
+++ b/SRC/claqps.f
@@ -108,7 +108,7 @@
       COMPLEX            AKK
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEMM, CGEMV, CLARFP, CSWAP
+      EXTERNAL           CGEMM, CGEMV, CLARFG, CSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CONJG, MAX, MIN, NINT, REAL, SQRT
@@ -162,9 +162,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL CLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL CLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL CLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL CLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )
diff --git a/SRC/clatrz.f b/SRC/clatrz.f
index da5e42d4..68d22a12 100644
--- a/SRC/clatrz.f
+++ b/SRC/clatrz.f
@@ -91,7 +91,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARFP, CLARZ
+      EXTERNAL           CLACGV, CLARFG, CLARZ
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG
@@ -116,7 +116,7 @@
 *
          CALL CLACGV( L, A( I, N-L+1 ), LDA )
          ALPHA = CONJG( A( I, I ) )
-         CALL CLARFP( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL CLARFG( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
          TAU( I ) = CONJG( TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
diff --git a/SRC/ctzrqf.f b/SRC/ctzrqf.f
index fb2bc230..38a38a7a 100644
--- a/SRC/ctzrqf.f
+++ b/SRC/ctzrqf.f
@@ -98,7 +98,7 @@
       INTRINSIC          CONJG, MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CAXPY, CCOPY, CGEMV, CGERC, CLACGV, CLARFP,
+      EXTERNAL           CAXPY, CCOPY, CGEMV, CGERC, CLACGV, CLARFG,
      $                   XERBLA
 *     ..
 *     .. Executable Statements ..
@@ -136,7 +136,7 @@
             A( K, K ) = CONJG( A( K, K ) )
             CALL CLACGV( N-M, A( K, M1 ), LDA )
             ALPHA = A( K, K )
-            CALL CLARFP( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
+            CALL CLARFG( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
             A( K, K ) = ALPHA
             TAU( K ) = CONJG( TAU( K ) )
 *
diff --git a/SRC/dgelq2.f b/SRC/dgelq2.f
index 96ef32a8..0eb8ca8e 100644
--- a/SRC/dgelq2.f
+++ b/SRC/dgelq2.f
@@ -74,7 +74,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i,i+1:n)
 *
-         CALL DLARFP( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
+         CALL DLARFG( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *
diff --git a/SRC/dgeql2.f b/SRC/dgeql2.f
index a684249c..02756c94 100644
--- a/SRC/dgeql2.f
+++ b/SRC/dgeql2.f
@@ -76,7 +76,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(1:m-k+i-1,n-k+i)
 *
-         CALL DLARFP( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
+         CALL DLARFG( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i,1:n-k+i-1) from the left
diff --git a/SRC/dgeqpf.f b/SRC/dgeqpf.f
index 4f783506..2c86e5a7 100644
--- a/SRC/dgeqpf.f
+++ b/SRC/dgeqpf.f
@@ -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
diff --git a/SRC/dgeqr2.f b/SRC/dgeqr2.f
index 78c3aac9..a1b9f234 100644
--- a/SRC/dgeqr2.f
+++ b/SRC/dgeqr2.f
@@ -74,7 +74,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL DLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL DLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *
diff --git a/SRC/dgerq2.f b/SRC/dgerq2.f
index e6357336..4bf23406 100644
--- a/SRC/dgerq2.f
+++ b/SRC/dgerq2.f
@@ -76,7 +76,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(m-k+i,1:n-k+i-1)
 *
-         CALL DLARFP( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
+         CALL DLARFG( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right
diff --git a/SRC/dlaqp2.f b/SRC/dlaqp2.f
index df455d20..b579fb4c 100644
--- a/SRC/dlaqp2.f
+++ b/SRC/dlaqp2.f
@@ -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,10 +125,10 @@
 *        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
diff --git a/SRC/dlaqps.f b/SRC/dlaqps.f
index 9e00472c..ddb72223 100644
--- a/SRC/dlaqps.f
+++ b/SRC/dlaqps.f
@@ -104,7 +104,7 @@
       DOUBLE PRECISION   AKK, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEMM, DGEMV, DLARFP, DSWAP
+      EXTERNAL           DGEMM, DGEMV, DLARFG, DSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, MAX, MIN, NINT, SQRT
@@ -152,9 +152,9 @@
 *        Generate elementary reflector H(k).
 *
          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
-            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
 *
          AKK = A( RK, K )
diff --git a/SRC/dlatrz.f b/SRC/dlatrz.f
index ca749356..4d4d489e 100644
--- a/SRC/dlatrz.f
+++ b/SRC/dlatrz.f
@@ -90,7 +90,7 @@
       INTEGER            I
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARFP, DLARZ
+      EXTERNAL           DLARFG, DLARZ
 *     ..
 *     .. Executable Statements ..
 *
@@ -112,7 +112,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        [ A(i,i) A(i,n-l+1:n) ]
 *
-         CALL DLARFP( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL DLARFG( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
 *
diff --git a/SRC/dtzrqf.f b/SRC/dtzrqf.f
index 9777036f..c1b425b4 100644
--- a/SRC/dtzrqf.f
+++ b/SRC/dtzrqf.f
@@ -96,7 +96,7 @@
       INTRINSIC          MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DAXPY, DCOPY, DGEMV, DGER, DLARFP, XERBLA
+      EXTERNAL           DAXPY, DCOPY, DGEMV, DGER, DLARFG, XERBLA
 *     ..
 *     .. Executable Statements ..
 *
@@ -130,7 +130,7 @@
 *           Use a Householder reflection to zero the kth row of A.
 *           First set up the reflection.
 *
-            CALL DLARFP( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
+            CALL DLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
 *
             IF( ( TAU( K ).NE.ZERO ) .AND. ( K.GT.1 ) ) THEN
 *
diff --git a/SRC/sgelq2.f b/SRC/sgelq2.f
index 16db9bb9..b4519a84 100644
--- a/SRC/sgelq2.f
+++ b/SRC/sgelq2.f
@@ -74,7 +74,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i,i+1:n)
 *
-         CALL SLARFP( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
+         CALL SLARFG( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *
diff --git a/SRC/sgeql2.f b/SRC/sgeql2.f
index 6d3bb751..6e64d814 100644
--- a/SRC/sgeql2.f
+++ b/SRC/sgeql2.f
@@ -76,7 +76,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(1:m-k+i-1,n-k+i)
 *
-         CALL SLARFP( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
+         CALL SLARFG( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i,1:n-k+i-1) from the left
diff --git a/SRC/sgeqpf.f b/SRC/sgeqpf.f
index 17816307..b32f8c81 100644
--- a/SRC/sgeqpf.f
+++ b/SRC/sgeqpf.f
@@ -92,7 +92,7 @@
       REAL               AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEQR2, SLARF, SLARFP, SORM2R, SSWAP, XERBLA
+      EXTERNAL           SGEQR2, SLARF, SLARFG, SORM2R, SSWAP, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -182,9 +182,9 @@
 *           Generate elementary reflector H(i)
 *
             IF( I.LT.M ) THEN
-               CALL SLARFP( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
+               CALL SLARFG( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
             ELSE
-               CALL SLARFP( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
+               CALL SLARFG( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
             END IF
 *
             IF( I.LT.N ) THEN
diff --git a/SRC/sgeqr2.f b/SRC/sgeqr2.f
index 4353cd63..c70fb981 100644
--- a/SRC/sgeqr2.f
+++ b/SRC/sgeqr2.f
@@ -74,7 +74,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL SLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL SLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *
diff --git a/SRC/sgerq2.f b/SRC/sgerq2.f
index e5279094..b8f11be0 100644
--- a/SRC/sgerq2.f
+++ b/SRC/sgerq2.f
@@ -76,7 +76,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(m-k+i,1:n-k+i-1)
 *
-         CALL SLARFP( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
+         CALL SLARFG( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right
diff --git a/SRC/slaqp2.f b/SRC/slaqp2.f
index 2f51ecde..10d5790d 100644
--- a/SRC/slaqp2.f
+++ b/SRC/slaqp2.f
@@ -88,7 +88,7 @@
       REAL               AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, SSWAP
+      EXTERNAL           SLARF, SLARFG, SSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -125,10 +125,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL SLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL SLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL SLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL SLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN
diff --git a/SRC/slaqps.f b/SRC/slaqps.f
index 3ff36978..343e16d6 100644
--- a/SRC/slaqps.f
+++ b/SRC/slaqps.f
@@ -104,7 +104,7 @@
       REAL               AKK, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEMM, SGEMV, SLARFP, SSWAP
+      EXTERNAL           SGEMM, SGEMV, SLARFG, SSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, NINT, REAL, SQRT
@@ -152,9 +152,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL SLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL SLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL SLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL SLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )
diff --git a/SRC/slatrz.f b/SRC/slatrz.f
index 3b0893a2..1adefa79 100644
--- a/SRC/slatrz.f
+++ b/SRC/slatrz.f
@@ -90,7 +90,7 @@
       INTEGER            I
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARFP, SLARZ
+      EXTERNAL           SLARFG, SLARZ
 *     ..
 *     .. Executable Statements ..
 *
@@ -112,7 +112,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        [ A(i,i) A(i,n-l+1:n) ]
 *
-         CALL SLARFP( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL SLARFG( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
 *
diff --git a/SRC/stzrqf.f b/SRC/stzrqf.f
index 30b4e206..a64965be 100644
--- a/SRC/stzrqf.f
+++ b/SRC/stzrqf.f
@@ -96,7 +96,7 @@
       INTRINSIC          MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SAXPY, SCOPY, SGEMV, SGER, SLARFP, XERBLA
+      EXTERNAL           SAXPY, SCOPY, SGEMV, SGER, SLARFG, XERBLA
 *     ..
 *     .. Executable Statements ..
 *
@@ -130,7 +130,7 @@
 *           Use a Householder reflection to zero the kth row of A.
 *           First set up the reflection.
 *
-            CALL SLARFP( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
+            CALL SLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
 *
             IF( ( TAU( K ).NE.ZERO ) .AND. ( K.GT.1 ) ) THEN
 *
diff --git a/SRC/zgelq2.f b/SRC/zgelq2.f
index 04ad7df5..8b91ff6e 100644
--- a/SRC/zgelq2.f
+++ b/SRC/zgelq2.f
@@ -74,7 +74,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -104,7 +104,7 @@
 *
          CALL ZLACGV( N-I+1, A( I, I ), LDA )
          ALPHA = A( I, I )
-         CALL ZLARFP( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
+         CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *
diff --git a/SRC/zgeql2.f b/SRC/zgeql2.f
index b8fb8636..3ae6d725 100644
--- a/SRC/zgeql2.f
+++ b/SRC/zgeql2.f
@@ -76,7 +76,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG, MAX, MIN
@@ -106,7 +106,7 @@
 *        A(1:m-k+i-1,n-k+i)
 *
          ALPHA = A( M-K+I, N-K+I )
-         CALL ZLARFP( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
+         CALL ZLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
 *
 *        Apply H(i)' to A(1:m-k+i,1:n-k+i-1) from the left
 *
diff --git a/SRC/zgeqpf.f b/SRC/zgeqpf.f
index 24000890..51e93c3f 100644
--- a/SRC/zgeqpf.f
+++ b/SRC/zgeqpf.f
@@ -96,7 +96,7 @@
       COMPLEX*16         AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZGEQR2, ZLARF, ZLARFP, ZSWAP, ZUNM2R
+      EXTERNAL           XERBLA, ZGEQR2, ZLARF, ZLARFG, ZSWAP, ZUNM2R
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DCMPLX, DCONJG, MAX, MIN, SQRT
@@ -186,7 +186,7 @@
 *           Generate elementary reflector H(i)
 *
             AII = A( I, I )
-            CALL ZLARFP( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
+            CALL ZLARFG( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
      $                   TAU( I ) )
             A( I, I ) = AII
 *
diff --git a/SRC/zgeqr2.f b/SRC/zgeqr2.f
index 2924f8c9..76dcba4d 100644
--- a/SRC/zgeqr2.f
+++ b/SRC/zgeqr2.f
@@ -74,7 +74,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG, MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL ZLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL ZLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *
diff --git a/SRC/zgerq2.f b/SRC/zgerq2.f
index 7fcccc4e..989330bd 100644
--- a/SRC/zgerq2.f
+++ b/SRC/zgerq2.f
@@ -76,7 +76,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -107,7 +107,7 @@
 *
          CALL ZLACGV( N-K+I, A( M-K+I, 1 ), LDA )
          ALPHA = A( M-K+I, N-K+I )
-         CALL ZLARFP( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) )
+         CALL ZLARFG( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right
 *
diff --git a/SRC/zlaqp2.f b/SRC/zlaqp2.f
index 2354a001..85a6223d 100644
--- a/SRC/zlaqp2.f
+++ b/SRC/zlaqp2.f
@@ -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,10 +128,10 @@
 *        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
diff --git a/SRC/zlaqps.f b/SRC/zlaqps.f
index fdff5f34..571af2ff 100644
--- a/SRC/zlaqps.f
+++ b/SRC/zlaqps.f
@@ -103,7 +103,7 @@
       COMPLEX*16         AKK
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZGEMM, ZGEMV, ZLARFP, ZSWAP
+      EXTERNAL           ZGEMM, ZGEMV, ZLARFG, ZSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, DCONJG, MAX, MIN, NINT, SQRT
@@ -157,9 +157,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL ZLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL ZLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL ZLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL ZLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )
diff --git a/SRC/zlatrz.f b/SRC/zlatrz.f
index b2b837da..06f021f2 100644
--- a/SRC/zlatrz.f
+++ b/SRC/zlatrz.f
@@ -91,7 +91,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZLACGV, ZLARFP, ZLARZ
+      EXTERNAL           ZLACGV, ZLARFG, ZLARZ
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG
@@ -116,7 +116,7 @@
 *
          CALL ZLACGV( L, A( I, N-L+1 ), LDA )
          ALPHA = DCONJG( A( I, I ) )
-         CALL ZLARFP( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL ZLARFG( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
          TAU( I ) = DCONJG( TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
diff --git a/SRC/ztzrqf.f b/SRC/ztzrqf.f
index 4571ca6c..1744ac6a 100644
--- a/SRC/ztzrqf.f
+++ b/SRC/ztzrqf.f
@@ -99,7 +99,7 @@
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZAXPY, ZCOPY, ZGEMV, ZGERC, ZLACGV,
-     $                   ZLARFP
+     $                   ZLARFG
 *     ..
 *     .. Executable Statements ..
 *
@@ -136,7 +136,7 @@
             A( K, K ) = DCONJG( A( K, K ) )
             CALL ZLACGV( N-M, A( K, M1 ), LDA )
             ALPHA = A( K, K )
-            CALL ZLARFP( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
+            CALL ZLARFG( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
             A( K, K ) = ALPHA
             TAU( K ) = DCONJG( TAU( K ) )
 *