179   for (CeedInt i = 0; i < 3; i++) {                           // Jacobian matrices for 3 directions  in ConvectiveFluxJacobian_Euler()  local
181       dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j];  in ConvectiveFluxJacobian_Euler()
183         dF[i][0][k + 1]     = ((i == k) ? 1. : 0.);  in ConvectiveFluxJacobian_Euler()
184 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (…  in ConvectiveFluxJacobian_Euler()
185 …dF[i][4][k + 1]     = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.…  in ConvectiveFluxJacobian_Euler()
187       dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.);  in ConvectiveFluxJacobian_Euler()
189     dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho);  in ConvectiveFluxJacobian_Euler()
190     dF[i][4][4] = u[i] * gamma;  in ConvectiveFluxJacobian_Euler()
211   for (CeedInt i = 0; i < 3; i++) {  in Tau_spatial()  local
213     CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i]));  in Tau_spatial()
215     CeedScalar fastest_wave = fabs(u[i]) + sound_speed;  in Tau_spatial()
216     Tau_x[i]                = c_tau * h / fastest_wave;  in Tau_spatial()
229   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ICsEuler()  local
230     const CeedScalar x[]  = {X[0][i], X[1][i], X[2][i]};  in ICsEuler()
234     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsEuler()
274   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Euler()  local
277     const CeedScalar rho      = q[0][i];  in Euler()
278     const CeedScalar u[3]     = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in Euler()
279     const CeedScalar E        = q[4][i];  in Euler()
280     const CeedScalar drho[3]  = {dq[0][0][i], dq[1][0][i], dq[2][0][i]};  in Euler()
282         {dq[0][1][i], dq[1][1][i], dq[2][1][i]},  in Euler()
283         {dq[0][2][i], dq[1][2][i], dq[2][2][i]},  in Euler()
284         {dq[0][3][i], dq[1][3][i], dq[2][3][i]}  in Euler()
286     const CeedScalar dE[3] = {dq[0][4][i], dq[1][4][i], dq[2][4][i]};  in Euler()
288     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in Euler()
311       v[j][i] = 0.;  in Euler()
312       for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.;  in Euler()
317 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd…  in Euler()
322 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u…  in Euler()
328 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[…  in Euler()
370 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in Euler()
395   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in IFunction_Euler()  local
398     const CeedScalar rho      = q[0][i];  in IFunction_Euler()
399     const CeedScalar u[3]     = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in IFunction_Euler()
400     const CeedScalar E        = q[4][i];  in IFunction_Euler()
401     const CeedScalar drho[3]  = {dq[0][0][i], dq[1][0][i], dq[2][0][i]};  in IFunction_Euler()
403         {dq[0][1][i], dq[1][1][i], dq[2][1][i]},  in IFunction_Euler()
404         {dq[0][2][i], dq[1][2][i], dq[2][2][i]},  in IFunction_Euler()
405         {dq[0][3][i], dq[1][3][i], dq[2][3][i]}  in IFunction_Euler()
407     const CeedScalar dE[3] = {dq[0][4][i], dq[1][4][i], dq[2][4][i]};  in IFunction_Euler()
409     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in IFunction_Euler()
431       v[j][i] = 0.;  in IFunction_Euler()
432       for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.;  in IFunction_Euler()
435     for (CeedInt j = 0; j < 5; j++) v[j][i] += wdetJ * q_dot[j][i];  in IFunction_Euler()
439 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd…  in IFunction_Euler()
444 …dv[k][j + 1][i] -= wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u…  in IFunction_Euler()
450 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[…  in IFunction_Euler()
475     for (CeedInt j = 0; j < 5; j++) strong_res[j] = q_dot[j][i] + strong_conv[j];  in IFunction_Euler()
496 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in IFunction_Euler()
507 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in IFunction_Euler()
534     for (CeedInt i = 0; i < 3; i++) mean_velocity[i] = 0.;  in TravelingVortex_Inflow()  local
541   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in TravelingVortex_Inflow()  local
543     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, normal);  in TravelingVortex_Inflow()
550     for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.;  in TravelingVortex_Inflow()
562       v[0][i] -= wdetJb * rho_inlet * face_normal;  in TravelingVortex_Inflow()
565 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_inlet * face_normal * mean_velocity[j…  in TravelingVortex_Inflow()
568       v[4][i] -= wdetJb * face_normal * (E_inlet + P_inlet);  in TravelingVortex_Inflow()
591   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Euler_Outflow()  local
594     const CeedScalar rho  = q[0][i];  in Euler_Outflow()
595     const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in Euler_Outflow()
596     const CeedScalar E    = q[4][i];  in Euler_Outflow()
599     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, normal);  in Euler_Outflow()
606     for (CeedInt j = 0; j < 5; j++) v[j][i] = 0;  in Euler_Outflow()
615       v[0][i] -= wdetJb * rho * u_normal;  in Euler_Outflow()
618 …   for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P);  in Euler_Outflow()
621       v[4][i] -= wdetJb * u_normal * (E + P);  in Euler_Outflow()