Home
last modified time | relevance | path

Searched refs:fshift (Results 1 – 17 of 17) sorted by relevance

/petsc/src/ksp/pc/impls/sor/
H A Dsor.c11 PetscReal fshift; member
33 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApply_SOR()
47 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApplyTranspose_SOR()
59 PetscCall(MatSOR(pc->pmat, b, jac->omega, stype, jac->fshift, its * jac->its, jac->lits, y)); in PCApplyRichardson_SOR()
76 …eal("-pc_sor_diagonal_shift", "Add to the diagonal entries", "", jac->fshift, &jac->fshift, NULL)); in PCSetFromOptions_SOR()
444 jac->fshift = 0.0; in PCCreate_SOR()
/petsc/src/mat/impls/sbaij/seq/
H A Drelax.h122 …qSBAIJ_ushort(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_SeqSBAIJ_ushort() argument
124 PetscErrorCode MatSOR_SeqSBAIJ(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, P… in MatSOR_SeqSBAIJ_ushort()
142 …if (fshift == -1.0) fshift = 0.0; /* negative fshift indicates do not error on zero diagonal; this… in MatSOR_SeqSBAIJ_ushort()
168 d = fshift + aa[ai[i]]; in MatSOR_SeqSBAIJ_ushort()
H A Dsbaij.c677 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqSBAIJ() local
688 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqSBAIJ()
690 if (fshift) { in MatAssemblyEnd_SeqSBAIJ()
694 PetscCall(PetscArraymove(ip - fshift, ip, N)); in MatAssemblyEnd_SeqSBAIJ()
695 PetscCall(PetscArraymove(ap - bs2 * fshift, ap, bs2 * N)); in MatAssemblyEnd_SeqSBAIJ()
700 fshift += imax[mbs - 1] - ailen[mbs - 1]; in MatAssemblyEnd_SeqSBAIJ()
707 …fshift || a->nounused != -1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Unused space detected in matrix: %"… in MatAssemblyEnd_SeqSBAIJ()
709 …t_FMT " unneeded, %" PetscInt_FMT " used\n", m, A->rmap->N, A->rmap->bs, fshift * bs2, a->nz * bs2… in MatAssemblyEnd_SeqSBAIJ()
715 A->info.nz_unneeded = (PetscReal)fshift * bs2; in MatAssemblyEnd_SeqSBAIJ()
1143 PetscInt fshift = 0, fshift_prev = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax, j, k; in MatEliminateZeros_SeqSBAIJ() local
[all …]
/petsc/src/mat/impls/blockmat/seq/
H A Dblockmat.c18 …Mat_Symmetric(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_BlockMat_Symmetric() argument
36 PetscCheck(!fshift, PETSC_COMM_SELF, PETSC_ERR_SUP, "No support yet for fshift"); in MatSOR_BlockMat_Symmetric()
120 …tSOR_BlockMat(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_BlockMat() argument
138 PetscCheck(!fshift, PETSC_COMM_SELF, PETSC_ERR_SUP, "No support yet for fshift"); in MatSOR_BlockMat()
599 PetscInt fshift = 0, i, j, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_BlockMat() local
609 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_BlockMat()
611 if (fshift) { in MatAssemblyEnd_BlockMat()
616 ip[j - fshift] = ip[j]; in MatAssemblyEnd_BlockMat()
617 ap[j - fshift] = ap[j]; in MatAssemblyEnd_BlockMat()
623 fshift += imax[m - 1] - ailen[m - 1]; in MatAssemblyEnd_BlockMat()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daij.c779 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()
[all …]
H A Daij.h159 PetscScalar fshift, omega; /* last used omega and fshift */ member
H A Dinode.c2420 …c PetscErrorCode MatInvertDiagonalForSOR_SeqAIJ_Inode(Mat A, PetscScalar omega, PetscScalar fshift) in MatInvertDiagonalForSOR_SeqAIJ_Inode() argument
2493 …_SeqAIJ_Inode(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_SeqAIJ_Inode() argument
2508 …PetscCheck(fshift == 0.0, PETSC_COMM_SELF, PETSC_ERR_SUP, "No support for fshift != 0.0; use -mat_… in MatSOR_SeqAIJ_Inode()
2509 PetscCall(MatInvertDiagonalForSOR_SeqAIJ_Inode(A, omega, fshift)); in MatSOR_SeqAIJ_Inode()
/petsc/src/mat/impls/aij/seq/aijsell/
H A Daijsell.c187 …OR_SeqAIJSELL(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_SeqAIJSELL() argument
193 PetscCall(MatSOR_SeqSELL(aijsell->S, bb, omega, flag, fshift, its, lits, xx)); in MatSOR_SeqAIJSELL()
/petsc/src/mat/impls/sell/seq/
H A Dsell.c831 static PetscErrorCode MatInvertDiagonalForSOR_SeqSELL(Mat A, PetscScalar omega, PetscScalar fshift) in MatInvertDiagonalForSOR_SeqSELL() argument
841 …if (a->idiagState == ((PetscObject)A)->state && a->omega == omega && a->fshift == fshift) PetscFun… in MatInvertDiagonalForSOR_SeqSELL()
852 if (omega == 1.0 && PetscRealPart(fshift) <= 0.0) { in MatInvertDiagonalForSOR_SeqSELL()
856 …PetscCheck(PetscRealPart(fshift), PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Zero diagonal on row %" … in MatInvertDiagonalForSOR_SeqSELL()
868 idiag[i] = omega / (fshift + val[diag[i]]); in MatInvertDiagonalForSOR_SeqSELL()
874 a->fshift = fshift; in MatInvertDiagonalForSOR_SeqSELL()
1190 PetscInt fshift = 1; in MatView_SeqSELL_ASCII() local
1202 …" PetscInt_FMT " %" PetscInt_FMT " %g %g\n", i + fshift, a->colidx[shift + a->sliceheight * j] + f… in MatView_SeqSELL_ASCII()
1204 … "%" PetscInt_FMT " %" PetscInt_FMT " %g\n", i + fshift, a->colidx[shift + a->sliceheight * j] + f… in MatView_SeqSELL_ASCII()
1687 PetscErrorCode MatSOR_SeqSELL(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, Pe… in MatSOR_SeqSELL() argument
[all …]
H A Dsell.h70 PetscScalar fshift, omega; /* last used omega and fshift */ member
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c217 …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()
2259 …fshift) 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()
[all …]
/petsc/src/mat/impls/sell/mpi/
H A Dmpisell.c1722 PetscErrorCode MatSOR_MPISELL(Mat matin, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift… in MatSOR_MPISELL() argument
1729 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1737 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1750 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_SYMMETRIC_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1754 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1766 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_FORWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1770 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISELL()
1782 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_BACKWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPISELL()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c2590 …_MPISBAIJ(Mat matin, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_MPISBAIJ() argument
2603 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISBAIJ()
2609 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, lits, xx)); in MatSOR_MPISBAIJ()
2645 … PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_SYMMETRIC_SWEEP, fshift, lits, lits, xx)); in MatSOR_MPISBAIJ()
2649 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISBAIJ()
2651 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPISBAIJ()
2666 … bb, omega, (MatSORType)(SOR_ZERO_INITIAL_GUESS | SOR_LOCAL_BACKWARD_SWEEP), fshift, lits, 1, xx)); in MatSOR_MPISBAIJ()
2716 …c1a, omega, (MatSORType)(SOR_ZERO_INITIAL_GUESS | SOR_LOCAL_FORWARD_SWEEP), fshift, lits, 1, xx1)); in MatSOR_MPISBAIJ()
/petsc/src/mat/impls/cdiagonal/
H A Dcdiagonal.c257 …ntDiagonal(Mat matin, Vec x, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_ConstantDiagonal() argument
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c2167 …R_MPIBAIJ(Mat matin, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_MPIBAIJ() argument
2174 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2182 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2195 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_SYMMETRIC_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2199 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2211 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_FORWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2215 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
2227 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_BACKWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIBAIJ()
/petsc/src/mat/impls/kaij/
H A Dkaij.c855 …atSOR_SeqKAIJ(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_SeqKAIJ() argument
869 PetscCheck(!fshift, PETSC_COMM_SELF, PETSC_ERR_SUP, "No support for diagonal shift"); in MatSOR_SeqKAIJ()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiaij.c1399 …OR_MPIAIJ(Mat matin, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, PetscInt its, Pet… in MatSOR_MPIAIJ() argument
1407 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1415 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1428 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_SYMMETRIC_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1432 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1444 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_FORWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1448 PetscCall((*mat->A->ops->sor)(mat->A, bb, omega, flag, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1460 PetscCall((*mat->A->ops->sor)(mat->A, bb1, omega, SOR_BACKWARD_SWEEP, fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1466 … bb, omega, (MatSORType)(SOR_ZERO_INITIAL_GUESS | SOR_LOCAL_BACKWARD_SWEEP), fshift, lits, 1, xx)); in MatSOR_MPIAIJ()
1485 …bb1, omega, (MatSORType)(SOR_ZERO_INITIAL_GUESS | SOR_LOCAL_FORWARD_SWEEP), fshift, lits, 1, xx1)); in MatSOR_MPIAIJ()