Lines Matching refs:CeedScalar
41 …R void InterpolateProfile(const CeedScalar wall_dist, CeedScalar ubar[3], CeedScalar cij[6], CeedS… in InterpolateProfile()
44 const CeedScalar *prof_wd = &stg_ctx->data[stg_ctx->offsets.wall_dist]; in InterpolateProfile()
45 const CeedScalar *prof_eps = &stg_ctx->data[stg_ctx->offsets.eps]; in InterpolateProfile()
46 const CeedScalar *prof_lt = &stg_ctx->data[stg_ctx->offsets.lt]; in InterpolateProfile()
47 const CeedScalar *prof_ubar = &stg_ctx->data[stg_ctx->offsets.ubar]; in InterpolateProfile()
48 const CeedScalar *prof_cij = &stg_ctx->data[stg_ctx->offsets.cij]; in InterpolateProfile()
59 CeedScalar coeff = (wall_dist - prof_wd[idx - 1]) / (prof_wd[idx] - prof_wd[idx - 1]); in InterpolateProfile()
100 CEED_QFUNCTION_HELPER CeedScalar Calc_qn(const CeedScalar kappa, const CeedScalar dkappa, const Cee… in Calc_qn()
101 const CeedScalar ke, const CeedScalar Ektot_inv) { in Calc_qn()
102 …const CeedScalar feta_x_fcut = exp(-Square(12 * kappa / keta) - Cube(4 * Max(kappa - 0.9 * kcut, 0… in Calc_qn()
107 …HELPER void SpectrumConstants(const CeedScalar wall_dist, const CeedScalar eps, const CeedScalar l… in SpectrumConstants()
108 … const CeedScalar nu, CeedScalar *hmax, CeedScalar *ke, CeedScalar *keta, CeedScalar *kcut) { in SpectrumConstants()
128 …TION_HELPER void CalcSpectrum(const CeedScalar wall_dist, const CeedScalar eps, const CeedScalar l… in CalcSpectrum()
129 … const CeedScalar nu, CeedScalar qn[], const StgShur14Context stg_ctx) { in CalcSpectrum()
131 const CeedScalar *kappa = &stg_ctx->data[stg_ctx->offsets.kappa]; in CalcSpectrum()
132 CeedScalar hmax, ke, keta, kcut, Ektot = 0.0; in CalcSpectrum()
137 const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1]; in CalcSpectrum()
157 …TION_HELPER void StgShur14Calc(const CeedScalar X[3], const CeedScalar t, const CeedScalar ubar[3]… in StgShur14Calc()
158 … const CeedScalar qn[], CeedScalar u[3], const StgShur14Context stg_ctx) { in StgShur14Calc()
160 const CeedScalar *kappa = &stg_ctx->data[stg_ctx->offsets.kappa]; in StgShur14Calc()
161 const CeedScalar *phi = &stg_ctx->data[stg_ctx->offsets.phi]; in StgShur14Calc()
162 const CeedScalar *sigma = &stg_ctx->data[stg_ctx->offsets.sigma]; in StgShur14Calc()
163 const CeedScalar *d = &stg_ctx->data[stg_ctx->offsets.d]; in StgShur14Calc()
164 CeedScalar xdotd, vp[3] = {0.}; in StgShur14Calc()
165 CeedScalar xhat[] = {0., X[1], X[2]}; in StgShur14Calc()
171 const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]); in StgShur14Calc()
198 …oid StgShur14Calc_PrecompEktot(const CeedScalar X[3], const CeedScalar t, const CeedScalar ubar[3]… in StgShur14Calc_PrecompEktot()
199 … const CeedScalar Ektot, const CeedScalar h_node_sep[3], const CeedScalar wall_dist, in StgShur14Calc_PrecompEktot()
200 … const CeedScalar eps, const CeedScalar lt, const CeedScalar nu, CeedScalar u[3], in StgShur14Calc_PrecompEktot()
203 const CeedScalar *kappa = &stg_ctx->data[stg_ctx->offsets.kappa]; in StgShur14Calc_PrecompEktot()
204 const CeedScalar *phi = &stg_ctx->data[stg_ctx->offsets.phi]; in StgShur14Calc_PrecompEktot()
205 const CeedScalar *sigma = &stg_ctx->data[stg_ctx->offsets.sigma]; in StgShur14Calc_PrecompEktot()
206 const CeedScalar *d = &stg_ctx->data[stg_ctx->offsets.d]; in StgShur14Calc_PrecompEktot()
207 CeedScalar hmax, ke, keta, kcut; in StgShur14Calc_PrecompEktot()
209 CeedScalar xdotd, vp[3] = {0.}; in StgShur14Calc_PrecompEktot()
210 CeedScalar xhat[] = {0., X[1], X[2]}; in StgShur14Calc_PrecompEktot()
216 const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]); in StgShur14Calc_PrecompEktot()
217 const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1]; in StgShur14Calc_PrecompEktot()
218 const CeedScalar qn = Calc_qn(kappa[n], dkappa, keta, kcut, ke, Ektot); in StgShur14Calc_PrecompEktot()
233 CEED_QFUNCTION(StgShur14Preprocess)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in StgShur14Preprocess()
234 const CeedScalar(*dXdx_q)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0]; in StgShur14Preprocess()
235 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in StgShur14Preprocess()
237 CeedScalar(*stg_data) = (CeedScalar(*))out[0]; in StgShur14Preprocess()
239 CeedScalar ubar[3], cij[6], eps, lt; in StgShur14Preprocess()
241 const CeedScalar dx = stg_ctx->dx; in StgShur14Preprocess()
242 const CeedScalar mu = stg_ctx->newtonian_ctx.mu; in StgShur14Preprocess()
243 const CeedScalar theta0 = stg_ctx->theta0; in StgShur14Preprocess()
244 const CeedScalar P0 = stg_ctx->P0; in StgShur14Preprocess()
245 const CeedScalar Rd = GasConstant(&stg_ctx->newtonian_ctx); in StgShur14Preprocess()
246 const CeedScalar rho = P0 / (Rd * theta0); in StgShur14Preprocess()
247 const CeedScalar nu = mu / rho; in StgShur14Preprocess()
250 const CeedScalar *kappa = &stg_ctx->data[stg_ctx->offsets.kappa]; in StgShur14Preprocess()
251 CeedScalar hmax, ke, keta, kcut; in StgShur14Preprocess()
254 const CeedScalar wall_dist = x[1][i]; in StgShur14Preprocess()
255 const CeedScalar dXdx[2][3] = { in StgShur14Preprocess()
260 CeedScalar h_node_sep[3]; in StgShur14Preprocess()
269 CeedScalar Ek_tot = 0; in StgShur14Preprocess()
271 const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1]; in StgShur14Preprocess()
281 CEED_QFUNCTION(ICsStg)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { in ICsStg()
282 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsStg()
283 const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[1]; in ICsStg()
284 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsStg()
288 CeedScalar qn[STG_NMODES_MAX], u[3], ubar[3], cij[6], eps, lt; in ICsStg()
289 const CeedScalar dx = stg_ctx->dx; in ICsStg()
290 const CeedScalar time = stg_ctx->time; in ICsStg()
291 const CeedScalar theta0 = stg_ctx->theta0; in ICsStg()
292 const CeedScalar P0 = stg_ctx->P0; in ICsStg()
293 const CeedScalar rho = P0 / (GasConstant(gas) * theta0); in ICsStg()
294 const CeedScalar nu = gas->mu / rho; in ICsStg()
297 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in ICsStg()
298 CeedScalar dXdx[3][3]; in ICsStg()
300 CeedScalar h_node_sep[3]; in ICsStg()
313 CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5] = {0.}; in ICsStg()
329 CEED_QFUNCTION(StgShur14Inflow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *cons… in StgShur14Inflow()
330 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in StgShur14Inflow()
331 const CeedScalar(*q_data_sur) = in[2]; in StgShur14Inflow()
332 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in StgShur14Inflow()
334 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in StgShur14Inflow()
335 CeedScalar(*jac_data_sur) = out[1]; in StgShur14Inflow()
338 CeedScalar qn[STG_NMODES_MAX], u[3], ubar[3], cij[6], eps, lt; in StgShur14Inflow()
342 const CeedScalar dx = stg_ctx->dx; in StgShur14Inflow()
343 const CeedScalar mu = stg_ctx->newtonian_ctx.mu; in StgShur14Inflow()
344 const CeedScalar time = stg_ctx->time; in StgShur14Inflow()
345 const CeedScalar theta0 = stg_ctx->theta0; in StgShur14Inflow()
346 const CeedScalar P0 = stg_ctx->P0; in StgShur14Inflow()
347 const CeedScalar cv = stg_ctx->newtonian_ctx.cv; in StgShur14Inflow()
348 const CeedScalar Rd = GasConstant(&stg_ctx->newtonian_ctx); in StgShur14Inflow()
349 const CeedScalar gamma = HeatCapacityRatio(&stg_ctx->newtonian_ctx); in StgShur14Inflow()
352 const CeedScalar rho = prescribe_T ? q[0][i] : P0 / (Rd * theta0); in StgShur14Inflow()
353 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in StgShur14Inflow()
354 CeedScalar wdetJb, dXdx[2][3], norm[3]; in StgShur14Inflow()
358 CeedScalar h_node_sep[3]; in StgShur14Inflow()
371 const CeedScalar E_kinetic = .5 * rho * Dot3(u, u); in StgShur14Inflow()
372 CeedScalar E_internal, P; in StgShur14Inflow()
383 const CeedScalar E = E_internal + E_kinetic; in StgShur14Inflow()
386 const CeedScalar u_normal = Dot3(norm, u); in StgShur14Inflow()
402 const CeedScalar U[] = {rho, u[0], u[1], u[2], E}, kmstress[6] = {0.}; in StgShur14Inflow()
409 …CTION(StgShur14Inflow_Jacobian)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in StgShur14Inflow_Jacobian()
410 const CeedScalar(*dq)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in StgShur14Inflow_Jacobian()
411 const CeedScalar(*q_data_sur)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in StgShur14Inflow_Jacobian()
412 const CeedScalar(*jac_data_sur)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[4]; in StgShur14Inflow_Jacobian()
413 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in StgShur14Inflow_Jacobian()
417 const CeedScalar cv = stg_ctx->newtonian_ctx.cv; in StgShur14Inflow_Jacobian()
418 const CeedScalar Rd = GasConstant(&stg_ctx->newtonian_ctx); in StgShur14Inflow_Jacobian()
419 const CeedScalar gamma = HeatCapacityRatio(&stg_ctx->newtonian_ctx); in StgShur14Inflow_Jacobian()
421 const CeedScalar theta0 = stg_ctx->theta0; in StgShur14Inflow_Jacobian()
430 const CeedScalar wdetJb = (implicit ? -1. : 1.) * q_data_sur[0][i]; in StgShur14Inflow_Jacobian()
433 CeedScalar velocity[3]; in StgShur14Inflow_Jacobian()
438 CeedScalar drho, dE, dP; in StgShur14Inflow_Jacobian()
442 CeedScalar dE_internal = drho * cv * theta0; in StgShur14Inflow_Jacobian()
443 CeedScalar dE_kinetic = .5 * drho * Dot3(velocity, velocity); in StgShur14Inflow_Jacobian()
451 const CeedScalar norm[3] = {q_data_sur[1][i], q_data_sur[2][i], q_data_sur[3][i]}; in StgShur14Inflow_Jacobian()
453 const CeedScalar u_normal = Dot3(norm, velocity); in StgShur14Inflow_Jacobian()
468 …NCTION(StgShur14InflowStrongQF)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in StgShur14InflowStrongQF()
469 const CeedScalar(*dXdx_q)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0]; in StgShur14InflowStrongQF()
470 const CeedScalar(*coords)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in StgShur14InflowStrongQF()
471 const CeedScalar(*scale) = (const CeedScalar(*))in[2]; in StgShur14InflowStrongQF()
472 const CeedScalar(*inv_Ektotal) = (const CeedScalar(*))in[3]; in StgShur14InflowStrongQF()
473 CeedScalar(*bcval)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in StgShur14InflowStrongQF()
477 CeedScalar u[3], ubar[3], cij[6], eps, lt; in StgShur14InflowStrongQF()
479 const CeedScalar dx = stg_ctx->dx; in StgShur14InflowStrongQF()
480 const CeedScalar time = stg_ctx->time; in StgShur14InflowStrongQF()
481 const CeedScalar theta0 = stg_ctx->theta0; in StgShur14InflowStrongQF()
482 const CeedScalar P0 = stg_ctx->P0; in StgShur14InflowStrongQF()
483 const CeedScalar rho = P0 / (GasConstant(gas) * theta0); in StgShur14InflowStrongQF()
484 const CeedScalar nu = gas->mu / rho; in StgShur14InflowStrongQF()
487 const CeedScalar x[] = {coords[0][i], coords[1][i], coords[2][i]}; in StgShur14InflowStrongQF()
488 const CeedScalar dXdx[2][3] = { in StgShur14InflowStrongQF()
493 CeedScalar h_node_sep[3]; in StgShur14InflowStrongQF()
503 CeedScalar qn[STG_NMODES_MAX]; in StgShur14InflowStrongQF()
511 CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5] = {0.}; in StgShur14InflowStrongQF()