Lines Matching refs:i
182 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsAdvection() local
183 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsAdvection()
187 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection()
200 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsAdvection2d() local
201 const CeedScalar x[] = {X[0][i], X[1][i]}; in ICsAdvection2d()
205 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection2d()
210 CEED_QFUNCTION_HELPER void StatePhysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, N… in StatePhysicalGradientFromReference_ND() argument
218 …dqi[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXd… in StatePhysicalGradientFromReference_ND()
230 …dqi[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXd… in StatePhysicalGradientFromReference_ND()
231 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
280 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in IFunction_AdvectionGeneric() local
281 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in IFunction_AdvectionGeneric()
285 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in IFunction_AdvectionGeneric()
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
292 for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0; // No Change in density or momentum in IFunction_AdvectionGeneric()
293 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in IFunction_AdvectionGeneric()
306 v[4][i] = wdetJ * q_dot[4][i]; // transient part (ALWAYS) in IFunction_AdvectionGeneric()
308 v[4][i] += wdetJ * strong_conv; in IFunction_AdvectionGeneric()
310 for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = -wdetJ * s.U.E_total * uX[j]; in IFunction_AdvectionGeneric()
316 for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total; in IFunction_AdvectionGeneric() local
318 for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] -= wdetJ * Fe_dXdx[k]; in IFunction_AdvectionGeneric()
327 grad_v[j][4][i] += wdetJ * TauS * uX[j] * strong_conv; in IFunction_AdvectionGeneric()
330 … CeedScalar divFdiff_i = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? divFdiff[i] : 0.; in IFunction_AdvectionGeneric()
331 grad_v[j][4][i] += wdetJ * TauS * uX[j] * (q_dot[4][i] + strong_conv + divFdiff_i); in IFunction_AdvectionGeneric()
358 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in MassFunction_AdvectionGeneric() local
359 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MassFunction_AdvectionGeneric()
362 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in MassFunction_AdvectionGeneric()
365 for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0; // No Change in density or momentum in MassFunction_AdvectionGeneric()
366 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in MassFunction_AdvectionGeneric()
370 v[4][i] = wdetJ * q_dot[4][i]; in MassFunction_AdvectionGeneric()
380 grad_v[j][4][i] = 0; in MassFunction_AdvectionGeneric()
383 grad_v[j][4][i] = wdetJ * TauS * q_dot[4][i] * uX[j]; in MassFunction_AdvectionGeneric()
415 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in RHSFunction_AdvectionGeneric() local
416 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in RHSFunction_AdvectionGeneric()
420 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in RHSFunction_AdvectionGeneric()
422 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in RHSFunction_AdvectionGeneric()
427 for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0; // No Change in density or momentum in RHSFunction_AdvectionGeneric()
428 v[f][i] = 0.; in RHSFunction_AdvectionGeneric()
443 v[4][i] = -wdetJ * strong_conv; in RHSFunction_AdvectionGeneric()
444 for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = 0; in RHSFunction_AdvectionGeneric()
446 for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = wdetJ * s.U.E_total * uX[j]; in RHSFunction_AdvectionGeneric()
447 v[4][i] = 0.; in RHSFunction_AdvectionGeneric()
453 for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total; in RHSFunction_AdvectionGeneric() local
455 for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] += wdetJ * Fe_dXdx[k]; in RHSFunction_AdvectionGeneric()
465 … CeedScalar divFdiff_i = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? divFdiff[i] : 0.; in RHSFunction_AdvectionGeneric()
466 grad_v[j][4][i] -= wdetJ * TauS * (strong_conv + divFdiff_i) * uX[j]; in RHSFunction_AdvectionGeneric()
506 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Advection_InOutFlowGeneric() local
507 const CeedScalar rho = q[0][i]; in Advection_InOutFlowGeneric()
508 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Advection_InOutFlowGeneric()
509 const CeedScalar E = q[4][i]; in Advection_InOutFlowGeneric()
512 QdataBoundaryUnpack_ND(dim, Q, i, q_data_sur, &wdetJb, NULL, normal); in Advection_InOutFlowGeneric()
519 v[j][i] = 0; in Advection_InOutFlowGeneric()
523 v[4][i] = -(1 - strong_form) * wdetJb * E * u_normal; in Advection_InOutFlowGeneric()
525 v[4][i] = -(1 - strong_form) * wdetJb * E_wind * u_normal; in Advection_InOutFlowGeneric()
548 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DivDiffusiveFluxVolumeRHS_AdvDif_Generic() local
551 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
555 GradUnpackND(Q, i, 5, dim, Grad_q, Grad_qn); in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
563 for (CeedInt k = 0; k < dim; k++) Grad_v[k][i] = -wdetJ * F_diff_dXdx[k]; in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
585 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic() local
588 QdataBoundaryGradientUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx, normal); in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
592 GradUnpackND(Q, i, 5, dim, Grad_q, Grad_qn); in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
598 v[i] = wdetJ * DotN(F_diff, normal, dim); in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
620 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DiffusiveFluxRHS_AdvDif_Generic() local
623 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in DiffusiveFluxRHS_AdvDif_Generic()
627 GradUnpackND(Q, i, 5, dim, Grad_q, Grad_qn); in DiffusiveFluxRHS_AdvDif_Generic()
632 for (CeedInt k = 0; k < dim; k++) v[k][i] = wdetJ * F_diff[k]; in DiffusiveFluxRHS_AdvDif_Generic()