dm/interface/dm.c

78   v->defaultConstraint.section = NULL;
742   PetscCall(PetscSectionDestroy(&(*dm)->defaultConstraint.section));
906   PetscCall(PetscOptionsBool("-dm_ignore_perm_output", "Ignore the local section permutation on output", "DMGetOutputDM", dm->ignorePermOutput, &dm->ignorePermOutput, NULL));
1115     PetscSection section, sectionGlobal;
1117     PetscCall(DMGetLocalSection(dm, &section));
1118     if (section) {
1124       PetscCall(PetscSectionGetChart(section, &pStart, &pEnd));
1125       PetscCall(PetscSectionGetStorageSize(section, &n));
1131         PetscCall(PetscSectionGetDof(section, p, &dof));
1132         PetscCall(PetscSectionGetConstraintDof(section, p, &cdof));
1133         PetscCall(PetscSectionGetConstraintIndices(section, p, &cdofs));
1950   PetscSection section, sectionGlobal;
1966   PetscCall(DMGetLocalSection(dm, &section));
1967   if (section) {
1972     PetscCall(PetscSectionGetNumFields(section, &nF));
1977       PetscCall(PetscSectionGetFieldComponents(section, f, &fieldNc[f]));
1987           PetscCall(PetscSectionGetFieldDof(section, p, f, &fdof));
1988           PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof));
2011           PetscCall(PetscSectionGetFieldDof(section, p, f, &fdof));
2012           PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof));
2023         PetscCall(PetscSectionGetFieldName(section, f, &fieldName));
2119     PetscSection section;
2122     PetscCall(DMGetLocalSection(dm, &section));
2123     if (section) PetscCall(PetscSectionGetNumFields(section, &numFields));
2124     if (section && numFields && dm->ops->createsubdm) {
2134           PetscCall(PetscSectionGetFieldName(section, f, &fieldName));
2798   PetscSection section, cSec;
2811     PetscCall(DMGetLocalSection(dm, &section));
2821         PetscCall(VecSetValuesSection(l, section, p, vals, INSERT_ALL_VALUES));
3017     PetscSection section, cSec;
3025     PetscCall(DMGetLocalSection(dm, &section));
3033         PetscCall(VecGetValuesSection(l, section, p, &vals));
4339 . section - The `PetscSection`
4342 . -dm_petscsection_view - View the section created by the `DM`
4351 PetscErrorCode DMGetLocalSection(DM dm, PetscSection *section)
4355   PetscAssertPointer(section, 2);
4380   *section = dm->localSection;
4389 - section - The `PetscSection`
4398 PetscErrorCode DMSetLocalSection(DM dm, PetscSection section)
4405   if (section) PetscValidHeaderSpecific(section, PETSC_SECTION_CLASSID, 2);
4406   PetscCall(PetscObjectReference((PetscObject)section));
4408   dm->localSection = section;
4409   if (section) PetscCall(PetscSectionGetNumFields(dm->localSection, &numFields));
4421   /* The global section and the SectionSF will be rebuilt
4467 + section - The `PetscSection` describing the range of the constraint matrix: relates rows of the constraint matrix to dofs of the default section.  Returns `NULL` if there are no local constraints.
4468 . mat     - The `Mat` that interpolates local constraints: its width should be the layout size of the default section.  Returns `NULL` if there are no local constraints.
4478 PetscErrorCode DMGetDefaultConstraints(DM dm, PetscSection *section, Mat *mat, Vec *bias)
4482   if (!dm->defaultConstraint.section && !dm->defaultConstraint.mat && dm->ops->createdefaultconstraints) PetscUseTypeMethod(dm, createdefaultconstraints);
4483   if (section) *section = dm->defaultConstraint.section;
4496 . section - The `PetscSection` describing the range of the constraint matrix: relates rows of the constraint matrix to dofs of the default section.  Must have a local communicator (`PETSC_COMM_SELF` or derivative).
4497 . mat     - The `Mat` that interpolates local constraints: its width should be the layout size of the default section:  `NULL` indicates no constraints.  Must have a local communicator (`PETSC_COMM_SELF` or derivative).
4511 PetscErrorCode DMSetDefaultConstraints(DM dm, PetscSection section, Mat mat, Vec bias)
4517   if (section) {
4518     PetscValidHeaderSpecific(section, PETSC_SECTION_CLASSID, 2);
4519     PetscCallMPI(MPI_Comm_compare(PETSC_COMM_SELF, PetscObjectComm((PetscObject)section), &result));
4520     PetscCheck(result == MPI_CONGRUENT || result == MPI_IDENT, PETSC_COMM_SELF, PETSC_ERR_ARG_NOTSAMECOMM, "constraint section must have local communicator");
4532   PetscCall(PetscObjectReference((PetscObject)section));
4533   PetscCall(PetscSectionDestroy(&dm->defaultConstraint.section));
4534   dm->defaultConstraint.section = section;
4644 . section - The `PetscSection`
4653 PetscErrorCode DMGetGlobalSection(DM dm, PetscSection *section)
4657   PetscAssertPointer(section, 2);
4671   *section = dm->globalSection;
4680 - section - The PetscSection, or `NULL`
4689 PetscErrorCode DMSetGlobalSection(DM dm, PetscSection section)
4693   if (section) PetscValidHeaderSpecific(section, PETSC_SECTION_CLASSID, 2);
4694   PetscCall(PetscObjectReference((PetscObject)section));
4696   dm->globalSection = section;
4698   if (section) PetscCall(DMDefaultSectionCheckConsistency_Internal(dm, dm->localSection, section));
4735     PetscSection section, gSection;
4737     PetscCall(DMGetLocalSection(dm, &section));
4738     if (section) {
4740       PetscCall(DMCreateSectionSF(dm, section, gSection));
4806   DMGetPointSF - Get the `PetscSF` encoding the parallel section point overlap for the `DM`.
4833   DMSetPointSF - Set the `PetscSF` encoding the parallel section point overlap for the `DM`.
6189   DMUseTensorOrder - Use a tensor product closure ordering for the default section
6528   PetscSection section;
6536   PetscCall(DMGetLocalSection(dm, &section));
6537   PetscCall(PetscSectionHasConstraints(section, &hasConstraints));
6538   PetscCall(PetscSectionGetPermutation(section, &perm));
6553     PetscCall(PetscSectionClone(section, &newSection));
6559     PetscCall(PetscSectionCreateGlobalSection(section, sf, usePerm, PETSC_TRUE, PETSC_FALSE, &gsection));
7228 // PetscClangLinter pragma ignore: -fdoc-section-header-unknown
7349 // PetscClangLinter pragma ignore: -fdoc-section-header-unknown
8091   PetscCheck(!dm->localSection, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_WRONGSTATE, "Cannot add boundary to DM after creating local section");
8877   of the same topology. This implies that the section (field data) on one
9658   DMReorderSectionSetDefault - Set flag indicating whether the local section should be reordered by default
9679   DMReorderSectionGetDefault - Get flag indicating whether the local section should be reordered by default
9704   DMReorderSectionSetType - Set the type of local section reordering
9725   DMReorderSectionGetType - Get the reordering type for the local section