Lines Matching refs:mdiag
1217 PetscCall(PetscFree3(a->idiag, a->mdiag, a->ssor_work)); in MatReset_SeqAIJ()
1698 PetscInt *mdiag = NULL; in MatShift_SeqAIJ() local
1700 PetscCall(PetscCalloc1(A->rmap->n, &mdiag)); in MatShift_SeqAIJ()
1703 mdiag[i] = 1; in MatShift_SeqAIJ()
1713 for (PetscInt i = 0; i < PetscMin(A->rmap->n, A->cmap->n); i++) a->imax[i] += mdiag[i]; in MatShift_SeqAIJ()
1718 …PetscCall(MatSetValues(A, 1, &i, a->imax[i] - mdiag[i], &oldj[oldi[i]], &olda[oldi[i]], ADD_VALUES… in MatShift_SeqAIJ()
1726 PetscCall(PetscFree(mdiag)); in MatShift_SeqAIJ()
1814 PetscScalar *idiag, *mdiag; in MatInvertDiagonalForSOR_SeqAIJ() local
1822 if (!a->idiag) PetscCall(PetscMalloc3(m, &a->idiag, m, &a->mdiag, m, &a->ssor_work)); in MatInvertDiagonalForSOR_SeqAIJ()
1824 mdiag = a->mdiag; in MatInvertDiagonalForSOR_SeqAIJ()
1829 mdiag[i] = v[diag[i]]; in MatInvertDiagonalForSOR_SeqAIJ()
1830 if (!PetscAbsScalar(mdiag[i])) { /* zero diagonal */ in MatInvertDiagonalForSOR_SeqAIJ()
1842 mdiag[i] = v[diag[i]]; in MatInvertDiagonalForSOR_SeqAIJ()
1858 const MatScalar *v, *idiag = NULL, *mdiag, *aa; in MatSOR_SeqAIJ() local
1873 mdiag = a->mdiag; in MatSOR_SeqAIJ()
1883 d = fshift + mdiag[i]; in MatSOR_SeqAIJ()
2002 x[i] = (1. - omega) * x[i] + (sum + mdiag[i] * x[i]) * idiag[i]; in MatSOR_SeqAIJ()
4801 b->mdiag = NULL; in MatCreate_SeqAIJ()