Merge pull request #4131 from martin-frbg/lapack878

Fix computation of UPLO in LAPACKE_?larfb (Reference-LAPACK PR 878)
This commit is contained in:
Martin Kroeker 2023-07-10 01:50:16 +02:00 committed by GitHub
commit bcebe9b4c9
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
8 changed files with 8 additions and 8 deletions

View File

@ -58,7 +58,7 @@ lapack_int LAPACKE_clarfb( int matrix_layout, char side, char trans, char direct
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
if( ( col && k > nrows_v ) || ( !col && k > ncols_v ) ) {
LAPACKE_xerbla( "LAPACKE_clarfb", -8 );

View File

@ -60,7 +60,7 @@ lapack_int LAPACKE_clarfb_work( int matrix_layout, char side, char trans,
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
ldc_t = MAX(1,m);
ldt_t = MAX(1,k);

View File

@ -57,7 +57,7 @@ lapack_int LAPACKE_dlarfb( int matrix_layout, char side, char trans, char direct
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
if( ( col && k > nrows_v ) || ( !col && k > ncols_v ) ) {
LAPACKE_xerbla( "LAPACKE_dlarfb", -8 );

View File

@ -59,7 +59,7 @@ lapack_int LAPACKE_dlarfb_work( int matrix_layout, char side, char trans,
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
ldc_t = MAX(1,m);
ldt_t = MAX(1,k);

View File

@ -57,7 +57,7 @@ lapack_int LAPACKE_slarfb( int matrix_layout, char side, char trans, char direct
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
if( ( col && k > nrows_v ) || ( !col && k > ncols_v ) ) {
LAPACKE_xerbla( "LAPACKE_slarfb", -8 );

View File

@ -59,7 +59,7 @@ lapack_int LAPACKE_slarfb_work( int matrix_layout, char side, char trans,
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
ldc_t = MAX(1,m);
ldt_t = MAX(1,k);

View File

@ -58,7 +58,7 @@ lapack_int LAPACKE_zlarfb( int matrix_layout, char side, char trans, char direct
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
if( ( col && k > nrows_v ) || ( !col && k > ncols_v ) ) {
LAPACKE_xerbla( "LAPACKE_zlarfb", -8 );

View File

@ -60,7 +60,7 @@ lapack_int LAPACKE_zlarfb_work( int matrix_layout, char side, char trans,
nrows_v = ( col && left ) ? m : ( ( col && !left ) ? n : ( !col ? k : 1) );
ncols_v = ( !col && left ) ? m : ( ( !col && !left ) ? n : ( col ? k : 1 ) );
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
ldc_t = MAX(1,m);
ldt_t = MAX(1,k);