477   Physics_Euler *eu   = (Physics_Euler *)phys->data;  in PhysicsSolution_Euler()  local
486   gamma = eu->gamma;  in PhysicsSolution_Euler()
488   if (eu->type == EULER_IV_SHOCK || eu->type == EULER_SS_SHOCK) {  in PhysicsSolution_Euler()
495     if (x[0] < 0.0 + x[1] * eu->itana) {  in PhysicsSolution_Euler()
498         amach     = eu->amach;  in PhysicsSolution_Euler()
511       uu->r = eu->rhoR;  in PhysicsSolution_Euler()
514   } else if (eu->type == EULER_SHOCK_TUBE) {  in PhysicsSolution_Euler()
523   } else if (eu->type == EULER_LINEAR_WAVE) {  in PhysicsSolution_Euler()
525   } else SETERRQ(mod->comm, PETSC_ERR_SUP, "Unknown type %d", eu->type);  in PhysicsSolution_Euler()
541   Physics_Euler   *eu   = (Physics_Euler *)phys->data;  in PhysicsBoundary_Euler_Wall()  local
552   if (eu->type == EULER_LINEAR_WAVE) { /* debug */  in PhysicsBoundary_Euler_Wall()
563   Physics_Euler   *eu   = (Physics_Euler *)phys->data;  in PhysicsFunctional_Euler()  local
568   f[eu->monitor.Density]  = x->r;  in PhysicsFunctional_Euler()
569   f[eu->monitor.Momentum] = NormDIM(x->ru);  in PhysicsFunctional_Euler()
570   f[eu->monitor.Energy]   = x->E;  in PhysicsFunctional_Euler()
571   f[eu->monitor.Speed]    = NormDIM(x->ru) / x->r;  in PhysicsFunctional_Euler()
572   PetscCall(Pressure_PG(eu->gamma, x, &p));  in PhysicsFunctional_Euler()
573   f[eu->monitor.Pressure] = p;  in PhysicsFunctional_Euler()
579   Physics_Euler *eu = (Physics_Euler *)phys->data;  in SetUpBC_Euler()  local
584   if (eu->type == EULER_LINEAR_WAVE) {  in SetUpBC_Euler()
598   Physics_Euler *eu;  in CreateQFunctionContext_Euler()  local
601   PetscCall(PetscCalloc1(1, &eu));  in CreateQFunctionContext_Euler()
603   eu->gamma = in->gamma;  in CreateQFunctionContext_Euler()
606 …tscCallCEED(CeedQFunctionContextSetData(*qfCtx, CEED_MEM_HOST, CEED_USE_POINTER, sizeof(*eu), eu));  in CreateQFunctionContext_Euler()
632   Physics_Euler *eu;  in PhysicsCreate_Euler()  local
637   PetscCall(PetscNew(&eu));  in PhysicsCreate_Euler()
638   phys->data     = eu;  in PhysicsCreate_Euler()
654     eu->gamma = 1.4;  in PhysicsCreate_Euler()
655     eu->amach = 2.02;  in PhysicsCreate_Euler()
656     eu->rhoR  = 3.0;  in PhysicsCreate_Euler()
657     eu->itana = 0.57735026918963; /* angle of Euler self similar (SS) shock */  in PhysicsCreate_Euler()
661 …  PetscCall(PetscOptionsReal("-eu_gamma", "Heat capacity ratio", "", eu->gamma, &eu->gamma, NULL));  in PhysicsCreate_Euler()
662     PetscCall(PetscOptionsReal("-eu_amach", "Shock speed (Mach)", "", eu->amach, &eu->amach, NULL));  in PhysicsCreate_Euler()
663 …PetscCall(PetscOptionsReal("-eu_rho2", "Density right of discontinuity", "", eu->rhoR, &eu->rhoR, …  in PhysicsCreate_Euler()
666     eu->itana = 1. / PetscTanReal(alpha * PETSC_PI / 180.0);  in PhysicsCreate_Euler()
671       eu->type = EULER_LINEAR_WAVE;  in PhysicsCreate_Euler()
677         eu->type = EULER_IV_SHOCK;  in PhysicsCreate_Euler()
682           eu->type = EULER_SS_SHOCK;  in PhysicsCreate_Euler()
687           eu->type = EULER_SHOCK_TUBE;  in PhysicsCreate_Euler()
696 …PetscCall(ModelFunctionalRegister(mod, "Speed", &eu->monitor.Speed, PhysicsFunctional_Euler, phys)…  in PhysicsCreate_Euler()
697 …PetscCall(ModelFunctionalRegister(mod, "Energy", &eu->monitor.Energy, PhysicsFunctional_Euler, phy…  in PhysicsCreate_Euler()
698 …PetscCall(ModelFunctionalRegister(mod, "Density", &eu->monitor.Density, PhysicsFunctional_Euler, p…  in PhysicsCreate_Euler()
699 …PetscCall(ModelFunctionalRegister(mod, "Momentum", &eu->monitor.Momentum, PhysicsFunctional_Euler,…  in PhysicsCreate_Euler()
700 …PetscCall(ModelFunctionalRegister(mod, "Pressure", &eu->monitor.Pressure, PhysicsFunctional_Euler,…  in PhysicsCreate_Euler()