Home
last modified time | relevance | path

Searched refs:s (Results 1 – 25 of 31) sorted by relevance

12

/honee/qfunctions/
H A Dnewtonian_state.h24 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U() argument
25 U[0] = s.density; in UnpackState_U()
26 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
27 U[4] = s.E_total; in UnpackState_U()
30 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y() argument
31 Y[0] = s.pressure; in UnpackState_Y()
32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
33 Y[4] = s.temperature; in UnpackState_Y()
36 CEED_QFUNCTION_HELPER void UnpackState_V(StateEntropy s, CeedScalar V[5]) { in UnpackState_V() argument
37 V[0] = s.S_density; in UnpackState_V()
[all …]
H A Dnewtonian.h24 State s = StateFromPrimitive(gas, context->reference); in ICsNewtonianIG() local
25 StateToQ(gas, s, q, state_var); in ICsNewtonianIG()
56 const State s = StateFromQ(gas, qi, state_var); in MassFunction_Newtonian() local
70 Tau_diagPrim(context->tau_coeffs, gas, s, dXdx, context->dt, Tau_d); in MassFunction_Newtonian()
71 Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff, stab); in MassFunction_Newtonian()
88 … InternalDampingLayer_Residual(const NewtonianIGProperties gas, const State s, const CeedScalar si… in InternalDampingLayer_Residual() argument
91 State damp_s = StateFromY_fwd(gas, s, damp_Y); in InternalDampingLayer_Residual()
114 CEED_QFUNCTION_HELPER void InternalDampingLayer_IFunction_Integrand(const State s, const NewtonianI… in InternalDampingLayer_IFunction_Integrand() argument
118 CeedScalar damp_state[5] = {s.Y.pressure - pressure, 0, 0, 0, 0}, idl_residual[5] = {0.}; in InternalDampingLayer_IFunction_Integrand()
119 InternalDampingLayer_Residual(gas, s, sigma_, damp_state, idl_residual); in InternalDampingLayer_IFunction_Integrand()
[all …]
H A Dsgs_dd_training.h44 const State s = StateFromQ(gas, qi, state_var); in ComputeSGS_DDAnisotropicTrainingDataNodal() local
48 …ComputeSgsDDInputs(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, eigenvectors, inputs, &g… in ComputeSGS_DDAnisotropicTrainingDataNodal()
52 …v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
53 …v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
54 …v[2 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] - Square(s.Y.velocity[2])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
55 …v[3 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] - s.Y.velocity[1] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
56 …v[4 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] - s.Y.velocity[0] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
57 …v[5 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XY][i] - s.Y.velocity[0] * s.Y.velocity[1]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
H A Dstabilization.h24 CEED_QFUNCTION_HELPER void StabilizationMatrix(const NewtonianIGProperties gas, const State s, cons… in StabilizationMatrix() argument
31 State ds = StateFromY_fwd(gas, s, dY); in StabilizationMatrix()
32 FluxInviscid_fwd(gas, s, ds, dF); in StabilizationMatrix()
40 …n(StabilizationType stab_type, const NewtonianIGProperties gas, const State s, const CeedScalar Ta… in Stabilization() argument
49 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
52 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
56 StabilizationMatrix(gas, s, Tau_d, strong_residual, stab); in Stabilization()
66 …u_diagPrim(TauDiagCoefficients tau_coeffs, NewtonianIGProperties gas, State s, const CeedScalar dX… in Tau_diagPrim() argument
76 const CeedScalar rho = s.U.density; in Tau_diagPrim()
89 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
[all …]
H A Dbc_slip.h22 State s = StateFromQ(gas, qi, state_var); in Slip() local
29 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip()
30 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip()
31 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip()
34 StateConservative flux = RiemannFlux_HLLC(gas, s, s_reflect, normal); in Slip()
75 State s = StateFromQ(gas, qi, state_var); in Slip_Jacobian() local
76 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Slip_Jacobian()
79 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip_Jacobian()
80 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip_Jacobian()
81 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip_Jacobian()
[all …]
H A Ddifferential_filter.h41 const State s = StateFromQ(gas, qi, state_var); in DifferentialFilter_RHS() local
43 v0[DIFF_FILTER_PRESSURE][i] = wdetJ * s.Y.pressure; in DifferentialFilter_RHS()
44 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
45 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
46 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
47 v0[DIFF_FILTER_TEMPERATURE][i] = wdetJ * s.Y.temperature; in DifferentialFilter_RHS()
48 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
51 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
[all …]
H A Dblasius.h114 State s = BlasiusSolution(context, x, x0, x_inflow, S_infty.U.density, &t12); in ICsBlasius() local
117 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsBlasius()
146 State s = BlasiusSolution(context, x, x0, context->x_inflow, rho_0, &t12); in Blasius_Inflow() local
153 s.U.density = s_int.U.density; in Blasius_Inflow()
154 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
156 s.U.E_total = s_int.U.E_total; in Blasius_Inflow()
157 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
161 FluxInviscid(gas, s, Flux_inviscid); in Blasius_Inflow()
201 State s = BlasiusSolution(context, x, x0, 0, rho_0, &t12); in Blasius_Inflow_Jacobian() local
209 CeedScalar dE_kinetic = .5 * drho * Dot3(s.Y.velocity, s.Y.velocity); in Blasius_Inflow_Jacobian()
[all …]
H A Dmonitor_totalkineticenergy.h18 const State s = StateFromQ(gas, qi, state_var); in MonitorTotalKineticEnergy() local
23 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in MonitorTotalKineticEnergy()
25 v[0][i] = wdetJ * 0.5 * s.U.density * Dot3(s.Y.velocity, s.Y.velocity); in MonitorTotalKineticEnergy()
32 v[3][i] = wdetJ * s.Y.pressure * div_u; in MonitorTotalKineticEnergy()
H A Dadvection.h210 …ysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s, in StatePhysicalGradientFromReference_ND() argument
220 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
233 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
240 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd… in Tau() argument
245 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
254 MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj); in Tau()
255 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * Squa… in Tau()
282 const State s = StateFromU(gas, qi); in IFunction_AdvectionGeneric() local
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
303 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in IFunction_AdvectionGeneric()
[all …]
H A Dbc_freestream.h27 const State s = StateFromQ(gas, qi, state_var); in Freestream() local
36 flux = RiemannFlux_HLL(gas, s, context->S_infty, normal); in Freestream()
39 flux = RiemannFlux_HLLC(gas, s, context->S_infty, normal); in Freestream()
96 State s = StateFromQ(gas, qi, state_var); in Freestream_Jacobian() local
97 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Freestream_Jacobian()
102 dflux = RiemannFlux_HLL_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
105 dflux = RiemannFlux_HLLC_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
H A Dbc_outflow.h198 State s = StateFromQ(gas, qi, state_var); in PressureOutflow() local
199 s.Y.pressure = outflow->pressure; in PressureOutflow()
206 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in PressureOutflow()
212 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in PressureOutflow()
215 FluxInviscid(gas, s, F_inviscid); in PressureOutflow()
268 State s = StateFromQ(gas, qi, state_var); in PressureOutflow_Jacobian() local
269 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in PressureOutflow_Jacobian()
270 s.Y.pressure = outflow->pressure; in PressureOutflow_Jacobian()
274 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_dq, dXdx, grad_ds); in PressureOutflow_Jacobian()
281 ViscousEnergyFlux_fwd(gas, s.Y, ds.Y, grad_ds, stress, dstress, dFe); in PressureOutflow_Jacobian()
[all …]
H A Dmonitor_cfl.h19 const State s = StateFromQ(gas, qi, state_var); in MonitorCFL() local
30 v[i] = CalculateCFL_2D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
40 v[i] = CalculateCFL_3D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
H A Dsgs_dd_model.h87 const State s = StateFromQ(gas, qi, state_var); in ComputeSgsDDNodal_Fused() local
90 ComputeSgsDD_Fused(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, km_sgs, sgsdd_ctx); in ComputeSgsDDNodal_Fused()
131 const State s = StateFromQ(gas, qi, state_var); in ComputeSgsDDNodal_Sequential_Inputs() local
134 …ComputeSgsDDInputs(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, eigenvectors, model_inpu… in ComputeSgsDDNodal_Sequential_Inputs()
225 const State s = StateFromQ(newt_ctx->gas, qi, state_var); in IFunction_NodalSgs() local
232 FluxSubgridStress(s.Y, km_sgs_i, Flux); in IFunction_NodalSgs()
H A Dvelocity_gradient_projection.h23 const State s = StateFromQ(gas, qi, state_var); in VelocityGradientProjectionRHS() local
25 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, (CeedScalar *)Grad_q, dXdx, grad_s); in VelocityGradientProjectionRHS()
H A Dtaylorgreen.h46 State s = StateFromY(gas, Y); in ICsTaylorGreen() local
47 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsTaylorGreen()
H A Ddensitycurrent.h135 State s = Exact_DC(3, 0., x, 5, ctx); in ICsDC() local
137 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsDC()
H A Dchannel.h67 State s = Exact_Channel(3, 0., x, 5, ctx); in ICsChannel() local
69 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsChannel()
H A Dstg_shur14.h317 State s = StateFromY(gas, Y); in ICsStg() local
318 StateToQ(gas, s, q, context->state_var); in ICsStg()
510 State s = StateFromY(gas, Y); in StgShur14InflowStrongQF() local
511 StateToQ(gas, s, q, stg_ctx->newt_ctx.state_var); in StgShur14InflowStrongQF()
/honee/qfunctions/spanstats/
H A Dturbulence.h54 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection() local
56 v[TURB_MEAN_DENSITY][i] = wdetJ * s.U.density; in ChildStatsCollection()
57 v[TURB_MEAN_PRESSURE][i] = wdetJ * s.Y.pressure; in ChildStatsCollection()
58 v[TURB_MEAN_PRESSURE_SQUARED][i] = wdetJ * Square(s.Y.pressure); in ChildStatsCollection()
59 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
60 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
61 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
62 v[TURB_MEAN_DENSITY_TEMPERATURE][i] = wdetJ * s.U.density * s.Y.temperature; in ChildStatsCollection()
63 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
64 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
[all …]
H A Dcflpe.h30 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection_CflPe() local
44 cfl = CalculateCFL_2D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
45 Pe = CalculatePe_2D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
57 cfl = CalculateCFL_3D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
58 Pe = CalculatePe_3D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
/honee/tests/
H A Dsmartsim_regression_framework.py29 s = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
30 s.bind(('', 0))
31 addr = s.getsockname()
32 s.close()
/honee/
H A DREADME.md34 HONEE requires libCEED's `main` development branch, which can be [cloned from Github](https://githu…
45 HONEE requires PETSc's `main` development branch, which can be [cloned from GitLab](https://gitlab.…
84 If this is not true, you can add a `--download-mpich` flag to PETSc's configuration command to have…
H A DCONTRIBUTING.md33 ## [Developer's Certificate of Origin 1.1](https://developercertificate.org/)
43 …itely and may be redistributed consistent with this project or the open source license(s) involved.
/honee/doc/
H A Druntime_options.md166 - Run unit tests of Riemann problem solver and it's Jacobian
182 - View PETSc `TS` solver configuration before it begins it's solve
461 - `m/s^2`
471 - `Pa s`
479 … `entropy` ($\frac{\gamma - s}{\gamma - 1} - \frac{\rho}{P} (e - c_v T),\ \frac{\rho}{P} \bm{u},\ …
H A Dexamples.md6 …014`. It is intended to test non-reflecting/Riemann boundary conditions. It's primarily intended f…
171 - `1/s`
207 - `m/s`
318 - `m/s`
410 - `m/s`
477 This test problem is based on Sod's Shock Tube (from{cite}`sodshocktubewiki`), a canonical test cas…
573 Some require specific advection profiles, other's can be used with multiple.
582 …This is similar to the bubble ICs, but uses a cosine wave to define the bubble and it's radius is …
705 - `1/s`

12