280 …                          CeedScalar grad_v_i[5][3], CeedScalar kmstress[6], CeedScalar Tau_d[3]) {  in IFunction_Newtonian_Integrand()
292   AXPY(-1, (CeedScalar *)F_total, (CeedScalar *)grad_v_i, 15);  in IFunction_Newtonian_Integrand()
304   AXPY(1, (CeedScalar *)stab, (CeedScalar *)grad_v_i, 15);  in IFunction_Newtonian_Integrand()
339     CeedScalar v_i[5] = {0.}, grad_v_i[5][3] = {{0.}}, kmstress[6], Tau_d[3], sigma = 0;  in IFunction_Newtonian()  local
340 …n_Integrand(s, grad_s, s_dot, divFdiff_i, x_i, gas, context, dXdx, v_i, grad_v_i, kmstress, Tau_d);  in IFunction_Newtonian()
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()
390 …                         const CeedScalar Tau_d[3], CeedScalar v_i[5], CeedScalar grad_v_i[5][3]) {  in IJacobian_Newtonian_Integrand()
403   AXPY(-1, (CeedScalar *)dF_total, (CeedScalar *)grad_v_i, 15);  in IJacobian_Newtonian_Integrand()
418   AXPY(1, (CeedScalar *)dstab, (CeedScalar *)grad_v_i, 15);  in IJacobian_Newtonian_Integrand()
447     CeedScalar v_i[5] = {0.}, grad_v_i[5][3] = {{0.}};  in IJacobian_Newtonian()  local
448     IJacobian_Newtonian_Integrand(s, ds, grad_ds, gas, context, kmstress, Tau_d, v_i, grad_v_i);  in IJacobian_Newtonian()
458 …0; k < 3; k++) grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1] * dXdx[k][1…  in IJacobian_Newtonian()