210 …ysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s,  in StatePhysicalGradientFromReference_ND()  argument
220         grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var);  in StatePhysicalGradientFromReference_ND()
233         grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var);  in StatePhysicalGradientFromReference_ND()
240 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd…  in Tau()  argument
245       MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX);  in Tau()
254       MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj);  in Tau()
255 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * Squa…  in Tau()
282     const State      s     = StateFromU(gas, qi);  in IFunction_AdvectionGeneric()  local
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad…  in IFunction_AdvectionGeneric()
303     MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX);  in IFunction_AdvectionGeneric()
304     CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim);  in IFunction_AdvectionGeneric()
310       for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = -wdetJ * s.U.E_total * uX[j];  in IFunction_AdvectionGeneric()
321     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in IFunction_AdvectionGeneric()
360     const State      s     = StateFromU(gas, qi);  in MassFunction_AdvectionGeneric()  local
374     MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX);  in MassFunction_AdvectionGeneric()
375     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in MassFunction_AdvectionGeneric()
417     const State      s     = StateFromU(gas, qi);  in RHSFunction_AdvectionGeneric()  local
422 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad…  in RHSFunction_AdvectionGeneric()
437     CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim);  in RHSFunction_AdvectionGeneric()
440     MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX);  in RHSFunction_AdvectionGeneric()
446       for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = wdetJ * s.U.E_total * uX[j];  in RHSFunction_AdvectionGeneric()
458     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in RHSFunction_AdvectionGeneric()