Patch from Lawrence Mulholland (NAG) on Wednesday May 6th 2009.

See forum: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1408.
Thanks Lawrence.

diff --git a/SRC/chbgst.f b/SRC/chbgst.f
index 491600b..2d4b793 100644 (file)
--- a/SRC/chbgst.f
+++ b/SRC/chbgst.f
@@ -507,7 +507,7 @@
   230    CONTINUE
 *
          IF( KB.GT.1 ) THEN
-            DO 240 J = N - 1, I2 + KA, -1
+            DO 240 J = N - 1, J2 + KA, -1
                RWORK( J-M ) = RWORK( J-KA-M )
                WORK( J-M ) = WORK( J-KA-M )
   240       CONTINUE
@@ -770,7 +770,7 @@
   460    CONTINUE
 *
          IF( KB.GT.1 ) THEN
-            DO 470 J = N - 1, I2 + KA, -1
+            DO 470 J = N - 1, J2 + KA, -1
                RWORK( J-M ) = RWORK( J-KA-M )
                WORK( J-M ) = WORK( J-KA-M )
   470       CONTINUE

diff --git a/SRC/zhbgst.f b/SRC/zhbgst.f
index 01a295e..088a80d 100644 (file)
--- a/SRC/zhbgst.f
+++ b/SRC/zhbgst.f
@@ -507,7 +507,7 @@
   230    CONTINUE
 *
          IF( KB.GT.1 ) THEN
-            DO 240 J = N - 1, I2 + KA, -1
+            DO 240 J = N - 1, J2 + KA, -1
                RWORK( J-M ) = RWORK( J-KA-M )
                WORK( J-M ) = WORK( J-KA-M )
   240       CONTINUE
@@ -771,7 +771,7 @@
   460    CONTINUE
 *
          IF( KB.GT.1 ) THEN
-            DO 470 J = N - 1, I2 + KA, -1
+            DO 470 J = N - 1, J2 + KA, -1
                RWORK( J-M ) = RWORK( J-KA-M )
                WORK( J-M ) = WORK( J-KA-M )
   470       CONTINUE