213 PetscErrorCode NS_NEWTONIAN_IG(ProblemData problem, DM dm, void *ctx, SimpleBC bc) {  in NS_NEWTONIAN_IG()  argument
233   problem->dim                          = 3;  in NS_NEWTONIAN_IG()
234   problem->jac_data_size_sur            = 11;  in NS_NEWTONIAN_IG()
235   problem->compute_exact_solution_error = PETSC_FALSE;  in NS_NEWTONIAN_IG()
236   problem->print_info                   = PRINT_NEWTONIAN;  in NS_NEWTONIAN_IG()
237   problem->uses_newtonian               = PETSC_TRUE;  in NS_NEWTONIAN_IG()
283       problem->ics.qfunction                       = ICsNewtonianIG_Conserv;  in NS_NEWTONIAN_IG()
284       problem->ics.qfunction_loc                   = ICsNewtonianIG_Conserv_loc;  in NS_NEWTONIAN_IG()
285       problem->apply_vol_rhs.qfunction             = RHSFunction_Newtonian;  in NS_NEWTONIAN_IG()
286       problem->apply_vol_rhs.qfunction_loc         = RHSFunction_Newtonian_loc;  in NS_NEWTONIAN_IG()
287       problem->apply_vol_ifunction.qfunction       = IFunction_Newtonian_Conserv;  in NS_NEWTONIAN_IG()
288       problem->apply_vol_ifunction.qfunction_loc   = IFunction_Newtonian_Conserv_loc;  in NS_NEWTONIAN_IG()
289       problem->apply_vol_ijacobian.qfunction       = IJacobian_Newtonian_Conserv;  in NS_NEWTONIAN_IG()
290       problem->apply_vol_ijacobian.qfunction_loc   = IJacobian_Newtonian_Conserv_loc;  in NS_NEWTONIAN_IG()
291       problem->apply_inflow.qfunction              = BoundaryIntegral_Conserv;  in NS_NEWTONIAN_IG()
292       problem->apply_inflow.qfunction_loc          = BoundaryIntegral_Conserv_loc;  in NS_NEWTONIAN_IG()
293       problem->apply_inflow_jacobian.qfunction     = BoundaryIntegral_Jacobian_Conserv;  in NS_NEWTONIAN_IG()
294       problem->apply_inflow_jacobian.qfunction_loc = BoundaryIntegral_Jacobian_Conserv_loc;  in NS_NEWTONIAN_IG()
297       problem->ics.qfunction                       = ICsNewtonianIG_Prim;  in NS_NEWTONIAN_IG()
298       problem->ics.qfunction_loc                   = ICsNewtonianIG_Prim_loc;  in NS_NEWTONIAN_IG()
299       problem->apply_vol_ifunction.qfunction       = IFunction_Newtonian_Prim;  in NS_NEWTONIAN_IG()
300       problem->apply_vol_ifunction.qfunction_loc   = IFunction_Newtonian_Prim_loc;  in NS_NEWTONIAN_IG()
301       problem->apply_vol_ijacobian.qfunction       = IJacobian_Newtonian_Prim;  in NS_NEWTONIAN_IG()
302       problem->apply_vol_ijacobian.qfunction_loc   = IJacobian_Newtonian_Prim_loc;  in NS_NEWTONIAN_IG()
303       problem->apply_inflow.qfunction              = BoundaryIntegral_Prim;  in NS_NEWTONIAN_IG()
304       problem->apply_inflow.qfunction_loc          = BoundaryIntegral_Prim_loc;  in NS_NEWTONIAN_IG()
305       problem->apply_inflow_jacobian.qfunction     = BoundaryIntegral_Jacobian_Prim;  in NS_NEWTONIAN_IG()
306       problem->apply_inflow_jacobian.qfunction_loc = BoundaryIntegral_Jacobian_Prim_loc;  in NS_NEWTONIAN_IG()
309       problem->ics.qfunction                       = ICsNewtonianIG_Entropy;  in NS_NEWTONIAN_IG()
310       problem->ics.qfunction_loc                   = ICsNewtonianIG_Entropy_loc;  in NS_NEWTONIAN_IG()
311       problem->apply_vol_ifunction.qfunction       = IFunction_Newtonian_Entropy;  in NS_NEWTONIAN_IG()
312       problem->apply_vol_ifunction.qfunction_loc   = IFunction_Newtonian_Entropy_loc;  in NS_NEWTONIAN_IG()
313       problem->apply_vol_ijacobian.qfunction       = IJacobian_Newtonian_Entropy;  in NS_NEWTONIAN_IG()
314       problem->apply_vol_ijacobian.qfunction_loc   = IJacobian_Newtonian_Entropy_loc;  in NS_NEWTONIAN_IG()
315       problem->apply_inflow.qfunction              = BoundaryIntegral_Entropy;  in NS_NEWTONIAN_IG()
316       problem->apply_inflow.qfunction_loc          = BoundaryIntegral_Entropy_loc;  in NS_NEWTONIAN_IG()
317       problem->apply_inflow_jacobian.qfunction     = BoundaryIntegral_Jacobian_Entropy;  in NS_NEWTONIAN_IG()
318       problem->apply_inflow_jacobian.qfunction_loc = BoundaryIntegral_Jacobian_Entropy_loc;  in NS_NEWTONIAN_IG()
329   PetscInt dim = problem->dim;  in NS_NEWTONIAN_IG()
374 …if (stab == STAB_SUPG && !implicit) problem->create_mass_operator = CreateKSPMassOperator_Newtonia…  in NS_NEWTONIAN_IG()
407   problem->dm_scale = meter;  in NS_NEWTONIAN_IG()
443   PetscCallCeed(ceed, CeedQFunctionContextCreate(user->ceed, &problem->ics.qfunction_context));  in NS_NEWTONIAN_IG()
445 …CeedQFunctionContextSetData(problem->ics.qfunction_context, CEED_MEM_HOST, CEED_USE_POINTER, sizeo…  in NS_NEWTONIAN_IG()
446 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(problem->ics.qfunction_context, CEED_MEM_HO…  in NS_NEWTONIAN_IG()
447 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(problem->ics.qfunction_context, "evaluation…  in NS_NEWTONIAN_IG()
462   problem->apply_vol_rhs.qfunction_context = newtonian_ig_context;  in NS_NEWTONIAN_IG()
463 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_vol_if…  in NS_NEWTONIAN_IG()
464 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_vol_ij…  in NS_NEWTONIAN_IG()
465 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_inflow…  in NS_NEWTONIAN_IG()
466 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_inflow…  in NS_NEWTONIAN_IG()
468 …if (bc->num_freestream > 0) PetscCall(FreestreamBCSetup(problem, dm, ctx, newtonian_ig_ctx, &refer…  in NS_NEWTONIAN_IG()
469 …if (bc->num_outflow > 0) PetscCall(OutflowBCSetup(problem, dm, ctx, newtonian_ig_ctx, &reference));  in NS_NEWTONIAN_IG()
470   if (bc->num_slip > 0) PetscCall(SlipBCSetup(problem, dm, ctx, newtonian_ig_context));  in NS_NEWTONIAN_IG()
478 PetscErrorCode PRINT_NEWTONIAN(User user, ProblemData problem, AppCtx app_ctx) {  in PRINT_NEWTONIAN()  argument
484 …PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->apply_vol_rhs.qfunction_context, CEED_MEM…  in PRINT_NEWTONIAN()
490 …PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->apply_vol_rhs.qfunction_context, &new…  in PRINT_NEWTONIAN()