| /petsc/src/ksp/pc/impls/gamg/ |
| H A D | util.c | 23 PetscErrorCode PCGAMGGetDataWithGhosts(Mat Gmat, PetscInt data_sz, PetscReal data_in[], PetscInt *a… in PCGAMGGetDataWithGhosts() argument
33 PetscValidHeaderSpecific(Gmat, MAT_CLASSID, 1); in PCGAMGGetDataWithGhosts()
34 mpimat = (Mat_MPIAIJ *)Gmat->data; in PCGAMGGetDataWithGhosts()
35 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Gmat, MATMPIAIJ, &isMPIAIJ)); in PCGAMGGetDataWithGhosts()
36 PetscCall(MatGetOwnershipRange(Gmat, &my0, &Iend)); in PCGAMGGetDataWithGhosts()
41 PetscCall(MatCreateVecs(Gmat, &tmp_crds, NULL)); in PCGAMGGetDataWithGhosts()
|
| H A D | geo.c | 535 Mat Gmat = *a_Gmat; in PCGAMGCoarsen_GEO() local
542 PetscCall(MatGetOwnershipRange(Gmat, &Istart, &Iend)); in PCGAMGCoarsen_GEO()
550 PetscCall(MatGetRow(Gmat, Ii, &ncols, NULL, NULL)); in PCGAMGCoarsen_GEO()
556 PetscCall(MatRestoreRow(Gmat, Ii, &ncols, NULL, NULL)); in PCGAMGCoarsen_GEO()
592 PetscCall(MatCoarsenSetAdjacency(crs, Gmat)); in PCGAMGCoarsen_GEO()
608 Mat Prol, Gmat; in PCGAMGProlongator_GEO() local
627 PetscCall(PetscCDGetMat(agg_lists, &Gmat)); // get auxiliary matrix for ghost edges in PCGAMGProlongator_GEO()
635 PetscCall(MatGetRow(Gmat, lid + my0, &ncols, NULL, NULL)); in PCGAMGProlongator_GEO()
637 PetscCall(MatRestoreRow(Gmat, lid + my0, &ncols, NULL, NULL)); in PCGAMGProlongator_GEO()
671 …PetscCall(getGIDsOnSquareGraph(pc, nLocalSelected, clid_flid, Gmat, &selected_2, &Gmat2, &crsGID)); in PCGAMGProlongator_GEO()
|
| H A D | agg.c | 724 Mat tGmat, Gmat = *a_Gmat; in PCGAMGCreateGraph_AGG() local
734 PetscCall(PetscObjectGetComm((PetscObject)Gmat, &comm)); in PCGAMGCreateGraph_AGG()
735 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Gmat, MATSEQAIJ, &isseqaij)); in PCGAMGCreateGraph_AGG()
736 PetscCall(MatGetType(Gmat, &jtype)); in PCGAMGCreateGraph_AGG()
745 PetscCall(MatGetSize(Gmat, &MM, &NN)); in PCGAMGCreateGraph_AGG()
746 PetscCall(MatGetOwnershipRange(Gmat, &Istart, &Iend)); in PCGAMGCreateGraph_AGG()
750 a = Gmat; in PCGAMGCreateGraph_AGG()
753 Mat_MPIAIJ *d = (Mat_MPIAIJ *)Gmat->data; in PCGAMGCreateGraph_AGG()
804 … PetscCall(MatPropagateSymmetryOptions(Gmat, tGmat)); /* Normal Mat options are not relevant ? */ in PCGAMGCreateGraph_AGG()
807 PetscCall(MatDestroy(&Gmat)); in PCGAMGCreateGraph_AGG()
[all …]
|
| H A D | gamg.c | 726 Mat Gmat, mat; in PCSetUp_GAMG() local
730 PetscCall(PCGAMGCreateGraph(pc, Aarr[level], &Gmat)); in PCSetUp_GAMG()
731 …PetscCall(pc_gamg->ops->coarsen(pc, &Gmat, &agg_lists)); // Gmat may have ghosts for QR aggregates… in PCSetUp_GAMG()
733 if (!mat) PetscCall(PetscCDSetMat(agg_lists, Gmat)); in PCSetUp_GAMG()
762 … if (mat == Gmat) PetscCall(PetscCDClearMat(agg_lists)); // take the Mat away from the list (yuck) in PCSetUp_GAMG()
772 PetscCall(MatCoarsenSetAdjacency(pc_gamg->asm_crs, Gmat)); in PCSetUp_GAMG()
791 … if (mat == Gmat) PetscCall(PetscCDClearMat(agg_lists)); // take the Mat away from the list (yuck) in PCSetUp_GAMG()
792 PetscCall(MatDestroy(&Gmat)); in PCSetUp_GAMG()
|
| /petsc/src/mat/graphops/coarsen/impls/mis/ |
| H A D | mis.c | 23 static PetscErrorCode MatCoarsenApply_MIS_private(IS perm, Mat Gmat, PetscBool strict_aggs, PetscCo… in MatCoarsenApply_MIS_private() argument
33 const PetscInt nloc = Gmat->rmap->n; in MatCoarsenApply_MIS_private()
40 PetscCall(PetscObjectGetComm((PetscObject)Gmat, &comm)); in MatCoarsenApply_MIS_private()
44 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Gmat, MATMPIAIJ, &isMPI)); in MatCoarsenApply_MIS_private()
46 mpimat = (Mat_MPIAIJ *)Gmat->data; in MatCoarsenApply_MIS_private()
50 …CompressedRow(mpimat->B, matB->nonzerorowcnt, &matB->compressedrow, matB->i, Gmat->rmap->n, -1.0)); in MatCoarsenApply_MIS_private()
52 matA = (Mat_SeqAIJ *)Gmat->data; in MatCoarsenApply_MIS_private()
53 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Gmat, MATSEQAIJ, &isAIJ)); in MatCoarsenApply_MIS_private()
56 PetscCall(MatGetOwnershipRange(Gmat, &my0, &Iend)); in MatCoarsenApply_MIS_private()
63 PetscCall(PetscSFCreate(PetscObjectComm((PetscObject)Gmat), &sf)); in MatCoarsenApply_MIS_private()
[all …]
|
| /petsc/src/mat/graphops/coarsen/impls/misk/ |
| H A D | misk.c | 46 static PetscErrorCode MatCoarsenApply_MISK_private(IS perm, const PetscInt misk, Mat Gmat, PetscCoa… in MatCoarsenApply_MISK_private() argument
56 PetscValidHeaderSpecific(Gmat, MAT_CLASSID, 3); in MatCoarsenApply_MISK_private()
59 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Gmat, MATMPIAIJ, &isMPI)); in MatCoarsenApply_MISK_private()
60 PetscCall(PetscObjectGetComm((PetscObject)Gmat, &comm)); in MatCoarsenApply_MISK_private()
63 PetscCall(PetscInfo(Gmat, "misk %" PetscInt_FMT "\n", misk)); in MatCoarsenApply_MISK_private()
65 if (misk > 1) PetscCall(MatDuplicate(Gmat, MAT_COPY_VALUES, &cMat)); in MatCoarsenApply_MISK_private()
66 else cMat = Gmat; in MatCoarsenApply_MISK_private()
221 …PetscCall(PetscInfo(Gmat, "\t removed %" PetscInt_FMT " of %" PetscInt_FMT " vertices. %" PetscIn… in MatCoarsenApply_MISK_private()
250 PetscCall(MatGetType(Gmat, &jtype)); in MatCoarsenApply_MISK_private()
282 } else if (cMat != Gmat) PetscCall(MatDestroy(&cMat)); in MatCoarsenApply_MISK_private()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 7737 Mat Gmat; in MatCreateGraph_Simple_AIJ() local
7758 PetscCall(MatCreate(comm, &Gmat)); in MatCreateGraph_Simple_AIJ()
7759 PetscCall(MatSetType(Gmat, jtype)); in MatCreateGraph_Simple_AIJ()
7760 PetscCall(MatSetSizes(Gmat, nloc, nloc, PETSC_DETERMINE, PETSC_DETERMINE)); in MatCreateGraph_Simple_AIJ()
7761 PetscCall(MatSetBlockSizes(Gmat, 1, 1)); in MatCreateGraph_Simple_AIJ()
7805 PetscCall(MatSeqAIJSetPreallocation(Gmat, 0, d_nnz)); in MatCreateGraph_Simple_AIJ()
7806 PetscCall(MatMPIAIJSetPreallocation(Gmat, 0, d_nnz, 0, o_nnz)); in MatCreateGraph_Simple_AIJ()
7840 PetscCall(MatSetValues(Gmat, 1, &grow, n / bs, AJ, AA, ADD_VALUES)); in MatCreateGraph_Simple_AIJ()
7883 PetscCall(MatSetValues(Gmat, 1, &grow, nc, AJ, AA, ADD_VALUES)); in MatCreateGraph_Simple_AIJ()
7886 PetscCall(MatAssemblyBegin(Gmat, MAT_FINAL_ASSEMBLY)); in MatCreateGraph_Simple_AIJ()
[all …]
|