85   CeedScalar rho, P, u[3] = {0.};  in Exact_ShockTube()  local
98   q[1] = rho * u[0];  in Exact_ShockTube()
99   q[2] = rho * u[1];  in Exact_ShockTube()
100   q[3] = rho * u[2];  in Exact_ShockTube()
101   q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.;  in Exact_ShockTube()
109 …cobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar…  in ConvectiveFluxJacobian_Euler()
111   CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2];  // Velocity square  in ConvectiveFluxJacobian_Euler()
114       dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j];  in ConvectiveFluxJacobian_Euler()
117 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (…  in ConvectiveFluxJacobian_Euler()
118 …     = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.) * u[i] * u[k];  in ConvectiveFluxJacobian_Euler()
122     dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho);  in ConvectiveFluxJacobian_Euler()
123     dF[i][4][4] = u[i] * gamma;  in ConvectiveFluxJacobian_Euler()
156 …_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const CeedScalar u[3], const CeedScalar…  in Tau_spatial()
162     CeedScalar fastest_wave = fabs(u[i]) + sound_speed;  in Tau_spatial()
227     const CeedScalar u[3]     = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in EulerShockTube()  local
246         du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho;  in EulerShockTube()
256 …const CeedScalar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E…  in EulerShockTube()
268 … j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rh…  in EulerShockTube()
273 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[…  in EulerShockTube()
274                                     (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]);  in EulerShockTube()
279 …nt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[j][1] + u[…  in EulerShockTube()
317     ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma);  in EulerShockTube()
339     Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau);  in EulerShockTube()