| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex6.c | 127 PetscInt closureSize, ki, kj, f, c, foff = 0; in ReadData2D() local
138 …, "%2.0f", (double)PetscRealPart(closure[(kj * (user->k[f] + 1) + ki) * user->Nc[f] + c + foff]))); in ReadData2D()
144 foff += PetscSqr(user->k[f] + 1); in ReadData2D()
160 PetscInt closureSize, ki, kj, kk, f, c, foff = 0; in ReadData3D() local
172 …alPart(closure[((kk * (user->k[f] + 1) + kj) * (user->k[f] + 1) + ki) * user->Nc[f] + c + foff]))); in ReadData3D()
180 foff += PetscSqr(user->k[f] + 1); in ReadData3D()
|
| /petsc/src/vec/is/section/interface/ |
| H A D | section.c | 1312 PetscCount foff; in PetscSectionSetUp() local
1320 for (f = 0, foff = s->atlasOff[q]; f < s->numFields; ++f) { in PetscSectionSetUp()
1323 PetscCall(PetscIntCast(foff, &sf->atlasOff[q])); in PetscSectionSetUp()
1324 foff += sf->atlasDof[q]; in PetscSectionSetUp()
1482 PetscInt foff; in PetscSectionCreateGlobalSection() local
1570 for (f = 0, foff = off; f < numFields; ++f) { in PetscSectionCreateGlobalSection()
1575 PetscCall(PetscSectionSetFieldOffset(gs, p, f, foff)); in PetscSectionCreateGlobalSection()
1577 foff = off < 0 ? foff - (dof - cdof) : foff + (dof - cdof); in PetscSectionCreateGlobalSection()
1953 PetscInt off, foff; in PetscSectionGetFieldPointOffset() local
1960 PetscCall(PetscSectionGetOffset(s->field[field], point, &foff)); in PetscSectionGetFieldPointOffset()
[all …]
|
| /petsc/src/snes/utils/dm/ |
| H A D | dminterpolatesnes.c | 893 PetscInt coff = 0, foff = 0, clSize; in DMInterpolationEvaluate() local
917 …for (PetscInt f = 0; f < Nb; ++f) interpolant[p * ctx->dof + coff + fc] += xa[foff + f] * basis[(0… in DMInterpolationEvaluate()
920 foff += Nb; in DMInterpolationEvaluate()
929 … for (PetscInt fc = 0; fc < Nc; ++fc) interpolant[p * ctx->dof + coff + fc] = xa[foff + fc]; in DMInterpolationEvaluate()
931 foff += Nc; in DMInterpolationEvaluate()
936 …PetscCheck(foff == clSize, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Total FE/FV space size %" PetscInt_FM… in DMInterpolationEvaluate()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexinterpolate.c | 540 PetscInt numFaces, foff = 0; in DMPlexInterpolateFaces_Internal() local
550 for (PetscInt cf = 0; cf < numFaces; foff += faceSizes[cf], ++cf) { in DMPlexInterpolateFaces_Internal()
553 const PetscInt *face = &faces[foff]; in DMPlexInterpolateFaces_Internal()
591 PetscInt numFaces, foff = 0; in DMPlexInterpolateFaces_Internal() local
600 for (PetscInt cf = 0; cf < numFaces; foff += faceSizes[cf], ++cf) { in DMPlexInterpolateFaces_Internal()
603 const PetscInt *face = &faces[foff]; in DMPlexInterpolateFaces_Internal()
727 PetscInt coneSize, numFaces, foff = 0, poff = 0; in DMPlexInterpolateFaces_Internal() local
741 for (PetscInt cf = 0; cf < numFaces; foff += faceSizes[cf], ++cf) { in DMPlexInterpolateFaces_Internal()
745 const PetscInt *face = &faces[foff]; in DMPlexInterpolateFaces_Internal()
|
| H A D | plexsection.c | 373 PetscInt numConst = 0, foff = 0; in DMPlexCreateSectionBCIndices() local
383 for (d = 0; d < fcdof; ++d) indices[numConst + d] = find[d] + foff; in DMPlexCreateSectionBCIndices()
385 foff += fdof; in DMPlexCreateSectionBCIndices()
|
| H A D | plexfem.c | 2203 PetscInt foff = 0; in DMPlexComputeClementInterpolant() local
2233 PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[foff], &qgeom, q, interpolant)); in DMPlexComputeClementInterpolant()
2234 for (fc = 0; fc < fNc; ++fc) val[foff + fc] += interpolant[fc] * wt; in DMPlexComputeClementInterpolant()
2237 foff += Nb; in DMPlexComputeClementInterpolant()
2549 PetscInt foff; in DMPlexComputeIntegral_Internal() local
2553 PetscCall(PetscDSGetFieldOffset(prob, f, &foff)); in DMPlexComputeIntegral_Internal()
2560 …obj_func(dim, Nf, NfAux, uOff, uOff_x, &u[totDim * c + foff], NULL, u_x, aOff, NULL, PetscSafePoin… in DMPlexComputeIntegral_Internal()
3906 PetscInt fdofIn, foff, foffIn; in DMPlexGetHybridCellFields() local
3912 PetscCall(PetscDSGetFieldOffsetCohesive(ds, f, &foff)); in DMPlexGetHybridCellFields()
3914 for (PetscInt i = 0; i < fdofIn; ++i) ul[foffIn + i] = xf[foff + i]; in DMPlexGetHybridCellFields()
[all …]
|
| H A D | plex.c | 6538 PetscInt fdof, foff, b; in DMPlexVecGetClosure_Fields_Static() local
6544 PetscCall(PetscSectionGetFieldOffset(section, point, f, &foff)); in DMPlexVecGetClosure_Fields_Static()
6545 varr = &vArray[foff]; in DMPlexVecGetClosure_Fields_Static()
7029 PetscInt fdof, foff, fcdof, foffset = *offset; in updatePointFields_private() local
7036 PetscCall(PetscSectionGetFieldOffset(section, point, f, &foff)); in updatePointFields_private()
7037 a = &array[foff]; in updatePointFields_private()
7099 PetscInt fdof, foff, fcdof, foffset = *offset; in updatePointFieldsBC_private() local
7108 PetscCall(PetscSectionGetFieldOffset(section, point, f, &foff)); in updatePointFieldsBC_private()
7109 a = &array[foff]; in updatePointFieldsBC_private()
7700 PetscInt numFields, foff, f; in DMPlexGetIndicesPointFields_Internal() local
[all …]
|
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 205 PetscInt fdof, foff, d; in SNESMonitorFields() local
208 PetscCall(PetscSectionGetFieldOffset(s, p, f, &foff)); in SNESMonitorFields()
209 for (d = 0; d < fdof; ++d) lnorms[f] += PetscRealPart(PetscSqr(r[foff + d])); in SNESMonitorFields()
|
| /petsc/src/ksp/pc/impls/patch/ |
| H A D | pcpatch.c | 2510 PetscInt fdof, d, foff; in PCSetUp_PATCH() local
2513 PetscCall(PetscSectionGetFieldOffset(s, p, f, &foff)); in PCSetUp_PATCH()
2514 for (d = 0; d < fdof; ++d, ++cdoff) cellDofs[f][cdoff] = foff + d; in PCSetUp_PATCH()
|