Lines Matching refs:q
64 …eneric(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, CeedScalar q[], void *ctx) { in Exact_AdvectionGeneric() argument
80 q[0] = 1.; in Exact_AdvectionGeneric()
81 q[1] = -(y - center[1]); in Exact_AdvectionGeneric()
82 q[2] = (x - center[0]); in Exact_AdvectionGeneric()
83 q[3] = 0; in Exact_AdvectionGeneric()
86 q[0] = 1.; in Exact_AdvectionGeneric()
87 q[1] = wind[0]; in Exact_AdvectionGeneric()
88 q[2] = wind[1]; in Exact_AdvectionGeneric()
89 q[3] = dim == 2 ? 0. : wind[2]; in Exact_AdvectionGeneric()
92 q[0] = 1.; in Exact_AdvectionGeneric()
93 q[1] = y / ly; in Exact_AdvectionGeneric()
94 q[2] = 0.; in Exact_AdvectionGeneric()
95 q[3] = 0.; in Exact_AdvectionGeneric()
106 q[4] = r <= rc ? (1. - r / rc) : 0.; in Exact_AdvectionGeneric()
110 q[4] = ((r <= rc) && (y < center[1])) ? (1. - r / rc) : 0.; in Exact_AdvectionGeneric()
114 … q[4] = ((r <= rc) && (y < center[1])) ? (1. - r / rc) * fmin(1.0, (center[1] - y) / 1.25) : 0.; in Exact_AdvectionGeneric()
117 q[4] = r <= rc ? .5 + .5 * cos(r * M_PI / rc) : 0; in Exact_AdvectionGeneric()
126 q[4] = r > half_width ? 0. : cos(2 * M_PI * r / half_width + M_PI) + 1.; in Exact_AdvectionGeneric()
136 q[4] = cross_product[2] > boundary_threshold ? 0 : 1; in Exact_AdvectionGeneric()
142 q[4] = 0; in Exact_AdvectionGeneric()
150 q[4] = sin(theta); in Exact_AdvectionGeneric()
153 q[4] = sin(theta) > 100 * CEED_EPSILON ? 1 : -1; in Exact_AdvectionGeneric()
162 q[4] = 1; // inflow and top boundary in Exact_AdvectionGeneric()
164 q[4] = 0; // lower wall in Exact_AdvectionGeneric()
167 if (y > bl_height) q[4] = 1; in Exact_AdvectionGeneric()
168 else q[4] = y / bl_height; in Exact_AdvectionGeneric()
184 CeedScalar q[5] = {0.}; in ICsAdvection() local
186 Exact_AdvectionGeneric(3, 0., x, 5, q, ctx); in ICsAdvection()
187 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection()
202 CeedScalar q[5] = {0.}; in ICsAdvection2d() local
204 Exact_AdvectionGeneric(2, context->time, x, 5, q, ctx); in ICsAdvection2d()
205 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection2d()
269 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; 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()
349 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; 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()
405 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; 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()
497 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; 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()