153 …                          const CeedScalar qn[], CeedScalar u[3], const StgShur14Context stg_ctx) {  in StgShur14Calc()
173   u[0] = ubar[0] + cij[0] * vp[0];  in StgShur14Calc()
174   u[1] = ubar[1] + cij[3] * vp[0] + cij[1] * vp[1];  in StgShur14Calc()
175   u[2] = ubar[2] + cij[4] * vp[0] + cij[5] * vp[1] + cij[2] * vp[2];  in StgShur14Calc()
195 …                   const CeedScalar eps, const CeedScalar lt, const CeedScalar nu, CeedScalar u[3],  in StgShur14Calc_PrecompEktot()
220   u[0] = ubar[0] + cij[0] * vp[0];  in StgShur14Calc_PrecompEktot()
221   u[1] = ubar[1] + cij[3] * vp[0] + cij[1] * vp[1];  in StgShur14Calc_PrecompEktot()
222   u[2] = ubar[2] + cij[4] * vp[0] + cij[5] * vp[1] + cij[2] * vp[2];  in StgShur14Calc_PrecompEktot()
291   CeedScalar                     qn[STG_NMODES_MAX], u[3], ubar[3], cij[6], eps, lt;  in ICsStg()  local
311       StgShur14Calc(x_i, time, ubar, cij, qn, u, stg_ctx);  in ICsStg()
313       for (CeedInt j = 0; j < 3; j++) u[j] = ubar[j];  in ICsStg()
316     CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5];  in ICsStg()
342   CeedScalar                  qn[STG_NMODES_MAX], u[3], ubar[3], cij[6], eps, lt;  in StgShur14Inflow()  local
370       StgShur14Calc(x, time, ubar, cij, qn, u, stg_ctx);  in StgShur14Inflow()
372       for (CeedInt j = 0; j < 3; j++) u[j] = ubar[j];  in StgShur14Inflow()
375     const CeedScalar E_kinetic = .5 * rho * Dot3(u, u);  in StgShur14Inflow()
390     const CeedScalar u_normal = Dot3(normal, u);  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()
406     const CeedScalar U[] = {rho, u[0], u[1], u[2], E}, kmstress[6] = {0.};  in StgShur14Inflow()
482   CeedScalar                  u[3], ubar[3], cij[6], eps, lt;  in StgShur14InflowStrongQF()  local
499 …14Calc_PrecompEktot(x, time, ubar, cij, inv_Ektotal[i], h_node_sep, x[1], eps, lt, nu, u, stg_ctx);  in StgShur14InflowStrongQF()
503         StgShur14Calc(x, time, ubar, cij, qn, u, stg_ctx);  in StgShur14InflowStrongQF()
506       for (CeedInt j = 0; j < 3; j++) u[j] = ubar[j];  in StgShur14InflowStrongQF()
509     CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5];  in StgShur14InflowStrongQF()