26     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsNewtonianIG()  local
74     for (CeedInt j = 0; j < 5; j++) {  in MassFunction_Newtonian()  local
76 …Grad_v[k][j][i] = wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[k…  in MassFunction_Newtonian()
206     for (int j = 0; j < 5; j++) U[j] = q[j][i];  in RHSFunction_Newtonian()  local
226     for (CeedInt j = 0; j < 5; j++) {  in RHSFunction_Newtonian()  local
227 …nt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1]…  in RHSFunction_Newtonian()
231     for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j];  in RHSFunction_Newtonian()  local
238       for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j];  in RHSFunction_Newtonian()  local
243       for (int j = 1; j < 5; j++) divFdiff_i[j] = divFdiff[j - 1][i];  in RHSFunction_Newtonian()  local
250     for (CeedInt j = 0; j < 5; j++) {  in RHSFunction_Newtonian()  local
251 …(CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in RHSFunction_Newtonian()
298   for (CeedInt j = 0; j < 5; j++) v_i[j] += U_dot[j] - body_force[j];  in IFunction_Newtonian_Integrand()  local
337       for (int j = 1; j < 5; j++) divFdiff_i[j] = divFdiff[j - 1][i];  in IFunction_Newtonian()  local
345     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j];  in IFunction_Newtonian()  local
346     for (CeedInt j = 0; j < 5; j++) {  in IFunction_Newtonian()  local
348 …grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1] * dXdx[k][1] + grad_v_i[j]…  in IFunction_Newtonian()
409   for (CeedInt j = 0; j < 5; j++) {  in IJacobian_Newtonian_Integrand()  local
410     dU_dot[j] = context->ijacobian_time_shift * dU[j];  in IJacobian_Newtonian_Integrand()
411     v_i[j]    = dU_dot[j] - dbody_force[j];  in IJacobian_Newtonian_Integrand()
453       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * v_i[j];  in IJacobian_Newtonian()  local
456     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j];  in IJacobian_Newtonian()  local
457     for (int j = 0; j < 5; j++) {  in IJacobian_Newtonian()  local
458 …int k = 0; k < 3; k++) grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1] * d…  in IJacobian_Newtonian()
513     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in BoundaryIntegral()  local
558     for (int j = 0; j < 5; j++) dqi[j] = dq[j][i];  in BoundaryIntegral_Jacobian()  local
579     for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j];  in BoundaryIntegral_Jacobian()  local
628     for (CeedInt j = 1; j < 5; j++) {  // Continuity has no diffusive flux, therefore skip  in DivDiffusiveFluxVolumeRHS_NS()  local
630         Grad_v[k][j - 1][i] = -wdetJ * Dot3(dXdx[k], Fdiff[j]);  in DivDiffusiveFluxVolumeRHS_NS()
682     for (CeedInt j = 1; j < 5; j++) v[j - 1][i] = wdetJ * Fdiff[j];  in DivDiffusiveFluxBoundaryRHS_NS()  local
730     for (CeedInt j = 1; j < 5; j++) {  // Continuity has no diffusive flux, therefore skip  in DiffusiveFluxRHS_NS()  local
732         v[(j - 1) * 3 + k][i] = wdetJ * Fdiff[j][k];  in DiffusiveFluxRHS_NS()