276   CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1];  in IFunction_AdvectionGeneric()  local
292       for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0;  // No Change in density or momentum  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()
318       for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] -= wdetJ * Fe_dXdx[k];  in IFunction_AdvectionGeneric()
327           grad_v[j][4][i] += wdetJ * TauS * uX[j] * strong_conv;  in IFunction_AdvectionGeneric()
331           grad_v[j][4][i] += wdetJ * TauS * uX[j] * (q_dot[4][i] + strong_conv + divFdiff_i);  in IFunction_AdvectionGeneric()
353   CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1];  in MassFunction_AdvectionGeneric()  local
365       for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0;  // No Change in density or momentum  in MassFunction_AdvectionGeneric()
380           grad_v[j][4][i] = 0;  in MassFunction_AdvectionGeneric()
383           grad_v[j][4][i] = wdetJ * TauS * q_dot[4][i] * uX[j];  in MassFunction_AdvectionGeneric()
411   CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1];  in RHSFunction_AdvectionGeneric()  local
427       for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0;  // No Change in density or momentum  in RHSFunction_AdvectionGeneric()
444       for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = 0;  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()
455       for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] += wdetJ * Fe_dXdx[k];  in RHSFunction_AdvectionGeneric()
466           grad_v[j][4][i] -= wdetJ * TauS * (strong_conv + divFdiff_i) * uX[j];  in RHSFunction_AdvectionGeneric()