| /petsc/src/dm/dt/fv/interface/ceed/ |
| H A D | fvceed.c | 18 PetscErrorCode PetscFVSetCeed(PetscFV fv, Ceed ceed) in PetscFVSetCeed() argument 21 PetscValidHeaderSpecific(fv, PETSCFV_CLASSID, 1); in PetscFVSetCeed() 22 if (fv->ceed == ceed) PetscFunctionReturn(PETSC_SUCCESS); in PetscFVSetCeed() 23 PetscCallCEED(CeedReferenceCopy(ceed, &fv->ceed)); in PetscFVSetCeed() 45 PetscErrorCode PetscFVGetCeedBasis(PetscFV fv, CeedBasis *basis) in PetscFVGetCeedBasis() argument 51 PetscValidHeaderSpecific(fv, PETSCFV_CLASSID, 1); in PetscFVGetCeedBasis() 53 if (!fv->ceedBasis && fv->ceed) { in PetscFVGetCeedBasis() 54 PetscCall(PetscFVGetSpatialDimension(fv, &dim)); in PetscFVGetCeedBasis() 55 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in PetscFVGetCeedBasis() 56 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscFVGetCeedBasis() [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex11.c | 11 PetscBool fv; // Use an FV discretization, instead of FE member 64 options->fv = PETSC_FALSE; in ProcessOptions() 73 …PetscCall(PetscOptionsBool("-fv", "Use FV instead of FE", __FILE__, options->fv, &options->fv, NUL… in ProcessOptions() 113 PetscFV fv; in CreateDiscretization() local 121 if (user->fv) { in CreateDiscretization() 122 PetscCall(PetscFVCreate(PETSC_COMM_SELF, &fv)); in CreateDiscretization() 123 PetscCall(PetscObjectSetName((PetscObject)fv, "fv")); in CreateDiscretization() 124 PetscCall(PetscFVSetNumComponents(fv, 1)); in CreateDiscretization() 125 PetscCall(PetscFVSetSpatialDimension(fv, dim)); in CreateDiscretization() 126 PetscCall(PetscFVCreateDualSpace(fv, ct)); in CreateDiscretization() [all …]
|
| H A D | ex2.c | 473 PetscFV fv = (PetscFV)obj; in InterpolateGradient() local 475 PetscCall(PetscFVGetQuadrature(fv, &quad)); in InterpolateGradient() 476 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in InterpolateGradient()
|
| /petsc/src/snes/interface/ |
| H A D | snesob.c | 145 PetscScalar fv, xv; in SNESObjectiveComputeFunctionDefaultFD() local 196 fv = (-ob2 + 6. * ob1 - 3. * ob - 2. * ob3) / (6. * dx); in SNESObjectiveComputeFunctionDefaultFD() 197 if (PetscAbsScalar(fv) > eps) { in SNESObjectiveComputeFunctionDefaultFD() 198 PetscCall(VecSetValues(F, 1, &i, &fv, INSERT_VALUES)); in SNESObjectiveComputeFunctionDefaultFD() 200 fv = 0.; in SNESObjectiveComputeFunctionDefaultFD() 201 PetscCall(VecSetValues(F, 1, &i, &fv, INSERT_VALUES)); in SNESObjectiveComputeFunctionDefaultFD()
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex49.c | 57 PetscFV fv; in SetupDiscretization() local 59 PetscCall(PetscFVCreate(PETSC_COMM_SELF, &fv)); in SetupDiscretization() 60 PetscCall(PetscFVSetType(fv, PETSCFVLEASTSQUARES)); in SetupDiscretization() 61 PetscCall(PetscFVSetNumComponents(fv, dim)); in SetupDiscretization() 62 PetscCall(PetscFVSetSpatialDimension(fv, dim)); in SetupDiscretization() 63 PetscCall(PetscFVSetFromOptions(fv)); in SetupDiscretization() 64 PetscCall(PetscFVSetUp(fv)); in SetupDiscretization() 65 PetscCall(PetscObjectSetName((PetscObject)fv, "vector")); in SetupDiscretization() 66 PetscCall(DMSetField(dm, 0, NULL, (PetscObject)fv)); in SetupDiscretization() 67 PetscCall(PetscFVDestroy(&fv)); in SetupDiscretization()
|
| H A D | ex3.c | 599 PetscFV fv; in TestFVGrad() local 618 PetscCall(PetscFVCreate(comm, &fv)); in TestFVGrad() 619 PetscCall(PetscFVSetType(fv, PETSCFVLEASTSQUARES)); in TestFVGrad() 620 PetscCall(PetscFVSetNumComponents(fv, user->numComponents)); in TestFVGrad() 621 PetscCall(PetscFVSetSpatialDimension(fv, dim)); in TestFVGrad() 622 PetscCall(PetscFVSetFromOptions(fv)); in TestFVGrad() 623 PetscCall(PetscFVSetUp(fv)); in TestFVGrad() 640 PetscCall(DMSetField(dmfv, 0, NULL, (PetscObject)fv)); in TestFVGrad() 644 PetscCall(DMPlexGetGradientDM(dmfv, fv, &dmgrad)); in TestFVGrad() 721 PetscCall(PetscFVDestroy(&fv)); in TestFVGrad()
|
| /petsc/src/ts/utils/libceed/ |
| H A D | dmplextsceed.c | 9 PetscFV fv; in DMPlexTSComputeRHSFunctionFVMCEED() local 19 PetscCall(DMGetField(dm, 0, NULL, (PetscObject *)&fv)); in DMPlexTSComputeRHSFunctionFVMCEED() 20 PetscCall(DMPlexInsertBoundaryValuesFVM(dm, fv, locX, time, NULL)); in DMPlexTSComputeRHSFunctionFVMCEED()
|
| /petsc/src/dm/impls/forest/tests/ |
| H A D | ex2.c | 190 PetscFV fv; in main() local 198 PetscCall(PetscFVCreate(PETSC_COMM_SELF, &fv)); in main() 199 PetscCall(PetscFVSetSpatialDimension(fv, dim)); in main() 200 PetscCall(PetscFVSetType(fv, PETSCFVLEASTSQUARES)); in main() 201 PetscCall(PetscFVSetNumComponents(fv, Nf)); in main() 204 PetscCall(PetscFVSetLimiter(fv, limiter)); in main() 206 PetscCall(PetscFVSetFromOptions(fv)); in main() 207 PetscCall(DMSetField(base, 0, NULL, (PetscObject)fv)); in main() 208 PetscCall(PetscFVDestroy(&fv)); in main()
|
| /petsc/src/dm/dt/fv/interface/ |
| H A D | fv.c | 1755 PetscErrorCode PetscFVClone(PetscFV fv, PetscFV *fvNew) in PetscFVClone() argument 1763 PetscCall(PetscFVGetDualSpace(fv, &Q)); in PetscFVClone() 1764 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscFVClone() 1767 PetscCall(PetscFVCreate(PetscObjectComm((PetscObject)fv), fvNew)); in PetscFVClone() 1769 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in PetscFVClone() 1771 PetscCall(PetscFVGetSpatialDimension(fv, &cdim)); in PetscFVClone() 1796 PetscErrorCode PetscFVRefine(PetscFV fv, PetscFV *fvRef) in PetscFVRefine() argument 1808 PetscCall(PetscFVGetDualSpace(fv, &Q)); in PetscFVRefine() 1809 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscFVRefine() 1813 PetscCall(DMRefine(K, PetscObjectComm((PetscObject)fv), &Kref)); in PetscFVRefine() [all …]
|
| /petsc/src/snes/tests/ |
| H A D | ex2.c | 229 PetscFV fv; in CreateDiscretization() local 232 PetscCall(PetscFVCreate(PetscObjectComm((PetscObject)dm), &fv)); in CreateDiscretization() 233 PetscCall(PetscObjectSetName((PetscObject)fv, "phi")); in CreateDiscretization() 234 PetscCall(PetscFVSetFromOptions(fv)); in CreateDiscretization() 235 PetscCall(PetscFVSetNumComponents(fv, Nc)); in CreateDiscretization() 237 PetscCall(PetscFVSetSpatialDimension(fv, cdim)); in CreateDiscretization() 238 PetscCall(DMSetField(dm, 0, NULL, (PetscObject)fv)); in CreateDiscretization() 239 PetscCall(PetscFVDestroy(&fv)); in CreateDiscretization()
|
| /petsc/include/petsc/private/ |
| H A D | petscfvimpl.h | 98 static inline PetscErrorCode PetscFVInterpolate_Static(PetscFV fv, const PetscScalar x[], PetscInt … in PetscFVInterpolate_Static() argument 103 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in PetscFVInterpolate_Static()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexvtu.c | 146 static PetscErrorCode DMGetFieldIfFV_Private(DM dm, PetscInt field, PetscFV *fv) in DMGetFieldIfFV_Private() argument 153 *fv = NULL; in DMGetFieldIfFV_Private() 158 if (fClass == PETSCFV_CLASSID) *fv = (PetscFV)f; in DMGetFieldIfFV_Private() 291 PetscFV fv = NULL; in DMPlexVTKWriteAll_VTU() local 300 PetscCall(DMGetFieldIfFV_Private(dmX, field, &fv)); in DMPlexVTKWriteAll_VTU() 311 if (fv) { in DMPlexVTKWriteAll_VTU() 312 PetscCall(PetscFVGetComponentName(fv, j, &compName)); in DMPlexVTKWriteAll_VTU() 333 if (fv) PetscCall(PetscFVGetComponentName(fv, j, &compName)); in DMPlexVTKWriteAll_VTU() 559 PetscFV fv = NULL; in DMPlexVTKWriteAll_VTU() local 565 PetscCall(DMGetFieldIfFV_Private(dmX, field, &fv)); in DMPlexVTKWriteAll_VTU() [all …]
|
| H A D | plexfem.c | 1064 PetscFV fv; in DMPlexInsertBoundaryValuesRiemann() local 1066 PetscCall(PetscDSGetDiscretization(prob, field, (PetscObject *)&fv)); in DMPlexInsertBoundaryValuesRiemann() 1069 PetscCall(PetscFVGetNumComponents(fv, &pdim)); in DMPlexInsertBoundaryValuesRiemann() 1389 PetscErrorCode DMPlexInsertBoundaryValuesFVM(DM dm, PetscFV fv, Vec locX, PetscReal time, Vec *locG… in DMPlexInsertBoundaryValuesFVM() argument 1396 PetscValidHeaderSpecific(fv, PETSCFV_CLASSID, 2); in DMPlexInsertBoundaryValuesFVM() 1404 PetscCall(DMPlexGetGradientDM(dm, fv, &dmGrad)); in DMPlexInsertBoundaryValuesFVM() 1411 …PetscCall(DMPlexReconstructGradients_Internal(dm, fv, fStart, fEnd, faceGeometryFVM, cellGeometryF… in DMPlexInsertBoundaryValuesFVM() 1495 PetscFV fv = (PetscFV)obj; in DMPlexComputeL2DiffLocal() local 1497 PetscCall(PetscFVGetQuadrature(fv, &quad)); in DMPlexComputeL2DiffLocal() 1498 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in DMPlexComputeL2DiffLocal() [all …]
|
| H A D | plexceed.c | 198 PetscFV fv; in DMPlexGetLocalOffsetsSupport() local 226 PetscCall(PetscDSGetDiscretization(ds, ds_field, (PetscObject *)&fv)); in DMPlexGetLocalOffsetsSupport() 227 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in DMPlexGetLocalOffsetsSupport()
|
| H A D | plexgeometry.c | 2737 PetscInt fv[4] = {0, 1, 2, 3}; in DMPlexComputeGeometryFVM_2D_Internal() local 2745 fv[2] = 3; in DMPlexComputeGeometryFVM_2D_Internal() 2746 fv[3] = 2; in DMPlexComputeGeometryFVM_2D_Internal() 2761 e0[d] = PetscRealPart(coords[cdim * fv[p + 1] + d]) - origin[d]; in DMPlexComputeGeometryFVM_2D_Internal() 2762 e1[d] = PetscRealPart(coords[cdim * fv[p + 2] + d]) - origin[d]; in DMPlexComputeGeometryFVM_2D_Internal() 2772 …m; d++) c[d] += a * PetscRealPart(origin[d] + coords[cdim * fv[p + 1] + d] + coords[cdim * fv[p + … in DMPlexComputeGeometryFVM_2D_Internal() 2856 PetscInt fv[4] = {0, 1, 2, 3}; in DMPlexComputeGeometryFVM_3D_Internal() local 2860 fv[2] = 3; in DMPlexComputeGeometryFVM_3D_Internal() 2861 fv[3] = 2; in DMPlexComputeGeometryFVM_3D_Internal() 2866 coordsTmp[0 * dim + d] = PetscRealPart(coords[faces[fOff + fv[0]] * dim + d]) - origin[d]; in DMPlexComputeGeometryFVM_3D_Internal() [all …]
|
| H A D | plexproject.c | 867 PetscFV fv; in DMProjectLocal_Generic_Plex() local 871 PetscCall(PetscDSGetDiscretization(ds, f, (PetscObject *)&fv)); in DMProjectLocal_Generic_Plex() 872 PetscCall(PetscFVGetDualSpace(fv, &cellsp[f])); in DMProjectLocal_Generic_Plex() 888 PetscFV fv; in DMProjectLocal_Generic_Plex() local 890 PetscCall(PetscDSGetDiscretization(dsIn, f, (PetscObject *)&fv)); in DMProjectLocal_Generic_Plex() 891 PetscCall(PetscFVGetDualSpace(fv, &cellspIn[f])); in DMProjectLocal_Generic_Plex()
|
| H A D | plexsection.c | 600 PetscFV fv = (PetscFV)dm->fields[f].disc; in DMCreateLocalSection_Plex() local 602 PetscCall(PetscFVGetNumComponents(fv, &numComp[f])); in DMCreateLocalSection_Plex()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex18.c | 751 PetscFV fv; in SetupDiscretization() local 765 PetscCall(PetscFVCreate(PetscObjectComm((PetscObject)dm), &fv)); in SetupDiscretization() 766 PetscCall(PetscObjectSetName((PetscObject)fv, "porosity")); in SetupDiscretization() 767 PetscCall(PetscFVSetFromOptions(fv)); in SetupDiscretization() 768 PetscCall(PetscFVSetNumComponents(fv, 1)); in SetupDiscretization() 769 PetscCall(PetscFVSetSpatialDimension(fv, dim)); in SetupDiscretization() 771 PetscCall(PetscFVSetQuadrature(fv, q)); in SetupDiscretization() 774 if (user->useFV) PetscCall(DMSetField(dm, 1, NULL, (PetscObject)fv)); in SetupDiscretization() 788 PetscCall(PetscFVDestroy(&fv)); in SetupDiscretization() 884 Vec *fv; in MonitorFunctionals() local [all …]
|
| /petsc/src/tao/unconstrained/impls/bmrm/ |
| H A D | bmrm.c | 24 PetscReal fr, fbest, fv, fc, fv0; in solve() local 138 fv = 0.0; in solve() 142 fv += xplus[i] * (0.5 * tplus[i] + f[i]); in solve() 146 if ((innerIter == 1 && fv >= fv0) || (innerIter > 1 && fv >= fr)) { in solve() 147 fv = 0.0; in solve() 151 fv += xplus[i] * (0.5 * tplus[i] + f[i]); in solve() 164 if (fv < fbest) { in solve() 165 fbest = fv; in solve() 166 fc = fv; in solve() 169 fc = (fc > fv ? fc : fv); in solve() [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex55k.kokkos.cxx | 45 PetscScalarKokkosOffsetView2D fv; in FormFunctionLocalVec() local 58 PetscCall(DMDAVecGetKokkosOffsetViewWrite(info->da, f, &fv)); in FormFunctionLocalVec() 69 fv(j, i) = 2.0 * (hydhx + hxdhy) * (xv(j, i) - mms_solution); in FormFunctionLocalVec() 105 fv(j, i) = uxx + uyy - hx * hy * (lambda * PetscExpScalar(u) + mms_forcing); in FormFunctionLocalVec() 110 PetscCall(DMDAVecRestoreKokkosOffsetViewWrite(info->da, f, &fv)); in FormFunctionLocalVec()
|
| H A D | ex19.c | 381 PetscScalar fu, fv, fomega, ftemp; in NonlinearGS() local 503 fv = 0.0; in NonlinearGS() 531 fv = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy - bjiv; in NonlinearGS() 577 yv = fv / dfvdv; in NonlinearGS() 619 pfnorm = PetscRealPart(fu * fu + fv * fv + fomega * fomega + ftemp * ftemp); in NonlinearGS()
|
| /petsc/src/tao/bound/utils/ |
| H A D | isutil.c | 29 PetscReal *fv, *rv; in TaoVecGetSubVec() local 77 PetscCall(VecGetArray(vfull, &fv)); in TaoVecGetSubVec() 81 for (i = 0; i < nlocal; ++i) rv[s[i] - flow] = fv[s[i] - flow]; in TaoVecGetSubVec() 83 PetscCall(VecRestoreArray(vfull, &fv)); in TaoVecGetSubVec()
|
| /petsc/src/dm/interface/ |
| H A D | dmceed.c | 328 PetscFV fv; in DMCeedCreateFVM_Internal() local 338 PetscCall(PetscDSGetDiscretization(ds, 0, (PetscObject *)&fv)); in DMCeedCreateFVM_Internal() 339 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in DMCeedCreateFVM_Internal()
|
| H A D | dmi.c | 255 PetscFV fv, fvNew; in DMSelectFields_Private() local 259 PetscCall(DMGetField(dm, field, NULL, (PetscObject *)&fv)); in DMSelectFields_Private() 260 PetscCall(PetscFVClone(fv, &fvNew)); in DMSelectFields_Private()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtds.c | 399 PetscFV fv = (PetscFV)obj; in PetscDSSetUp() local 401 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscDSSetUp() 437 PetscFV fv = (PetscFV)obj; in PetscDSSetUp() local 439 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscDSSetUp() 441 PetscCall(PetscFVSetQuadrature(fv, maxQuad[dim])); in PetscDSSetUp() 484 PetscFV fv = (PetscFV)obj; in PetscDSSetUp() local 486 PetscCall(PetscFVGetQuadrature(fv, &q)); in PetscDSSetUp() 487 PetscCall(PetscFVGetNumComponents(fv, &Nc)); in PetscDSSetUp() 489 PetscCall(PetscFVGetCellTabulation(fv, &prob->T[f])); in PetscDSSetUp()
|