42   for (int i = 0; i < 4; i++) {  in ChebyshevEval()  local
44     f[i] = table[i][0] * Tf[0] + table[i][1] * Tf[1] + table[i][2] * Tf[2];  in ChebyshevEval()
46   for (int i = 3; i < N; i++) {  in ChebyshevEval()  local
48     table[0][i % 3] = 2 * x * table[0][(i - 1) % 3] - table[0][(i - 2) % 3];  in ChebyshevEval()
52       table[j][i % 3] = i * (2 * table[j - 1][(i - 1) % 3] + table[j][(i - 2) % 3] / (i - 2));  in ChebyshevEval()
55       f[j] += table[j][i % 3] * Tf[i];  in ChebyshevEval()
58   for (int i = 1; i < 4; i++) {  in ChebyshevEval()  local
60     for (int j = 0; j < i; j++) f[i] *= dX_deta;  in ChebyshevEval()
112   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ICsBlasius()  local
113     const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]};  in ICsBlasius()
118     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsBlasius()
138   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Blasius_Inflow()  local
140     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, normal);  in Blasius_Inflow()
144     const CeedScalar x[3] = {X[0][i], X[1][i], 0.};  in Blasius_Inflow()
148     for (CeedInt j = 0; j < 5; j++) qi[j] = q[j][i];  in Blasius_Inflow()
171     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in Blasius_Inflow()
193   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Blasius_Inflow_Jacobian()  local
195     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, normal);  in Blasius_Inflow_Jacobian()
199     const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]};  in Blasius_Inflow_Jacobian()
207       drho                   = dq[0][i];  in Blasius_Inflow_Jacobian()
215       dE   = dq[4][i];  in Blasius_Inflow_Jacobian()
221     v[0][i] = -wdetJb * drho * u_normal;  in Blasius_Inflow_Jacobian()
223       v[j + 1][i] = -wdetJb * (drho * u_normal * s.Y.velocity[j] + normal[j] * dP);  in Blasius_Inflow_Jacobian()
225     v[4][i] = -wdetJb * u_normal * (dE + dP);  in Blasius_Inflow_Jacobian()