17 CEED_QFUNCTION_HELPER void dYFromTau(CeedScalar Y[5], CeedScalar Tau_d[3], CeedScalar dY[5]) {  in dYFromTau()
28 …d StabilizationMatrix(NewtonianIdealGasContext gas, State s, CeedScalar Tau_d[3], CeedScalar stron…  in StabilizationMatrix()
29                                                CeedScalar stab[5][3]) {  in StabilizationMatrix()
30   CeedScalar        dY[5];  in StabilizationMatrix()
40     CeedScalar dF_i[5];  in StabilizationMatrix()
46 … Stabilization(NewtonianIdealGasContext gas, State s, CeedScalar Tau_d[3], State ds[3], CeedScalar…  in Stabilization()
47                                          const CeedScalar body_force[5], CeedScalar stab[5][3]) {  in Stabilization()
49   CeedScalar strong_residual[5] = {0};  in Stabilization()
71 …au_diagPrim(NewtonianIdealGasContext gas, State s, const CeedScalar dXdx[3][3], const CeedScalar d…  in Tau_diagPrim()
72                                         CeedScalar Tau_d[3]) {  in Tau_diagPrim()
74   const CeedScalar Ctau_t = gas->Ctau_t;  in Tau_diagPrim()
75   const CeedScalar Ctau_v = gas->Ctau_v;  in Tau_diagPrim()
76   const CeedScalar Ctau_C = gas->Ctau_C;  in Tau_diagPrim()
77   const CeedScalar Ctau_M = gas->Ctau_M;  in Tau_diagPrim()
78   const CeedScalar Ctau_E = gas->Ctau_E;  in Tau_diagPrim()
79   const CeedScalar cv     = gas->cv;  in Tau_diagPrim()
80   const CeedScalar mu     = gas->mu;  in Tau_diagPrim()
81   const CeedScalar rho    = s.U.density;  in Tau_diagPrim()
83   CeedScalar tau;  in Tau_diagPrim()
84   CeedScalar dts;  in Tau_diagPrim()
85   CeedScalar fact;  in Tau_diagPrim()
87   CeedScalar gijd_mat[3][3] = {{0.}}, velocity_term;  in Tau_diagPrim()
93     CeedScalar gij_uj[3] = {0.};  in Tau_diagPrim()
98 …Square(dts) + velocity_term) + Ctau_v * Square(mu) * DotN((CeedScalar *)gijd_mat, (CeedScalar *)gi…  in Tau_diagPrim()