65   const SetupContextAdv context = (SetupContextAdv)ctx;  in Exact_AdvectionGeneric()  local
66   const CeedScalar      rc      = context->rc;  in Exact_AdvectionGeneric()
67   const CeedScalar      lx      = context->lx;  in Exact_AdvectionGeneric()
68   const CeedScalar      ly      = context->ly;  in Exact_AdvectionGeneric()
69   const CeedScalar      lz      = dim == 2 ? 0. : context->lz;  in Exact_AdvectionGeneric()
70   const CeedScalar     *wind    = context->wind;  in Exact_AdvectionGeneric()
78   switch (context->wind_type) {  in Exact_AdvectionGeneric()
99   switch (context->initial_condition_type) {  in Exact_AdvectionGeneric()
103       switch (context->bubble_continuity_type) {  in Exact_AdvectionGeneric()
125       CeedScalar half_width = context->lx / 2;  in Exact_AdvectionGeneric()
131       CeedScalar       inflow_to_point[3] = {x - context->lx / 2, y, 0};  in Exact_AdvectionGeneric()
139 …(x > context->lx - boundary_threshold && wind[0] > boundary_threshold) ||  // outflow at +x bounda…  in Exact_AdvectionGeneric()
140 …(y > context->ly - boundary_threshold && wind[1] > boundary_threshold)     // outflow at +y bounda…  in Exact_AdvectionGeneric()
147       CeedScalar theta = context->wave_frequency * DotN(X, wind, dim) + context->wave_phase;  in Exact_AdvectionGeneric()
148       switch (context->wave_type) {  in Exact_AdvectionGeneric()
166         CeedScalar bl_height = ly * context->bl_height_factor;  in Exact_AdvectionGeneric()
198   const SetupContextAdv context    = (SetupContextAdv)ctx;  in ICsAdvection2d()  local
204     Exact_AdvectionGeneric(2, context->time, x, 5, q, ctx);  in ICsAdvection2d()
240 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd…  in Tau()  argument
241   switch (context->stabilization_tau) {  in Tau()
246       return context->CtauS / sqrt(DotN(uX, uX, dim));  in Tau()
255 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * Squa…  in Tau()
256 … Square(context->diffusion_coeff) * DotN(gijd_mat, gijd_mat, dim * dim) * Square(context->Ctau_d));  in Tau()
267   AdvectionContext context = (AdvectionContext)ctx;  in IFunction_AdvectionGeneric()  local
273 …const CeedScalar(*divFdiff)          = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? in[5]…  in IFunction_AdvectionGeneric()
307     if (context->strong_form) {     // Strong Galerkin convection term: v div(E u)  in IFunction_AdvectionGeneric()
316       for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total;  in IFunction_AdvectionGeneric()
321     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in IFunction_AdvectionGeneric()
323       switch (context->stabilization) {  in IFunction_AdvectionGeneric()
330 …    CeedScalar divFdiff_i = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? divFdiff[i] : 0.;  in IFunction_AdvectionGeneric()
355   AdvectionContext      context = (AdvectionContext)ctx;  in MassFunction_AdvectionGeneric()  local
375     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in MassFunction_AdvectionGeneric()
377       switch (context->stabilization) {  in MassFunction_AdvectionGeneric()
403   AdvectionContext context = (AdvectionContext)ctx;  in RHSFunction_AdvectionGeneric()  local
408 …const CeedScalar(*divFdiff)      = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? in[4] : N…  in RHSFunction_AdvectionGeneric()
442     if (context->strong_form) {  // Strong Galerkin convection term: v div(E u)  in RHSFunction_AdvectionGeneric()
453       for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total;  in RHSFunction_AdvectionGeneric()
458     const CeedScalar TauS = Tau(context, s, dXdx, dim);  in RHSFunction_AdvectionGeneric()
460       switch (context->stabilization) {  in RHSFunction_AdvectionGeneric()
465 …    CeedScalar divFdiff_i = context->divFdiff_method != DIV_DIFF_FLUX_PROJ_NONE ? divFdiff[i] : 0.;  in RHSFunction_AdvectionGeneric()
501   AdvectionContext context     = (AdvectionContext)ctx;  in Advection_InOutFlowGeneric()  local
502   const CeedScalar E_wind      = context->E_wind;  in Advection_InOutFlowGeneric()
503   const CeedScalar strong_form = context->strong_form;  in Advection_InOutFlowGeneric()
504   const bool       is_implicit = context->implicit;  in Advection_InOutFlowGeneric()
546   AdvectionContext context = (AdvectionContext)ctx;  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()  local
558       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
583   AdvectionContext context = (AdvectionContext)ctx;  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()  local
595       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
618   AdvectionContext context = (AdvectionContext)ctx;  in DiffusiveFluxRHS_AdvDif_Generic()  local
630       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DiffusiveFluxRHS_AdvDif_Generic()