Lines Matching refs:i
22 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsNewtonianIG() local
26 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsNewtonianIG()
53 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in MassFunction_Newtonian() local
54 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MassFunction_Newtonian()
55 const CeedScalar qi_dot[5] = {q_dot[0][i], q_dot[1][i], q_dot[2][i], q_dot[3][i], q_dot[4][i]}; in MassFunction_Newtonian()
59 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MassFunction_Newtonian()
63 v[f][i] = wdetJ * qi_dot[f]; in MassFunction_Newtonian()
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()
95 for (int i = 0; i < 5; i++) damp_residual[i] += U[i]; in InternalDampingLayer_Residual() local
203 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in RHSFunction_Newtonian() local
205 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in RHSFunction_Newtonian()
206 for (int j = 0; j < 5; j++) U[j] = q[j][i]; in RHSFunction_Newtonian()
207 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in RHSFunction_Newtonian()
211 StatePhysicalGradientFromReference(Q, i, gas, s, STATEVAR_CONSERVATIVE, Grad_q, dXdx, grad_s); in RHSFunction_Newtonian()
227 …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()
231 for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j]; in RHSFunction_Newtonian()
238 for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j]; in RHSFunction_Newtonian()
243 for (int j = 1; j < 5; j++) divFdiff_i[j] = divFdiff[j - 1][i]; in RHSFunction_Newtonian()
251 …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()
324 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in IFunction_Newtonian() local
325 const CeedScalar q_i[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in IFunction_Newtonian()
326 const CeedScalar q_i_dot[5] = {q_dot[0][i], q_dot[1][i], q_dot[2][i], q_dot[3][i], q_dot[4][i]}; in IFunction_Newtonian()
327 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in IFunction_Newtonian()
332 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in IFunction_Newtonian()
334 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, grad_q, dXdx, grad_s); in IFunction_Newtonian()
337 for (int j = 1; j < 5; j++) divFdiff_i[j] = divFdiff[j - 1][i]; in IFunction_Newtonian()
345 for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j]; 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()
351 StoredValuesPack(Q, i, 0, 5, q_i, jac_data); in IFunction_Newtonian()
352 StoredValuesPack(Q, i, 5, 6, kmstress, jac_data); in IFunction_Newtonian()
353 StoredValuesPack(Q, i, 11, 3, Tau_d, jac_data); in IFunction_Newtonian()
354 if (context->idl_enable) StoredValuesPack(Q, i, 14, 1, &sigma, jac_data); in IFunction_Newtonian()
433 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in IJacobian_Newtonian() local
434 const CeedScalar dq_i[5] = {dq[0][i], dq[1][i], dq[2][i], dq[3][i], dq[4][i]}; in IJacobian_Newtonian()
436 StoredValuesUnpack(Q, i, 0, 5, jac_data, qi); in IJacobian_Newtonian()
437 StoredValuesUnpack(Q, i, 5, 6, jac_data, kmstress); in IJacobian_Newtonian()
438 StoredValuesUnpack(Q, i, 11, 3, jac_data, Tau_d); in IJacobian_Newtonian()
443 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in IJacobian_Newtonian()
445 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, grad_dq, dXdx, grad_ds); in IJacobian_Newtonian()
451 StoredValuesUnpack(Q, i, 14, 1, jac_data, &sigma); 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()
458 …for (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… in IJacobian_Newtonian()
490 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in BoundaryIntegral() local
491 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in BoundaryIntegral()
495 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in BoundaryIntegral()
499 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in BoundaryIntegral()
513 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in BoundaryIntegral()
516 StoredValuesPack(Q, i, 0, 5, qi, jac_data_sur); in BoundaryIntegral()
517 StoredValuesPack(Q, i, 5, 6, kmstress, jac_data_sur); in BoundaryIntegral()
550 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in BoundaryIntegral_Jacobian() local
552 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in BoundaryIntegral_Jacobian()
556 StoredValuesUnpack(Q, i, 0, 5, jac_data_sur, qi); in BoundaryIntegral_Jacobian()
557 StoredValuesUnpack(Q, i, 5, 6, jac_data_sur, kmstress); in BoundaryIntegral_Jacobian()
558 for (int j = 0; j < 5; j++) dqi[j] = dq[j][i]; in BoundaryIntegral_Jacobian()
564 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_dq, dXdx, grad_ds); in BoundaryIntegral_Jacobian()
579 for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in BoundaryIntegral_Jacobian()
608 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DivDiffusiveFluxVolumeRHS_NS() local
609 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in DivDiffusiveFluxVolumeRHS_NS()
614 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in DivDiffusiveFluxVolumeRHS_NS()
619 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in DivDiffusiveFluxVolumeRHS_NS()
630 Grad_v[k][j - 1][i] = -wdetJ * Dot3(dXdx[k], Fdiff[j]); in DivDiffusiveFluxVolumeRHS_NS()
661 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DivDiffusiveFluxBoundaryRHS_NS() local
662 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in DivDiffusiveFluxBoundaryRHS_NS()
667 QdataBoundaryGradientUnpack_3D(Q, i, q_data, &wdetJ, dXdx, normal); in DivDiffusiveFluxBoundaryRHS_NS()
672 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in DivDiffusiveFluxBoundaryRHS_NS()
682 for (CeedInt j = 1; j < 5; j++) v[j - 1][i] = wdetJ * Fdiff[j]; in DivDiffusiveFluxBoundaryRHS_NS()
710 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DiffusiveFluxRHS_NS() local
711 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in DiffusiveFluxRHS_NS()
716 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in DiffusiveFluxRHS_NS()
721 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in DiffusiveFluxRHS_NS()
732 v[(j - 1) * 3 + k][i] = wdetJ * Fdiff[j][k]; in DiffusiveFluxRHS_NS()