Refs #723. Avoid out of boundary for getf2.
This commit is contained in:
parent
f5df444ceb
commit
13f0f8c10e
|
@ -94,6 +94,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
|
||||
ipiv[j + offset] = jp + offset;
|
||||
jp--;
|
||||
temp1 = *(b + jp);
|
||||
|
|
|
@ -98,6 +98,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
|
||||
ipiv[j + offset] = jp + offset;
|
||||
jp--;
|
||||
|
||||
|
|
Loading…
Reference in New Issue