Lines Matching refs:fhat
396 static PetscErrorCode TestL2ProjectionParticlesToField(DM dm, DM sw, Vec fhat, AppCtx *user) in TestL2ProjectionParticlesToField() argument
399 Vec fields[1] = {fhat}; in TestL2ProjectionParticlesToField()
408 PetscCall(computeFEMMoments(dm, fhat, fmoments, user)); in TestL2ProjectionParticlesToField()
417 static PetscErrorCode TestL2ProjectionFieldToParticles(DM dm, DM sw, Vec fhat, AppCtx *user) in TestL2ProjectionFieldToParticles() argument
420 Vec fields[1] = {fhat}; in TestL2ProjectionFieldToParticles()
429 PetscCall(computeFEMMoments(dm, fhat, fmoments, user)); in TestL2ProjectionFieldToParticles()
554 Vec f, rhs, fhat, grad; /* Particle field f, \int phi_i f, FEM field */ in TestFieldGradientProjection() local
563 PetscCall(DMGetGlobalVector(dm, &fhat)); in TestFieldGradientProjection()
580 PetscCall(DMPlexSNESComputeJacobianFEM(dm, fhat, M, M, user)); in TestFieldGradientProjection()
582 PetscCall(InterpolateGradient(dm, fhat, grad)); in TestFieldGradientProjection()
588 PetscCall(PetscObjectSetName((PetscObject)fhat, "fhat")); in TestFieldGradientProjection()
589 PetscCall(VecViewFromOptions(fhat, NULL, "-fhat_view")); in TestFieldGradientProjection()
602 PetscCall(DMRestoreGlobalVector(dm, &fhat)); in TestFieldGradientProjection()
661 Vec fhat; in main() local
664 PetscCall(DMGetGlobalVector(dm, &fhat)); in main()
665 PetscCall(TestL2ProjectionParticlesToField(dm, sw, fhat, &user)); in main()
666 PetscCall(TestL2ProjectionFieldToParticles(dm, sw, fhat, &user)); in main()
668 PetscCall(DMRestoreGlobalVector(dm, &fhat)); in main()