Lines Matching refs:idiag
836 PetscScalar *idiag, *mdiag; in MatInvertDiagonalForSOR_SeqSELL() local
845 if (!a->idiag) { in MatInvertDiagonalForSOR_SeqSELL()
846 PetscCall(PetscMalloc3(m, &a->idiag, m, &a->mdiag, m, &a->ssor_work)); in MatInvertDiagonalForSOR_SeqSELL()
850 idiag = a->idiag; in MatInvertDiagonalForSOR_SeqSELL()
862 idiag[i] = 1.0 / val[diag[i]]; in MatInvertDiagonalForSOR_SeqSELL()
868 idiag[i] = omega / (fshift + 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
1702 idiag = a->idiag; in MatSOR_SeqSELL()
1720 x[i] = sum * idiag[i]; in MatSOR_SeqSELL()
1732 x[i] = sum * idiag[i]; in MatSOR_SeqSELL()
1734 x[i] = (1. - omega) * x[i] + sum * idiag[i]; /* omega in idiag */ in MatSOR_SeqSELL()
1753 x[i] = (1. - omega) * x[i] + sum * idiag[i]; /* omega in idiag */ in MatSOR_SeqSELL()
1766 x[i] = (1. - omega) * x[i] + (sum + mdiag[i] * x[i]) * idiag[i]; in MatSOR_SeqSELL()
1770 x[i] = (1. - omega) * x[i] + sum * idiag[i]; /* omega in idiag */ in MatSOR_SeqSELL()
2168 b->idiag = NULL; in MatCreate_SeqSELL()
2265 c->idiag = NULL; in MatDuplicateNoCreate_SeqSELL()