116 …                                                            CeedScalar v_i[5], CeedScalar *sigma) {  in InternalDampingLayer_IFunction_Integrand()
120   AXPY(1, idl_residual, v_i, 5);  in InternalDampingLayer_IFunction_Integrand()
136                                                                     CeedScalar v_i[5]) {  in InternalDampingLayer_IJacobian_Integrand()
139   AXPY(1, idl_residual, v_i, 5);  in InternalDampingLayer_IJacobian_Integrand()
279 …            const NewtonianIdealGasContext context, const CeedScalar dXdx[3][3], CeedScalar v_i[5],  in IFunction_Newtonian_Integrand()
298   for (CeedInt j = 0; j < 5; j++) v_i[j] += U_dot[j] - body_force[j];  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 …IFunction_Newtonian_Integrand(s, grad_s, s_dot, divFdiff_i, x_i, gas, context, dXdx, v_i, grad_v_i…  in IFunction_Newtonian()
343                                                v_i, &sigma);  in IFunction_Newtonian()
345     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * 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()
411     v_i[j]    = dU_dot[j] - dbody_force[j];  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()
452       InternalDampingLayer_IJacobian_Integrand(s, ds, gas, sigma, v_i);  in IJacobian_Newtonian()
453       for (int j = 0; j < 5; j++) v[j][i] += wdetJ * v_i[j];  in IJacobian_Newtonian()
456     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j];  in IJacobian_Newtonian()