Lines Matching refs:CeedScalar
23 CeedScalar theta0;
24 CeedScalar thetaC;
25 CeedScalar P0;
26 CeedScalar N;
27 CeedScalar cv;
28 CeedScalar cp;
29 CeedScalar time;
30 CeedScalar mid_point;
31 CeedScalar P_high;
32 CeedScalar rho_high;
33 CeedScalar P_low;
34 CeedScalar rho_low;
39 CeedScalar Cyzb;
40 CeedScalar Byzb;
41 CeedScalar c_tau;
77 …N_HELPER CeedInt Exact_ShockTube(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, C… in Exact_ShockTube()
80 const CeedScalar mid_point = context->mid_point; // Midpoint of the domain in Exact_ShockTube()
81 const CeedScalar P_high = context->P_high; // Driver section pressure in Exact_ShockTube()
82 const CeedScalar rho_high = context->rho_high; // Driver section density in Exact_ShockTube()
83 const CeedScalar P_low = context->P_low; // Driven section pressure in Exact_ShockTube()
84 const CeedScalar rho_low = context->rho_low; // Driven section density in Exact_ShockTube()
87 const CeedScalar gamma = 1.4; // ratio of specific heats in Exact_ShockTube()
88 const CeedScalar x = X[0]; // Coordinates in Exact_ShockTube()
90 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube()
114 …void ConvectiveFluxJacobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u… in ConvectiveFluxJacobian_Euler()
115 const CeedScalar gamma) { in ConvectiveFluxJacobian_Euler()
116 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
139 CEED_QFUNCTION_HELPER CeedScalar Covariant_length_along_vector(CeedScalar vec[3], const CeedScalar … in Covariant_length_along_vector()
140 CeedScalar vec_norm = sqrt(vec[0] * vec[0] + vec[1] * vec[1] + vec[2] * vec[2]); in Covariant_length_along_vector()
141 CeedScalar vec_dot_jacobian[3] = {0.0}; in Covariant_length_along_vector()
147 CeedScalar norm_vec_dot_jacobian = in Covariant_length_along_vector()
149 CeedScalar h = 2.0 * vec_norm / norm_vec_dot_jacobian; in Covariant_length_along_vector()
168 …UNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const CeedScalar… in Tau_spatial()
169 const CeedScalar c_tau) { in Tau_spatial()
172 …CeedScalar h = 2 / sqrt(dXdx[0][i] * dXdx[0][i] + dXdx[1][i] * dXdx[1][i] + dXdx[2][i] * dXdx[2][i… in Tau_spatial()
174 CeedScalar fastest_wave = fabs(u[i]) + sound_speed; in Tau_spatial()
182 CEED_QFUNCTION(ICsShockTube)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *… in ICsShockTube()
183 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsShockTube()
184 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsShockTube()
187 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsShockTube()
188 CeedScalar q[5]; in ICsShockTube()
221 CEED_QFUNCTION(EulerShockTube)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const… in EulerShockTube()
222 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in EulerShockTube()
223 const CeedScalar(*dq)[5][CEED_Q_VLA] = (const CeedScalar(*)[5][CEED_Q_VLA])in[1]; in EulerShockTube()
224 const CeedScalar(*q_data) = in[2]; in EulerShockTube()
225 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in EulerShockTube()
226 CeedScalar(*dv)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1]; in EulerShockTube()
228 const CeedScalar gamma = 1.4; in EulerShockTube()
231 const CeedScalar Cyzb = context->Cyzb; in EulerShockTube()
232 const CeedScalar Byzb = context->Byzb; in EulerShockTube()
233 const CeedScalar c_tau = context->c_tau; in EulerShockTube()
238 const CeedScalar rho = q[0][i]; in EulerShockTube()
239 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in EulerShockTube()
240 const CeedScalar E = q[4][i]; in EulerShockTube()
241 const CeedScalar drho[3] = {dq[0][0][i], dq[1][0][i], dq[2][0][i]}; in EulerShockTube()
242 const CeedScalar dU[3][3] = { in EulerShockTube()
247 const CeedScalar dE[3] = {dq[0][4][i], dq[1][4][i], dq[2][4][i]}; in EulerShockTube()
248 CeedScalar wdetJ, dXdx[3][3]; in EulerShockTube()
251 CeedScalar du[3][3] = {{0}}; in EulerShockTube()
252 CeedScalar drhodx[3] = {0}; in EulerShockTube()
253 CeedScalar dEdx[3] = {0}; in EulerShockTube()
254 CeedScalar dUdx[3][3] = {{0}}; in EulerShockTube()
255 CeedScalar dXdxdXdxT[3][3] = {{0}}; in EulerShockTube()
268 …const CeedScalar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E… in EulerShockTube()
295 CeedScalar drho_norm = 0.0; // magnitude of the density gradient in EulerShockTube()
296 CeedScalar j_vec[3] = {0.0}; // unit vector aligned with the density gradient in EulerShockTube()
297 CeedScalar h_shock = 0.0; // element lengthscale in EulerShockTube()
298 CeedScalar acoustic_vel = 0.0; // characteristic velocity, acoustic speed in EulerShockTube()
299 CeedScalar tau_shock = 0.0; // timescale in EulerShockTube()
300 CeedScalar nu_shock = 0.0; // artificial diffusion in EulerShockTube()
328 CeedScalar jacob_F_conv[3][5][5] = {{{0.}}}; in EulerShockTube()
332 CeedScalar dqdx[5][3]; in EulerShockTube()
340 CeedScalar strong_conv[5] = {0}; in EulerShockTube()
349 const CeedScalar sound_speed = sqrt(gamma * P / rho); in EulerShockTube()
350 CeedScalar Tau_x[3] = {0.}; in EulerShockTube()
353 CeedScalar stab[5][3] = {0}; in EulerShockTube()