75     CeedScalar       wdetJ, dXdx[3][3];  in MassFunction_Newtonian()  local
76     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in MassFunction_Newtonian()
80       v[f][i] = wdetJ * qi_dot[f];  in MassFunction_Newtonian()
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()
163     CeedScalar       U[5], wdetJ, dXdx[3][3];  in RHSFunction_Newtonian()  local
166     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in RHSFunction_Newtonian()
186 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * …  in RHSFunction_Newtonian()
190     for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j];  in RHSFunction_Newtonian()
196       for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j];  in RHSFunction_Newtonian()
205 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] *…  in RHSFunction_Newtonian()
238     CeedScalar wdetJ, dXdx[3][3];  in IFunction_Newtonian()  local
239     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in IFunction_Newtonian()
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()
270     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * (U_dot[j] - body_force[j]);  in IFunction_Newtonian()
276       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * idl_residual[j];  in IFunction_Newtonian()
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()
321     CeedScalar wdetJ, dXdx[3][3];  in IJacobian_Newtonian()  local
322     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in IJacobian_Newtonian()
352 …for (int k = 0; k < 3; k++) Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * dFlux[j][0] + dXdx[k][1] * dF…  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()
365       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * idl_residual[j];  in IJacobian_Newtonian()
374 …for (int k = 0; k < 3; k++) Grad_v[k][j][i] += wdetJ * (dstab[j][0] * dXdx[k][0] + dstab[j][1] * d…  in IJacobian_Newtonian()