23   Ceed             ceed = honee->ceed;  in PRINT_ADVECTION()  local
30   PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->ics.qfctx, CEED_MEM_HOST, &setup_ctx));  in PRINT_ADVECTION()
31 …PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->apply_vol_rhs.qfctx, CEED_MEM_HOST, &adve…  in PRINT_ADVECTION()
68   PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->ics.qfctx, &setup_ctx));  in PRINT_ADVECTION()
69 …PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->apply_vol_rhs.qfctx, &advection_ctx));  in PRINT_ADVECTION()
77   Ceed                 ceed = honee->ceed;  in CreateKSPMassOperator_AdvectionStabilized()  local
93     PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_rhs_ctx->op, &sub_ops));  in CreateKSPMassOperator_AdvectionStabilized()
94     PetscCallCeed(ceed, CeedOperatorGetFieldByName(sub_ops[sub_op_index], "q", &op_field));  in CreateKSPMassOperator_AdvectionStabilized()
95     PetscCallCeed(ceed, CeedOperatorFieldGetData(op_field, NULL, &elem_restr_q, &basis_q, NULL));  in CreateKSPMassOperator_AdvectionStabilized()
96     PetscCallCeed(ceed, CeedOperatorGetFieldByName(sub_ops[sub_op_index], "qdata", &op_field));  in CreateKSPMassOperator_AdvectionStabilized()
97     PetscCallCeed(ceed, CeedOperatorFieldGetData(op_field, NULL, &elem_restr_qd, NULL, &q_data));  in CreateKSPMassOperator_AdvectionStabilized()
99     PetscCallCeed(ceed, CeedOperatorGetContext(sub_ops[sub_op_index], &qfctx));  in CreateKSPMassOperator_AdvectionStabilized()
102   PetscCallCeed(ceed, CeedElemRestrictionGetNumComponents(elem_restr_q, &num_comp_q));  in CreateKSPMassOperator_AdvectionStabilized()
103   PetscCallCeed(ceed, CeedElemRestrictionGetNumComponents(elem_restr_qd, &q_data_size));  in CreateKSPMassOperator_AdvectionStabilized()
107 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, MassFunction_Advection2D, MassFunction_Ad…  in CreateKSPMassOperator_AdvectionStabilized()
110 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, MassFunction_Advection, MassFunction_Adve…  in CreateKSPMassOperator_AdvectionStabilized()
114   PetscCallCeed(ceed, CeedQFunctionSetContext(qf_mass, qfctx));  in CreateKSPMassOperator_AdvectionStabilized()
115   PetscCallCeed(ceed, CeedQFunctionSetUserFlopsEstimate(qf_mass, 0));  in CreateKSPMassOperator_AdvectionStabilized()
116   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "q_dot", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_AdvectionStabilized()
117   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "q", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_AdvectionStabilized()
118   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "qdata", q_data_size, CEED_EVAL_NONE));  in CreateKSPMassOperator_AdvectionStabilized()
119   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "v", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_AdvectionStabilized()
120   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "Grad_v", 5 * dim, CEED_EVAL_GRAD));  in CreateKSPMassOperator_AdvectionStabilized()
122   PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_mass, NULL, NULL, op_mass));  in CreateKSPMassOperator_AdvectionStabilized()
123 …PetscCallCeed(ceed, CeedOperatorSetName(*op_mass, "RHS Mass Operator, Advection-Diffusion Stabiliz…  in CreateKSPMassOperator_AdvectionStabilized()
124 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "q_dot", elem_restr_q, basis_q, CEED_VECTOR_ACT…  in CreateKSPMassOperator_AdvectionStabilized()
125   PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "q", elem_restr_q, basis_q, honee->q_ceed));  in CreateKSPMassOperator_AdvectionStabilized()
126 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "qdata", elem_restr_qd, CEED_BASIS_NONE, q_data…  in CreateKSPMassOperator_AdvectionStabilized()
127 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "v", elem_restr_q, basis_q, CEED_VECTOR_ACTIVE)…  in CreateKSPMassOperator_AdvectionStabilized()
128 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "Grad_v", elem_restr_q, basis_q, CEED_VECTOR_AC…  in CreateKSPMassOperator_AdvectionStabilized()
130   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in CreateKSPMassOperator_AdvectionStabilized()
131   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qd));  in CreateKSPMassOperator_AdvectionStabilized()
132   PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in CreateKSPMassOperator_AdvectionStabilized()
133   PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in CreateKSPMassOperator_AdvectionStabilized()
134   PetscCallCeed(ceed, CeedQFunctionContextDestroy(&qfctx));  in CreateKSPMassOperator_AdvectionStabilized()
135   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_mass));  in CreateKSPMassOperator_AdvectionStabilized()
147   Ceed                 ceed       = honee->ceed;  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()  local
161 …if (honee->op_ifunction) PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_ifunction, …  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
162     else PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_rhs_ctx->op, &sub_ops));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
163     PetscCallCeed(ceed, CeedOperatorGetContext(sub_ops[sub_op_index], &advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
165   PetscCallCeed(ceed, CeedOperatorCreateComposite(ceed, op_rhs));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
174 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
175 …PetscCall(DMPlexCeedBasisCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, 0, 0, &ba…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
177 …PetscCall(QDataGet(ceed, projection->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, &elem_restr_qd, &…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
180 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DivDiffusiveFluxVolumeRHS_AdvDif_2D, DivD…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
184 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DivDiffusiveFluxVolumeRHS_AdvDif_3D, DivD…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
190     PetscCallCeed(ceed, CeedQFunctionSetContext(qf_rhs_volume, advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
191 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_volume, "Grad_q", num_comp_q * dim, CEED_EVAL_GRA…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
192     PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_volume, "qdata", q_data_size, CEED_EVAL_NONE));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
193 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_volume, "diffusive flux RHS", projection->num_co…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
195     PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_rhs_volume, NULL, NULL, &op_rhs_volume));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
196 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_volume, "Grad_q", elem_restr_q, basis_q, CEED_VECT…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
197 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_volume, "qdata", elem_restr_qd, CEED_BASIS_NONE, q…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
198 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_volume, "diffusive flux RHS", elem_restr_diff_flux…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
200     PetscCallCeed(ceed, CeedOperatorCompositeAddSub(*op_rhs, op_rhs_volume));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
202     PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
203     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qd));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
204     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
205     PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
206     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_diff_flux_volume));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
207     PetscCallCeed(ceed, CeedBasisDestroy(&basis_diff_flux));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
208     PetscCallCeed(ceed, CeedOperatorDestroy(&op_rhs_volume));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
209     PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_rhs_volume));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
221 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DivDiffusiveFluxBoundaryRHS_AdvDif_2D, Di…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
225 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DivDiffusiveFluxBoundaryRHS_AdvDif_3D, Di…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
231     PetscCallCeed(ceed, CeedQFunctionSetContext(qf_rhs_boundary, advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
232 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_boundary, "Grad_q", num_comp_q * dim, CEED_EVAL_G…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
233 … PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_boundary, "qdata", q_data_size, CEED_EVAL_NONE));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
234 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_boundary, "diffusive flux RHS", projection->num_…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
258 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, honee->dm, face_orientation_label, orientation, he…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
259 …PetscCall(DMPlexCeedBasisCellToFaceCreate(ceed, honee->dm, face_orientation_label, orientation, or…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
260 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, projection->dm, face_orientation_label, orientatio…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
262 …PetscCall(DMPlexCeedBasisCreate(ceed, projection->dm, face_orientation_label, orientation, height_…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
263 …PetscCall(QDataBoundaryGradientGet(ceed, honee->dm, face_orientation_label, orientation, &elem_res…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
265 …      PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_rhs_boundary, NULL, NULL, &op_rhs_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
266 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_boundary, "Grad_q", elem_restr_q, basis_q, CEED_VE…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
267 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_boundary, "qdata", elem_restr_qdata, CEED_BASIS_NO…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
268 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_boundary, "diffusive flux RHS", elem_restr_diff_fl…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
271         PetscCallCeed(ceed, CeedOperatorCompositeAddSub(*op_rhs, op_rhs_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
273         PetscCallCeed(ceed, CeedOperatorDestroy(&op_rhs_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
274         PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qdata));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
275         PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
276         PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_diff_flux_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
277         PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
278         PetscCallCeed(ceed, CeedBasisDestroy(&basis_diff_flux_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
279         PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
284     PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_rhs_boundary));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
287   PetscCallCeed(ceed, CeedQFunctionContextDestroy(&advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
299   Ceed                 ceed       = honee->ceed;  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()  local
318 …if (honee->op_ifunction) PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_ifunction, …  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
319     else PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_rhs_ctx->op, &sub_ops));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
320     PetscCallCeed(ceed, CeedOperatorGetContext(sub_ops[sub_op_index], &advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
322 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
323 …PetscCall(DMPlexCeedBasisCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, 0, 0, &ba…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
324 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, projection->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_V…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
325 …PetscCall(DMPlexCeedBasisCreate(ceed, projection->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, heig…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
326 …PetscCall(QDataGet(ceed, projection->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, &elem_restr_qd, &…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
330 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DiffusiveFluxRHS_AdvDif_2D, DiffusiveFlux…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
333 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DiffusiveFluxRHS_AdvDif_3D, DiffusiveFlux…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
338   PetscCallCeed(ceed, CeedQFunctionSetContext(qf_rhs, advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
339   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "Grad_q", num_comp_q * dim, CEED_EVAL_GRAD));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
340   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "qdata", q_data_size, CEED_EVAL_NONE));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
341 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs, "F_diff RHS", projection->num_comp, CEED_EVAL_I…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
343   PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_rhs, NULL, NULL, op_rhs));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
344 …PetscCallCeed(ceed, CeedOperatorSetField(*op_rhs, "Grad_q", elem_restr_q, basis_q, CEED_VECTOR_ACT…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
345 …PetscCallCeed(ceed, CeedOperatorSetField(*op_rhs, "qdata", elem_restr_qd, CEED_BASIS_NONE, q_data)…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
346 …PetscCallCeed(ceed, CeedOperatorSetField(*op_rhs, "F_diff RHS", elem_restr_diff_flux, basis_diff_f…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
348   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_rhs));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
349   PetscCallCeed(ceed, CeedQFunctionContextDestroy(&advection_qfctx));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
350   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
351   PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
352   PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
353   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qd));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
354   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_diff_flux));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
355   PetscCallCeed(ceed, CeedBasisDestroy(&basis_diff_flux));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
390   Ceed                       ceed  = honee->ceed;  in NS_ADVECTION()  local
533   PetscCallCeed(ceed, CeedQFunctionContextCreate(honee->ceed, &ics_qfctx));  in NS_ADVECTION()
534 …PetscCallCeed(ceed, CeedQFunctionContextSetData(ics_qfctx, CEED_MEM_HOST, CEED_USE_POINTER, sizeof…  in NS_ADVECTION()
535 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(ics_qfctx, CEED_MEM_HOST, FreeContextPetsc)…  in NS_ADVECTION()
537   PetscCallCeed(ceed, CeedQFunctionContextCreate(honee->ceed, &advection_qfctx));  in NS_ADVECTION()
538 …PetscCallCeed(ceed, CeedQFunctionContextSetData(advection_qfctx, CEED_MEM_HOST, CEED_USE_POINTER, …  in NS_ADVECTION()
539 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(advection_qfctx, CEED_MEM_HOST, FreeContext…  in NS_ADVECTION()
540 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(advection_qfctx, "timestep size", offsetof(…  in NS_ADVECTION()
564 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(advection_qfctx, &problem->apply_vol_ifuncti…  in NS_ADVECTION()
606       PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(advection_qfctx, &honee_bc->qfctx));  in NS_ADVECTION()