| /petsc/src/mat/interface/ftn-custom/ |
| H A D | zmatrixf.c | 35 …createsubmatrices_(Mat *mat, PetscInt *n, IS *isrow, IS *iscol, MatReuse *scall, F90Array1d *ptr, … in matcreatesubmatrices_() argument 39 if (*scall == MAT_INITIAL_MATRIX) { in matcreatesubmatrices_() 40 *ierr = MatCreateSubMatrices(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatrices_() 44 *ierr = MatCreateSubMatrices(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatrices_() 48 …atesubmatricesmpi_(Mat *mat, PetscInt *n, IS *isrow, IS *iscol, MatReuse *scall, F90Array1d *ptr, … in matcreatesubmatricesmpi_() argument 52 if (*scall == MAT_INITIAL_MATRIX) { in matcreatesubmatricesmpi_() 53 *ierr = MatCreateSubMatricesMPI(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatricesmpi_() 59 *ierr = MatCreateSubMatricesMPI(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatricesmpi_()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpb_aij.c | 3 PetscErrorCode MatGetMultiProcBlock_MPIAIJ(Mat mat, MPI_Comm subComm, MatReuse scall, Mat *subMat) in MatGetMultiProcBlock_MPIAIJ() argument 17 if (scall == MAT_INITIAL_MATRIX) { in MatGetMultiProcBlock_MPIAIJ() 47 if (scall == MAT_INITIAL_MATRIX) { in MatGetMultiProcBlock_MPIAIJ() 89 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscFree(nnz)); in MatGetMultiProcBlock_MPIAIJ()
|
| H A D | mpiov.c | 1083 PetscErrorCode MatCreateSubMatrix_MPIAIJ_All(Mat A, MatCreateSubMatrixOption flag, MatReuse scall, … in MatCreateSubMatrix_MPIAIJ_All() argument 1095 if (scall == MAT_INITIAL_MATRIX) { in MatCreateSubMatrix_MPIAIJ_All() 1213 …Local(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, PetscBool allcolu… in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() argument 1248 PetscValidLogicalCollectiveEnum(C, scall, 5); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() 1268 if (scall == MAT_INITIAL_MATRIX) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() 1844 if (scall == MAT_INITIAL_MATRIX || !iscolsorted) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() 1907 if (scall == MAT_INITIAL_MATRIX) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() 1919 …_SingleIS(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_MPIAIJ_SingleIS() argument 1926 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscCalloc1(2, submat)); in MatCreateSubMatrices_MPIAIJ_SingleIS() 1933 …PetscCall(MatCreateSubMatrices_MPIAIJ_SingleIS_Local(C, ismax, isrow, iscol, scall, allcolumns, *s… in MatCreateSubMatrices_MPIAIJ_SingleIS() [all …]
|
| H A D | mpiaij.c | 4570 …teMPIMatConcatenateSeqMat_MPIAIJ(MPI_Comm comm, Mat inmat, PetscInt n, MatReuse scall, Mat *outmat) in MatCreateMPIMatConcatenateSeqMat_MPIAIJ() argument 4579 if (scall == MAT_INITIAL_MATRIX) { /* symbolic phase */ in MatCreateMPIMatConcatenateSeqMat_MPIAIJ() 5028 …eateMPIAIJSumSeqAIJ(MPI_Comm comm, Mat seqmat, PetscInt m, PetscInt n, MatReuse scall, Mat *mpimat) in MatCreateMPIAIJSumSeqAIJ() argument 5036 if (scall == MAT_INITIAL_MATRIX) { in MatCreateMPIAIJSumSeqAIJ() 5045 …if (scall == MAT_INITIAL_MATRIX) PetscCall(MatCreateMPIAIJSumSeqAIJSymbolic(comm, seqmat, m, n, mp… in MatCreateMPIAIJSumSeqAIJ() 5120 PetscErrorCode MatMPIAIJGetLocalMat(Mat A, MatReuse scall, Mat *A_loc) in MatMPIAIJGetLocalMat() argument 5137 if (scall == MAT_INITIAL_MATRIX) { in MatMPIAIJGetLocalMat() 5140 } else if (scall == MAT_REUSE_MATRIX) { in MatMPIAIJGetLocalMat() 5157 if (scall == MAT_INITIAL_MATRIX) { in MatMPIAIJGetLocalMat() 5194 } else if (scall == MAT_REUSE_MATRIX) { in MatMPIAIJGetLocalMat() [all …]
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpb_baij.c | 3 PetscErrorCode MatGetMultiProcBlock_MPIBAIJ(Mat mat, MPI_Comm subComm, MatReuse scall, Mat *subMat) in MatGetMultiProcBlock_MPIBAIJ() argument 18 if (scall == MAT_INITIAL_MATRIX) { in MatGetMultiProcBlock_MPIBAIJ() 52 if (scall == MAT_INITIAL_MATRIX) { in MatGetMultiProcBlock_MPIBAIJ() 108 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscFree(nnz)); in MatGetMultiProcBlock_MPIBAIJ()
|
| H A D | baijov.c | 524 …s_MPIBAIJ(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_MPIBAIJ() argument 563 if (scall == MAT_INITIAL_MATRIX) { in MatCreateSubMatrices_MPIBAIJ() 587 …ubMatrices_MPIBAIJ_local(C, max_no, isrow_block + pos, iscol_block + pos, scall, *submat + pos, sy… in MatCreateSubMatrices_MPIBAIJ() 589 if (scall == MAT_INITIAL_MATRIX) { /* submat[pos] is a dummy matrix */ in MatCreateSubMatrices_MPIBAIJ() 597 if (scall == MAT_INITIAL_MATRIX && ismax) { in MatCreateSubMatrices_MPIBAIJ() 613 …local(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, Mat *submats, Pet… argument 686 if (scall == MAT_REUSE_MATRIX) {
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mmdense.c | 25 …_MPIDense(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_MPIDense() argument 31 if (scall != MAT_REUSE_MATRIX) PetscCall(PetscCalloc1(ismax + 1, submat)); in MatCreateSubMatrices_MPIDense() 44 …PetscCall(MatCreateSubMatrices_MPIDense_Local(C, max_no, isrow + pos, iscol + pos, scall, *submat … in MatCreateSubMatrices_MPIDense() 50 …ense_Local(Mat C, PetscInt ismax, const IS isrow[], const IS iscol[], MatReuse scall, Mat *submats) in MatCreateSubMatrices_MPIDense_Local() argument 274 if (scall == MAT_REUSE_MATRIX) { in MatCreateSubMatrices_MPIDense_Local()
|
| /petsc/src/tao/matrix/ |
| H A D | submatfree.c | 222 PetscErrorCode MatCreateSubMatrices_SMF(Mat A, PetscInt n, IS *irow, IS *icol, MatReuse scall, Mat … in MatCreateSubMatrices_SMF() argument 227 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscCalloc1(n + 1, B)); in MatCreateSubMatrices_SMF() 229 for (i = 0; i < n; i++) PetscCall(MatCreateSubMatrix_SMF(A, irow[i], icol[i], scall, &(*B)[i])); in MatCreateSubMatrices_SMF()
|
| H A D | adamat.c | 233 …cErrorCode MatCreateSubMatrices_ADA(Mat A, PetscInt n, IS *irow, IS *icol, MatReuse scall, Mat **B) in MatCreateSubMatrices_ADA() argument 238 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscCalloc1(n + 1, B)); in MatCreateSubMatrices_ADA() 239 for (i = 0; i < n; i++) PetscCall(MatCreateSubMatrix_ADA(A, irow[i], icol[i], scall, &(*B)[i])); in MatCreateSubMatrices_ADA()
|
| /petsc/src/ksp/pc/impls/fieldsplit/ |
| H A D | fieldsplit.c | 731 MatReuse scall; in PCSetUp_FieldSplit() local 737 scall = MAT_INITIAL_MATRIX; in PCSetUp_FieldSplit() 738 } else scall = MAT_REUSE_MATRIX; in PCSetUp_FieldSplit() 745 …if (!pmat) PetscCall(MatCreateSubMatrix(pc->pmat, ilink->is, ilink->is_col, scall, &jac->pmat[i])); in PCSetUp_FieldSplit() 759 MatReuse scall; in PCSetUp_FieldSplit() local 765 scall = MAT_INITIAL_MATRIX; in PCSetUp_FieldSplit() 766 } else scall = MAT_REUSE_MATRIX; in PCSetUp_FieldSplit() 769 PetscCall(MatCreateSubMatrix(pc->mat, ilink->is, ilink->is_col, scall, &jac->mat[i])); in PCSetUp_FieldSplit() 802 MatReuse scall; in PCSetUp_FieldSplit() local 806 scall = MAT_INITIAL_MATRIX; in PCSetUp_FieldSplit() [all …]
|
| /petsc/src/mat/impls/adj/mpi/ |
| H A D | mpiadj.c | 133 …at, PetscInt n, const IS irow[], const IS icol[], PetscBool subcomm, MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_MPIAdj_Private() argument 189 if (scall == MAT_INITIAL_MATRIX) { in MatCreateSubMatrices_MPIAdj_Private() 209 …cesMPI_MPIAdj(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatricesMPI_MPIAdj() argument 213 PetscCall(MatCreateSubMatrices_MPIAdj_Private(mat, n, irow, icol, PETSC_TRUE, scall, submat)); in MatCreateSubMatricesMPI_MPIAdj() 217 …trices_MPIAdj(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_MPIAdj() argument 221 PetscCall(MatCreateSubMatrices_MPIAdj_Private(mat, n, irow, icol, PETSC_FALSE, scall, submat)); in MatCreateSubMatrices_MPIAdj()
|
| /petsc/src/ksp/pc/impls/asm/ |
| H A D | asm.c | 158 MatReuse scall = MAT_REUSE_MATRIX; in PCSetUp_ASM() local 274 scall = MAT_INITIAL_MATRIX; in PCSetUp_ASM() 282 scall = MAT_INITIAL_MATRIX; in PCSetUp_ASM() 287 if (scall == MAT_REUSE_MATRIX && osm->sub_mat_type) { in PCSetUp_ASM() 290 scall = MAT_INITIAL_MATRIX; in PCSetUp_ASM() 296 PetscCall(MatCreateSubMatrices(pc->pmat, osm->n_local_true, osm->is, osm->is, scall, &osm->pmat)); in PCSetUp_ASM() 297 if (scall == MAT_INITIAL_MATRIX) { in PCSetUp_ASM() 393 PetscCall(MatCreateSubMatrices(pc->pmat, osm->n_local_true, osm->is, cis, scall, &osm->lmats)); in PCSetUp_ASM()
|
| /petsc/src/mat/impls/normal/ |
| H A D | normmh.c | 9 …rmalHermitian(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_NormalHermitian() argument 22 if (scall != MAT_REUSE_MATRIX) PetscCall(PetscCalloc1(n, submat)); in MatCreateSubMatrices_NormalHermitian()
|
| H A D | normm.c | 26 …trices_Normal(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_Normal() argument 39 if (scall != MAT_REUSE_MATRIX) PetscCall(PetscCalloc1(n, submat)); in MatCreateSubMatrices_Normal()
|
| /petsc/src/mat/interface/ |
| H A D | matrix.c | 7258 …teSubMatrices(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices() argument 7273 if (n && scall == MAT_REUSE_MATRIX) { in MatCreateSubMatrices() 7281 PetscUseTypeMethod(mat, createsubmatrices, n, irow, icol, scall, submat); in MatCreateSubMatrices() 7319 …ubMatricesMPI(Mat mat, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatricesMPI() argument 7334 if (n && scall == MAT_REUSE_MATRIX) { in MatCreateSubMatricesMPI() 7343 PetscUseTypeMethod(mat, createsubmatricesmpi, n, irow, icol, scall, submat); in MatCreateSubMatricesMPI() 10117 PetscErrorCode MatPtAP(Mat A, Mat P, MatReuse scall, PetscReal fill, Mat *C) in MatPtAP() argument 10120 if (scall == MAT_REUSE_MATRIX) MatCheckProduct(*C, 5); in MatPtAP() 10121 …PetscCheck(scall != MAT_INPLACE_MATRIX, PetscObjectComm((PetscObject)A), PETSC_ERR_SUP, "Inplace p… in MatPtAP() 10123 if (scall == MAT_INITIAL_MATRIX) { in MatPtAP() [all …]
|
| /petsc/src/ksp/pc/impls/bjacobi/ |
| H A D | bjacobi.c | 1022 MatReuse scall; in PCSetUp_BJacobi_Multiblock() local 1040 scall = MAT_INITIAL_MATRIX; in PCSetUp_BJacobi_Multiblock() 1117 scall = MAT_INITIAL_MATRIX; in PCSetUp_BJacobi_Multiblock() 1118 } else scall = MAT_REUSE_MATRIX; in PCSetUp_BJacobi_Multiblock() 1121 PetscCall(MatCreateSubMatrices(pmat, n_local, bjac->is, bjac->is, scall, &bjac->pmat)); in PCSetUp_BJacobi_Multiblock() 1124 PetscCall(MatCreateSubMatrices(mat, n_local, bjac->is, bjac->is, scall, &bjac->mat)); in PCSetUp_BJacobi_Multiblock()
|
| /petsc/src/mat/impls/cdiagonal/ |
| H A D | cdiagonal.c | 88 …onstantDiagonal(Mat A, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_ConstantDiagonal() argument 94 PetscCall(MatCreateSubMatrices(B, n, irow, icol, scall, submat)); in MatCreateSubMatrices_ConstantDiagonal()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matrart.c | 331 PetscErrorCode MatRARt_SeqAIJ_SeqAIJ(Mat A, Mat R, MatReuse scall, PetscReal fill, Mat *C) in MatRARt_SeqAIJ_SeqAIJ() argument 337 if (scall == MAT_INITIAL_MATRIX) { in MatRARt_SeqAIJ_SeqAIJ()
|
| /petsc/src/ksp/pc/impls/gasm/ |
| H A D | gasm.c | 304 MatReuse scall = MAT_REUSE_MATRIX; in PCSetUp_GASM() local 539 scall = MAT_INITIAL_MATRIX; in PCSetUp_GASM() 546 scall = MAT_INITIAL_MATRIX; in PCSetUp_GASM() 549 …PetscCall(MatCreateSubMatrix(pc->pmat, osm->permutationIS, osm->permutationIS, scall, &osm->permut… in PCSetUp_GASM() 560 PetscCall(MatCreateSubMatricesMPI(pc->pmat, osm->n, osm->ois, osm->ois, scall, &osm->pmat)); in PCSetUp_GASM() 562 PetscCall(MatCreateSubMatrices(pc->pmat, osm->n, osm->ois, osm->ois, scall, &osm->pmat)); in PCSetUp_GASM() 564 if (scall == MAT_INITIAL_MATRIX) { in PCSetUp_GASM()
|
| /petsc/src/ksp/pc/impls/ml/ |
| H A D | ml.c | 169 static PetscErrorCode MatConvert_MPIAIJ_ML(Mat A, MatType newtype, MatReuse scall, Mat *Aloc) in MatConvert_MPIAIJ_ML() argument 182 if (scall == MAT_INITIAL_MATRIX) { in MatConvert_MPIAIJ_ML() 218 } else if (scall == MAT_REUSE_MATRIX) { in MatConvert_MPIAIJ_ML() 231 …e SETERRQ(PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid MatReuse %d", (int)scall); in MatConvert_MPIAIJ_ML()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 105 static PetscErrorCode MatCreateSubMatrix_SeqSBAIJ_Private(Mat A, IS isrow, IS iscol, MatReuse scall… in MatCreateSubMatrix_SeqSBAIJ_Private() argument 138 if (scall == MAT_REUSE_MATRIX) { in MatCreateSubMatrix_SeqSBAIJ_Private() 227 PetscErrorCode MatCreateSubMatrix_SeqSBAIJ(Mat A, IS isrow, IS iscol, MatReuse scall, Mat *B) in MatCreateSubMatrix_SeqSBAIJ() argument 251 if (sym) PetscCall(MatCreateSubMatrix_SeqSBAIJ_Private(A, is1, is2, scall, B, sym)); in MatCreateSubMatrix_SeqSBAIJ() 257 …PetscCheck(scall != MAT_INPLACE_MATRIX, PETSC_COMM_SELF, PETSC_ERR_SUP, "MAT_INPLACE_MATRIX not su… in MatCreateSubMatrix_SeqSBAIJ() 258 if (scall == MAT_REUSE_MATRIX) PetscCall(MatCopy(C[0], *B, SAME_NONZERO_PATTERN)); in MatCreateSubMatrix_SeqSBAIJ() 275 …SubMatrices_SeqSBAIJ(Mat A, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *B[]) in MatCreateSubMatrices_SeqSBAIJ() argument 280 if (scall == MAT_INITIAL_MATRIX) PetscCall(PetscCalloc1(n, B)); in MatCreateSubMatrices_SeqSBAIJ() 282 …for (i = 0; i < n; i++) PetscCall(MatCreateSubMatrix_SeqSBAIJ(A, irow[i], icol[i], scall, &(*B)[i]… in MatCreateSubMatrices_SeqSBAIJ()
|
| /petsc/src/mat/impls/aij/mpi/mpihipsparse/ |
| H A D | mpiaijhipsparse.hip.cxx | 240 static PetscErrorCode MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE(Mat A, MatReuse scall, IS *glob, Ma… in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE() argument 247 PetscCall(MatSeqAIJHIPSPARSEMergeMats(Ad, Ao, scall, A_loc)); in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE()
|
| /petsc/src/mat/impls/aij/mpi/mpicusparse/ |
| H A D | mpiaijcusparse.cu | 240 static PetscErrorCode MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE(Mat A, MatReuse scall, IS *glob, Mat… in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE() argument 247 PetscCall(MatSeqAIJCUSPARSEMergeMats(Ad, Ao, scall, A_loc)); in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE()
|
| /petsc/src/mat/impls/htool/ |
| H A D | htool.cxx | 158 …bMatrices_Htool(Mat A, PetscInt n, const IS irow[], const IS icol[], MatReuse scall, Mat *submat[]) in MatCreateSubMatrices_Htool() argument 171 if (scall != MAT_REUSE_MATRIX) PetscCall(PetscCalloc1(n, submat)); in MatCreateSubMatrices_Htool() 177 …if (scall != MAT_REUSE_MATRIX) PetscCall(MatCreateDense(PETSC_COMM_SELF, nrow, m, nrow, m, nullptr… in MatCreateSubMatrices_Htool()
|
| /petsc/src/mat/impls/blockmat/seq/ |
| H A D | blockmat.c | 537 static PetscErrorCode MatCreateSubMatrix_BlockMat(Mat A, IS isrow, IS iscol, MatReuse scall, Mat *B) in MatCreateSubMatrix_BlockMat() argument 559 if (scall == MAT_REUSE_MATRIX) { in MatCreateSubMatrix_BlockMat()
|