Home
last modified time | relevance | path

Searched refs:mdiag (Results 1 – 4 of 4) sorted by relevance

/petsc/src/mat/impls/sell/seq/
H A Dsell.c836 PetscScalar *idiag, *mdiag; in MatInvertDiagonalForSOR_SeqSELL() local
846 PetscCall(PetscMalloc3(m, &a->idiag, m, &a->mdiag, m, &a->ssor_work)); in MatInvertDiagonalForSOR_SeqSELL()
849 mdiag = a->mdiag; in MatInvertDiagonalForSOR_SeqSELL()
854 mdiag[i] = val[diag[i]]; in MatInvertDiagonalForSOR_SeqSELL()
855 if (!PetscAbsScalar(mdiag[i])) { /* zero diagonal */ in MatInvertDiagonalForSOR_SeqSELL()
867 mdiag[i] = val[diag[i]]; in MatInvertDiagonalForSOR_SeqSELL()
899 PetscCall(PetscFree3(a->idiag, a->mdiag, a->ssor_work)); in MatDestroy_SeqSELL()
1691 const MatScalar *idiag = NULL, *mdiag; in MatSOR_SeqSELL() local
1703 mdiag = a->mdiag; in MatSOR_SeqSELL()
1766 x[i] = (1. - omega) * x[i] + (sum + mdiag[i] * x[i]) * idiag[i]; in MatSOR_SeqSELL()
[all …]
H A Dsell.h67 PetscScalar *mdiag, *idiag; /* diagonal values, inverse of diagonal entries */ member
/petsc/src/mat/impls/aij/seq/
H A Daij.c1217 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()
[all …]
H A Daij.h156 PetscScalar *mdiag, *idiag; /* diagonal values, inverse of diagonal entries */ member