Lines Matching refs:adiag
751 const PetscInt *adiag; in MatSolve_SeqAIJ_inplace() local
756 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolve_SeqAIJ_inplace()
773 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_inplace()
781 v = aa + adiag[i] + 1; in MatSolve_SeqAIJ_inplace()
782 vi = aj + adiag[i] + 1; in MatSolve_SeqAIJ_inplace()
783 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_inplace()
786 x[*c--] = tmp[i] = sum * aa[adiag[i]]; in MatSolve_SeqAIJ_inplace()
807 const PetscInt *adiag; in MatMatSolve_SeqAIJ_inplace() local
817 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatMatSolve_SeqAIJ_inplace()
834 nz = adiag[i] - ai[i]; in MatMatSolve_SeqAIJ_inplace()
841 v = aa + adiag[i] + 1; in MatMatSolve_SeqAIJ_inplace()
842 vi = aj + adiag[i] + 1; in MatMatSolve_SeqAIJ_inplace()
843 nz = ai[i + 1] - adiag[i] - 1; in MatMatSolve_SeqAIJ_inplace()
846 x[c[i]] = tmp[i] = sum * aa[adiag[i]]; in MatMatSolve_SeqAIJ_inplace()
865 const PetscInt *adiag; in MatMatSolve_SeqAIJ() local
880 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatMatSolve_SeqAIJ()
905 v = aa + adiag[i + 1] + 1; in MatMatSolve_SeqAIJ()
906 vi = aj + adiag[i + 1] + 1; in MatMatSolve_SeqAIJ()
907 nz = adiag[i] - adiag[i + 1] - 1; in MatMatSolve_SeqAIJ()
929 const PetscInt *adiag = a->diag; in MatMatSolveTranspose_SeqAIJ() local
944 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatMatSolveTranspose_SeqAIJ()
960 v = aa + adiag[i + 1] + 1; in MatMatSolveTranspose_SeqAIJ()
961 vi = aj + adiag[i + 1] + 1; in MatMatSolveTranspose_SeqAIJ()
962 nz = adiag[i] - adiag[i + 1] - 1; in MatMatSolveTranspose_SeqAIJ()
996 const PetscInt *r, *c, *rout, *cout, *adiag; in MatSolve_SeqAIJ_InplaceWithPerm() local
1006 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolve_SeqAIJ_InplaceWithPerm()
1024 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_InplaceWithPerm()
1033 v = aa + adiag[i] + 1; in MatSolve_SeqAIJ_InplaceWithPerm()
1034 vi = aj + adiag[i] + 1; in MatSolve_SeqAIJ_InplaceWithPerm()
1035 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_InplaceWithPerm()
1038 x[*c--] = tmp[row] = sum * aa[adiag[i]]; in MatSolve_SeqAIJ_InplaceWithPerm()
1054 const PetscInt *ai = a->i, *aj = a->j, *adiag; in MatSolve_SeqAIJ_NaturalOrdering_inplace() local
1067 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1073 fortransolveaij_(&n, x, ai, aj, adiag, aa, b); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1081 nz = adiag[i] - ai_i; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1089 adiag_i = adiag[i]; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1111 const PetscInt *rout, *cout, *r, *c, *vi, *ai = a->i, *aj = a->j, *adiag; in MatSolveAdd_SeqAIJ_inplace() local
1119 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolveAdd_SeqAIJ_inplace()
1135 nz = adiag[i] - ai[i]; in MatSolveAdd_SeqAIJ_inplace()
1143 v = aa + adiag[i] + 1; in MatSolveAdd_SeqAIJ_inplace()
1144 vi = aj + adiag[i] + 1; in MatSolveAdd_SeqAIJ_inplace()
1145 nz = ai[i + 1] - adiag[i] - 1; in MatSolveAdd_SeqAIJ_inplace()
1148 tmp[i] = sum * aa[adiag[i]]; in MatSolveAdd_SeqAIJ_inplace()
1167 const PetscInt *rout, *cout, *r, *c, *vi, *ai = a->i, *aj = a->j, *adiag; in MatSolveAdd_SeqAIJ() local
1175 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolveAdd_SeqAIJ()
1200 v = aa + adiag[n - 1]; in MatSolveAdd_SeqAIJ()
1201 vi = aj + adiag[n - 1]; in MatSolveAdd_SeqAIJ()
1203 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveAdd_SeqAIJ()
1283 const PetscInt *rout, *cout, *r, *c, *adiag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqAIJ() local
1291 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolveTranspose_SeqAIJ()
1307 v = aa + adiag[i + 1] + 1; in MatSolveTranspose_SeqAIJ()
1308 vi = aj + adiag[i + 1] + 1; in MatSolveTranspose_SeqAIJ()
1309 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTranspose_SeqAIJ()
1401 const PetscInt *rout, *cout, *r, *c, *adiag, *ai = a->i, *aj = a->j, *vi; in MatSolveTransposeAdd_SeqAIJ() local
1410 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatSolveTransposeAdd_SeqAIJ()
1426 v = aa + adiag[i + 1] + 1; in MatSolveTransposeAdd_SeqAIJ()
1427 vi = aj + adiag[i + 1] + 1; in MatSolveTransposeAdd_SeqAIJ()
1428 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTransposeAdd_SeqAIJ()
1476 const PetscInt n = A->rmap->n, *ai = a->i, *aj, *adiag; in MatILUFactorSymbolic_SeqAIJ_ilu0() local
1482 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1504 nz = adiag[i] - ai[i]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1513 nz = ai[i + 1] - adiag[i] - 1; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1514 aj = a->j + adiag[i] + 1; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1719 const PetscInt *adiag; in MatCholeskyFactorNumeric_SeqAIJ() local
1722 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatCholeskyFactorNumeric_SeqAIJ()
1731 d = aa[adiag[i]]; in MatCholeskyFactorNumeric_SeqAIJ()
1886 const PetscInt *adiag; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1889 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, NULL)); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1898 d = aa[adiag[i]]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2057 const PetscInt *rip, *riip, *adiag; in MatICCFactorSymbolic_SeqAIJ() local
2070 PetscCall(MatGetDiagonalMarkers_SeqAIJ(A, &adiag, &diagDense)); in MatICCFactorSymbolic_SeqAIJ()
2082 ncols = ai[i + 1] - adiag[i]; in MatICCFactorSymbolic_SeqAIJ()
2089 aj = a->j + adiag[i] + 1; /* 1st entry of U(i,:) without diagonal */ in MatICCFactorSymbolic_SeqAIJ()
2090 ncols = ai[i + 1] - adiag[i] - 1; in MatICCFactorSymbolic_SeqAIJ()
2423 const PetscInt *ai = a->i, *aj = a->j, *adiag = a->diag, *vi; in MatSolve_SeqAIJ_NaturalOrdering() local
2451 v = aa + adiag[i + 1] + 1; in MatSolve_SeqAIJ_NaturalOrdering()
2452 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqAIJ_NaturalOrdering()
2453 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqAIJ_NaturalOrdering()
2470 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag, nz; in MatSolve_SeqAIJ() local
2504 v = aa + adiag[i + 1] + 1; in MatSolve_SeqAIJ()
2505 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqAIJ()
2506 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqAIJ()