diff --git a/src/mat/impls/sbaij/seq/sbaij.c b/src/mat/impls/sbaij/seq/sbaij.c index c256ca1..93bddf4 100644 --- a/src/mat/impls/sbaij/seq/sbaij.c +++ b/src/mat/impls/sbaij/seq/sbaij.c @@ -87,12 +87,12 @@ static PetscErrorCode MatGetRowIJ_SeqSBAIJ(Mat A,PetscInt oshift,PetscBool symme ierr = PetscMalloc2((n+1)*bs,ia,nz*bs,ja);CHKERRQ(ierr); for (i=0; ii[i]*bs+j+oshift; + (*ia)[i*bs+j] = a->i[i]*bs+j+oshift; } } for (i=0; ij[i]*bs+j+oshift; + (*ja)[i*bs+j] = a->j[i]*bs+j+oshift; } } } else { /* blockcompressed */