Home
last modified time | relevance | path

Searched refs:x (Results 1 – 25 of 26) sorted by relevance

12

/honee/qfunctions/
H A Dtaylorgreen.h36 CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsTaylorGreen() local
38 ScaleN(x, 2 * M_PI / context->lx, 3); in ICsTaylorGreen()
40 …ference.pressure + (density0 * Square(V0) / 16) * (cos(2 * x[0]) + cos(2 * x[1])) * (cos(2 * x[2] … 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()
H A Dblasius.h33 CEED_QFUNCTION_HELPER void ChebyshevEval(int N, const double *Tf, double x, double eta_max, double … in ChebyshevEval() argument
37 {1, x, 2 * x * x - 1}, in ChebyshevEval()
38 {0, 1, 4 * x }, in ChebyshevEval()
48 table[0][i % 3] = 2 * x * table[0][(i - 1) % 3] - table[0][(i - 2) % 3]; in ChebyshevEval()
67 State CEED_QFUNCTION_HELPER(BlasiusSolution)(const BlasiusContext blasius, const CeedScalar x[3], c… in BlasiusSolution()
74 CeedScalar eta = x[1] * sqrt(U_infty / (nu * (x0 + x[0] - x_inflow))); in BlasiusSolution()
82 *t12 = mu * U_infty * f[2] * sqrt(U_infty / (nu * (x0 + x[0] - x_inflow))); in BlasiusSolution()
86 Y[2] = 0.5 * sqrt(nu * U_infty / (x0 + x[0] - x_inflow)) * (eta * f[1] - f[0]); in BlasiusSolution()
113 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in ICsBlasius() local
114 State s = BlasiusSolution(context, x, x0, x_inflow, S_infty.U.density, &t12); in ICsBlasius()
[all …]
H A Dutils.h23 CEED_QFUNCTION_HELPER CeedScalar Square(CeedScalar x) { return x * x; } in Square() argument
24 CEED_QFUNCTION_HELPER CeedScalar Cube(CeedScalar x) { return x * x * x; } in Cube() argument
37 … void CopyN(const CeedScalar *x, CeedScalar *y, const CeedInt N) { CeedPragmaSIMD for (CeedInt i =… in CopyN() argument
50 CEED_QFUNCTION_HELPER void AXPY(CeedScalar alpha, const CeedScalar *x, CeedScalar *y, CeedInt N) { in AXPY() argument
51 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) y[i] += alpha * x[i]; in AXPY()
81 CEED_QFUNCTION_HELPER void MatVecNM(const CeedScalar *A, const CeedScalar *x, const CeedInt N, cons… in MatVecNM() argument
85 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M); in MatVecNM()
88 … i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[i] += A[j * M + i] * x[j]; } in MatVecNM()
94 CEED_QFUNCTION_HELPER void MatVec3(const CeedScalar A[3][3], const CeedScalar x[3], const CeedTrans… in MatVec3()
95 MatVecNM((const CeedScalar *)A, (const CeedScalar *)x, 3, 3, transpose_A, (CeedScalar *)b); in MatVec3()
[all …]
H A Dchannel.h40 const CeedScalar x[3] = {0, X[1], X[2]}; in Exact_Channel() local
44 const CeedScalar theta = theta0 * (1 + (Pr * Ec / 3) * (1 - Square(Square((x[1] - center) / H)))); in Exact_Channel()
47 Y[1] = umax * (1 - Square((x[1] - center) / H)); in Exact_Channel()
66 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsChannel() local
67 State s = Exact_Channel(3, 0., x, 5, ctx); in ICsChannel()
99 const CeedScalar x[3] = {0, X[1][i], X[2][i]}; in Channel_Inflow() local
102 State s_exact = Exact_Channel(3, 0., x, 5, ctx); in Channel_Inflow()
H A Ddifferential_filter.h78 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in DifferentialFilter_LHS_N() local
87 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in DifferentialFilter_LHS_N()
185 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in DifferentialFilter_MMS_IC() local
189 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in DifferentialFilter_MMS_IC()
H A Ddensitycurrent.h97 const CeedScalar x = X[0]; in Exact_DC() local
102 CeedScalar rr[3] = {x - center[0], y - center[1], z - center[2]}; in Exact_DC()
134 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsDC() local
135 State s = Exact_DC(3, 0., x, 5, ctx); in ICsDC()
H A Dgaussianwave.h39 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in IC_GaussianWave() local
40 const CeedScalar x0 = x[0] - xc; in IC_GaussianWave()
41 const CeedScalar y0 = x[1] - yc; in IC_GaussianWave()
H A Dstg_shur14.h244 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in StgShur14Preprocess() local
262 const CeedScalar wall_dist = x[1][i]; in StgShur14Preprocess()
284 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsStg() local
300 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in ICsStg()
357 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in StgShur14Inflow() local
370 StgShur14Calc(x, time, ubar, cij, qn, u, stg_ctx); in StgShur14Inflow()
491 const CeedScalar x[] = {coords[0][i], coords[1][i], coords[2][i]}; in StgShur14InflowStrongQF() local
499 …StgShur14Calc_PrecompEktot(x, time, ubar, cij, inv_Ektotal[i], h_node_sep, x[1], eps, lt, nu, u, s… in StgShur14InflowStrongQF()
503 StgShur14Calc(x, time, ubar, cij, qn, u, stg_ctx); in StgShur14InflowStrongQF()
H A Dadvection.h76 const CeedScalar x = X[0], y = X[1], z = dim == 2 ? 0. : X[2]; in Exact_AdvectionGeneric() local
82 q[2] = (x - center[0]); in Exact_AdvectionGeneric()
102 CeedScalar r = sqrt(Square(x - x0[0]) + Square(y - x0[1]) + Square(z - x0[2])); in Exact_AdvectionGeneric()
124 CeedScalar r = sqrt(Square(x - center[0]) + Square(y - center[1])); in Exact_AdvectionGeneric()
131 CeedScalar inflow_to_point[3] = {x - context->lx / 2, y, 0}; in Exact_AdvectionGeneric()
137 …if ((x < boundary_threshold && wind[0] < boundary_threshold) || // outflow at -x bo… in Exact_AdvectionGeneric()
139 …(x > context->lx - boundary_threshold && wind[0] > boundary_threshold) || // outflow at +x bounda… in Exact_AdvectionGeneric()
161 if ((x < boundary_threshold) || (y > ly - boundary_threshold)) { in Exact_AdvectionGeneric()
183 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsAdvection() local
186 Exact_AdvectionGeneric(3, 0., x, 5, q, ctx); in ICsAdvection()
[all …]
H A Dbc_outflow.h18 CEED_QFUNCTION_HELPER CeedScalar Softplus(CeedScalar x, CeedScalar width) { in Softplus() argument
19 if (x > 40 * width) return x; in Softplus()
20 return width * log1p(exp(x / width)); in Softplus()
23 CEED_QFUNCTION_HELPER CeedScalar Softplus_fwd(CeedScalar x, CeedScalar dx, CeedScalar width) { in Softplus_fwd() argument
24 if (x > 40 * width) return 1; in Softplus_fwd()
25 const CeedScalar t = exp(x / width); in Softplus_fwd()
H A Dshocktube.h83 const CeedScalar x = X[0]; // Coordinates in Exact_ShockTube() local
88 if (x <= mid_point + 200 * CEED_EPSILON) { in Exact_ShockTube()
175 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsShockTube() local
178 Exact_ShockTube(3, 0., x, 5, q, ctx); in ICsShockTube()
H A Dnewtonian.h194 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in RHSFunction_Newtonian() local
205 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in RHSFunction_Newtonian()
316 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[4]; in IFunction_Newtonian() local
327 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in IFunction_Newtonian()
H A Dsgs_dd_model.h33 CEED_QFUNCTION_HELPER void LeakyReLU(CeedScalar *x, const CeedScalar alpha, const CeedInt N) { in LeakyReLU() argument
34 for (CeedInt i = 0; i < N; i++) x[i] *= (x[i] < 0 ? alpha : 1.); in LeakyReLU()
H A Deulervortex.h69 const CeedScalar x = X[0], y = X[1]; // Coordinates in Exact_Euler() local
74 const CeedScalar x0 = x - xc; in Exact_Euler()
230 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsEuler() local
233 Exact_Euler(3, context->curr_time, x, 5, q, ctx); in ICsEuler()
/honee/include/
H A Dpetsc-ceed-utils.h126 PetscScalar *x; in VecPetscToCeed() local
129 PetscCall(VecGetArrayAndMemType(X_petsc, &x, mem_type)); in VecPetscToCeed()
130 …eturnCeed(x_ceed), CeedVectorSetArray(x_ceed, MemTypePetscToCeed(*mem_type), CEED_USE_POINTER, x)); in VecPetscToCeed()
146 PetscScalar *x; in VecCeedToPetsc() local
149 …lCeed(CeedVectorReturnCeed(x_ceed), CeedVectorTakeArray(x_ceed, MemTypePetscToCeed(mem_type), &x)); in VecCeedToPetsc()
150 PetscCall(VecRestoreArrayAndMemType(X_petsc, &x)); in VecCeedToPetsc()
166 PetscScalar *x; in VecReadPetscToCeed() local
169 PetscCall(VecGetArrayReadAndMemType(X_petsc, (const PetscScalar **)&x, mem_type)); in VecReadPetscToCeed()
170 …eturnCeed(x_ceed), CeedVectorSetArray(x_ceed, MemTypePetscToCeed(*mem_type), CEED_USE_POINTER, x)); in VecReadPetscToCeed()
186 PetscScalar *x; in VecReadCeedToPetsc() local
[all …]
H A Dnavierstokes.h361 PetscErrorCode KSPPostSolve_Honee(KSP ksp, Vec rhs, Vec x, void *ctx);
/honee/qfunctions/spanstats/
H A Dturbulence.h99 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in ChildStatsCollectionMMSTest() local
107 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in ChildStatsCollectionMMSTest()
119 const CeedScalar(*x)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in ChildStatsCollectionMMSTest_Error() local
124 const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]}; in ChildStatsCollectionMMSTest_Error()
/honee/doc/
H A Dexamples.md12 \rho &= \rho_\infty\left(1+A\exp\left(\frac{-(\bar{x}^2 + \bar{y}^2)}{2\sigma^2}\right)\right) \\
14 E &= \frac{p_\infty}{\gamma -1}\left(1+A\exp\left(\frac{-(\bar{x}^2 + \bar{y}^2)}{2\sigma^2}\right)…
18 …e amplitude and width of the perturbation, respectively, and $(\bar{x}, \bar{y}) = (x-x_e, y-y_e)$…
21 The boundary conditions are freestream in the x and y directions. When using an HLL (Harten, Lax, v…
63 A cylinder with diameter $D=1$ is centered at $(0,0)$ in a computational domain $-4.5 \leq x \leq 1…
119 …is defined in terms of the Exner pressure, $\pi(\bm{x},t)$, and potential temperature, $\theta(\bm…
122 …rac{P_0}{( c_p - c_v)\theta(\bm{x},t)} \pi(\bm{x},t)^{\frac{c_v}{ c_p - c_v}} \, , \\ e &= c_v \th…
246 Algorithmically, a base node distribution is defined at the inlet (assumed to be $\min(x)$) and the…
253 The node locations used exactly at the inlet (assumed to be $\min(x)$) and linearly stretched/squee…
385 u &= V'_x + V \sin(\hat x) \cos(\hat y) \sin(\hat z) \\
[all …]
H A Dauxiliary.md13 \langle \phi \rangle(x,y) = \frac{1}{L_z + (T_f - T_0)}\int_0^{L_z} \int_{T_0}^{T_f} \phi(x, y, z, …
22 The function $\langle \phi \rangle (x,y)$ is represented on a 2-D finite element grid, taken from t…
26 Define a function space on the parent grid as $\mathcal{V}_p^\mathrm{parent} = \{ \bm v(\bm x) \in …
34 \langle \phi \rangle_z(x,y,t) = \frac{1}{L_z} \int_0^{L_z} \phi(x, y, z, t) \mathrm{d}z
41 \bm M u_N = \int_0^{L_x} \int_0^{L_y} u \psi^\mathrm{parent}_N \mathrm{d}y \mathrm{d}x
47 …ft [\frac{1}{L_z} \int_0^{L_z} \phi(x,y,z,t) \mathrm{d}z \right ] \psi^\mathrm{parent}_N(x,y) \mat…
53 \bm M [\langle \phi \rangle_z]_N = \frac{1}{L_z} \int_\Omega \phi(x,y,z,t) \psi^\mathrm{parent}_N(x…
58 This assumption means quadrature points in the full domain have the same $(x,y)$ coordinate locatio…
73 \bm M [\langle \phi \rangle]_N = \frac{1}{L_z + (T_f - T_0)} \int_\Omega \int_{T_0}^{T_f} \phi(x,y,…
75 where the integral $\int_{T_0}^{T_f} \phi(x,y,z,t) \mathrm{d}t$ is calculated on a running basis.
[all …]
H A Dtheory.md79 \bm{q}_N (\bm{x},t)^{(e)} = \sum_{k=1}^{P}\psi_k (\bm{x})\bm{q}_k^{(e)}
91 with $\mathcal{V}_p = \{ \bm v(\bm x) \in H^{1}(\Omega_e) \,|\, \bm v(\bm x_e(\bm X)) \in P_p(\bm{I…
228 …can be pulled back to the reference element as $\bm u_{\bm X} = \nabla_{\bm x}\bm X \cdot \bm u$, …
229 …boundary layer element of dimension $(1, \epsilon)$, for which $\nabla_{\bm x} \bm X = \bigl(\begi…
232 …of a unit vector $\hat{\bm n}$ is given by $\lVert \bigl(\nabla_{\bm X} \bm x\bigr)^T \hat{\bm n} …
233 While $\nabla_{\bm X} \bm x$ is readily computable, its inverse $\nabla_{\bm x} \bm X$ is needed di…
284 where $\bm g = \nabla_{\bm x} \bm{X}^T \cdot \nabla_{\bm x} \bm{X}$ is the metric tensor and $\Vert…
316 …ector in direction $i$, and $\nabla_{x_i} = \hat{\bm n}_i \cdot \nabla_{\bm x}$ is the derivative …
515 \bm{u}(\bm{x}, t) = \bm{\overline{u}}(\bm{x}) + \bm{C}(\bm{x}) \cdot \bm{v}'
520 …= 2 \sqrt{3/2} \sum^N_{n=1} \sqrt{q^n(\bm{x})} \bm{\sigma}^n \cos(\kappa^n \bm{d}^n \cdot \bm{\hat…
[all …]
H A Druntime_options.md265 - Use symmetry boundary conditions, for the x component, on this list of faces
296 - +x
304 - -x
331 - +x
335 - -x
/honee/examples/postprocess/
H A Dvortexshedding.py32 x="Time",
38 x="Time",
/honee/tests/createPyTorchModel/
H A Dupdate_weights.py34 def forward(self, x): argument
35 return self.net(x)
/honee/src/
H A Dhonee-ksp.c9 PetscErrorCode KSPPostSolve_Honee(KSP ksp, Vec rhs, Vec x, void *ctx) { in KSPPostSolve_Honee() argument
H A Dmat-ceed.c1436 const PetscScalar *x; in MatCeedContextCreate() local
1440 PetscCall(VecGetArrayReadAndMemType(X, &x, &(*ctx)->mem_type)); in MatCeedContextCreate()
1441 PetscCall(VecRestoreArrayReadAndMemType(X, &x)); in MatCeedContextCreate()

12