Skip to content

Lapacke x larfb nancheck #742

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 9 commits into from
Oct 31, 2022
Merged

Lapacke x larfb nancheck #742

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
a7e0724
Used trapezoidal NaN check in LAPACKE_?larfb
STCE2023 Oct 30, 2022
2da7000
Fixed bug in trapezoidal NaN check.
STCE2023 Oct 30, 2022
c43d27a
Paranthesis formatting.
STCE2023 Oct 30, 2022
70c4244
Paranthesis formatting.
STCE2023 Oct 30, 2022
6b5c5bb
Implemented transpose function for trapezoidal matrices.
STCE2023 Oct 30, 2022
c91f13e
Fixed rows and cols order in ?larfb
STCE2023 Oct 30, 2022
87cb1d1
Added trapezoidal transpositions to LAPACKE_?larfb_work.
STCE2023 Oct 30, 2022
de9eb94
Updated lapacke_utils header
STCE2023 Oct 30, 2022
6ef9ae9
Updated build system.
STCE2023 Oct 30, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Paranthesis formatting.
  • Loading branch information
@STCE2023
STCE2023 committed Oct 30, 2022
commit 70c4244912d50dcbb1fa10686514632910da734a
12 changes: 6 additions & 6 deletions LAPACKE/utils/lapacke_ctz_nancheck.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,23 +106,23 @@ lapack_logical LAPACKE_ctz_nancheck( int matrix_layout, char direct, char uplo,
lapack_int tri_offset = 0;
lapack_int tri_n = MIN(m,n);
lapack_int rect_offset = -1;
lapack_int rect_m = (m > n) ? m - n : m;
lapack_int rect_n = (n > m) ? n - m : n;
lapack_int rect_m = ( m > n ) ? m - n : m;
lapack_int rect_n = ( n > m ) ? n - m : n;

/* Fix offsets depending on the shape of the matrix */
if( front ) {
if( lower && m > n) {
if( lower && m > n ) {
rect_offset = tri_n * ( !colmaj ? lda : 1 );
} else if( !lower && n > m) {
} else if( !lower && n > m ) {
rect_offset = tri_n * ( colmaj ? lda : 1 );
}
} else {
if( m > n) {
if( m > n ) {
tri_offset = rect_m * ( !colmaj ? lda : 1 );
if( !lower ) {
rect_offset = 0;
}
} else if( n > m) {
} else if( n > m ) {
tri_offset = rect_n * ( colmaj ? lda : 1 );
if( lower ) {
rect_offset = 0;
Expand Down
12 changes: 6 additions & 6 deletions LAPACKE/utils/lapacke_dtz_nancheck.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -105,23 +105,23 @@ lapack_logical LAPACKE_dtz_nancheck( int matrix_layout, char direct, char uplo,
lapack_int tri_offset = 0;
lapack_int tri_n = MIN(m,n);
lapack_int rect_offset = -1;
lapack_int rect_m = (m > n) ? m - n : m;
lapack_int rect_n = (n > m) ? n - m : n;
lapack_int rect_m = ( m > n ) ? m - n : m;
lapack_int rect_n = ( n > m ) ? n - m : n;

/* Fix offsets depending on the shape of the matrix */
if( front ) {
if( lower && m > n) {
if( lower && m > n ) {
rect_offset = tri_n * ( !colmaj ? lda : 1 );
} else if( !lower && n > m) {
} else if( !lower && n > m ) {
rect_offset = tri_n * ( colmaj ? lda : 1 );
}
} else {
if( m > n) {
if( m > n ) {
tri_offset = rect_m * ( !colmaj ? lda : 1 );
if( !lower ) {
rect_offset = 0;
}
} else if( n > m) {
} else if( n > m ) {
tri_offset = rect_n * ( colmaj ? lda : 1 );
if( lower ) {
rect_offset = 0;
Expand Down
12 changes: 6 additions & 6 deletions LAPACKE/utils/lapacke_stz_nancheck.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -105,23 +105,23 @@ lapack_logical LAPACKE_stz_nancheck( int matrix_layout, char direct, char uplo,
lapack_int tri_offset = 0;
lapack_int tri_n = MIN(m,n);
lapack_int rect_offset = -1;
lapack_int rect_m = (m > n) ? m - n : m;
lapack_int rect_n = (n > m) ? n - m : n;
lapack_int rect_m = ( m > n ) ? m - n : m;
lapack_int rect_n = ( n > m ) ? n - m : n;

/* Fix offsets depending on the shape of the matrix */
if( front ) {
if( lower && m > n) {
if( lower && m > n ) {
rect_offset = tri_n * ( !colmaj ? lda : 1 );
} else if( !lower && n > m) {
} else if( !lower && n > m ) {
rect_offset = tri_n * ( colmaj ? lda : 1 );
}
} else {
if( m > n) {
if( m > n ) {
tri_offset = rect_m * ( !colmaj ? lda : 1 );
if( !lower ) {
rect_offset = 0;
}
} else if( n > m) {
} else if( n > m ) {
tri_offset = rect_n * ( colmaj ? lda : 1 );
if( lower ) {
rect_offset = 0;
Expand Down
12 changes: 6 additions & 6 deletions LAPACKE/utils/lapacke_ztz_nancheck.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,23 +106,23 @@ lapack_logical LAPACKE_ztz_nancheck( int matrix_layout, char direct, char uplo,
lapack_int tri_offset = 0;
lapack_int tri_n = MIN(m,n);
lapack_int rect_offset = -1;
lapack_int rect_m = (m > n) ? m - n : m;
lapack_int rect_n = (n > m) ? n - m : n;
lapack_int rect_m = ( m > n ) ? m - n : m;
lapack_int rect_n = ( n > m ) ? n - m : n;

/* Fix offsets depending on the shape of the matrix */
if( front ) {
if( lower && m > n) {
if( lower && m > n ) {
rect_offset = tri_n * ( !colmaj ? lda : 1 );
} else if( !lower && n > m) {
} else if( !lower && n > m ) {
rect_offset = tri_n * ( colmaj ? lda : 1 );
}
} else {
if( m > n) {
if( m > n ) {
tri_offset = rect_m * ( !colmaj ? lda : 1 );
if( !lower ) {
rect_offset = 0;
}
} else if( n > m) {
} else if( n > m ) {
tri_offset = rect_n * ( colmaj ? lda : 1 );
if( lower ) {
rect_offset = 0;
Expand Down