284     CeedScalar wdetJ, dXdx[9];  in IFunction_AdvectionGeneric()  local
285     QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx);  in IFunction_AdvectionGeneric()
293       v[f][i] = wdetJ * q_dot[f][i];                          // K Mass/transient term  in IFunction_AdvectionGeneric()
306     v[4][i] = wdetJ * q_dot[4][i];  // transient part (ALWAYS)  in IFunction_AdvectionGeneric()
308       v[4][i] += wdetJ * strong_conv;  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()
361     CeedScalar       wdetJ, dXdx[9];  in MassFunction_AdvectionGeneric()  local
362     QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx);  in MassFunction_AdvectionGeneric()
366       v[f][i] = wdetJ * q_dot[f][i];                          // K Mass/transient term  in MassFunction_AdvectionGeneric()
370     v[4][i] = wdetJ * q_dot[4][i];  in MassFunction_AdvectionGeneric()
383           grad_v[j][4][i] = wdetJ * TauS * q_dot[4][i] * uX[j];  in MassFunction_AdvectionGeneric()
419     CeedScalar wdetJ, dXdx[9];  in RHSFunction_AdvectionGeneric()  local
420     QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx);  in RHSFunction_AdvectionGeneric()
443       v[4][i] = -wdetJ * strong_conv;  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()
549     CeedScalar wdetJ, dXdx[9], F_diff[3] = {0.};  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()  local
551     QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx);  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
563     for (CeedInt k = 0; k < dim; k++) Grad_v[k][i] = -wdetJ * F_diff_dXdx[k];  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
586     CeedScalar wdetJ, normal[3], dXdx[9], F_diff[3] = {0.};  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()  local
588     QdataBoundaryGradientUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx, normal);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
598     v[i] = wdetJ * DotN(F_diff, normal, dim);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
621     CeedScalar wdetJ, dXdx[9], F_diff[3] = {0.};  in DiffusiveFluxRHS_AdvDif_Generic()  local
623     QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx);  in DiffusiveFluxRHS_AdvDif_Generic()
632     for (CeedInt k = 0; k < dim; k++) v[k][i] = wdetJ * F_diff[k];  in DiffusiveFluxRHS_AdvDif_Generic()