23     CeedScalar q[5];  in ICsNewtonianIG()  local
25     StateToQ(gas, s, q, state_var);  in ICsNewtonianIG()
26     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsNewtonianIG()
45   const CeedScalar(*q)[CEED_Q_VLA]     = (const CeedScalar(*)[CEED_Q_VLA])in[1];  in MassFunction_Newtonian()  local
54     const CeedScalar qi[5]     = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in MassFunction_Newtonian()
191   const CeedScalar(*q)[CEED_Q_VLA]        = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in RHSFunction_Newtonian()  local
206     for (int j = 0; j < 5; j++) U[j] = q[j][i];  in RHSFunction_Newtonian()
312   const CeedScalar(*q)[CEED_Q_VLA]        = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in IFunction_Newtonian()  local
325     const CeedScalar q_i[5]     = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in IFunction_Newtonian()
481   const CeedScalar(*q)[CEED_Q_VLA]       = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in BoundaryIntegral()  local
491     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in BoundaryIntegral()
599   const CeedScalar(*q)[CEED_Q_VLA]   = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in DivDiffusiveFluxVolumeRHS_NS()  local
609     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in DivDiffusiveFluxVolumeRHS_NS()
652   const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in DivDiffusiveFluxBoundaryRHS_NS()  local
662     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in DivDiffusiveFluxBoundaryRHS_NS()
701   const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in DiffusiveFluxRHS_NS()  local
711     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in DiffusiveFluxRHS_NS()