Lines Matching refs:fshift

217 …atSOR_SeqBAIJ(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet…  in MatSOR_SeqBAIJ()  argument
231 PetscCheck(!fshift, PETSC_COMM_SELF, PETSC_ERR_SUP, "No support for diagonal shift"); in MatSOR_SeqBAIJ()
2217 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqBAIJ() local
2229 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqBAIJ()
2231 if (fshift) { in MatAssemblyEnd_SeqBAIJ()
2235 PetscCall(PetscArraymove(ip - fshift, ip, N)); in MatAssemblyEnd_SeqBAIJ()
2236 if (!A->structure_only) PetscCall(PetscArraymove(ap - bs2 * fshift, ap, bs2 * N)); in MatAssemblyEnd_SeqBAIJ()
2241 fshift += imax[mbs - 1] - ailen[mbs - 1]; in MatAssemblyEnd_SeqBAIJ()
2259fshift) PetscCheck(a->nounused != -1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Unused space detected in m… in MatAssemblyEnd_SeqBAIJ()
2260 …t_FMT " unneeded, %" PetscInt_FMT " used\n", m, A->cmap->n, A->rmap->bs, fshift * bs2, a->nz * bs2… in MatAssemblyEnd_SeqBAIJ()
2266 A->info.nz_unneeded = (PetscReal)fshift * bs2; in MatAssemblyEnd_SeqBAIJ()
2924 PetscInt fshift = 0, fshift_prev = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax, j, k; in MatEliminateZeros_SeqBAIJ() local
2940 if (zero && (aj[k] != i - 1 || !keep)) fshift++; in MatEliminateZeros_SeqBAIJ()
2943 aj[k - fshift] = aj[k]; in MatEliminateZeros_SeqBAIJ()
2944 PetscCall(PetscArraymove(ap - bs2 * fshift, ap, bs2)); in MatEliminateZeros_SeqBAIJ()
2948 fshift_prev = fshift; in MatEliminateZeros_SeqBAIJ()
2949 ailen[i - 1] = imax[i - 1] = ai[i] - fshift - ai[i - 1]; in MatEliminateZeros_SeqBAIJ()
2950 a->nonzerorowcnt += ((ai[i] - fshift - ai[i - 1]) > 0); in MatEliminateZeros_SeqBAIJ()
2953 if (fshift) { in MatEliminateZeros_SeqBAIJ()
2955 ai[mbs] -= fshift; in MatEliminateZeros_SeqBAIJ()
2958 …eliminated: %" PetscInt_FMT "; nonzeros left: %" PetscInt_FMT "\n", m, A->cmap->n, fshift, a->nz)); in MatEliminateZeros_SeqBAIJ()
2960 A->info.nz_unneeded += (PetscReal)fshift; in MatEliminateZeros_SeqBAIJ()