25   CeedScalar                       theta0;
26   CeedScalar                       thetaC;
27   CeedScalar                       P0;
28   CeedScalar                       N;
29   CeedScalar                       rc;
30   CeedScalar                       center[3];
31   CeedScalar                       dc_axis[3];
86 CEED_QFUNCTION_HELPER State Exact_DC(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf…  in Exact_DC()
89   const CeedScalar            theta0  = context->theta0;  in Exact_DC()
90   const CeedScalar            thetaC  = context->thetaC;  in Exact_DC()
91   const CeedScalar            P0      = context->P0;  in Exact_DC()
92   const CeedScalar            N       = context->N;  in Exact_DC()
93   const CeedScalar            rc      = context->rc;  in Exact_DC()
94   const CeedScalar           *center  = context->center;  in Exact_DC()
95   const CeedScalar           *dc_axis = context->dc_axis;  in Exact_DC()
97   const CeedScalar            cp      = gas->cp;  in Exact_DC()
98   const CeedScalar            cv      = gas->cv;  in Exact_DC()
99   const CeedScalar            Rd      = cp - cv;  in Exact_DC()
100   const CeedScalar           *g_vec   = gas->g;  in Exact_DC()
101   const CeedScalar            g       = -g_vec[2];  in Exact_DC()
105   const CeedScalar x = X[0];  in Exact_DC()
106   const CeedScalar y = X[1];  in Exact_DC()
107   const CeedScalar z = X[2];  in Exact_DC()
110   CeedScalar rr[3] = {x - center[0], y - center[1], z - center[2]};  in Exact_DC()
113   const CeedScalar r           = sqrt(Dot3(rr, rr));  in Exact_DC()
114   const CeedScalar delta_theta = r <= rc ? thetaC * (1. + cos(M_PI * r / rc)) / 2. : 0.;  in Exact_DC()
115   const CeedScalar theta       = theta0 * exp(Square(N) * z / g) + delta_theta;  in Exact_DC()
118   const CeedScalar Pi = 1. + Square(g) * (exp(-Square(N) * z / g) - 1.) / (cp * theta0 * Square(N));  in Exact_DC()
121   CeedScalar Y[5] = {0.};  in Exact_DC()
134 CEED_QFUNCTION(ICsDC)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) {  in ICsDC()
135   const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0];  in ICsDC()
136   CeedScalar(*q0)[CEED_Q_VLA]      = (CeedScalar(*)[CEED_Q_VLA])out[0];  in ICsDC()
142     const CeedScalar x[]  = {X[0][i], X[1][i], X[2][i]};  in ICsDC()
144     CeedScalar       q[5] = {0};  in ICsDC()