Home
last modified time | relevance | path

Searched refs:Y (Results 1 – 25 of 27) sorted by relevance

12

/honee/qfunctions/
H A Dnewtonian_state.h21 StatePrimitive Y; member
30 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y()
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()
53 CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy()
54 CeedScalar e_internal = gas.cv * s.Y.temperature; in TotalSpecificEnthalpy()
55 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
59 CeedScalar de_kinetic = Dot3(ds.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy_fwd()
60 CeedScalar de_internal = gas.cv * ds.Y.temperature; in TotalSpecificEnthalpy_fwd()
[all …]
H A Dblasius.h73 CeedScalar U_infty = Norm3(S_infty.Y.velocity); in BlasiusSolution()
84 CeedScalar Y[5]; in BlasiusSolution() local
85 Y[1] = U_infty * f[1]; in BlasiusSolution()
86 Y[2] = 0.5 * sqrt(nu * U_infty / (x0 + x[0] - x_inflow)) * (eta * f[1] - f[0]); in BlasiusSolution()
87 Y[3] = 0.; in BlasiusSolution()
88 Y[4] = S_infty.Y.temperature * h[0]; in BlasiusSolution()
89 Y[0] = rho_infty / h[0] * Rd * Y[4]; in BlasiusSolution()
90 return StateFromY(blasius->newt_ctx.gas, Y); in BlasiusSolution()
108 const CeedScalar U_infty = Norm3(S_infty.Y.velocity); in ICsBlasius()
135 const CeedScalar U_infty = Norm3(S_infty.Y.velocity); in Blasius_Inflow()
[all …]
H A Driemann_solver.h163 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL()
164 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL()
198 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL_fwd()
199 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL_fwd()
200 CeedScalar du_left = Dot3(dleft.Y.velocity, normal); in RiemannFlux_HLL_fwd()
201 CeedScalar du_right = Dot3(dright.Y.velocity, normal); in RiemannFlux_HLL_fwd()
223 side.Y.velocity[0] + (s_star - u_side) * normal[0], in RiemannFlux_HLLC_Star()
224 side.Y.velocity[1] + (s_star - u_side) * normal[1], in RiemannFlux_HLLC_Star()
225 side.Y.velocity[2] + (s_star - u_side) * normal[2], in RiemannFlux_HLLC_Star()
228 + (s_star - u_side) * (s_star + side.Y.pressure / denom) in RiemannFlux_HLLC_Star()
[all …]
H A Dchannel.h45 CeedScalar Y[5] = {0.}; in Exact_Channel() local
46 Y[0] = P0; in Exact_Channel()
47 Y[1] = umax * (1 - Square((x[1] - center) / H)); in Exact_Channel()
48 Y[2] = 0.; in Exact_Channel()
49 Y[3] = 0.; in Exact_Channel()
50 Y[4] = theta; in Exact_Channel()
52 return StateFromY(gas, Y); in Exact_Channel()
110 const CeedScalar P = s_inside.Y.pressure; in Channel_Inflow()
113 const CeedScalar e_internal = gas.cv * s_exact.Y.temperature; in Channel_Inflow()
115 const CeedScalar E_kinetic = .5 * rho_in * Dot3(s_exact.Y.velocity, s_exact.Y.velocity); in Channel_Inflow()
[all …]
H A Dsgs_dd_training.h52 …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 Ddifferential_filter.h43 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()
52 v1[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[2]; in DifferentialFilter_RHS()
[all …]
H A Dbc_slip.h29 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()
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()
85 const CeedScalar dvel_normal = Dot3(ds.Y.velocity, normal); in Slip_Jacobian()
86 … for (CeedInt j = 0; j < 3; j++) dvel_reflect[j] = ds.Y.velocity[j] - 2. * normal[j] * dvel_normal; in Slip_Jacobian()
87 …const CeedScalar dY_reflect[5] = {ds.Y.pressure, dvel_reflect[0], dvel_reflect[1], dvel_reflect[2]… in Slip_Jacobian()
H A Dtaylorgreen.h37 CeedScalar q[5], Y[5]; in ICsTaylorGreen() local
40 …Y[0] = reference.pressure + (density0 * Square(V0) / 16) * (cos(2 * x[0]) + cos(2 * x[1])) * (cos(… in ICsTaylorGreen()
41 Y[1] = context->u[0] + V0 * sin(x[0]) * cos(x[1]) * cos(x[2]); in ICsTaylorGreen()
42 Y[2] = context->u[1] - V0 * cos(x[0]) * sin(x[1]) * cos(x[2]); in ICsTaylorGreen()
43 Y[3] = context->u[2]; in ICsTaylorGreen()
44 Y[4] = reference.temperature; in ICsTaylorGreen()
46 State s = StateFromY(gas, Y); in ICsTaylorGreen()
H A Dstabilization.h13 CEED_QFUNCTION_HELPER void dYFromTau(const CeedScalar Y[5], const CeedScalar Tau_d[3], CeedScalar d… in dYFromTau()
14 dY[0] = Tau_d[0] * Y[0]; in dYFromTau()
15 dY[1] = Tau_d[1] * Y[1]; in dYFromTau()
16 dY[2] = Tau_d[1] * Y[2]; in dYFromTau()
17 dY[3] = Tau_d[1] * Y[3]; in dYFromTau()
18 dY[4] = Tau_d[2] * Y[4]; in dYFromTau()
89 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
90 velocity_term = Dot3(s.Y.velocity, gij_uj); in Tau_diagPrim()
H A Ddensitycurrent.h113 CeedScalar Y[5] = {0.}; in Exact_DC() local
114 Y[0] = P0 * pow(Pi, cp / Rd); in Exact_DC()
115 Y[1] = 0.0; in Exact_DC()
116 Y[2] = 0.0; in Exact_DC()
117 Y[3] = 0.0; in Exact_DC()
118 Y[4] = Pi * theta; in Exact_DC()
120 return StateFromY(gas, Y); in Exact_DC()
H A Dgaussianwave.h42 …const CeedScalar e_kinetic = 0.5 * S_infty.U.density * Dot3(S_infty.Y.velocity, S_infty.Y.velocity… in IC_GaussianWave()
47 U[1] = S_infty.Y.velocity[0] * U[0]; in IC_GaussianWave()
48 U[2] = S_infty.Y.velocity[1] * U[0]; in IC_GaussianWave()
49 U[3] = S_infty.Y.velocity[2] * U[0]; in IC_GaussianWave()
50 U[4] = S_infty.Y.pressure / (gamma - 1) * perturbation + e_kinetic; in IC_GaussianWave()
H A Dbc_outflow.h55 StatePrimitive y_ext = s_int.Y; in RiemannOutflow()
72 ViscousEnergyFlux(gas, s_int.Y, grad_s, stress, Fe); in RiemannOutflow()
129 StatePrimitive y_ext = s_int.Y, dy_ext = ds_int.Y; in RiemannOutflow_Jacobian()
134 const CeedScalar u_normal = Dot3(s_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
135 const CeedScalar du_normal = Dot3(ds_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
154 ViscousEnergyFlux_fwd(gas, s_int.Y, ds_int.Y, grad_ds, stress, dstress, dFe); in RiemannOutflow_Jacobian()
199 s.Y.pressure = outflow->pressure; in PressureOutflow()
212 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in PressureOutflow()
270 s.Y.pressure = outflow->pressure; in PressureOutflow_Jacobian()
271 ds.Y.pressure = 0.; in PressureOutflow_Jacobian()
[all …]
H A Dmonitor_totalkineticenergy.h25 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.h245 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()
299 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
303 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in IFunction_AdvectionGeneric()
304 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
374 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in MassFunction_AdvectionGeneric()
434 div_u += grad_s[k].Y.velocity[j]; in RHSFunction_AdvectionGeneric()
437 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in RHSFunction_AdvectionGeneric()
440 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in RHSFunction_AdvectionGeneric()
H A Dnewtonian.h118 CeedScalar damp_state[5] = {s.Y.pressure - pressure, 0, 0, 0, 0}, idl_residual[5] = {0.}; in InternalDampingLayer_IFunction_Integrand()
137 CeedScalar damp_state[5] = {ds.Y.pressure, 0, 0, 0, 0}, idl_residual[5] = {0.}; in InternalDampingLayer_IJacobian_Integrand()
217 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in RHSFunction_Newtonian()
236 …CeedScalar damp_state[5] = {s.Y.pressure - idl_pressure, 0, 0, 0, 0}, idl_residual[5] = {0.}; in RHSFunction_Newtonian()
289 ViscousEnergyFlux(gas, s.Y, grad_s, stress, F_visc_energy); in IFunction_Newtonian_Integrand()
400 ViscousEnergyFlux_fwd(gas, s.Y, ds.Y, grad_ds, stress, dstress, dF_visc_energy); in IJacobian_Newtonian_Integrand()
505 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in BoundaryIntegral()
571 ViscousEnergyFlux_fwd(gas, s.Y, ds.Y, grad_ds, stress, dstress, dFe); in BoundaryIntegral_Jacobian()
623 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in DivDiffusiveFluxVolumeRHS_NS()
676 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in DivDiffusiveFluxBoundaryRHS_NS()
[all …]
H A Dsgs_dd_model.h203 CEED_QFUNCTION_HELPER int FluxSubgridStress(const StatePrimitive Y, const CeedScalar km_sgs[6], Cee… in FluxSubgridStress() argument
210 Flux[4][j] = Y.velocity[0] * sgs[0][j] + Y.velocity[1] * sgs[1][j] + Y.velocity[2] * sgs[2][j]; in FluxSubgridStress()
232 FluxSubgridStress(s.Y, km_sgs_i, Flux); in IFunction_NodalSgs()
H A Dmonitor_cfl.h30 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()
/honee/qfunctions/spanstats/
H A Dturbulence.h57 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()
65 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Z][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
69 v[TURB_MEAN_MOMENTUMFLUX_XX][i] = wdetJ * s.U.momentum[0] * s.Y.velocity[0]; in ChildStatsCollection()
[all …]
H A Dcflpe.h44 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/src/
H A Dpetsc_ops.c208 …Core(Vec X, Vec X_loc, CeedVector x_ceed, CeedVector y_ceed, Vec Y_loc, Vec Y, OperatorApplyContex… in ApplyCeedOperator_Core() argument
219 PetscCall(PetscLogEventBegin(HONEE_CeedOperatorApply, X, Y, 0, 0)); in ApplyCeedOperator_Core()
224 PetscCall(PetscLogEventEnd(HONEE_CeedOperatorApply, X, Y, 0, 0)); in ApplyCeedOperator_Core()
229 if (Y) PetscCall(DMLocalToGlobal(ctx->dm_y, Y_loc, ADD_VALUES, Y)); in ApplyCeedOperator_Core()
233 PetscErrorCode ApplyCeedOperatorGlobalToGlobal(Vec X, Vec Y, OperatorApplyContext ctx) { in ApplyCeedOperatorGlobalToGlobal() argument
237 PetscCall(VecZeroEntries(Y)); in ApplyCeedOperatorGlobalToGlobal()
243 PetscCall(ApplyCeedOperator_Core(X, X_loc, ctx->x_ceed, ctx->y_ceed, Y_loc, Y, ctx, false)); in ApplyCeedOperatorGlobalToGlobal()
251 PetscErrorCode ApplyCeedOperatorLocalToGlobal(Vec X_loc, Vec Y, OperatorApplyContext ctx) { in ApplyCeedOperatorLocalToGlobal() argument
255 PetscCall(VecZeroEntries(Y)); in ApplyCeedOperatorLocalToGlobal()
260 PetscCall(ApplyCeedOperator_Core(NULL, X_loc, ctx->x_ceed, ctx->y_ceed, Y_loc, Y, ctx, false)); in ApplyCeedOperatorLocalToGlobal()
H A Dmat-ceed.c388 Vec Y; in MatCreateCeed() local
390 PetscCall(DMGetGlobalVector(dm_y, &Y)); in MatCreateCeed()
391 PetscCall(VecGetSize(Y, &Y_g_size)); in MatCreateCeed()
392 PetscCall(VecGetLocalSize(Y, &Y_l_size)); in MatCreateCeed()
393 PetscCall(DMRestoreGlobalVector(dm_y, &Y)); in MatCreateCeed()
1630 PetscErrorCode MatMult_Ceed(Mat A, Vec X, Vec Y) { in MatMult_Ceed() argument
1635 PetscCall(PetscLogEventBegin(ctx->log_event_mult, A, X, Y, NULL)); in MatMult_Ceed()
1654 PetscCall(PetscLogEventBegin(ctx->log_event_ceed_mult, A, X, Y, NULL)); in MatMult_Ceed()
1661 PetscCall(PetscLogEventEnd(ctx->log_event_ceed_mult, A, X, Y, NULL)); in MatMult_Ceed()
1668 PetscCall(VecZeroEntries(Y)); in MatMult_Ceed()
[all …]
/honee/include/
H A Dmat-ceed-impl.h40 PETSC_CEED_EXTERN PetscErrorCode MatMult_Ceed(Mat A, Vec X, Vec Y);
41 PETSC_CEED_EXTERN PetscErrorCode MatMultTranspose_Ceed(Mat A, Vec Y, Vec X);
H A Dpetsc_ops.h28 PetscErrorCode ApplyCeedOperatorGlobalToGlobal(Vec X, Vec Y, OperatorApplyContext ctx);
30 PetscErrorCode ApplyCeedOperatorLocalToGlobal(Vec X_loc, Vec Y, OperatorApplyContext ctx);
/honee/problems/
H A Dblasius.c23 CeedScalar U_infty = Norm3(S_infty.Y.velocity); in CompressibleBlasiusResidual()
27 …PetscScalar Ma = Mach(gas, S_infty.Y.temperature, U_infty), Pr = Prandtl(gas), gamma = HeatCapacit… in CompressibleBlasiusResidual()
61 r[N] = h[0] - blasius->T_wall / S_infty.Y.temperature; in CompressibleBlasiusResidual()
221 PetscCall(SetupStg(comm, dm, problem, honee, weakT, S_infty.Y.temperature, S_infty.Y.pressure)); in NS_BLASIUS()
H A Dbc_freestream.c357 CeedScalar u_left0 = Dot3(left0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
358 CeedScalar u_right0 = Dot3(right0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
373 CeedScalar du_left = Dot3(dleft0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
374 CeedScalar du_right = Dot3(dright0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
391 CeedScalar u_left1 = Dot3(left1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
392 CeedScalar u_right1 = Dot3(right1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()

12