12 static PetscErrorCode VelocityGradientProjectionCreateDM(NodalProjectionData grad_velo_proj, Honee …  in VelocityGradientProjectionCreateDM()  argument
16   grad_velo_proj->num_comp = 9;  // 9 velocity gradient  in VelocityGradientProjectionCreateDM()
18   PetscCall(DMClone(honee->dm, &grad_velo_proj->dm));  in VelocityGradientProjectionCreateDM()
19   PetscCall(DMSetMatrixPreallocateSkip(grad_velo_proj->dm, PETSC_TRUE));  in VelocityGradientProjectionCreateDM()
20   PetscCall(PetscObjectSetName((PetscObject)grad_velo_proj->dm, "Velocity Gradient Projection"));  in VelocityGradientProjectionCreateDM()
22 …TRUE, PETSC_TRUE, honee->app_ctx->degree, 1, honee->app_ctx->q_extra, 1, &grad_velo_proj->num_comp,  in VelocityGradientProjectionCreateDM()
23                                grad_velo_proj->dm));  in VelocityGradientProjectionCreateDM()
25   PetscCall(DMGetLocalSection(grad_velo_proj->dm, &section));  in VelocityGradientProjectionCreateDM()
41   NodalProjectionData grad_velo_proj;  in VelocityGradientProjectionSetup()  local
49   PetscCall(PetscNew(&grad_velo_proj));  in VelocityGradientProjectionSetup()
50   PetscCall(VelocityGradientProjectionCreateDM(grad_velo_proj, honee, honee->app_ctx->degree));  in VelocityGradientProjectionSetup()
53   PetscCall(DMGetDimension(grad_velo_proj->dm, &dim));  in VelocityGradientProjectionSetup()
55 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, grad_velo_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFAU…  in VelocityGradientProjectionSetup()
57 …PetscCall(DMPlexCeedBasisCreate(ceed, grad_velo_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, …  in VelocityGradientProjectionSetup()
58 …PetscCall(QDataGet(ceed, grad_velo_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, &elem_restr_q…  in VelocityGradientProjectionSetup()
90 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_assemble, "velocity gradient", grad_velo_proj->n…  in VelocityGradientProjectionSetup()
99 …ratorApplyContextCreate(honee->dm, grad_velo_proj->dm, ceed, op_rhs_assemble, NULL, NULL, NULL, NU…  in VelocityGradientProjectionSetup()
113     PetscCall(HoneeMassQFunctionCreate(ceed, grad_velo_proj->num_comp, q_data_size, &qf_mass));  in VelocityGradientProjectionSetup()
119     PetscCall(MatCreateCeed(grad_velo_proj->dm, grad_velo_proj->dm, op_mass, NULL, &mat_mass));  in VelocityGradientProjectionSetup()
122     PetscCall(KSPCreate(PetscObjectComm((PetscObject)grad_velo_proj->dm), &grad_velo_proj->ksp));  in VelocityGradientProjectionSetup()
123     PetscCall(KSPSetOptionsPrefix(grad_velo_proj->ksp, "velocity_gradient_projection_"));  in VelocityGradientProjectionSetup()
126       PetscCall(KSPGetPC(grad_velo_proj->ksp, &pc));  in VelocityGradientProjectionSetup()
129       PetscCall(KSPSetType(grad_velo_proj->ksp, KSPPREONLY));  in VelocityGradientProjectionSetup()
131     PetscCall(KSPSetFromOptions_WithMatCeed(grad_velo_proj->ksp, mat_mass));  in VelocityGradientProjectionSetup()
138   *pgrad_velo_proj = grad_velo_proj;  in VelocityGradientProjectionSetup()
147 PetscErrorCode VelocityGradientProjectionApply(NodalProjectionData grad_velo_proj, Vec Q_loc, Vec V…  in VelocityGradientProjectionApply()  argument
148   OperatorApplyContext l2_rhs_ctx = grad_velo_proj->l2_rhs_ctx;  in VelocityGradientProjectionApply()
154   PetscCall(KSPSolve(grad_velo_proj->ksp, VelocityGradient, VelocityGradient));  in VelocityGradientProjectionApply()