From c3cf0a404a76b23097c708991167058c29ac199a Mon Sep 17 00:00:00 2001 From: james Date: Fri, 22 Mar 2013 20:54:57 +0000 Subject: [PATCH] corrected bug found by Tyrone Rees posted on email list: "I think there's a bug in the routine dlaqp2.f in version 3.4.2 of LAPACK (and possibly earlier). Line 220 of the source is IF( I.LE.N ) THEN whereas for the other versions the same line reads IF( I.LT.N ) THEN This difference appears to make the routine dgeqp3 fall over when asked to do a factorization of the last two columns of a 10x5 matrix. The NAG fortran compiler gives the error: Subscript out of range for assumed-size array A - Access to element 22 but actual argument has only 20 elements Program terminated by fatal error The single precision version runs without a hiccup. Best wishes, Tyrone" --- SRC/dlaqp2.f | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SRC/dlaqp2.f b/SRC/dlaqp2.f index 0c04436..4d41be0 100644 --- a/SRC/dlaqp2.f +++ b/SRC/dlaqp2.f @@ -217,7 +217,7 @@ CALL DLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) ) END IF * - IF( I.LE.N ) THEN + IF( I.LT.N ) THEN * * Apply H(i)**T to A(offset+i:m,i+1:n) from the left. * -- 2.7.4