19 static PetscErrorCode PRINT_NEWTONIAN(Honee honee, ProblemData problem, AppCtx app_ctx) {  in PRINT_NEWTONIAN()  argument
25 …PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->apply_vol_rhs.qfctx, CEED_MEM_HOST, &newt…  in PRINT_NEWTONIAN()
27                         "  Problem:\n"  in PRINT_NEWTONIAN()
28                         "    Problem Name                       : %s\n"  in PRINT_NEWTONIAN()
31   PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->apply_vol_rhs.qfctx, &newt_ctx));  in PRINT_NEWTONIAN()
370 PetscErrorCode NS_NEWTONIAN_IG(ProblemData problem, DM dm, void *ctx) {  in NS_NEWTONIAN_IG()  argument
412   PetscOptionsBegin(comm, NULL, "Options for Newtonian Ideal Gas based problem", NULL);  in NS_NEWTONIAN_IG()
531   // Set problem information  in NS_NEWTONIAN_IG()
532   problem->num_comps_jac_data = 14;  in NS_NEWTONIAN_IG()
533   if (newtonian_ig_ctx->idl_enable) problem->num_comps_jac_data += 1;  in NS_NEWTONIAN_IG()
534   problem->compute_exact_solution_error = PETSC_FALSE;  in NS_NEWTONIAN_IG()
535   problem->print_info                   = PRINT_NEWTONIAN;  in NS_NEWTONIAN_IG()
536   problem->num_components               = 5;  in NS_NEWTONIAN_IG()
537   PetscCall(PetscMalloc1(problem->num_components, &problem->component_names));  in NS_NEWTONIAN_IG()
545 …problem->ics                 = (HoneeQFSpec){.qf_func_ptr = ICsNewtonianIG_Conserv, .qf_loc = ICsN…  in NS_NEWTONIAN_IG()
546 …problem->apply_vol_rhs       = (HoneeQFSpec){.qf_func_ptr = RHSFunction_Newtonian, .qf_loc = RHSFu…  in NS_NEWTONIAN_IG()
547 …problem->apply_vol_ifunction = (HoneeQFSpec){.qf_func_ptr = IFunction_Newtonian_Conserv, .qf_loc =…  in NS_NEWTONIAN_IG()
548 …problem->apply_vol_ijacobian = (HoneeQFSpec){.qf_func_ptr = IJacobian_Newtonian_Conserv, .qf_loc =…  in NS_NEWTONIAN_IG()
549 … = 0; i < 5; i++) PetscCall(PetscStrallocpy(conserv_component_names[i], &problem->component_names[…  in NS_NEWTONIAN_IG()
552 …problem->ics                 = (HoneeQFSpec){.qf_func_ptr = ICsNewtonianIG_Prim, .qf_loc = ICsNewt…  in NS_NEWTONIAN_IG()
553 …problem->apply_vol_ifunction = (HoneeQFSpec){.qf_func_ptr = IFunction_Newtonian_Prim, .qf_loc = IF…  in NS_NEWTONIAN_IG()
554 …problem->apply_vol_ijacobian = (HoneeQFSpec){.qf_func_ptr = IJacobian_Newtonian_Prim, .qf_loc = IJ…  in NS_NEWTONIAN_IG()
555 …for (PetscInt i = 0; i < 5; i++) PetscCall(PetscStrallocpy(prim_component_names[i], &problem->comp…  in NS_NEWTONIAN_IG()
558 …problem->ics                 = (HoneeQFSpec){.qf_func_ptr = ICsNewtonianIG_Entropy, .qf_loc = ICsN…  in NS_NEWTONIAN_IG()
559 …problem->apply_vol_ifunction = (HoneeQFSpec){.qf_func_ptr = IFunction_Newtonian_Entropy, .qf_loc =…  in NS_NEWTONIAN_IG()
560 …problem->apply_vol_ijacobian = (HoneeQFSpec){.qf_func_ptr = IJacobian_Newtonian_Entropy, .qf_loc =…  in NS_NEWTONIAN_IG()
561 … = 0; i < 5; i++) PetscCall(PetscStrallocpy(entropy_component_names[i], &problem->component_names[…  in NS_NEWTONIAN_IG()
565   problem->ics.qfctx           = ics_qfctx;  in NS_NEWTONIAN_IG()
566   problem->apply_vol_rhs.qfctx = newtonian_ig_qfctx;  in NS_NEWTONIAN_IG()
567 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_qfctx, &problem->apply_vol_ifun…  in NS_NEWTONIAN_IG()
568 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_qfctx, &problem->apply_vol_ijac…  in NS_NEWTONIAN_IG()
570 …if (stab == STAB_SUPG && !implicit) problem->create_mass_operator = CreateKSPMassOperator_Newtonia…  in NS_NEWTONIAN_IG()
611   for (PetscCount b = 0; b < problem->num_bc_defs; b++) {  in NS_NEWTONIAN_IG()
612     BCDefinition bc_def = problem->bc_defs[b];  in NS_NEWTONIAN_IG()
617       PetscCall(SlipBCSetup(bc_def, problem, dm, ctx, newtonian_ig_qfctx));  in NS_NEWTONIAN_IG()
619       PetscCall(FreestreamBCSetup(bc_def, problem, dm, ctx, newtonian_ig_ctx, &reference));  in NS_NEWTONIAN_IG()
621       PetscCall(OutflowBCSetup(bc_def, problem, dm, ctx, newtonian_ig_ctx, &reference));  in NS_NEWTONIAN_IG()