Lines Matching refs:fshift
779 PetscInt nzd = 0, fshift = 1, *sptr; in MatView_SeqAIJ_ASCII() local
815 … (a->j[j] >= i) PetscCall(PetscViewerASCIIPrintf(viewer, " %" PetscInt_FMT " ", a->j[j] + fshift)); in MatView_SeqAIJ_ASCII()
871 PetscInt fshift = 1; in MatView_SeqAIJ_ASCII() local
882 …intf(viewer, "%" PetscInt_FMT " %" PetscInt_FMT " %g %g\n", i + fshift, a->j[j] + fshift, (double)… in MatView_SeqAIJ_ASCII()
884 …IPrintf(viewer, "%" PetscInt_FMT " %" PetscInt_FMT " %g\n", i + fshift, a->j[j] + fshift, (double)… in MatView_SeqAIJ_ASCII()
1100 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqAIJ() local
1116 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqAIJ()
1118 if (fshift) { in MatAssemblyEnd_SeqAIJ()
1122 PetscCall(PetscArraymove(ip - fshift, ip, N)); in MatAssemblyEnd_SeqAIJ()
1123 if (!A->structure_only) PetscCall(PetscArraymove(ap - fshift, ap, N)); in MatAssemblyEnd_SeqAIJ()
1128 fshift += imax[m - 1] - ailen[m - 1]; in MatAssemblyEnd_SeqAIJ()
1143 …fshift || a->nounused != -1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Unused space detected in matrix: %"… in MatAssemblyEnd_SeqAIJ()
1144 …orage space: %" PetscInt_FMT " unneeded, %" PetscInt_FMT " used\n", m, A->cmap->n, fshift, a->nz)); in MatAssemblyEnd_SeqAIJ()
1150 A->info.nz_unneeded = (PetscReal)fshift; in MatAssemblyEnd_SeqAIJ()
1809 static PetscErrorCode MatInvertDiagonalForSOR_SeqAIJ(Mat A, PetscScalar omega, PetscScalar fshift) in MatInvertDiagonalForSOR_SeqAIJ() argument
1819 …if (a->idiagState == ((PetscObject)A)->state && a->omega == omega && a->fshift == fshift) PetscFun… in MatInvertDiagonalForSOR_SeqAIJ()
1827 if (omega == 1.0 && PetscRealPart(fshift) <= 0.0) { in MatInvertDiagonalForSOR_SeqAIJ()
1831 …PetscCheck(PetscRealPart(fshift), PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Zero diagonal on row %" … in MatInvertDiagonalForSOR_SeqAIJ()
1843 idiag[i] = omega / (fshift + v[diag[i]]); in MatInvertDiagonalForSOR_SeqAIJ()
1850 a->fshift = fshift; in MatInvertDiagonalForSOR_SeqAIJ()
1854 PetscErrorCode MatSOR_SeqAIJ(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, Pet… in MatSOR_SeqAIJ() argument
1864 if (a->inode.use && a->inode.checked && omega == 1.0 && fshift == 0.0) { in MatSOR_SeqAIJ()
1865 PetscCall(MatSOR_SeqAIJ_Inode(A, bb, omega, flag, fshift, its, lits, xx)); in MatSOR_SeqAIJ()
1869 PetscCall(MatInvertDiagonalForSOR_SeqAIJ(A, omega, fshift)); in MatSOR_SeqAIJ()
1883 d = fshift + mdiag[i]; in MatSOR_SeqAIJ()
4804 b->fshift = 0.0; in MatCreate_SeqAIJ()
5315 PetscInt m = A->rmap->n, fshift = 0, fshift_prev = 0, i, k; in MatEliminateZeros_SeqAIJ() local
5324 if (aa[k] == 0 && (aj[k] != i - 1 || !keep)) fshift++; in MatEliminateZeros_SeqAIJ()
5327 aa[k - fshift] = aa[k]; in MatEliminateZeros_SeqAIJ()
5328 aj[k - fshift] = aj[k]; in MatEliminateZeros_SeqAIJ()
5332 fshift_prev = fshift; in MatEliminateZeros_SeqAIJ()
5334 ailen[i - 1] = imax[i - 1] = ai[i] - fshift - ai[i - 1]; in MatEliminateZeros_SeqAIJ()
5335 a->nonzerorowcnt += ((ai[i] - fshift - ai[i - 1]) > 0); in MatEliminateZeros_SeqAIJ()
5338 if (fshift) { in MatEliminateZeros_SeqAIJ()
5340 ai[m] -= fshift; in MatEliminateZeros_SeqAIJ()
5343 …eliminated: %" PetscInt_FMT "; nonzeros left: %" PetscInt_FMT "\n", m, A->cmap->n, fshift, a->nz)); in MatEliminateZeros_SeqAIJ()
5345 A->info.nz_unneeded += (PetscReal)fshift; in MatEliminateZeros_SeqAIJ()