112   for (CeedInt i = 0; i < 3; i++) {                           // Jacobian matrices for 3 directions  in ConvectiveFluxJacobian_Euler()  local
114       dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j];  in ConvectiveFluxJacobian_Euler()
116         dF[i][0][k + 1]     = ((i == k) ? 1. : 0.);  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 …dF[i][4][k + 1]     = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.…  in ConvectiveFluxJacobian_Euler()
120       dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.);  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()
158   for (CeedInt i = 0; i < 3; i++) {  in Tau_spatial()  local
160     CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i]));  in Tau_spatial()
162     CeedScalar fastest_wave = fabs(u[i]) + sound_speed;  in Tau_spatial()
163     Tau_x[i]                = c_tau * h / fastest_wave;  in Tau_spatial()
174   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ICsShockTube()  local
175     const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]};  in ICsShockTube()
180     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsShockTube()
223   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in EulerShockTube()  local
226     const CeedScalar rho      = q[0][i];  in EulerShockTube()
227     const CeedScalar u[3]     = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in EulerShockTube()
228     const CeedScalar E        = q[4][i];  in EulerShockTube()
229     const CeedScalar drho[3]  = {dq[0][0][i], dq[1][0][i], dq[2][0][i]};  in EulerShockTube()
231         {dq[0][1][i], dq[1][1][i], dq[2][1][i]},  in EulerShockTube()
232         {dq[0][2][i], dq[1][2][i], dq[2][2][i]},  in EulerShockTube()
233         {dq[0][3][i], dq[1][3][i], dq[2][3][i]}  in EulerShockTube()
235     const CeedScalar dE[3] = {dq[0][4][i], dq[1][4][i], dq[2][4][i]};  in EulerShockTube()
237     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in EulerShockTube()
262       v[j][i] = 0;  in EulerShockTube()
263       for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0;  in EulerShockTube()
268 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd…  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()
279 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[…  in EulerShockTube()
304       for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * nu_shock * drhodx[j];  in EulerShockTube()
307         for (CeedInt j = 0; j < 3; j++) dv[j][k][i] -= wdetJ * nu_shock * du[k][j];  in EulerShockTube()
310       for (CeedInt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * nu_shock * dEdx[j];  in EulerShockTube()
354 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in EulerShockTube()