Home
last modified time | relevance | path

Searched refs:q (Results 1 – 24 of 24) sorted by relevance

/honee/qfunctions/
H A Dadvection.h64 …eneric(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, CeedScalar q[], void *ctx) { in Exact_AdvectionGeneric() argument
80 q[0] = 1.; in Exact_AdvectionGeneric()
81 q[1] = -(y - center[1]); in Exact_AdvectionGeneric()
82 q[2] = (x - center[0]); in Exact_AdvectionGeneric()
83 q[3] = 0; in Exact_AdvectionGeneric()
86 q[0] = 1.; in Exact_AdvectionGeneric()
87 q[1] = wind[0]; in Exact_AdvectionGeneric()
88 q[2] = wind[1]; in Exact_AdvectionGeneric()
89 q[3] = dim == 2 ? 0. : wind[2]; in Exact_AdvectionGeneric()
92 q[0] = 1.; in Exact_AdvectionGeneric()
[all …]
H A Deulervortex.h58 …_Euler(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, CeedScalar q[], void *ctx) { in Exact_Euler() argument
96 q[0] = rho; in Exact_Euler()
97 q[1] = rho * u[0]; in Exact_Euler()
98 q[2] = rho * u[1]; in Exact_Euler()
99 q[3] = rho * u[2]; in Exact_Euler()
100 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
107 q[0] = rho; in Exact_Euler()
108 q[1] = rho * u[0]; in Exact_Euler()
109 q[2] = rho * u[1]; in Exact_Euler()
110 q[3] = rho * u[2]; in Exact_Euler()
[all …]
H A Dsgs_dd_model.h69 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ComputeSgsDDNodal_Fused() local
79 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ComputeSgsDDNodal_Fused()
112 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ComputeSgsDDNodal_Sequential_Inputs() local
123 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ComputeSgsDDNodal_Sequential_Inputs()
216 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in IFunction_NodalSgs() local
224 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in IFunction_NodalSgs()
H A Dchannel.h68 CeedScalar q[5]; in ICsChannel() local
69 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsChannel()
70 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsChannel()
79 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Channel_Inflow() local
108 for (CeedInt j = 0; j < 5; j++) q_inside[j] = q[j][i]; in Channel_Inflow()
141 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Channel_Outflow() local
153 const CeedScalar rho = q[0][i]; in Channel_Outflow()
154 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Channel_Outflow()
155 const CeedScalar E = q[4][i]; in Channel_Outflow()
H A Dnewtonian.h23 CeedScalar q[5]; in ICsNewtonianIG() local
25 StateToQ(gas, s, q, state_var); in ICsNewtonianIG()
26 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsNewtonianIG()
45 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in MassFunction_Newtonian() local
54 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MassFunction_Newtonian()
191 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in RHSFunction_Newtonian() local
206 for (int j = 0; j < 5; j++) U[j] = q[j][i]; in RHSFunction_Newtonian()
312 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in IFunction_Newtonian() local
325 const CeedScalar q_i[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in IFunction_Newtonian()
481 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in BoundaryIntegral() local
[all …]
H A Dvelocity_gradient_projection.h11 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in VelocityGradientProjectionRHS() local
19 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in VelocityGradientProjectionRHS()
H A Dshocktube.h72 …ckTube(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, CeedScalar q[], void *ctx) { in Exact_ShockTube() argument
97 q[0] = rho; in Exact_ShockTube()
98 q[1] = rho * u[0]; in Exact_ShockTube()
99 q[2] = rho * u[1]; in Exact_ShockTube()
100 q[3] = rho * u[2]; in Exact_ShockTube()
101 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
176 CeedScalar q[5]; in ICsShockTube() local
178 Exact_ShockTube(3, 0., x, 5, q, ctx); in ICsShockTube()
180 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsShockTube()
210 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in EulerShockTube() local
[all …]
H A Dmonitor_totalkineticenergy.h10 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in MonitorTotalKineticEnergy() local
17 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MonitorTotalKineticEnergy()
H A Dmonitor_cfl.h11 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in MonitorCFL() local
18 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MonitorCFL()
H A Dsgs_dd_training.h25 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ComputeSGS_DDAnisotropicTrainingDataNodal() local
36 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ComputeSGS_DDAnisotropicTrainingDataNodal()
H A Ddifferential_filter.h30 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in DifferentialFilter_RHS() local
39 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in DifferentialFilter_RHS()
75 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in DifferentialFilter_LHS_N() local
126 v[j][i] = wdetJ * q[j][i]; in DifferentialFilter_LHS_N()
162 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in DifferentialFilter_MMS_RHS() local
168 v[0][i] = wdetJ * q[0][i]; in DifferentialFilter_MMS_RHS()
H A Dbc_outflow.h39 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in RiemannOutflow() local
52 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in RiemannOutflow()
187 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in PressureOutflow() local
197 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in PressureOutflow()
H A Dbc_slip.h13 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Slip() local
21 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in Slip()
H A Dbc_freestream.h17 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Freestream() local
26 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in Freestream()
H A Dtaylorgreen.h37 CeedScalar q[5], Y[5]; in ICsTaylorGreen() local
47 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsTaylorGreen()
48 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsTaylorGreen()
H A Dstg_shur14.h316 CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5]; in ICsStg() local
318 StateToQ(gas, s, q, context->state_var); in ICsStg()
320 q0[j][i] = q[j]; in ICsStg()
333 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in StgShur14Inflow() local
356 const CeedScalar rho = prescribe_T ? q[0][i] : P0 / (Rd * theta0); in StgShur14Inflow()
383 E_internal = q[4][i] - E_kinetic; // uses prescribed rho and u, E from solution in StgShur14Inflow()
509 CeedScalar Y[5] = {P0, u[0], u[1], u[2], theta0}, q[5]; in StgShur14InflowStrongQF() local
511 StateToQ(gas, s, q, stg_ctx->newt_ctx.state_var); in StgShur14InflowStrongQF()
514 q[4] = 0.; // Don't set energy in StgShur14InflowStrongQF()
517 q[0] = 0; // Don't set pressure in StgShur14InflowStrongQF()
[all …]
H A Ddensitycurrent.h136 CeedScalar q[5]; in ICsDC() local
137 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsDC()
138 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsDC()
H A Dblasius.h115 CeedScalar q[5]; in ICsBlasius() local
117 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsBlasius()
118 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsBlasius()
126 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Blasius_Inflow() local
148 for (CeedInt j = 0; j < 5; j++) qi[j] = q[j][i]; in Blasius_Inflow()
/honee/doc/
H A Dtheory.md31 \frac{\partial \bm{q}}{\partial t} + \nabla \cdot \bm{F}(\bm{q}) -S(\bm{q}) = 0 \, ,
37 \bm{q} =
54 \bm{F}(\bm{q}) &=
65 S(\bm{q}) &=
79 \bm{q}_N (\bm{x},t)^{(e)} = \sum_{k=1}^{P}\psi_k (\bm{x})\bm{q}_k^{(e)}
88 …Omega} \bm v \cdot \left(\frac{\partial \bm{q}_N}{\partial t} + \nabla \cdot \bm{F}(\bm{q}_N) - \b…
97 \int_{\Omega} \bm v \cdot \left( \frac{\partial \bm{q}_N}{\partial t} - \bm{S}(\bm{q}_N) \right) \…
98 - \int_{\Omega} \nabla \bm v \!:\! \bm{F}(\bm{q}_N)\,dV & \\
120 \bm{q}_N^{n+1} = \bm{q}_N^n + \Delta t \sum_{i=1}^{s} b_i k_i \, ,
127 k_1 &= f(t^n, \bm{q}_N^n)\\
[all …]
H A Dexamples.md482 \int_{\Omega} \nu_{SHOCK} \nabla \bm v \!:\! \nabla \bm q dV
/honee/qfunctions/spanstats/
H A Dcflpe.h20 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ChildStatsCollection_CflPe() local
29 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ChildStatsCollection_CflPe()
H A Dturbulence.h42 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ChildStatsCollection() local
53 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ChildStatsCollection()
117 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ChildStatsCollectionMMSTest_Error() local
126 v[0][i] = wdetJ * Square(ChildStatsCollectionTest_Exact(x_i) - q[0][i]); in ChildStatsCollectionMMSTest_Error()
/honee/examples/meshes/
H A DMakefile3 ALL = $(foreach q,1,$(foreach n,08 12 20 40,cylinder-q$(q)-n$(n).msh))
/honee/src/
H A Ddm-utils.c295 PetscQuadrature *q) { in PetscDTUniformTensorQuadrature() argument
362 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, q)); in PetscDTUniformTensorQuadrature()
363 PetscCall(PetscQuadratureSetCellType(*q, ct)); in PetscDTUniformTensorQuadrature()
364 PetscCall(PetscQuadratureSetOrder(*q, num_points)); in PetscDTUniformTensorQuadrature()
365 PetscCall(PetscQuadratureSetData(*q, dim, num_comp, num_total_points, coords, weights)); in PetscDTUniformTensorQuadrature()
366 PetscCall(PetscObjectChangeTypeName((PetscObject)*q, "UniformTensor")); in PetscDTUniformTensorQuadrature()
615 for (CeedInt q = 0; q < Q; q++) { in ComputeFieldTabulationP2C() local
618 …(*interp)[q * P + p_ceed] = tabulation->T[0][((face * Q + q) * P * num_comp + p_petsc) * num_comp … in ComputeFieldTabulationP2C()
620 …(*grad)[(d * Q + q) * P + p_ceed] = tabulation->T[1][(((face * Q + q) * P * num_comp + p_petsc) * … in ComputeFieldTabulationP2C()