Bug report: array out-of-bounds reference in DLAQR5

[mat_cross]

In the section of the code "Special case: 2-by-2 reflection (if needed)" the logical test

Code: Select all

IF( BMP22 .AND. ( V( 1, M22 ).NE.ZERO ) )

must be split into a nested test

Code: Select all

IF (BMP22) THEN, IF (V(1,M22).NE.ZERO) THEN...

If .NOT. BMP22 then M22 can exceed the second extent of V (and recall that logical expressions can be evaluated in any order in Fortran).

[Julien Langou]

Hello Mathew, Thanks for the bug report, this is fixed in the svn and will be corrected in the next release. Cheers, Julien.
(See http://www.netlib.org/lapack/Errata/ -- bug 0061.)