43     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsNewtonianIG()  local
91     for (CeedInt j = 0; j < 5; j++) {  in MassFunction_Newtonian()  local
93 …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()
165     for (int j = 0; j < 5; j++) U[j] = q[j][i];  in RHSFunction_Newtonian()  local
185     for (CeedInt j = 0; j < 5; j++) {  in RHSFunction_Newtonian()  local
186 …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()
190     for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j];  in RHSFunction_Newtonian()  local
196       for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j];  in RHSFunction_Newtonian()  local
204     for (CeedInt j = 0; j < 5; j++) {  in RHSFunction_Newtonian()  local
205 …(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()
256     for (CeedInt j = 0; j < 5; j++) {  in IFunction_Newtonian()  local
258 …Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1] + dXdx[k][2] * Flux[…  in IFunction_Newtonian()
266     for (int j = 0; j < 5; j++) qi_dot[j] = q_dot[j][i];  in IFunction_Newtonian()  local
270     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * (U_dot[j] - body_force[j]);  in IFunction_Newtonian()  local
276       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * idl_residual[j];  in IFunction_Newtonian()  local
282     for (CeedInt j = 0; j < 5; j++) {  in IFunction_Newtonian()  local
284 …Grad_v[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[…  in IFunction_Newtonian()
331     for (int j = 0; j < 5; j++) dqi[j] = dq[j][i];  in IJacobian_Newtonian()  local
351     for (int j = 0; j < 5; j++) {  in IJacobian_Newtonian()  local
352 … k = 0; k < 3; k++) Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * dFlux[j][0] + dXdx[k][1] * dFlux[j][1…  in IJacobian_Newtonian()
358 …for (int j = 0; j < 5; j++) v[j][i] = wdetJ * (context->ijacobian_time_shift * dU[j] - dbody_force…  in IJacobian_Newtonian()  local
365       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * idl_residual[j];  in IJacobian_Newtonian()  local
370     for (CeedInt j = 0; j < 5; j++) U_dot[j] = context->ijacobian_time_shift * dU[j];  in IJacobian_Newtonian()  local
373     for (int j = 0; j < 5; j++) {  in IJacobian_Newtonian()  local
374 …or (int k = 0; k < 3; k++) Grad_v[k][j][i] += wdetJ * (dstab[j][0] * dXdx[k][0] + dstab[j][1] * dX…  in IJacobian_Newtonian()
428     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in BoundaryIntegral()  local
472     for (int j = 0; j < 5; j++) dqi[j] = dq[j][i];  in BoundaryIntegral_Jacobian()  local
493     for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j];  in BoundaryIntegral_Jacobian()  local