Array out-of-bounds reference in xLAQR5.
authorlangou <langou@users.noreply.github.com>
Tue, 17 Aug 2010 15:45:24 +0000 (15:45 +0000)
committerlangou <langou@users.noreply.github.com>
Tue, 17 Aug 2010 15:45:24 +0000 (15:45 +0000)
commit150cd6266c7d887da544f7f23bfd7fd610044f66
tree7308322b129bfc0ade8fc01794ffdb308aaa6c22
parent4ff24babc1c98305f30b50c4f363fdef3ac7f446
Array out-of-bounds reference in xLAQR5.

Bug report and fix from Mathew Cross (NAG) on Sat Aug 14 2010:

In the section of the code "Special case: 2-by-2 reflection (if needed)" the logical test
IF( BMP22 .AND. ( V( 1, M22 ).NE.ZERO ) )
must be split into a nested test
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).

See http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1949.
SRC/claqr5.f
SRC/dlaqr5.f
SRC/slaqr5.f
SRC/zlaqr5.f