Email from Tim Hopkins (University of Kent), March 13th (item #4)

Used A(1,1) and T(1,1) in call to *LARFG (around line 177)
to make arguments 2 and 5 scalars rather than 2d-arrays.

diff --git a/SRC/cgeqrt3.f b/SRC/cgeqrt3.f
index a5b55c1d32ae85abe9966112ecb140e51e368533..62b76e868dfd86edfcc28339bc8a9be621982ff0 100644 (file)
--- a/SRC/cgeqrt3.f
+++ b/SRC/cgeqrt3.f
@@ -177,7 +177,7 @@
 *
 *        Compute Householder transform when N=1
 *
-         CALL CLARFG( M, A, A( MIN( 2, M ), 1 ), 1, T )
+         CALL CLARFG( M, A(1,1), A( MIN( 2, M ), 1 ), 1, T(1,1) )
 *         
       ELSE
 *

diff --git a/SRC/dgeqrt3.f b/SRC/dgeqrt3.f
index c5f57a29f4183c5422c136aea5fff2a9ed4f5485..07b6d488d2eda4b5265e144760ff44f70711a7d2 100644 (file)
--- a/SRC/dgeqrt3.f
+++ b/SRC/dgeqrt3.f
@@ -177,7 +177,7 @@
 *
 *        Compute Householder transform when N=1
 *
-         CALL DLARFG( M, A, A( MIN( 2, M ), 1 ), 1, T )
+         CALL DLARFG( M, A(1,1), A( MIN( 2, M ), 1 ), 1, T(1,1) )
 *         
       ELSE
 *

diff --git a/SRC/sgeqrt3.f b/SRC/sgeqrt3.f
index 86a43f67de8f415e819a996c6bd5ec43902f5e44..872f8b402cd9df1aabf4e07e25958fd3894bf911 100644 (file)
--- a/SRC/sgeqrt3.f
+++ b/SRC/sgeqrt3.f
@@ -177,7 +177,7 @@
 *
 *        Compute Householder transform when N=1
 *
-         CALL SLARFG( M, A, A( MIN( 2, M ), 1 ), 1, T )
+         CALL SLARFG( M, A(1,1), A( MIN( 2, M ), 1 ), 1, T(1,1) )
 *         
       ELSE
 *

diff --git a/SRC/zgeqrt3.f b/SRC/zgeqrt3.f
index 8926b9980b8df86a600febf4459d38a8844d5cb1..f06f1c5305ab527597e97ea81d479bf0c4833b6f 100644 (file)
--- a/SRC/zgeqrt3.f
+++ b/SRC/zgeqrt3.f
@@ -177,7 +177,7 @@
 *
 *        Compute Householder transform when N=1
 *
-         CALL ZLARFG( M, A, A( MIN( 2, M ), 1 ), 1, T )
+         CALL ZLARFG( M, A(1,1), A( MIN( 2, M ), 1 ), 1, T(1,1) )
 *         
       ELSE
 *