73   const SetupContextAdv context = (SetupContextAdv)ctx;  in Exact_AdvectionGeneric()  local
74   const CeedScalar      rc      = context->rc;  in Exact_AdvectionGeneric()
75   const CeedScalar      lx      = context->lx;  in Exact_AdvectionGeneric()
76   const CeedScalar      ly      = context->ly;  in Exact_AdvectionGeneric()
77   const CeedScalar      lz      = dim == 2 ? 0. : context->lz;  in Exact_AdvectionGeneric()
78   const CeedScalar     *wind    = context->wind;  in Exact_AdvectionGeneric()
87   switch (context->initial_condition_type) {  in Exact_AdvectionGeneric()
99   switch (context->wind_type) {  in Exact_AdvectionGeneric()
116   switch (context->initial_condition_type) {  in Exact_AdvectionGeneric()
119       switch (context->bubble_continuity_type) {  in Exact_AdvectionGeneric()
138       CeedScalar half_width = context->lx / 2;  in Exact_AdvectionGeneric()
143       CeedScalar       inflow_to_point[3] = {x - context->lx / 2, y, 0};  in Exact_AdvectionGeneric()
151 …(x > context->lx - boundary_threshold && wind[0] > boundary_threshold) ||  // outflow at +x bounda…  in Exact_AdvectionGeneric()
152 …(y > context->ly - boundary_threshold && wind[1] > boundary_threshold)     // outflow at +y bounda…  in Exact_AdvectionGeneric()
184   const SetupContextAdv context    = (SetupContextAdv)ctx;  in ICsAdvection2d()  local
190     Exact_AdvectionGeneric(2, context->time, x, 5, q, ctx);  in ICsAdvection2d()
257 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd…  in Tau()  argument
258   switch (context->stabilization_tau) {  in Tau()
263       return context->CtauS / sqrt(DotN(uX, uX, dim));  in Tau()
270 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * cont…  in Tau()
290   AdvectionContext                 context   = (AdvectionContext)ctx;  in IFunction_AdvectionGeneric()  local
326     if (context->strong_form) {  // Strong Galerkin convection term: v div(E u)  in IFunction_AdvectionGeneric()
335       for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total;  in IFunction_AdvectionGeneric()
340     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in IFunction_AdvectionGeneric()
341     for (CeedInt j = 0; j < dim; j++) switch (context->stabilization) {  in IFunction_AdvectionGeneric()
373   AdvectionContext                 context    = (AdvectionContext)ctx;  in MassFunction_AdvectionGeneric()  local
394     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in MassFunction_AdvectionGeneric()
395     for (CeedInt j = 0; j < dim; j++) switch (context->stabilization) {  in MassFunction_AdvectionGeneric()
428   AdvectionContext                 context    = (AdvectionContext)ctx;  in RHSFunction_AdvectionGeneric()  local
459     if (context->strong_form) {  // Strong Galerkin convection term: v div(E u)  in RHSFunction_AdvectionGeneric()
470       for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total;  in RHSFunction_AdvectionGeneric()
475     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in RHSFunction_AdvectionGeneric()
476     for (CeedInt j = 0; j < dim; j++) switch (context->stabilization) {  in RHSFunction_AdvectionGeneric()
516   AdvectionContext context     = (AdvectionContext)ctx;  in Advection_InOutFlowGeneric()  local
517   const CeedScalar E_wind      = context->E_wind;  in Advection_InOutFlowGeneric()
518   const CeedScalar strong_form = context->strong_form;  in Advection_InOutFlowGeneric()
519   const bool       is_implicit = context->implicit;  in Advection_InOutFlowGeneric()