APPLYING INTEL PATCHES sent to Julie on Feb 19th 2016 by Dima from INTEL (dmitry...

[PATCH 29/42] Fix lapacke_?gb_nancheck - misuse of leading dimension as matrix size

diff --git a/LAPACKE/utils/lapacke_cgb_nancheck.c b/LAPACKE/utils/lapacke_cgb_nancheck.c
index 9446b355adce193e15726a09f51cab5e1b7924a4..400e806f2ad0d698b77958644dadeb2e4bfa4d90 100644 (file)
--- a/LAPACKE/utils/lapacke_cgb_nancheck.c
+++ b/LAPACKE/utils/lapacke_cgb_nancheck.c
@@ -46,14 +46,14 @@ lapack_logical LAPACKE_cgb_nancheck( int matrix_layout, lapack_int m,
 
     if( matrix_layout == LAPACK_COL_MAJOR ) {
         for( j = 0; j < n; j++ ) {
-            for( i = MAX( ku-j, 0 ); i < MIN3( ldab, m+ku-j, kl+ku+1 );
+            for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 );
                  i++ ) {
                 if( LAPACK_CISNAN( ab[i+(size_t)j*ldab] ) )
                      return (lapack_logical) 1;
             }
         }
     } else if ( matrix_layout == LAPACK_ROW_MAJOR ) {
-        for( j = 0; j < MIN( n, ldab ); j++ ) {
+        for( j = 0; j < n; j++ ) {
             for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 ); i++ ) {
                 if( LAPACK_CISNAN( ab[(size_t)i*ldab+j] ) )
                      return (lapack_logical) 1;

diff --git a/LAPACKE/utils/lapacke_dgb_nancheck.c b/LAPACKE/utils/lapacke_dgb_nancheck.c
index 5fa13e75dcd2bc93e74a20c76acd82557cb946a7..e94d488d73d38a2c41f5b8f0835158a951fb33ed 100644 (file)
--- a/LAPACKE/utils/lapacke_dgb_nancheck.c
+++ b/LAPACKE/utils/lapacke_dgb_nancheck.c
@@ -46,14 +46,14 @@ lapack_logical LAPACKE_dgb_nancheck( int matrix_layout, lapack_int m,
 
     if( matrix_layout == LAPACK_COL_MAJOR ) {
         for( j = 0; j < n; j++ ) {
-            for( i = MAX( ku-j, 0 ); i < MIN3( ldab, m+ku-j, kl+ku+1 );
+            for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 );
                  i++ ) {
                 if( LAPACK_DISNAN( ab[i+(size_t)j*ldab] ) )
                      return (lapack_logical) 1;
             }
         }
     } else if ( matrix_layout == LAPACK_ROW_MAJOR ) {
-        for( j = 0; j < MIN( n, ldab ); j++ ) {
+        for( j = 0; j < n; j++ ) {
             for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 ); i++ ) {
                 if( LAPACK_DISNAN( ab[(size_t)i*ldab+j] ) )
                      return (lapack_logical) 1;

diff --git a/LAPACKE/utils/lapacke_sgb_nancheck.c b/LAPACKE/utils/lapacke_sgb_nancheck.c
index dcebcf6b570db7fbca9a47d05f5076bd6773a9a4..e20cb2421c743c209f0be25732ba0a1a3e1e9144 100644 (file)
--- a/LAPACKE/utils/lapacke_sgb_nancheck.c
+++ b/LAPACKE/utils/lapacke_sgb_nancheck.c
@@ -46,14 +46,14 @@ lapack_logical LAPACKE_sgb_nancheck( int matrix_layout, lapack_int m,
 
     if( matrix_layout == LAPACK_COL_MAJOR ) {
         for( j = 0; j < n; j++ ) {
-            for( i = MAX( ku-j, 0 ); i < MIN3( ldab, m+ku-j, kl+ku+1 );
+            for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 );
                  i++ ) {
                 if( LAPACK_SISNAN( ab[i+(size_t)j*ldab] ) )
                      return (lapack_logical) 1;
             }
         }
     } else if ( matrix_layout == LAPACK_ROW_MAJOR ) {
-        for( j = 0; j < MIN( n, ldab ); j++ ) {
+        for( j = 0; j < n; j++ ) {
             for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 ); i++ ) {
                 if( LAPACK_SISNAN( ab[(size_t)i*ldab+j] ) )
                      return (lapack_logical) 1;

diff --git a/LAPACKE/utils/lapacke_zgb_nancheck.c b/LAPACKE/utils/lapacke_zgb_nancheck.c
index 019a73578b1911cb8a02fcb5b6828f2478ff9c7d..8f0ffc1d69b8d7ea0c3d41b5139f78b8f2845bcd 100644 (file)
--- a/LAPACKE/utils/lapacke_zgb_nancheck.c
+++ b/LAPACKE/utils/lapacke_zgb_nancheck.c
@@ -46,14 +46,14 @@ lapack_logical LAPACKE_zgb_nancheck( int matrix_layout, lapack_int m,
 
     if( matrix_layout == LAPACK_COL_MAJOR ) {
         for( j = 0; j < n; j++ ) {
-            for( i = MAX( ku-j, 0 ); i < MIN3( ldab, m+ku-j, kl+ku+1 );
+            for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 );
                  i++ ) {
                 if( LAPACK_ZISNAN( ab[i+(size_t)j*ldab] ) )
                      return (lapack_logical) 1;
             }
         }
     } else if ( matrix_layout == LAPACK_ROW_MAJOR ) {
-        for( j = 0; j < MIN( n, ldab ); j++ ) {
+        for( j = 0; j < n; j++ ) {
             for( i = MAX( ku-j, 0 ); i < MIN( m+ku-j, kl+ku+1 ); i++ ) {
                 if( LAPACK_ZISNAN( ab[(size_t)i*ldab+j] ) )
                      return (lapack_logical) 1;