21 static PetscErrorCode PRINT_ADVECTION(Honee honee, ProblemData problem, AppCtx app_ctx) {  in PRINT_ADVECTION()  argument
22   MPI_Comm         comm = honee->comm;  in PRINT_ADVECTION()
23   Ceed             ceed = honee->ceed;  in PRINT_ADVECTION()
29   PetscCall(DMGetDimension(honee->dm, &dim));  in PRINT_ADVECTION()
76 PetscErrorCode CreateKSPMassOperator_AdvectionStabilized(Honee honee, CeedOperator *op_mass) {  in CreateKSPMassOperator_AdvectionStabilized()  argument
77   Ceed                 ceed = honee->ceed;  in CreateKSPMassOperator_AdvectionStabilized()
87   PetscCall(DMGetDimension(honee->dm, &dim));  in CreateKSPMassOperator_AdvectionStabilized()
93     PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_rhs_ctx->op, &sub_ops));  in CreateKSPMassOperator_AdvectionStabilized()
125   PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "q", elem_restr_q, basis_q, honee->q_ceed));  in CreateKSPMassOperator_AdvectionStabilized()
146 static PetscErrorCode DivDiffFluxProjectionCreateRHS_Direct_AdvDif(Honee honee, DivDiffFluxProjecti…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()  argument
147   Ceed                 ceed       = honee->ceed;  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
155   PetscCall(DMGetFieldNumComps(honee->dm, 0, &num_comp_q));  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
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()
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()
188 …PetscCheck(qf_rhs_volume, honee->comm, PETSC_ERR_SUP, "%s not valid for DM of dimension %" PetscIn…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
230     PetscCall(QDataBoundaryGradientGetNumComponents(honee->dm, &q_data_size));  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()
263 …PetscCall(QDataBoundaryGradientGet(ceed, honee->dm, face_orientation_label, orientation, &elem_res…  in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
298 static PetscErrorCode DivDiffFluxProjectionCreateRHS_Indirect_AdvDif(Honee honee, DivDiffFluxProjec…  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()  argument
299   Ceed                 ceed       = honee->ceed;  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
312   PetscCall(DMGetFieldNumComps(honee->dm, 0, &num_comp_q));  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
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()
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()
336 …PetscCheck(qf_rhs, honee->comm, PETSC_ERR_SUP, "%s not valid for DM of dimension %" PetscInt_FMT, …  in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif()
388   Honee                      honee = *(Honee *)ctx;  in NS_ADVECTION()  local
389   MPI_Comm                   comm  = honee->comm;  in NS_ADVECTION()
390   Ceed                       ceed  = honee->ceed;  in NS_ADVECTION()
408   CeedScalar     Ctau_a           = PetscPowScalarInt(honee->app_ctx->degree, 2);  in NS_ADVECTION()
409   CeedScalar     Ctau_d           = PetscPowScalarInt(honee->app_ctx->degree, 4);  in NS_ADVECTION()
495   Units units = honee->units;  in NS_ADVECTION()
520   honee->phys->implicit            = implicit;  in NS_ADVECTION()
531   advection_ctx->divFdiff_method   = honee->app_ctx->divFdiffproj_method;  in NS_ADVECTION()
533   PetscCallCeed(ceed, CeedQFunctionContextCreate(honee->ceed, &ics_qfctx));  in NS_ADVECTION()
537   PetscCallCeed(ceed, CeedQFunctionContextCreate(honee->ceed, &advection_qfctx));  in NS_ADVECTION()
570 …PetscCall(DivDiffFluxProjectionCreate(honee, honee->app_ctx->divFdiffproj_method, 1, &honee->diff_…  in NS_ADVECTION()
571   if (honee->diff_flux_proj) {  in NS_ADVECTION()
572     DivDiffFluxProjectionData diff_flux_proj = honee->diff_flux_proj;  in NS_ADVECTION()
579     switch (honee->diff_flux_proj->method) {  in NS_ADVECTION()
591 …SETERRQ(PetscObjectComm((PetscObject)honee->dm), PETSC_ERR_ARG_WRONG, "Should not reach here with …  in NS_ADVECTION()
592                 DivDiffFluxProjectionMethods[honee->app_ctx->divFdiffproj_method]);  in NS_ADVECTION()
607       honee_bc->honee              = honee;  in NS_ADVECTION()