14 static PetscErrorCode CreateKSPMassOperator_Unstabilized(Honee honee, CeedOperator *op_mass) {  in CreateKSPMassOperator_Unstabilized()  argument
15   Ceed                ceed = honee->ceed;  in CreateKSPMassOperator_Unstabilized()
28     PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(honee->op_rhs_ctx->op, &sub_ops));  in CreateKSPMassOperator_Unstabilized()
54 static PetscErrorCode CreateKSPMass(Honee honee, ProblemData problem) {  in CreateKSPMass()  argument
55   Ceed         ceed = honee->ceed;  in CreateKSPMass()
56   DM           dm   = honee->dm;  in CreateKSPMass()
60   if (problem->create_mass_operator) PetscCall(problem->create_mass_operator(honee, &op_mass));  in CreateKSPMass()
61   else PetscCall(CreateKSPMassOperator_Unstabilized(honee, &op_mass));  in CreateKSPMass()
73     PetscCall(KSPCreate(comm, &honee->mass_ksp));  in CreateKSPMass()
74     PetscCall(KSPSetOptionsPrefix(honee->mass_ksp, "mass_"));  in CreateKSPMass()
75     PetscCall(PetscObjectSetName((PetscObject)honee->mass_ksp, "Explicit Mass"));  in CreateKSPMass()
78       PetscCall(KSPGetPC(honee->mass_ksp, &pc));  in CreateKSPMass()
81       PetscCall(KSPSetType(honee->mass_ksp, KSPPREONLY));  in CreateKSPMass()
83     PetscCall(KSPSetFromOptions_WithMatCeed(honee->mass_ksp, mat_mass));  in CreateKSPMass()
92 PetscErrorCode SetupLibceed(Ceed ceed, DM dm, Honee honee, AppCtx app_ctx, ProblemData problem) {  in SetupLibceed()  argument
122     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_q, &honee->q_ceed, NULL));  in SetupLibceed()
123     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_q, &honee->q_dot_ceed, NULL));  in SetupLibceed()
124     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_q, &honee->g_ceed, NULL));  in SetupLibceed()
131     PetscCheck(honee->diff_flux_proj, honee->comm, PETSC_ERR_ARG_WRONGSTATE,  in SetupLibceed()
133 …PetscCall(DivDiffFluxProjectionGetOperatorFieldData(honee->diff_flux_proj, &elem_restr_diff_flux, …  in SetupLibceed()
154 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_ics, "evaluation time", &honee->phys->ics_…  in SetupLibceed()
155 …tscCall(OperatorApplyContextCreate(NULL, dm, honee->ceed, op_ics, x_coord, NULL, NULL, honee->Q_lo…  in SetupLibceed()
173 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_vol, "div F_diff", honee->diff_flux_proj->num_dif…  in SetupLibceed()
203 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_ifunction_vol, "div F_diff", honee->diff_flux_proj->n…  in SetupLibceed()
211 …d(ceed, CeedOperatorSetField(op_ifunction_vol, "q dot", elem_restr_q, basis_q, honee->q_dot_ceed));  in SetupLibceed()
249   if (!honee->phys->implicit) {  // RHS  in SetupLibceed()
256 …peratorApplyContextCreate(dm, dm, ceed, op_rhs, honee->q_ceed, honee->g_ceed, honee->Q_loc, NULL, …  in SetupLibceed()
259 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_rhs, "solution time", &honee->phys->soluti…  in SetupLibceed()
260 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_rhs, "timestep size", &honee->phys->timest…  in SetupLibceed()
263     PetscCall(CreateKSPMass(honee, problem));  in SetupLibceed()
264 …PetscCheck(app_ctx->sgs_model_type == SGS_MODEL_NONE, honee->comm, PETSC_ERR_SUP, "SGS modeling no…  in SetupLibceed()
269     PetscCallCeed(ceed, CeedOperatorCreateComposite(ceed, &honee->op_ifunction));  in SetupLibceed()
270     PetscCallCeed(ceed, CeedOperatorCompositeAddSub(honee->op_ifunction, op_ifunction_vol));  in SetupLibceed()
275 …num_bc_defs; b++) PetscCall(BCDefinitionAddOperators(problem->bc_defs[b], honee->op_ifunction, op_…  in SetupLibceed()
278 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(honee->op_ifunction, "solution time", &honee-…  in SetupLibceed()
279 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(honee->op_ifunction, "timestep size", &honee-…  in SetupLibceed()
282       PetscCall(MatCreateCeed(honee->dm, honee->dm, op_ijacobian, NULL, &honee->mat_ijacobian));  in SetupLibceed()
283       PetscCall(MatCeedSetLocalVectors(honee->mat_ijacobian, honee->Q_dot_loc, NULL));  in SetupLibceed()
286 … if (app_ctx->sgs_model_type == SGS_MODEL_DATA_DRIVEN) PetscCall(SgsDDSetup(ceed, honee, problem));  in SetupLibceed()
289   if (problem->use_strong_bc_ceed) PetscCall(SetupStrongBC_Ceed(ceed, dm, honee, problem));  in SetupLibceed()
290   if (app_ctx->sgs_train_enable) PetscCall(SGS_DD_TrainingSetup(ceed, honee));  in SetupLibceed()
291 …thod != DIV_DIFF_FLUX_PROJ_NONE) PetscCall(DivDiffFluxProjectionSetup(honee, honee->diff_flux_proj…  in SetupLibceed()