| /honee/qfunctions/ |
| H A D | eulervortex.h | 81 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
91 P = rho * T; in Exact_Euler()
100 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
129 P = 1.; in Exact_Euler()
131 rho = P / (R * T); in Exact_Euler()
142 P = 1.; in Exact_Euler()
144 rho = P / (R * T); in Exact_Euler()
156 P = 1.; in Exact_Euler()
158 rho = P / (R * T); in Exact_Euler()
306 P = E_internal * (gamma - 1.); // P = pressure in Euler() local
[all …]
|
| H A D | shocktube.h | 85 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube() local
90 P = P_high; in Exact_ShockTube()
93 P = P_low; in Exact_ShockTube()
101 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
257 P = E_internal * (gamma - 1); // P = pressure in EulerShockTube() local
273 …1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[1] + (j ==… in EulerShockTube()
274 (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]); in EulerShockTube()
279 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… in EulerShockTube()
299 acoustic_vel = sqrt(gamma * P / rho); in EulerShockTube()
337 const CeedScalar sound_speed = sqrt(gamma * P / rho); in EulerShockTube()
|
| H A D | channel.h | 110 const CeedScalar P = s_inside.Y.pressure; in Channel_Inflow() local
114 const CeedScalar rho_in = P / ((gamma - 1) * e_internal); in Channel_Inflow()
129 …= 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_in * u_normal * s_exact.Y.velocity[j] + norm[j] * P); in Channel_Inflow()
132 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Inflow()
162 const CeedScalar P = context->P0; // pressure in Channel_Outflow() local
169 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Channel_Outflow()
172 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Outflow()
|
| H A D | stg_shur14.h | 376 CeedScalar E_internal, P; in StgShur14Inflow() local
381 P = rho * Rd * theta0; // interior rho with exterior T in StgShur14Inflow()
384 P = E_internal * (gamma - 1.); in StgShur14Inflow()
401 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P); in StgShur14Inflow()
404 v[4][i] -= wdetJb * u_normal * (E + P); in StgShur14Inflow()
|
| H A D | utils.h | 230 …t CeedScalar *mat_A, const CeedScalar *mat_B, CeedScalar *mat_C, CeedInt N, CeedInt M, CeedInt P) { in MatMatNM() argument
233 for (CeedInt k = 0; k < P; k++) mat_C[i * M + j] += mat_A[i * P + k] * mat_B[k * M + j]; in MatMatNM()
|
| /honee/src/ |
| H A D | dm-utils.c | 594 PetscInt field_offset = 0, num_comp, P, Q, dim; in ComputeFieldTabulationP2C() local
600 P = tabulation->Nb / num_comp; in ComputeFieldTabulationP2C()
604 PetscCall(PetscCalloc1(P * Q, interp)); in ComputeFieldTabulationP2C()
605 PetscCall(PetscCalloc1(P * Q * dim, grad)); in ComputeFieldTabulationP2C()
616 for (CeedInt p_ceed = 0; p_ceed < P; p_ceed++) { in ComputeFieldTabulationP2C()
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()
828 …PetscInt num_comp = basis_tabulation->Nc, P = basis_tabulation->Nb / num_comp, Q = basis_tabulatio… in DMPlexCeedBasisCreate() local
830 …PetscCallCeed(ceed, CeedBasisCreateH1(ceed, elem_topo, num_comp, P, Q, interp, grad, q_points, q_w… in DMPlexCeedBasisCreate()
999 PetscInt num_derivatives = 1, num_comp, P, Q = -1; in DMPlexCeedBasisCellToFaceCreate() local
[all …]
|
| /honee/doc/ |
| H A D | theory.md | 14 \frac{\partial \bm{U}}{\partial t} + \nabla \cdot \left( \frac{\bm{U} \otimes \bm{U}}{\rho} + P \bm…
15 \frac{\partial E}{\partial t} + \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\s…
20 …$k$ the thermal conductivity constant, $T$ represents the temperature, and $P$ the pressure, given…
23 P = \left( {c_p}/{c_v} -1\right) \left( E - {\bm{U}\cdot\bm{U}}/{(2 \rho)} \right) \, ,
57 {(\bm{U} \otimes \bm{U})}/{\rho} + P \bm{I}_3 \\
58 {(E + P)\bm{U}}/{\rho}
79 \bm{q}_N (\bm{x},t)^{(e)} = \sum_{k=1}^{P}\psi_k (\bm{x})\bm{q}_k^{(e)}
82 with $P=p+1$ the number of nodes in the element $e$.
173 Galerkin methods produce oscillations for transport-dominated problems (any time the cell Péclet nu…
219 …bm U + \bm U \otimes \diff\bm U)/\rho - (\bm U \otimes \bm U)/\rho^2 \diff\rho + \diff P \bm I_3 \\
[all …]
|
| H A D | runtime_options.md | 479 …$P, \bm{u}, T$), or `entropy` ($\frac{\gamma - s}{\gamma - 1} - \frac{\rho}{P} (e - c_v T),\ \frac…
|
| H A D | examples.md | 424 \frac{\partial \bm{U}}{\partial t} + \nabla \cdot \left( \frac{\bm{U} \otimes \bm{U}}{\rho} + P \bm…
425 \frac{\partial E}{\partial t} + \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} \right) &= 0 \, , \\
429 …}`zhang2011verification`, the mean flow for this problem is $\rho=1$, $P=1$, $T=P/\rho= 1$ (Specif…
436 There is no perturbation in the entropy $S=P/\rho^\gamma$ ($\delta S=0)$.
477 …c Vortex problem. The default initial conditions are $P=1$, $\rho=1$ for the driver section and $P…
649 - $P^2$
654 - $P^4$
|
| H A D | auxiliary.md | 140 …e prefix `cflpe` (e.g. `-ts_monitor_spanstats_cflpe`) obtains statistics for CFL and Péclet number.
|
| /honee/problems/ |
| H A D | newtonian.c | 741 const CeedScalar P = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in UnitTests_Newtonian() local
748 const CeedScalar entropy = log(P) - gamma * log(rho); in UnitTests_Newtonian()
749 const CeedScalar rho_div_p = rho / P; in UnitTests_Newtonian()
750 const CeedScalar Y0[5] = {P, u[0], u[1], u[2], T}; in UnitTests_Newtonian()
|