| /petsc/src/ts/impls/mimex/ |
| H A D | mimex.c | 102 PetscInt fdof, fcdof, d; in SNESTSFormFunction_Mimex() local
105 PetscCall(PetscSectionGetFieldConstraintDof(s, p, f, &fcdof)); in SNESTSFormFunction_Mimex()
108 for (d = 0; d < fdof - fcdof; ++d) axs[d] = a[d]; in SNESTSFormFunction_Mimex()
174 PetscInt fdof, fcdof, d; in TSStep_Mimex_Split() local
177 PetscCall(PetscSectionGetFieldConstraintDof(s, p, f, &fcdof)); in TSStep_Mimex_Split()
180 for (d = 0; d < fdof - fcdof; ++d) as[d] = au[d]; in TSStep_Mimex_Split()
196 PetscInt fdof, fcdof, d; in TSStep_Mimex_Split() local
199 PetscCall(PetscSectionGetFieldConstraintDof(s, p, f, &fcdof)); in TSStep_Mimex_Split()
202 for (d = 0; d < fdof - fcdof; ++d) as[d] += dt * au[d]; in TSStep_Mimex_Split()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexsection.c | 304 PetscInt fdof, fcdof, c, j; in DMPlexCreateSectionBCIndicesField() local
310 fcdof = fdof; in DMPlexCreateSectionBCIndicesField()
314 PetscCall(PetscSectionGetFieldConstraintDof(section, p, field, &fcdof)); in DMPlexCreateSectionBCIndicesField()
317 for (d = 0; d < fcdof; ++d) { in DMPlexCreateSectionBCIndicesField()
324 for (c = d; c < fcdof; ++c) indices[c] = -1; in DMPlexCreateSectionBCIndicesField()
325 fcdof = d; in DMPlexCreateSectionBCIndicesField()
327 PetscCall(PetscSectionSetFieldConstraintDof(section, p, field, fcdof)); in DMPlexCreateSectionBCIndicesField()
377 PetscInt fcdof, fdof; in DMPlexCreateSectionBCIndices() local
380 PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof)); in DMPlexCreateSectionBCIndices()
383 for (d = 0; d < fcdof; ++d) indices[numConst + d] = find[d] + foff; in DMPlexCreateSectionBCIndices()
[all …]
|
| H A D | plex.c | 7029 PetscInt fdof, foff, fcdof, foffset = *offset; in updatePointFields_private() local
7035 PetscCall(PetscSectionGetFieldConstraintDof(section, point, f, &fcdof)); in updatePointFields_private()
7038 if (!fcdof || setBC) { in updatePointFields_private()
7057 if ((cind < fcdof) && (b == fcdofs[cind])) { in updatePointFields_private()
7065 if ((cind < fcdof) && (b == fcdofs[cind])) { in updatePointFields_private()
7075 if ((cind < fcdof) && (b == fcdofs[cind])) { in updatePointFields_private()
7083 if ((cind < fcdof) && (b == fcdofs[cind])) { in updatePointFields_private()
7099 PetscInt fdof, foff, fcdof, foffset = *offset; in updatePointFieldsBC_private() local
7107 PetscCall(PetscSectionGetFieldConstraintDof(section, point, f, &fcdof)); in updatePointFieldsBC_private()
7110 if (fcdof) { in updatePointFieldsBC_private()
[all …]
|
| /petsc/src/dm/interface/ |
| H A D | dmi.c | 118 PetscInt fdof, fcdof, sfdof, sfcdof = 0; in PetscSectionSelectFields_Private() local
122 PetscCall(PetscSectionGetFieldConstraintDof(s, point, fields[f], &fcdof)); in PetscSectionSelectFields_Private()
141 pSubSize += fdof - fcdof; in PetscSectionSelectFields_Private()
176 PetscInt fdof, fcdof, poff = 0; in PetscSectionSelectFields_Private() local
181 PetscCall(PetscSectionGetFieldConstraintDof(s, point, f2, &fcdof)); in PetscSectionSelectFields_Private()
182 poff += fdof - fcdof; in PetscSectionSelectFields_Private()
185 PetscCall(PetscSectionGetFieldConstraintDof(s, point, fields[f], &fcdof)); in PetscSectionSelectFields_Private()
209 … for (PetscInt fc = 0; fc < fdof - fcdof; ++fc, ++subOff) subIndices[subOff] = goff + poff + fc; in PetscSectionSelectFields_Private()
|
| H A D | dm.c | 1985 PetscInt fdof, fcdof, fpdof; in DMCreateFieldIS() local
1988 PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof)); in DMCreateFieldIS()
1989 fpdof = fdof - fcdof; in DMCreateFieldIS()
2009 PetscInt fdof, fcdof, fc; in DMCreateFieldIS() local
2012 PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof)); in DMCreateFieldIS()
2013 … for (fc = 0; fc < fdof - fcdof; ++fc, ++fieldSizes[f]) fieldIndices[f][fieldSizes[f]] = goff++; in DMCreateFieldIS()
|
| /petsc/include/petsc/private/ |
| H A D | dmimpl.h | 501 PetscInt loff, lfoff, fdof, fcdof, ffcdof, f; in DMGetGlobalFieldOffset_Private() local
509 PetscCall(PetscSectionGetFieldConstraintDof(dm->localSection, point, field, &fcdof)); in DMGetGlobalFieldOffset_Private()
515 *end = *start < 0 ? *start - (fdof - fcdof) : *start + fdof - fcdof; in DMGetGlobalFieldOffset_Private()
526 const PetscInt fcdof = fs->bc ? fs->bc->atlasDof[point - fs->bc->pStart] : 0; in DMGetGlobalFieldOffset_Private() local
534 *end = *start < 0 ? *start - (fdof - fcdof) : *start + fdof - fcdof; in DMGetGlobalFieldOffset_Private()
|
| /petsc/src/vec/is/section/interface/ |
| H A D | section.c | 128 PetscInt dof, cdof, fcdof = 0; in PetscSectionCopy_Internal() local
140 if (force_constrained) fcdof = dof; in PetscSectionCopy_Internal()
141 else PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof)); in PetscSectionCopy_Internal()
142 if (fcdof) PetscCall(PetscSectionSetFieldConstraintDof(newSection, p, f, fcdof)); in PetscSectionCopy_Internal()
148 PetscInt off, cdof, fcdof = 0; in PetscSectionCopy_Internal() local
163 PetscCall(PetscSectionGetFieldConstraintDof(newSection, p, f, &fcdof)); in PetscSectionCopy_Internal()
164 if (fcdof) { in PetscSectionCopy_Internal()
|