diff options
| author | julie <julielangou@users.noreply.github.com> | 2015-07-18 23:53:25 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2015-07-18 23:53:25 +0000 |
| commit | 8e26f833b884e2b3feed3b4ae02b6715ae2b5ab4 (patch) | |
| tree | 3787f04cbfe23776a840dc7df51f7a98121fca49 | |
| parent | a002ad4689cfe9647630a8eac71e24277d823938 (diff) | |
| download | lapack-8e26f833b884e2b3feed3b4ae02b6715ae2b5ab4.tar.gz lapack-8e26f833b884e2b3feed3b4ae02b6715ae2b5ab4.tar.bz2 lapack-8e26f833b884e2b3feed3b4ae02b6715ae2b5ab4.zip | |
Commit bug reported by Alexander Kobotov (INTEL) on LAPACK Forum
Link: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4762&p=11439#p11439
Recent commit r1533 introduced out of bound access which could lead to a crash.
It results scaling of a vector instead of a single number.
Fix is simple:
zhgeqz.f:669
- CALL ZSCAL( ILAST, SIGNBC, H( ILAST, ILAST ), 1 )
+ CALL ZSCAL( 1, SIGNBC, H( ILAST, ILAST ), 1 )
Similar for chgeqz on line 668
| -rw-r--r-- | SRC/chgeqz.f | 2 | ||||
| -rw-r--r-- | SRC/zhgeqz.f | 2 |
2 files changed, 2 insertions, 2 deletions
diff --git a/SRC/chgeqz.f b/SRC/chgeqz.f index 7d2e810d..3cfcaab2 100644 --- a/SRC/chgeqz.f +++ b/SRC/chgeqz.f @@ -665,7 +665,7 @@ CALL CSCAL( ILAST+1-IFRSTM, SIGNBC, H( IFRSTM, ILAST ), $ 1 ) ELSE - CALL CSCAL( ILAST, SIGNBC, H( ILAST, ILAST ), 1 ) + CALL CSCAL( 1, SIGNBC, H( ILAST, ILAST ), 1 ) END IF IF( ILZ ) $ CALL CSCAL( N, SIGNBC, Z( 1, ILAST ), 1 ) diff --git a/SRC/zhgeqz.f b/SRC/zhgeqz.f index 3df67545..2699e1ed 100644 --- a/SRC/zhgeqz.f +++ b/SRC/zhgeqz.f @@ -666,7 +666,7 @@ CALL ZSCAL( ILAST+1-IFRSTM, SIGNBC, H( IFRSTM, ILAST ), $ 1 ) ELSE - CALL ZSCAL( ILAST, SIGNBC, H( ILAST, ILAST ), 1 ) + CALL ZSCAL( 1, SIGNBC, H( ILAST, ILAST ), 1 ) END IF IF( ILZ ) $ CALL ZSCAL( N, SIGNBC, Z( 1, ILAST ), 1 ) |