22   CeedScalar momentum[3];  member
33   for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i];  in UnpackState_U()
73   for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density;  in StatePrimitiveFromConservative()
85     dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density;  in StatePrimitiveFromConservative_fwd()
158   for (int i = 0; i < 3; i++) U.momentum[i] = U.density * Y.velocity[i];  in StateConservativeFromPrimitive()
170     dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dY.velocity[i];  in StateConservativeFromPrimitive_fwd()
185   const CeedScalar e_kinetic  = .5 * Dot3(U.momentum, U.momentum) / U.density;  in StateEntropyFromConservative()
191   for (int i = 0; i < 3; i++) V.S_momentum[i] = U.momentum[i] / p;  in StateEntropyFromConservative()
199   const CeedScalar e_kinetic   = .5 * Dot3(s.U.momentum, s.U.momentum) / s.U.density;  in StateEntropyFromConservative_fwd()
200 …const CeedScalar de_kinetic  = (Dot3(s.U.momentum, dU.momentum) - e_kinetic * dU.density) / s.U.de…  in StateEntropyFromConservative_fwd()
209     dV.S_momentum[i] = (dU.momentum[i] - s.U.momentum[i] * dp / p) / p;  in StateEntropyFromConservative_fwd()
224   for (int i = 0; i < 3; i++) U.momentum[i] = U.density * velocity[i];  in StateConservativeFromEntropy()
244 …for (int i = 0; i < 3; i++) dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dvelocit…  in StateConservativeFromEntropy_fwd()
274     for (int j = 0; j < 3; j++) R.momentum[j] += a[i] * X[i].momentum[j];  in StateConservativeMult()
284 …for (int i = 0; i < 3; i++) R.momentum[i] = a * X.momentum[i] + b * Y.momentum[i] + c * Z.momentum…  in StateConservativeAXPBYPCZ()
331   s.U.momentum[0] = U[1];  in StateFromU()
332   s.U.momentum[1] = U[2];  in StateFromU()
333   s.U.momentum[2] = U[3];  in StateFromU()
342   ds.U.momentum[0] = dU[1];  in StateFromU_fwd()
343   ds.U.momentum[1] = dU[2];  in StateFromU_fwd()
344   ds.U.momentum[2] = dU[3];  in StateFromU_fwd()
432     Flux[i].density = s.U.momentum[i];  in FluxInviscid()
433 …for (CeedInt j = 0; j < 3; j++) Flux[i].momentum[j] = s.U.momentum[i] * s.Y.velocity[j] + s.Y.pres…  in FluxInviscid()
440     dFlux[i].density = ds.U.momentum[i];  in FluxInviscid_fwd()
442 …dFlux[i].momentum[j] = ds.U.momentum[i] * s.Y.velocity[j] + s.U.momentum[i] * ds.Y.velocity[j] + d…  in FluxInviscid_fwd()
453     for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += Flux[i].momentum[j] * normal[i];  in FluxInviscidDotNormal()
464     for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += dFlux[i].momentum[j] * normal[i];  in FluxInviscidDotNormal_fwd()
484     for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j];  in FluxTotal()
495       Flux[k + 1] += (F_inviscid[j].momentum[k] - stress[k][j]) * normal[j];  in FluxTotal_Boundary()
504   for (CeedInt k = 0; k < 3; k++) Flux[k + 1] = F_inviscid_normal.momentum[k];  in FluxTotal_RiemannBoundary()