| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sro.c | 40 const PetscInt *rip,*riip; in MatReorderingSeqSBAIJ()
44 PetscCall(ISGetIndices(perm,&rip)); in MatReorderingSeqSBAIJ()
49 …for (i=0; i<mbs; i++) PetscCheck(rip[i] == riip[i],PETSC_COMM_SELF,PETSC_ERR_ARG_INCOMP,"Non-symme… in MatReorderingSeqSBAIJ()
83 if (rip[k] < rip[i]) aj[j] = i; in MatReorderingSeqSBAIJ()
127 PetscCall(ISRestoreIndices(perm,&rip)); in MatReorderingSeqSBAIJ()
|
| H A D | sbaijfact.c | 42 const PetscInt *rip, *ai, *aj; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
52 PetscCall(ISGetIndices(perm, &rip)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
88 jmin = ai[rip[k]] + 1; /* exclude diag[k] */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
89 jmax = ai[rip[k] + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
91 vj = rip[aj[j]]; /* col. value */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
176 PetscCall(ISRestoreIndices(perm, &rip)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
226 const PetscInt *rip, *ai = a->i, *aj = a->j; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
247 PetscCall(ISGetIndices(perm, &rip)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
332 PetscCall(ISRestoreIndices(perm, &rip)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
391 const PetscInt *rip, *ai, *aj; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() local
[all …]
|
| H A D | sbaijfact2.c | 2282 const PetscInt *rip, mbs = a->mbs, *ai, *aj; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2300 PetscCall(ISGetIndices(perm, &rip)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2323 nz = ai[rip[k] + 1] - ai[rip[k]]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2324 j = ai[rip[k]]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2326 vj = rip[aj[j++]]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2422 PetscCall(ISRestoreIndices(perm, &rip)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2475 const PetscInt *rip; in MatICCFactorSymbolic_SeqSBAIJ() local
2519 PetscCall(ISGetIndices(perm, &rip)); in MatICCFactorSymbolic_SeqSBAIJ()
2549 PetscCall(PetscIncompleteLLInit(ncols, cols, am, rip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ()
2614 PetscCall(ISRestoreIndices(perm, &rip)); in MatICCFactorSymbolic_SeqSBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/seq/bas/ |
| H A D | basfactor.c | 11 const PetscInt *rip, *riip, *adiag; in MatICCFactorSymbolic_SeqAIJ_Bas() local
40 PetscCall(ISGetIndices(perm, &rip)); in MatICCFactorSymbolic_SeqAIJ_Bas()
46 PetscCall(spbas_apply_reordering(&Pattern_0, rip, riip)); in MatICCFactorSymbolic_SeqAIJ_Bas()
99 const PetscInt *rip, *riip; in MatCholeskyFactorNumeric_SeqAIJ_Bas() local
119 PetscCall(ISGetIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ_Bas()
127 …ierr = (int)spbas_incomplete_cholesky(A, rip, riip, Pattern, droptol, shiftnz, &matrix_LT, &succes… in MatCholeskyFactorNumeric_SeqAIJ_Bas()
138 PetscCall(ISRestoreIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ_Bas()
|
| H A D | spbas_cholesky.h | 239 PetscErrorCode spbas_incomplete_cholesky(Mat A, const PetscInt *rip, const PetscInt *riip, spbas_ma… in spbas_incomplete_cholesky() argument
309 A_nnz = ai[rip[i] + 1] - ai[rip[i]]; in spbas_incomplete_cholesky()
310 A_icol = &aj[ai[rip[i]]]; in spbas_incomplete_cholesky()
311 A_val = &aa[ai[rip[i]]]; in spbas_incomplete_cholesky()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 752 const PetscInt *rip; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
775 PetscCall(ISGetIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
791 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
792 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
794 col = rip[aj[j]]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
870 PetscCall(ISRestoreIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
1023 const PetscInt *rip; in MatICCFactorSymbolic_SeqBAIJ() local
1050 PetscCall(ISGetIndices(perm, &rip)); in MatICCFactorSymbolic_SeqBAIJ()
1106 ncols = ai[rip[k] + 1] - ai[rip[k]]; in MatICCFactorSymbolic_SeqBAIJ()
1110 i = rip[*(aj + ai[rip[k]] + j)]; in MatICCFactorSymbolic_SeqBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 1709 const PetscInt *rip, *riip; in MatCholeskyFactorNumeric_SeqAIJ() local
1744 PetscCall(ISGetIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ()
1767 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqAIJ()
1768 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqAIJ()
1836 PetscCall(ISRestoreIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ()
1876 const PetscInt *rip, *riip; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1911 PetscCall(ISGetIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1931 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1932 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2003 PetscCall(ISRestoreIndices(ip, &rip)); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
[all …]
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 736 const PetscInt *irip, *rip; in MatSeqAIJHIPSPARSEICCAnalysisAndCopyToGPU() local
740 PetscCall(ISGetIndices(ip, &rip)); in MatSeqAIJHIPSPARSEICCAnalysisAndCopyToGPU()
743 hipsparseTriFactors->rpermIndices->assign(rip, rip + n); in MatSeqAIJHIPSPARSEICCAnalysisAndCopyToGPU()
747 PetscCall(ISRestoreIndices(ip, &rip)); in MatSeqAIJHIPSPARSEICCAnalysisAndCopyToGPU()
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 996 const PetscInt *irip, *rip; in MatSeqAIJCUSPARSEICCAnalysisAndCopyToGPU() local
1000 PetscCall(ISGetIndices(ip, &rip)); in MatSeqAIJCUSPARSEICCAnalysisAndCopyToGPU()
1002 cusparseTriFactors->rpermIndices->assign(rip, rip + n); in MatSeqAIJCUSPARSEICCAnalysisAndCopyToGPU()
1007 PetscCall(ISRestoreIndices(ip, &rip)); in MatSeqAIJCUSPARSEICCAnalysisAndCopyToGPU()
|