| /honee/qfunctions/ |
| H A D | eulervortex.h | 182 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
183 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); in ConvectiveFluxJacobian_Euler()
184 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
185 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
295 for (CeedInt k = 0; k < 3; k++) { in Euler() local
296 drhodx[j] += drho[k] * dXdx[k][j]; in Euler()
297 dEdx[j] += dE[k] * dXdx[k][j]; in Euler()
299 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in Euler()
300 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in Euler()
312 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler() local
[all …]
|
| H A D | shocktube.h | 115 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
116 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); in ConvectiveFluxJacobian_Euler()
117 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
118 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
245 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
246 du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho; in EulerShockTube()
247 drhodx[j] += drho[k] * dXdx[k][j]; in EulerShockTube()
248 dEdx[j] += dE[k] * dXdx[k][j]; in EulerShockTube()
250 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in EulerShockTube()
251 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in EulerShockTube()
[all …]
|
| H A D | setupgeo_helpers.h | 149 for (CeedInt k = 0; k < 2; k++) { in InvertBoundaryMappingJacobian_3D() local
150 for (CeedInt l = 0; l < 3; l++) dxdXTdxdX[j][k] += dxdX[l][j] * dxdX[l][k]; in InvertBoundaryMappingJacobian_3D()
165 for (CeedInt k = 0; k < 3; k++) { in InvertBoundaryMappingJacobian_3D() local
166 dXdx[j][k] = 0; in InvertBoundaryMappingJacobian_3D()
167 for (CeedInt l = 0; l < 2; l++) dXdx[j][k] += dxdXTdxdX_inv[l][j] * dxdX[k][l]; in InvertBoundaryMappingJacobian_3D()
|
| H A D | advection.h | 215 for (CeedInt k = 0; k < 2; k++) { in StatePhysicalGradientFromReference_ND() local
218 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * … in StatePhysicalGradientFromReference_ND()
220 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
227 for (CeedInt k = 0; k < 3; k++) { in StatePhysicalGradientFromReference_ND() local
230 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * … in StatePhysicalGradientFromReference_ND()
231 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
233 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
298 for (CeedInt k = 0; k < dim; k++) { in IFunction_AdvectionGeneric() local
299 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
318 for (CeedInt k = 0; k < dim; k++) grad_v[k][4][i] -= wdetJ * Fe_dXdx[k]; in IFunction_AdvectionGeneric() local
[all …]
|
| H A D | utils.h | 131 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[k * N + j]; in MatMatN() local
138 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[j * N + k]; in MatMatN() local
149 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[k * N + j]; in MatMatN() local
156 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[j * N + k]; in MatMatN() local
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() local
|
| H A D | newtonian.h | 75 for (CeedInt k = 0; k < 3; k++) { in MassFunction_Newtonian() local
76 …Grad_v[k][j][i] = wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[k… in MassFunction_Newtonian()
227 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * … in RHSFunction_Newtonian() local
251 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] *… in RHSFunction_Newtonian() local
347 for (CeedInt k = 0; k < 3; k++) in IFunction_Newtonian() local
348 …grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1] * dXdx[k][1] + grad_v_i[j]… in IFunction_Newtonian()
458 …for (int k = 0; k < 3; k++) grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1… in IJacobian_Newtonian() local
629 for (CeedInt k = 0; k < 3; k++) { in DivDiffusiveFluxVolumeRHS_NS() local
630 Grad_v[k][j - 1][i] = -wdetJ * Dot3(dXdx[k], Fdiff[j]); in DivDiffusiveFluxVolumeRHS_NS()
731 for (CeedInt k = 0; k < 3; k++) { in DiffusiveFluxRHS_NS() local
[all …]
|
| H A D | newtonian_state.h | 46 …_QFUNCTION_HELPER CeedScalar Prandtl(NewtonianIGProperties gas) { return gas.cp * gas.mu / gas.k; } in Prandtl()
481 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j]; in FluxTotal() local
491 for (CeedInt k = 0; k < 3; k++) { in FluxTotal_Boundary() local
492 Flux[k + 1] += (F_inviscid[j].momentum[k] - stress[k][j]) * normal[j]; in FluxTotal_Boundary()
501 for (CeedInt k = 0; k < 3; k++) Flux[k + 1] = F_inviscid_normal.momentum[k]; in FluxTotal_RiemannBoundary() local
504 for (CeedInt k = 0; k < 3; k++) { in FluxTotal_RiemannBoundary() local
505 Flux[k + 1] -= stress[k][j] * normal[j]; in FluxTotal_RiemannBoundary()
563 …s[0][i] - Y.velocity[1] * stress[1][i] - Y.velocity[2] * stress[2][i] - gas.k * grad_s[i].Y.temper… in ViscousEnergyFlux()
571 … Y.velocity[2] * dstress[2][i] - dY.velocity[2] * stress[2][i] - gas.k * grad_ds[i].Y.temperature; in ViscousEnergyFlux_fwd()
590 for (CeedInt k = 0; k < 5; k++) dqi[k] = dq[k][j]; in StatePhysicalGradientFromReference() local
[all …]
|
| H A D | differential_filter.h | 128 for (int k = 0; k < 3; k++) { in DifferentialFilter_LHS_N() local
129 …dq[k] = Grad_q[0 * N + j][i] * dXdx[0][k] + Grad_q[1 * N + j][i] * dXdx[1][k] + Grad_q[2 * N + j][… in DifferentialFilter_LHS_N()
133 for (int k = 0; k < 3; k++) { in DifferentialFilter_LHS_N() local
134 Grad_v[k * N + j][i] = wdetJ * dq_dXdx_a[k]; in DifferentialFilter_LHS_N()
|
| H A D | velocity_gradient_projection.h | 31 for (CeedInt k = 0; k < 3; k++) { in VelocityGradientProjectionRHS() local
32 v[j * 3 + k][i] = wdetJ * grad_velocity[j][k]; in VelocityGradientProjectionRHS()
|
| H A D | sgs_dd_model.h | 209 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = sgs[k][j]; in FluxSubgridStress() local
234 for (CeedInt k = 0; k < 3; k++) { in IFunction_NodalSgs() local
236 …Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1] + dXdx[k][2] * Flux[… in IFunction_NodalSgs()
|
| H A D | utils_eigensolver_jacobi.h | 254 for (CeedInt k = 0; k < N; k++) SwapScalar(&evec[i * N + k], &evec[i_max * N + k]); in SortRows() local
|
| H A D | channel.h | 34 const CeedScalar k = gas.k; in Exact_Channel() local
42 const CeedScalar Pr = mu / (cp * k); in Exact_Channel()
|
| H A D | newtonian_types.h | 23 CeedScalar k; member
|
| /honee/src/ |
| H A D | dm-utils.c | 347 for (PetscInt k = 0; k < num_points; k++) { in PetscDTUniformTensorQuadrature() local
348 coords[((i * num_points + j) * num_points + k) * dim + 0] = coords_1d[i]; in PetscDTUniformTensorQuadrature()
349 coords[((i * num_points + j) * num_points + k) * dim + 1] = coords_1d[j]; in PetscDTUniformTensorQuadrature()
350 coords[((i * num_points + j) * num_points + k) * dim + 2] = coords_1d[k]; in PetscDTUniformTensorQuadrature()
351 …for (PetscInt c = 0; c < num_comp; c++) weights[((i * num_points + j) * num_points + k) * num_comp… in PetscDTUniformTensorQuadrature()
|
| /honee/tests/junit-xml/junit_xml/ |
| H A D | __init__.py | 157 for k, v in self.properties.items():
158 attrs = {"name": decode(k, encoding), "value": decode(v, encoding)}
|
| /honee/problems/ |
| H A D | newtonian.c | 395 .k = 0.02638, in NS_NEWTONIAN_IG()
422 …alar("-k", "Thermal conductivity", NULL, newtonian_ig_ctx->gas.k, &newtonian_ig_ctx->gas.k, NULL)); in NS_NEWTONIAN_IG()
480 newtonian_ig_ctx->gas.k *= units->W_per_m_K; in NS_NEWTONIAN_IG()
|
| /honee/doc/ |
| H A D | theory.md | 15 … \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\sigma} - k \nabla T \right) - \…
20 … matrix, $\bm{b}$ is a body force vector (e.g., gravity vector $\bm{g}$), $k$ the thermal conduct…
63 - \bm{u} \cdot \bm{\sigma} - k \nabla T
79 \bm{q}_N (\bm{x},t)^{(e)} = \sum_{k=1}^{P}\psi_k (\bm{x})\bm{q}_k^{(e)}
595 kn[k^n];
641 | $\{\kappa^n\}_{n=1}^N$ | k^n | No | Yes |
|
| H A D | references.bib | 171 author = {Shakib, Farzin and Hughes, Thomas JR and Johan, Zden{\v{e}}k},
|
| H A D | runtime_options.md | 473 * - `-k`
|
| /honee/tests/ |
| H A D | junit_common.py | 370 for k, v in other.items():
545 for k, v in spec.key_values.items():
|