| /petsc/src/dm/tutorials/ |
| H A D | swarm_ex1.c | 16 DM dms; in ex1_1() local 26 PetscCall(DMCreate(PETSC_COMM_WORLD, &dms)); in ex1_1() 27 PetscCall(DMSetType(dms, DMSWARM)); in ex1_1() 28 PetscCall(PetscObjectSetName((PetscObject)dms, "Particles")); in ex1_1() 30 PetscCall(DMSwarmInitializeFieldRegister(dms)); in ex1_1() 31 PetscCall(DMSwarmRegisterPetscDatatypeField(dms, "viscosity", 1, PETSC_REAL)); in ex1_1() 32 PetscCall(DMSwarmRegisterPetscDatatypeField(dms, "strain", 1, PETSC_REAL)); in ex1_1() 33 PetscCall(DMSwarmFinalizeFieldRegister(dms)); in ex1_1() 34 PetscCall(DMSwarmSetLocalSizes(dms, 5 + rank, 4)); in ex1_1() 35 PetscCall(DMView(dms, PETSC_VIEWER_STDOUT_WORLD)); in ex1_1() [all …]
|
| H A D | swarm_ex2.c | 12 DM dms; in ex2_1() local 20 PetscCall(DMCreate(PETSC_COMM_WORLD, &dms)); in ex2_1() 21 PetscCall(DMSetType(dms, DMSWARM)); in ex2_1() 22 PetscCall(PetscObjectSetName((PetscObject)dms, "Particles")); in ex2_1() 23 PetscCall(DMSwarmInitializeFieldRegister(dms)); in ex2_1() 24 PetscCall(DMSwarmRegisterPetscDatatypeField(dms, "viscosity", 1, PETSC_REAL)); in ex2_1() 25 PetscCall(DMSwarmRegisterPetscDatatypeField(dms, "strain", 3, PETSC_REAL)); in ex2_1() 26 PetscCall(DMSwarmFinalizeFieldRegister(dms)); in ex2_1() 27 PetscCall(DMSwarmSetLocalSizes(dms, 5 + rank, 4)); in ex2_1() 28 PetscCall(DMView(dms, PETSC_VIEWER_STDOUT_WORLD)); in ex2_1() [all …]
|
| H A D | swarm_ex3.c | 99 PetscErrorCode SwarmViewGP(DM dms, const char prefix[]) in SwarmViewGP() argument 113 PetscCall(DMSwarmGetLocalSize(dms, &npoints)); in SwarmViewGP() 114 PetscCall(DMSwarmGetField(dms, DMSwarmPICField_coor, &bs, NULL, (void **)&array)); in SwarmViewGP() 115 PetscCall(DMSwarmGetField(dms, "itag", NULL, NULL, (void **)&iarray)); in SwarmViewGP() 117 PetscCall(DMSwarmRestoreField(dms, "itag", NULL, NULL, (void **)&iarray)); in SwarmViewGP() 118 PetscCall(DMSwarmRestoreField(dms, DMSwarmPICField_coor, &bs, NULL, (void **)&array)); in SwarmViewGP() 129 DM dms, dmcell, dmregular; in ex3_1() local 158 PetscCall(DMCreate(PETSC_COMM_WORLD, &dms)); in ex3_1() 159 PetscCall(DMSetType(dms, DMSWARM)); in ex3_1() 160 PetscCall(DMSetDimension(dms, 2)); in ex3_1() [all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex3.c | 10 DM dms[NDMS]; in main() local 17 …DECIDE, PETSC_DECIDE, PETSC_DECIDE, 2, 3, 4, 5, DMSTAG_STENCIL_BOX, 1, NULL, NULL, NULL, &dms[0])); in main() 18 PetscCall(DMSetUp(dms[0])); in main() 19 …DECIDE, PETSC_DECIDE, PETSC_DECIDE, 2, 3, 4, 5, DMSTAG_STENCIL_BOX, 1, NULL, NULL, NULL, &dms[1])); in main() 20 PetscCall(DMSetUp(dms[1])); in main() 23 PetscCall(DMStagCreateCompatibleDMStag(dms[0], 0, 1, 0, 1, &dms[2])); in main() 26 …DECIDE, PETSC_DECIDE, PETSC_DECIDE, 2, 3, 4, 5, DMSTAG_STENCIL_BOX, 2, NULL, NULL, NULL, &dms[3])); in main() 31 PetscCall(DMGetCompatibility(dms[i], dms[i], &compatible, &set)); in main() 32 …PetscCheck(set && compatible, PetscObjectComm((PetscObject)dms[i]), PETSC_ERR_PLIB, "DM %" PetscIn… in main() 38 PetscCall(DMGetCompatibility(dms[0], dms[i], &compatible, &set)); in main() [all …]
|
| /petsc/src/dm/interface/ftn-custom/ |
| H A D | zdmf.c | 17 PETSC_EXTERN void dmcreatesuperdm_(DM dms[], PetscInt *len, IS ***is, DM *superdm, PetscErrorCode *… in dmcreatesuperdm_() argument 19 *ierr = DMCreateSuperDM(dms, *len, *is, superdm); in dmcreatesuperdm_() 35 PETSC_EXTERN void dmcreatefielddecompositiongetisdm_(DM *dm, F90Array1d *iss, F90Array1d *dms, Pets… in dmcreatefielddecompositiongetisdm_() argument 41 if (iss && dms) { in dmcreatefielddecompositiongetisdm_() 45 } else if (dms) { in dmcreatefielddecompositiongetisdm_() 51 if (dms) *ierr = F90Array1dCreate(tdm, MPIU_FORTRANADDR, 1, n, dms PETSC_F90_2PTR_PARAM(ptrd2)); in dmcreatefielddecompositiongetisdm_() 54 PETSC_EXTERN void dmcreatefielddecompositionrestoreisdm_(DM *dm, F90Array1d *iss, F90Array1d *dms, … in dmcreatefielddecompositionrestoreisdm_() argument 70 if (dms) { in dmcreatefielddecompositionrestoreisdm_() 72 *ierr = F90Array1dAccess(dms, MPIU_FORTRANADDR, (void **)&tdm PETSC_F90_2PTR_PARAM(ptrd2)); in dmcreatefielddecompositionrestoreisdm_() 74 *ierr = F90Array1dDestroy(dms, MPIU_FORTRANADDR PETSC_F90_2PTR_PARAM(ptrd2)); in dmcreatefielddecompositionrestoreisdm_()
|
| /petsc/src/dm/ftn-mod/ |
| H A D | petscdmmod.F90 | 74 subroutine DMCreateFieldDecompositionGetISDM(dm, iss, dms, ierr) 80 DM, pointer :: dms(:) 86 subroutine DMCreateFieldDecompositionRestoreISDM(dm, iss, dms, ierr) 92 DM, pointer :: dms(:) 143 subroutine DMCreateFieldDecomposition(dm, n, names, iss, dms, ierr) argument 147 DM, pointer :: dms(:) 157 call DMCreateFieldDecompositionGetISDM(dm, iss, dms, ierr) 160 subroutine DMDestroyFieldDecomposition(dm, n, names, iss, dms, ierr) argument 164 DM, pointer :: dms(:) 171 call DMCreateFieldDecompositionRestoreISDM(dm, iss, dms, ierr)
|
| /petsc/src/snes/interface/ |
| H A D | snesj2.c | 69 DMSNES dms; in SNESComputeJacobianDefaultColor() local 91 PetscCall(DMGetDMSNES(dm, &dms)); in SNESComputeJacobianDefaultColor() 92 if (dms->ops->computemffunction) { in SNESComputeJacobianDefaultColor() 146 DMSNES dms; in SNESPruneJacobianColor() local 163 PetscCall(DMGetDMSNES(dm, &dms)); in SNESPruneJacobianColor()
|
| H A D | snesj.c | 62 DMSNES dms; in SNESComputeJacobianDefault() local 91 PetscCall(DMGetDMSNES(dm, &dms)); in SNESComputeJacobianDefault() 92 if (dms->ops->computemffunction) { in SNESComputeJacobianDefault() 129 if (dms->ops->computemffunction) { in SNESComputeJacobianDefault()
|
| /petsc/src/snes/mf/ |
| H A D | snesmfj.c | 86 DMSNES dms; in MatAssemblyEnd_SNESMF() local 95 PetscCall(DMGetDMSNES(dm, &dms)); in MatAssemblyEnd_SNESMF() 96 …if ((j->func == (PetscErrorCode (*)(void *, Vec, Vec))SNESComputeFunction) && !dms->ops->computemf… in MatAssemblyEnd_SNESMF() 287 DMSNES dms; in MatCreateSNESMF() local 290 PetscCall(DMGetDMSNES(dm, &dms)); in MatCreateSNESMF() 291 …PetscCall(MatMFFDSetFunction(*J, (PetscErrorCode (*)(void *, Vec, Vec))(dms->ops->computemffunctio… in MatCreateSNESMF()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | petscdmshell.pxi | 341 names, ises, dms = decomp(Dm, *args, **kargs) 348 elif dms is not None: 349 clen[0] = <PetscInt>len(dms) 365 if dmlist != NULL and dms is not None: 366 CHKERR(PetscMalloc(len(dms)*sizeof(PetscDM), dmlist)) 367 for i in range(len(dms)): 368 dmlist[0][i] = (<DM?>dms[i]).dm 369 CHKERR(PetscINCREF((<DM?>dms[i]).obj)) 387 names, innerises, outerises, dms = decomp(Dm, *args, **kargs) 396 elif dms is not None: [all …]
|
| H A D | DM.pyx | 2345 cdef list dms = [] 2354 dms.append(dm) 2357 dms.append(None) 2372 return (names, isets, dms)
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex73.c | 534 PetscErrorCode DMShellCreate_ShellDA(DM da, DM *dms) in DMShellCreate_ShellDA() argument 538 PetscCall(DMShellCreate(PetscObjectComm((PetscObject)da), dms)); in DMShellCreate_ShellDA() 539 PetscCall(DMShellSetContext(*dms, da)); in DMShellCreate_ShellDA() 540 PetscCall(DMShellSetCreateGlobalVector(*dms, DMCreateGlobalVector_ShellDA)); in DMShellCreate_ShellDA() 541 PetscCall(DMShellSetCreateLocalVector(*dms, DMCreateLocalVector_ShellDA)); in DMShellCreate_ShellDA() 542 PetscCall(DMShellSetCreateMatrix(*dms, DMCreateMatrix_ShellDA)); in DMShellCreate_ShellDA() 543 PetscCall(DMShellSetCoarsen(*dms, DMCoarsen_ShellDA)); in DMShellCreate_ShellDA() 544 PetscCall(DMShellSetCreateInterpolation(*dms, DMCreateInterpolation_ShellDA)); in DMShellCreate_ShellDA() 546 *dms = NULL; in DMShellCreate_ShellDA() 775 DM *dms; in test_hierarchy() local [all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_dmshell.py | 94 dms = [dm.clone() for _ in range(len(ises))] 95 return _, ises, dms 99 dms = [dm.clone() for _ in range(len(ises))] 100 return names, ises, dms 110 names, ises, dms = self.dm.createFieldDecomposition() 112 self.assertEqual(len(dms), len(ises)) 119 for dm in dms:
|
| /petsc/src/dm/interface/ |
| H A D | dmi.c | 448 PetscErrorCode DMCreateSectionSuperDM(DM dms[], PetscInt len, IS *is[], DM *superdm) in DMCreateSectionSuperDM() argument 456 PetscCall(PetscObjectGetComm((PetscObject)dms[0], &comm)); in DMCreateSectionSuperDM() 460 PetscCall(DMGetLocalSection(dms[i], §ions[i])); in DMCreateSectionSuperDM() 461 PetscCall(DMGetGlobalSection(dms[i], §ionGlobals[i])); in DMCreateSectionSuperDM() 520 if (len && dms[0]->probs) { in DMCreateSectionSuperDM() 526 PetscCall(DMGetField(dms[i], f, NULL, &disc)); in DMCreateSectionSuperDM() 538 if (dms[i]->nullspaceConstructors) { in DMCreateSectionSuperDM() 539 (*superdm)->nullspaceConstructors[supf] = dms[i]->nullspaceConstructors[f]; in DMCreateSectionSuperDM()
|
| H A D | dm.c | 2200 PetscErrorCode DMCreateSuperDM(DM dms[], PetscInt n, IS *is[], DM *superdm) in DMCreateSuperDM() argument 2205 PetscAssertPointer(dms, 1); in DMCreateSuperDM() 2206 for (i = 0; i < n; ++i) PetscValidHeaderSpecific(dms[i], DM_CLASSID, 1); in DMCreateSuperDM() 2211 DM dm = dms[0]; in DMCreateSuperDM() 2213 PetscCall((*dm->ops->createsuperdm)(dms, n, is, superdm)); in DMCreateSuperDM()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexdd.c | 47 …nDecomposition_Plex(DM dm, PetscInt *nsub, char ***names, IS **innerises, IS **outerises, DM **dms) in DMCreateDomainDecomposition_Plex() argument 67 if (dms) PetscCall(PetscCalloc1(n, dms)); in DMCreateDomainDecomposition_Plex() 246 if (dms) (*dms)[0] = odm; in DMCreateDomainDecomposition_Plex()
|
| H A D | plexcreate.c | 5408 DM *dms, coarseDM; in DMSetFromOptions_Plex() local 5412 PetscCall(PetscMalloc1(refine, &dms)); in DMSetFromOptions_Plex() 5413 PetscCall(DMRefineHierarchy(dm, refine, dms)); in DMSetFromOptions_Plex() 5415 PetscCall(DMPlexSwap_Static(dm, dms[refine - 1])); in DMSetFromOptions_Plex() 5417 PetscCall(DMSetCoarseDM(dm, dms[0])); in DMSetFromOptions_Plex() 5420 PetscCall(DMSetCoarseDM(dm, dms[refine - 2])); in DMSetFromOptions_Plex() 5422 PetscCall(DMSetCoarseDM(dms[0], dms[refine - 1])); in DMSetFromOptions_Plex() 5423 PetscCall(DMPlexSetRegularRefinement(dms[0], PETSC_TRUE)); in DMSetFromOptions_Plex() 5425 PetscCall(DMSetCoarseDM(dms[refine - 1], coarseDM)); in DMSetFromOptions_Plex() 5429 PetscCall(DMSetFromOptions_NonRefinement_Plex(dms[r], PetscOptionsObject)); in DMSetFromOptions_Plex() [all …]
|
| /petsc/src/ksp/pc/impls/mg/ |
| H A D | mg.c | 870 DM *dms; in PCSetUp_MG() local 996 PetscCall(PetscMalloc1(n, &dms)); in PCSetUp_MG() 997 dms[n - 1] = pc->dm; in PCSetUp_MG() 999 for (i = n - 2; i > -1; i--) PetscCall(DMCoarsen(dms[i + 1], MPI_COMM_NULL, &dms[i])); in PCSetUp_MG() 1003 PetscCall(KSPSetDM(mglevels[i]->smoothd, dms[i])); in PCSetUp_MG() 1007 PetscCall(KSPSetDM(mglevels[i]->smoothu, dms[i])); in PCSetUp_MG() 1012 PetscCall(KSPSetDM(mglevels[i]->cr, dms[i])); in PCSetUp_MG() 1017 PetscCall(DMCreateInterpolation(dms[i], dms[i + 1], &p, &rscale)); in PCSetUp_MG() 1023 PetscCall(DMHasCreateRestriction(dms[i], &dmhasrestrict)); in PCSetUp_MG() 1025 PetscCall(DMCreateRestriction(dms[i], dms[i + 1], &p)); in PCSetUp_MG() [all …]
|
| /petsc/src/dm/impls/composite/ |
| H A D | pack.c | 1067 DM *dms; in DMCreateFieldIS_Composite() local 1075 PetscCall(PetscMalloc1(nDM, &dms)); in DMCreateFieldIS_Composite() 1077 PetscCall(DMCompositeGetEntriesArray(dm, dms)); in DMCreateFieldIS_Composite() 1085 PetscCall(PetscObjectGetOptionsPrefix((PetscObject)dms[i], &splitname)); in DMCreateFieldIS_Composite() 1101 PetscCall(PetscFree(dms)); in DMCreateFieldIS_Composite() 1271 PetscErrorCode DMCompositeGetEntriesArray(DM dm, DM dms[]) in DMCompositeGetEntriesArray() argument 1283 for (next = com->next, i = 0; next; next = next->next, i++) dms[i] = next->dm; in DMCompositeGetEntriesArray() 1333 DM dm = (DM)odm, *dms; in DMSetUpGLVisViewer_Composite() local 1344 PetscCall(PetscMalloc1(n, &dms)); in DMSetUpGLVisViewer_Composite() 1345 PetscCall(DMCompositeGetEntriesArray(dm, dms)); in DMSetUpGLVisViewer_Composite() [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex13.c | 380 DM dmErr, dmErrAux, dms[2]; in ComputeAdjoint() local 404 dms[0] = dm; in ComputeAdjoint() 405 dms[1] = dm; in ComputeAdjoint() 406 PetscCall(DMCreateSuperDM(dms, 2, &subis, &dmErrAux)); in ComputeAdjoint() 410 PetscCall(DMGetLocalSection(dms[0], &sec)); in ComputeAdjoint() 412 PetscCall(DMGetLocalSection(dms[1], &sec)); in ComputeAdjoint()
|
| /petsc/doc/manual/ |
| H A D | dmplex.md | 480 group /topologies/plexA/dms 481 group /topologies/plexA/dms/dmA 482 dataset /topologies/plexA/dms/dmA/order 483 group /topologies/plexA/dms/dmA/section 484 dataset /topologies/plexA/dms/dmA/section/atlasDof 485 dataset /topologies/plexA/dms/dmA/section/atlasOff 486 group /topologies/plexA/dms/dmA/vecs 487 group /topologies/plexA/dms/dmA/vecs/vecA 488 dataset /topologies/plexA/dms/dmA/vecs/vecA/vecA
|
| /petsc/src/snes/impls/patch/ |
| H A D | snespatch.c | 373 PetscErrorCode SNESPatchSetDiscretisationInfo(SNES snes, PetscInt nsubspaces, DM *dms, PetscInt *bs… in SNESPatchSetDiscretisationInfo() argument 382 …PetscCall(PCPatchSetDiscretisationInfo(patch->pc, nsubspaces, dms, bs, nodesPerCell, cellNodeMap, … in SNESPatchSetDiscretisationInfo()
|
| /petsc/src/ksp/pc/impls/fieldsplit/ |
| H A D | fieldsplit.c | 412 DM *dms; in PCFieldSplitSetDefaults() local 416 PetscCall(DMCreateFieldDecomposition(pc->dm, &numFields, &fieldNames, &fields, &dms)); in PCFieldSplitSetDefaults() 449 for (j = 0; j < numFields; j++) PetscCall(DMDestroy(dms + j)); in PCFieldSplitSetDefaults() 450 PetscCall(PetscFree(dms)); in PCFieldSplitSetDefaults() 451 PetscCall(PetscMalloc1(i, &dms)); in PCFieldSplitSetDefaults() 452 for (j = 0; j < i; ++j) dms[j] = subdm[j]; in PCFieldSplitSetDefaults() 456 if (dms) { in PCFieldSplitSetDefaults() 461 PetscCall(PetscObjectSetOptionsPrefix((PetscObject)dms[i], prefix)); in PCFieldSplitSetDefaults() 462 PetscCall(KSPSetDM(ilink->ksp, dms[i])); in PCFieldSplitSetDefaults() 464 PetscCall(PetscObjectIncrementTabLevel((PetscObject)dms[i], (PetscObject)ilink->ksp, 0)); in PCFieldSplitSetDefaults() [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex30.c | 1612 DM *dms; in CreateMesh() local 1627 PetscCall(PetscCalloc1(refine, &dms)); in CreateMesh() 1628 if (isHierarchy) PetscCall(DMRefineHierarchy(*dm, refine, dms)); in CreateMesh() 1631 DM dmr = dms[r]; in CreateMesh() 1650 PetscCall(PetscFree(dms)); in CreateMesh()
|
| /petsc/src/ksp/pc/impls/patch/ |
| H A D | pcpatch.c | 504 PetscErrorCode PCPatchSetDiscretisationInfo(PC pc, PetscInt nsubspaces, DM *dms, PetscInt *bs, Pets… in PCPatchSetDiscretisationInfo() argument 526 PetscCall(DMGetLocalSection(dms[i], &patch->dofSection[i])); in PCPatchSetDiscretisationInfo() 528 PetscCall(DMGetSectionSF(dms[i], &sfs[i])); in PCPatchSetDiscretisationInfo()
|