| /petsc/include/petsc/private/ |
| H A D | petscfeimpl.h | 296 …atic(PetscFE fe, const PetscScalar x[], PetscFEGeom *fegeom, PetscInt q, PetscScalar interpolant[]) in PetscFEInterpolate_Static() argument
307 interpolant[fc] = 0.0; in PetscFEInterpolate_Static()
308 for (PetscInt f = 0; f < Nb; ++f) interpolant[fc] += x[f] * basis[(q * Nb + f) * Nc + fc]; in PetscFEInterpolate_Static()
311 PetscCall(PetscFEPushforward(fe, fegeom, 1, interpolant)); in PetscFEInterpolate_Static()
315 …tscTabulation T, const PetscScalar x[], PetscFEGeom *fegeom, PetscInt q, PetscScalar interpolant[]) in PetscFEInterpolateAtPoints_Static() argument
325 interpolant[fc] = 0.0; in PetscFEInterpolateAtPoints_Static()
326 for (f = 0; f < Nb; ++f) interpolant[fc] += x[f] * basis[(q * Nb + f) * Nc + fc]; in PetscFEInterpolateAtPoints_Static()
329 PetscCall(PetscFEPushforward(fe, fegeom, 1, interpolant)); in PetscFEInterpolateAtPoints_Static()
333 … fe, PetscInt k, const PetscScalar x[], PetscFEGeom *fegeom, PetscInt q, PetscScalar interpolant[]) in PetscFEInterpolateGradient_Static() argument
349 for (d = 0; d < cdim; ++d) interpolant[fc * cdim + d] = 0.0; in PetscFEInterpolateGradient_Static()
[all …]
|
| H A D | petscfvimpl.h | 98 …PetscFVInterpolate_Static(PetscFV fv, const PetscScalar x[], PetscInt q, PetscScalar interpolant[]) in PetscFVInterpolate_Static() argument
104 PetscCall(PetscArraycpy(interpolant, x, Nc)); in PetscFVInterpolate_Static()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 1465 PetscScalar *funcVal, *interpolant; in DMPlexComputeL2DiffLocal() local
1504 …PetscCall(PetscMalloc6(numComponents, &funcVal, numComponents, &interpolant, coordDim * (Nq + 1), … in DMPlexComputeL2DiffLocal()
1555 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1558 …ASSID) PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[fieldOffset], &qgeom, q, interpolant)); in DMPlexComputeL2DiffLocal()
1559 …TSCFV_CLASSID) PetscCall(PetscFVInterpolate_Static((PetscFV)obj, &x[fieldOffset], q, interpolant)); in DMPlexComputeL2DiffLocal()
1565 …(double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)Pe… in DMPlexComputeL2DiffLocal()
1566 elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]; in DMPlexComputeL2DiffLocal()
1576 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1589 PetscScalar *funcVal, *interpolant; in DMComputeL2GradientDiff_Plex() local
1621 ….J, coordDim * coordDim * Nq, &fegeom.invJ, numComponents * coordDim, &interpolant, Nq, &fegeom.de… in DMComputeL2GradientDiff_Plex()
[all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex2.c | 446 PetscScalar *interpolant, *gradsum; in InterpolateGradient() local
482 …oc6(coordDim * numComponents * 2, &gradsum, coordDim * numComponents, &interpolant, coordDim * Nq,… in InterpolateGradient()
516 …Call(PetscFEInterpolateGradient_Static((PetscFE)obj, 1, &x[fieldOffset], &fegeom, q, interpolant)); in InterpolateGradient()
521 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.… in InterpolateGradient()
544 PetscCall(PetscFree6(gradsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in InterpolateGradient()
|
| /petsc/src/snes/utils/dm/ |
| H A D | dminterpolatesnes.c | 884 PetscScalar *interpolant; in DMInterpolationEvaluate() local
889 PetscCall(VecGetArrayWrite(v, &interpolant)); in DMInterpolationEvaluate()
916 interpolant[p * ctx->dof + coff + fc] = 0.0; in DMInterpolationEvaluate()
917 …for (PetscInt f = 0; f < Nb; ++f) interpolant[p * ctx->dof + coff + fc] += xa[foff + f] * basis[(0… in DMInterpolationEvaluate()
929 … for (PetscInt fc = 0; fc < Nc; ++fc) interpolant[p * ctx->dof + coff + fc] = xa[foff + fc]; in DMInterpolationEvaluate()
939 PetscCall(VecRestoreArrayWrite(v, &interpolant)); in DMInterpolationEvaluate()
|
| H A D | dmadapt.c | 899 PetscScalar *x = NULL, *field, *gradient, *interpolant, *interpolantGrad; in DMAdaptorComputeErrorIndicator_Gradient() local
913 PetscCall(PetscMalloc2(Nc, &interpolant, cdim * Nc, &interpolantGrad)); in DMAdaptorComputeErrorIndicator_Gradient()
922 PetscCall(PetscArrayzero(interpolant, Nc)); in DMAdaptorComputeErrorIndicator_Gradient()
925 …scFEInterpolateFieldAndGradient_Static((PetscFE)obj, 1, x, &fegeom, q, interpolant, interpolantGra… in DMAdaptorComputeErrorIndicator_Gradient()
929 field[fc] += interpolant[fc] * wt * fegeom.detJ[q]; in DMAdaptorComputeErrorIndicator_Gradient()
935 PetscCall(PetscFree2(interpolant, interpolantGrad)); in DMAdaptorComputeErrorIndicator_Gradient()
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 1399 …ons from the coarse space in the fine space, keeping the energy of the interpolant about the same.…
1413 Now we would like the interpolant of the coarse representer to the fine grid to be as close as poss…
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | DMPlex.pyx | 3227 The output `Vec` which holds the Clement interpolant of the gradient.
|