16 PetscErrorCode VelocityGradientProjectionCreateDM(NodalProjectionData grad_velo_proj, User user, Pe…  in VelocityGradientProjectionCreateDM()  argument
20   grad_velo_proj->num_comp = 9;  // 9 velocity gradient  in VelocityGradientProjectionCreateDM()
22   PetscCall(DMClone(user->dm, &grad_velo_proj->dm));  in VelocityGradientProjectionCreateDM()
23   PetscCall(PetscObjectSetName((PetscObject)grad_velo_proj->dm, "Velocity Gradient Projection"));  in VelocityGradientProjectionCreateDM()
25 …C_TRUE, PETSC_TRUE, user->app_ctx->degree, 1, user->app_ctx->q_extra, 1, &grad_velo_proj->num_comp,  in VelocityGradientProjectionCreateDM()
26                                grad_velo_proj->dm));  in VelocityGradientProjectionCreateDM()
28   PetscCall(DMGetLocalSection(grad_velo_proj->dm, &section));  in VelocityGradientProjectionCreateDM()
44   NodalProjectionData grad_velo_proj;  in VelocityGradientProjectionSetup()  local
54   PetscCall(PetscNew(&grad_velo_proj));  in VelocityGradientProjectionSetup()
56   PetscCall(VelocityGradientProjectionCreateDM(grad_velo_proj, user, user->app_ctx->degree));  in VelocityGradientProjectionSetup()
59   PetscCall(DMGetDimension(grad_velo_proj->dm, &dim));  in VelocityGradientProjectionSetup()
63 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, grad_velo_proj->dm, domain_label, label_value, hei…  in VelocityGradientProjectionSetup()
65 …PetscCall(CreateBasisFromPlex(ceed, grad_velo_proj->dm, domain_label, label_value, height, dm_fiel…  in VelocityGradientProjectionSetup()
88 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_assemble, "velocity gradient", grad_velo_proj->n…  in VelocityGradientProjectionSetup()
97 …eratorApplyContextCreate(user->dm, grad_velo_proj->dm, ceed, op_rhs_assemble, NULL, NULL, NULL, NU…  in VelocityGradientProjectionSetup()
100   PetscCall(CreateMassQFunction(ceed, grad_velo_proj->num_comp, q_data_size, &qf_mass));  in VelocityGradientProjectionSetup()
108     MPI_Comm comm = PetscObjectComm((PetscObject)grad_velo_proj->dm);  in VelocityGradientProjectionSetup()
110     PetscCall(MatCreateCeed(grad_velo_proj->dm, grad_velo_proj->dm, op_mass, NULL, &mat_mass));  in VelocityGradientProjectionSetup()
112     PetscCall(KSPCreate(comm, &grad_velo_proj->ksp));  in VelocityGradientProjectionSetup()
113     PetscCall(KSPSetOptionsPrefix(grad_velo_proj->ksp, "velocity_gradient_projection_"));  in VelocityGradientProjectionSetup()
116       PetscCall(KSPGetPC(grad_velo_proj->ksp, &pc));  in VelocityGradientProjectionSetup()
119       PetscCall(KSPSetType(grad_velo_proj->ksp, KSPPREONLY));  in VelocityGradientProjectionSetup()
121     PetscCall(KSPSetFromOptions_WithMatCeed(grad_velo_proj->ksp, mat_mass));  in VelocityGradientProjectionSetup()
125   *pgrad_velo_proj = grad_velo_proj;  in VelocityGradientProjectionSetup()
136 PetscErrorCode VelocityGradientProjectionApply(NodalProjectionData grad_velo_proj, Vec Q_loc, Vec V…  in VelocityGradientProjectionApply()  argument
137   OperatorApplyContext l2_rhs_ctx = grad_velo_proj->l2_rhs_ctx;  in VelocityGradientProjectionApply()
143   PetscCall(KSPSolve(grad_velo_proj->ksp, VelocityGradient, VelocityGradient));  in VelocityGradientProjectionApply()