215       for (CeedInt k = 0; k < 2; k++) {  in StatePhysicalGradientFromReference_ND()  local
218 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * …  in StatePhysicalGradientFromReference_ND()
220         grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var);  in StatePhysicalGradientFromReference_ND()
227       for (CeedInt k = 0; k < 3; k++) {  in StatePhysicalGradientFromReference_ND()  local
230 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * …  in StatePhysicalGradientFromReference_ND()
231                    grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k];  in StatePhysicalGradientFromReference_ND()
233         grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var);  in StatePhysicalGradientFromReference_ND()
293       v[f][i] = wdetJ * q_dot[f][i];                          // K Mass/transient term  in IFunction_AdvectionGeneric()
298       for (CeedInt k = 0; k < dim; k++) {  in IFunction_AdvectionGeneric()  local
299         div_u += grad_s[k].Y.velocity[j];  in IFunction_AdvectionGeneric()
318       for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] -= wdetJ * Fe_dXdx[k];  in IFunction_AdvectionGeneric()  local
366       v[f][i] = wdetJ * q_dot[f][i];                          // K Mass/transient term  in MassFunction_AdvectionGeneric()
433       for (CeedInt k = 0; k < dim; k++) {  in RHSFunction_AdvectionGeneric()  local
434         div_u += grad_s[k].Y.velocity[j];  in RHSFunction_AdvectionGeneric()
455       for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] += wdetJ * Fe_dXdx[k];  in RHSFunction_AdvectionGeneric()  local
563     for (CeedInt k = 0; k < dim; k++) Grad_v[k][i] = -wdetJ * F_diff_dXdx[k];  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()  local
632     for (CeedInt k = 0; k < dim; k++) v[k][i] = wdetJ * F_diff[k];  in DiffusiveFluxRHS_AdvDif_Generic()  local