1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 // SPDX-License-Identifier: BSD-2-Clause
28 CEED_QFUNCTION_HELPER void StabilizationMatrix(NewtonianIdealGasContext gas, State s, CeedScalar Ta…  in StabilizationMatrix()  argument
37   State ds = StateFromY_fwd(gas, s, dY);  in StabilizationMatrix()
38   FluxInviscid_fwd(gas, s, ds, dF);  in StabilizationMatrix()
46 CEED_QFUNCTION_HELPER void Stabilization(NewtonianIdealGasContext gas, State s, CeedScalar Tau_d[3]…  in Stabilization()  argument
48   // -- Stabilization method: none (Galerkin), SU, or SUPG  in Stabilization()
50   switch (gas->stabilization) {  in Stabilization()
54       FluxInviscidStrong(gas, s, ds, strong_residual);  in Stabilization()
57       FluxInviscidStrong(gas, s, ds, strong_residual);  in Stabilization()
58       for (CeedInt j = 0; j < 5; j++) strong_residual[j] += U_dot[j] - body_force[j];  in Stabilization()
61   StabilizationMatrix(gas, s, Tau_d, strong_residual, stab);  in Stabilization()
69 //   Tau[i] = itau=0 which is diagonal-Shakib (3 values still but not spatial)
71 CEED_QFUNCTION_HELPER void Tau_diagPrim(NewtonianIdealGasContext gas, State s, const CeedScalar dXd…  in Tau_diagPrim()  argument
74   const CeedScalar Ctau_t = gas->Ctau_t;  in Tau_diagPrim()
75   const CeedScalar Ctau_v = gas->Ctau_v;  in Tau_diagPrim()
76   const CeedScalar Ctau_C = gas->Ctau_C;  in Tau_diagPrim()
77   const CeedScalar Ctau_M = gas->Ctau_M;  in Tau_diagPrim()
78   const CeedScalar Ctau_E = gas->Ctau_E;  in Tau_diagPrim()
79   const CeedScalar cv     = gas->cv;  in Tau_diagPrim()
80   const CeedScalar mu     = gas->mu;  in Tau_diagPrim()
81   const CeedScalar rho    = s.U.density;  in Tau_diagPrim()
94     MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj);  in Tau_diagPrim()
95     velocity_term = Dot3(s.Y.velocity, gij_uj);  in Tau_diagPrim()