Lines Matching refs:odm
6 static PetscErrorCode DMTransferMaterialParameters(DM dm, PetscSF sf, DM odm) in DMTransferMaterialParameters() argument
19 PetscCall(DMClone(odm, &odmAux)); in DMTransferMaterialParameters()
20 PetscCall(DMGetCoordinateDM(odm, &ocdm)); in DMTransferMaterialParameters()
41 PetscCall(DMSetAuxiliaryVec(odm, NULL, 0, 0, oA)); in DMTransferMaterialParameters()
49 DM odm; in DMCreateDomainDecomposition_Plex() local
74 …cCall(DMPlexDistributeOverlap_Internal(dm, ddovl + 1, PETSC_COMM_SELF, oname, &migrationSF, &odm)); in DMCreateDomainDecomposition_Plex()
75 if (!odm) PetscCall(DMClone(dm, &odm)); in DMCreateDomainDecomposition_Plex()
79 PetscCall(DMPlexSetMaxProjectionHeight(odm, mh)); in DMCreateDomainDecomposition_Plex()
86 PetscCall(DMGetDS(odm, &ds)); in DMCreateDomainDecomposition_Plex()
89 PetscCall(PetscSectionCreate(PetscObjectComm((PetscObject)odm), &tsec)); in DMCreateDomainDecomposition_Plex()
99 PetscCall(DMTransferMaterialParameters(dm, migrationSF, odm)); in DMCreateDomainDecomposition_Plex()
101 PetscCall(DMViewFromOptions(odm, (PetscObject)dm, "-dm_plex_dd_overlap_dm_view")); in DMCreateDomainDecomposition_Plex()
141 … PetscCall(PetscViewerVTKOpen(PetscObjectComm((PetscObject)odm), name, FILE_MODE_WRITE, &viewer)); in DMCreateDomainDecomposition_Plex()
142 PetscCall(DMView(odm, viewer)); in DMCreateDomainDecomposition_Plex()
163 PetscCall(PetscSectionCreate(PetscObjectComm((PetscObject)odm), &tsec)); in DMCreateDomainDecomposition_Plex()
175 PetscCall(DMGetLabel(odm, oname, &label)); in DMCreateDomainDecomposition_Plex()
177 PetscCall(DMPlexLabelComplete(odm, label)); in DMCreateDomainDecomposition_Plex()
178 PetscCall(DMGetLocalSection(odm, &tsec)); in DMCreateDomainDecomposition_Plex()
183 PetscCall(DMSetLocalSection(odm, sec)); in DMCreateDomainDecomposition_Plex()
185 PetscCall(DMRemoveLabel(odm, oname, NULL)); in DMCreateDomainDecomposition_Plex()
188 PetscCall(DMSetLocalSection(odm, sec)); in DMCreateDomainDecomposition_Plex()
192 PetscCall(DMGetSectionSF(odm, §ionSF)); in DMCreateDomainDecomposition_Plex()
193 PetscCall(DMGetLocalSection(odm, &olsec)); in DMCreateDomainDecomposition_Plex()
194 PetscCall(DMGetGlobalSection(odm, &ogsec)); in DMCreateDomainDecomposition_Plex()
213 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)odm), ni, lidxs, PETSC_OWN_POINTER, &li_is)… in DMCreateDomainDecomposition_Plex()
233 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)odm), c, lidxs, PETSC_OWN_POINTER, &ll_is)); in DMCreateDomainDecomposition_Plex()
236 PetscCall(PetscObjectCompose((PetscObject)odm, "__Plex_DD_IS_gi", (PetscObject)gi_is)); in DMCreateDomainDecomposition_Plex()
237 PetscCall(PetscObjectCompose((PetscObject)odm, "__Plex_DD_IS_li", (PetscObject)li_is)); in DMCreateDomainDecomposition_Plex()
238 PetscCall(PetscObjectCompose((PetscObject)odm, "__Plex_DD_IS_go", (PetscObject)go_is)); in DMCreateDomainDecomposition_Plex()
239 PetscCall(PetscObjectCompose((PetscObject)odm, "__Plex_DD_IS_gl", (PetscObject)gl_is)); in DMCreateDomainDecomposition_Plex()
240 PetscCall(PetscObjectCompose((PetscObject)odm, "__Plex_DD_IS_ll", (PetscObject)ll_is)); in DMCreateDomainDecomposition_Plex()
246 if (dms) (*dms)[0] = odm; in DMCreateDomainDecomposition_Plex()
247 else PetscCall(DMDestroy(&odm)); in DMCreateDomainDecomposition_Plex()
311 DM odm; in DMCreateNeumannOverlap_Plex() local
327 PetscCall(DMPlexDistributeOverlap(dm, 1, &sf, &odm)); in DMCreateNeumannOverlap_Plex()
328 if (!odm) { in DMCreateNeumannOverlap_Plex()
338 PetscCall(DMSetLocalSection(odm, osec)); in DMCreateNeumannOverlap_Plex()
339 PetscCall(DMCopyDisc(dm, odm)); in DMCreateNeumannOverlap_Plex()
341 PetscCall(DMPlexSetMaxProjectionHeight(odm, mh)); in DMCreateNeumannOverlap_Plex()
346 PetscCall(DMTransferMaterialParameters(dm, sf, odm)); in DMCreateNeumannOverlap_Plex()
350 PetscCall(PetscObjectSetName((PetscObject)odm, "OVL")); in DMCreateNeumannOverlap_Plex()
351 PetscCall(DMViewFromOptions(odm, NULL, "-dm_plex_view_neumann_overlap")); in DMCreateNeumannOverlap_Plex()
357 PetscCall(DMSetMatType(odm, MATIS)); in DMCreateNeumannOverlap_Plex()
358 PetscCall(DMCreateMatrix(odm, &pJ)); in DMCreateNeumannOverlap_Plex()
366 PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)odm), n, idxs, PETSC_COPY_VALUES, ovl)); in DMCreateNeumannOverlap_Plex()
369 PetscCall(DMDestroy(&odm)); in DMCreateNeumannOverlap_Plex()