46   for (CeedInt i = 0; i < nprofs; i++) {  in InterpolateProfile()  local
47     if (wall_dist < prof_wd[i]) {  in InterpolateProfile()
48       idx = i;  in InterpolateProfile()
165     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc()  local
171   for (CeedInt i = 0; i < 3; i++) vp[i] *= 2 * sqrt(1.5);  in StgShur14Calc()  local
210     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc_PrecompEktot()  local
218   for (CeedInt i = 0; i < 3; i++) vp[i] *= 2 * sqrt(1.5);  in StgShur14Calc_PrecompEktot()  local
261   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in StgShur14Preprocess()  local
262     const CeedScalar wall_dist = x[1][i];  in StgShur14Preprocess()
264     StoredValuesUnpack(Q, i, 0, 9, dXdx_q, (CeedScalar *)dXdx);  in StgShur14Preprocess()
277     stg_data[i] = Ek_tot != 0 ? 1 / Ek_tot : 0;  in StgShur14Preprocess()
299   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ICsStg()  local
300     const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]};  in ICsStg()
304     InvertMappingJacobian_3D(Q, i, J, dXdx, NULL);  in ICsStg()
320       q0[j][i] = q[j];  in ICsStg()
355   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in StgShur14Inflow()  local
356     const CeedScalar rho = prescribe_T ? q[0][i] : P0 / (Rd * theta0);  in StgShur14Inflow()
357     const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]};  in StgShur14Inflow()
359     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal);  in StgShur14Inflow()
367     InterpolateProfile(X[1][i], ubar, cij, &eps, &lt, stg_ctx);  in StgShur14Inflow()
369       CalcSpectrum(X[1][i], eps, lt, h_node_sep, mu / rho, qn, stg_ctx);  in StgShur14Inflow()
383       E_internal = q[4][i] - E_kinetic;  // uses prescribed rho and u, E from solution  in StgShur14Inflow()
394     for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.;  in StgShur14Inflow()
398     v[0][i] -= wdetJb * rho * u_normal;  in StgShur14Inflow()
401     for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P);  in StgShur14Inflow()
404     v[4][i] -= wdetJb * u_normal * (E + P);  in StgShur14Inflow()
407     StoredValuesPack(Q, i, 0, 5, U, jac_data_sur);  in StgShur14Inflow()
408     StoredValuesPack(Q, i, 5, 6, kmstress, jac_data_sur);  in StgShur14Inflow()
429   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in StgShur14Inflow_Jacobian()  local
435     const CeedScalar wdetJb = (implicit ? -1. : 1.) * q_data_sur[0][i];  in StgShur14Inflow_Jacobian()
439     for (CeedInt j = 0; j < 3; j++) velocity[j] = jac_data_sur[5 + j][i];  in StgShur14Inflow_Jacobian()
446       drho                   = dq[0][i];  in StgShur14Inflow_Jacobian()
453       dE   = dq[4][i];  in StgShur14Inflow_Jacobian()
456     const CeedScalar normal[3] = {q_data_sur[1][i], q_data_sur[2][i], q_data_sur[3][i]};  in StgShur14Inflow_Jacobian()
460     v[0][i] = -wdetJb * drho * u_normal;  in StgShur14Inflow_Jacobian()
461 …for (int j = 0; j < 3; j++) v[j + 1][i] = -wdetJb * (drho * u_normal * velocity[j] + normal[j] * d…  in StgShur14Inflow_Jacobian()
462     v[4][i] = -wdetJb * u_normal * (dE + dP);  in StgShur14Inflow_Jacobian()
490   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in StgShur14InflowStrongQF()  local
491     const CeedScalar x[] = {coords[0][i], coords[1][i], coords[2][i]};  in StgShur14InflowStrongQF()
493     StoredValuesUnpack(Q, i, 0, 9, dXdx_q, (CeedScalar *)dXdx);  in StgShur14InflowStrongQF()
496     InterpolateProfile(coords[1][i], ubar, cij, &eps, &lt, stg_ctx);  in StgShur14InflowStrongQF()
499 …StgShur14Calc_PrecompEktot(x, time, ubar, cij, inv_Ektotal[i], h_node_sep, x[1], eps, lt, nu, u, s…  in StgShur14InflowStrongQF()
502         CalcSpectrum(coords[1][i], eps, lt, h_node_sep, nu, qn, stg_ctx);  in StgShur14InflowStrongQF()
524       bcval[j][i] = scale[i] * q[j];  in StgShur14InflowStrongQF()