floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Merge pull request #4207 from martin-frbg/issue4174-2 

Clarify the comment on the out-of-bounds check in ?GETF2
This commit is contained in:
Martin Kroeker 2023-08-26 12:05:37 +02:00 committed by GitHub
parent 9c4ae4d4fb 3b6050ac04
commit 12d8f219d6
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
lapack/getf2/getf2_k.c
View File

@ -95,7 +95,7 @@ blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa,
GEMV_N(m - j, j, 0, dm1, a + j, lda, b, 1, b + j, 1, sb);
jp = j + IAMAX_K(m - j, b + j, 1);
if (jp>m) jp = m; //avoid out of boundary
if (jp>m) jp = m; //avoid out of boundary when the iamax kernel does not cope with NaN in input, see gh issue 723
ipiv[j + offset] = jp + offset;
jp--;
temp1 = *(b + jp);

2
lapack/getf2/zgetf2_k.c
View File

@ -99,7 +99,7 @@ blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa,
GEMV_N(m - j, j, 0, dm1, ZERO, a + j * 2, lda, b, 1, b + j * 2, 1, sb);
jp = j + IAMAX_K(m - j, b + j * 2, 1);
if (jp>m) jp = m; //avoid out of boundary
if (jp>m) jp = m; //avoid out of boundary when the iamax kernel does not cope with NaN in input, see gh issue 723
ipiv[j + offset] = jp + offset;
jp--;