Home
last modified time | relevance | path

Searched refs:X (Results 1 – 18 of 18) sorted by relevance

/honee/qfunctions/
H A Ddensitycurrent.h79 CEED_QFUNCTION_HELPER State Exact_DC(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf… in Exact_DC() argument
97 const CeedScalar x = X[0]; in Exact_DC()
98 const CeedScalar y = X[1]; in Exact_DC()
99 const CeedScalar z = X[2]; in Exact_DC()
127 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsDC() local
134 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsDC()
H A Dblasius.h23 …CeedScalar *X; // !< Chebyshev polynomial… member
75 CeedScalar X = 2 * (eta / blasius->eta_max) - 1.; in BlasiusSolution() local
79 ChebyshevEval(N, blasius->Tf_cheb, X, blasius->eta_max, f); in BlasiusSolution()
80 ChebyshevEval(N - 1, blasius->Th_cheb, X, blasius->eta_max, h); in BlasiusSolution()
97 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsBlasius() local
113 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in ICsBlasius()
128 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Blasius_Inflow() local
144 const CeedScalar x[3] = {X[0][i], X[1][i], 0.}; in Blasius_Inflow()
180 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Blasius_Inflow_Jacobian() local
199 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in Blasius_Inflow_Jacobian()
H A Dchannel.h24 CEED_QFUNCTION_HELPER State Exact_Channel(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedI… in Exact_Channel() argument
40 const CeedScalar x[3] = {0, X[1], X[2]}; in Exact_Channel()
59 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsChannel() local
66 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsChannel()
81 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Channel_Inflow() local
99 const CeedScalar x[3] = {0, X[1][i], X[2][i]}; in Channel_Inflow()
H A Dtaylorgreen.h24 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsTaylorGreen() local
36 CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsTaylorGreen()
H A Dgaussianwave.h21 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in IC_GaussianWave() local
39 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in IC_GaussianWave()
H A Dstg_shur14.h152 CEED_QFUNCTION_HELPER void StgShur14Calc(const CeedScalar X[3], const CeedScalar t, const CeedScala… in StgShur14Calc()
160 CeedScalar xhat[] = {0., X[1], X[2]}; in StgShur14Calc()
163 xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1); in StgShur14Calc()
193 CEED_QFUNCTION_HELPER void StgShur14Calc_PrecompEktot(const CeedScalar X[3], const CeedScalar t, co… in StgShur14Calc_PrecompEktot()
205 CeedScalar xhat[] = {0., X[1], X[2]}; in StgShur14Calc_PrecompEktot()
208 xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1); in StgShur14Calc_PrecompEktot()
335 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in StgShur14Inflow() local
357 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in StgShur14Inflow()
367 InterpolateProfile(X[1][i], ubar, cij, &eps, &lt, stg_ctx); in StgShur14Inflow()
369 CalcSpectrum(X[1][i], eps, lt, h_node_sep, mu / rho, qn, stg_ctx); in StgShur14Inflow()
H A Dadvection.h64 CEED_QFUNCTION_HELPER int Exact_AdvectionGeneric(CeedInt dim, CeedScalar time, const CeedScalar X[]… in Exact_AdvectionGeneric() argument
76 const CeedScalar x = X[0], y = X[1], z = dim == 2 ? 0. : X[2]; in Exact_AdvectionGeneric()
147 CeedScalar theta = context->wave_frequency * DotN(X, wind, dim) + context->wave_phase; in Exact_AdvectionGeneric()
179 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsAdvection() local
183 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsAdvection()
196 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsAdvection2d() local
201 const CeedScalar x[] = {X[0][i], X[1][i]}; in ICsAdvection2d()
H A Dshocktube.h72 CEED_QFUNCTION_HELPER CeedInt Exact_ShockTube(CeedInt dim, CeedScalar time, const CeedScalar X[], C… in Exact_ShockTube() argument
83 const CeedScalar x = X[0]; // Coordinates in Exact_ShockTube()
171 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsShockTube() local
175 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsShockTube()
H A Deulervortex.h58 CEED_QFUNCTION_HELPER int Exact_Euler(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt N… in Exact_Euler() argument
69 const CeedScalar x = X[0], y = X[1]; // Coordinates in Exact_Euler()
224 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsEuler() local
230 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsEuler()
H A Dnewtonian_state.h263 …eConservative StateConservativeMult(CeedInt n, const CeedScalar a[], const StateConservative X[]) { in StateConservativeMult() argument
266 R.density += a[i] * X[i].density; in StateConservativeMult()
267 for (int j = 0; j < 3; j++) R.momentum[j] += a[i] * X[i].momentum[j]; in StateConservativeMult()
268 R.E_total += a[i] * X[i].E_total; in StateConservativeMult()
273 CEED_QFUNCTION_HELPER StateConservative StateConservativeAXPBYPCZ(CeedScalar a, StateConservative X… in StateConservativeAXPBYPCZ() argument
276 R.density = a * X.density + b * Y.density + c * Z.density; in StateConservativeAXPBYPCZ()
277 …for (int i = 0; i < 3; i++) R.momentum[i] = a * X.momentum[i] + b * Y.momentum[i] + c * Z.momentum… in StateConservativeAXPBYPCZ()
278 R.E_total = a * X.E_total + b * Y.E_total + c * Z.E_total; in StateConservativeAXPBYPCZ()
/honee/problems/
H A Dblasius.c17 PetscErrorCode CompressibleBlasiusResidual(SNES snes, Vec X, Vec R, void *ctx) { in CompressibleBlasiusResidual() argument
29 PetscCall(VecGetArrayRead(X, &Tf)); in CompressibleBlasiusResidual()
43 ChebyshevEval(N, Tf, blasius->X[i], blasius->eta_max, f); in CompressibleBlasiusResidual()
44 ChebyshevEval(N - 1, Th, blasius->X[i], blasius->eta_max, h); in CompressibleBlasiusResidual()
68 PetscCall(VecRestoreArrayRead(X, &Tf)); in CompressibleBlasiusResidual()
83 PetscCall(PetscMalloc2(N - 3, &blasius->X, N - 3, &w)); in ComputeChebyshevCoefficients()
84 PetscCall(PetscDTGaussQuadrature(N - 3, -1., 1., blasius->X, w)); in ComputeChebyshevCoefficients()
107 PetscCall(PetscFree2(blasius->X, w)); in ComputeChebyshevCoefficients()
/honee/src/
H A Dpetsc_ops.c208 PetscErrorCode ApplyCeedOperator_Core(Vec X, Vec X_loc, CeedVector x_ceed, CeedVector y_ceed, Vec Y… in ApplyCeedOperator_Core() argument
214 if (X) PetscCall(DMGlobalToLocal(ctx->dm_x, X, INSERT_VALUES, X_loc)); in ApplyCeedOperator_Core()
219 PetscCall(PetscLogEventBegin(HONEE_CeedOperatorApply, X, Y, 0, 0)); in ApplyCeedOperator_Core()
224 PetscCall(PetscLogEventEnd(HONEE_CeedOperatorApply, X, Y, 0, 0)); in ApplyCeedOperator_Core()
233 PetscErrorCode ApplyCeedOperatorGlobalToGlobal(Vec X, Vec Y, OperatorApplyContext ctx) { in ApplyCeedOperatorGlobalToGlobal() argument
243 PetscCall(ApplyCeedOperator_Core(X, X_loc, ctx->x_ceed, ctx->y_ceed, Y_loc, Y, ctx, false)); in ApplyCeedOperatorGlobalToGlobal()
267 PetscErrorCode ApplyCeedOperatorGlobalToLocal(Vec X, Vec Y_loc, OperatorApplyContext ctx) { in ApplyCeedOperatorGlobalToLocal() argument
274 PetscCall(ApplyCeedOperator_Core(X, X_loc, ctx->x_ceed, ctx->y_ceed, Y_loc, NULL, ctx, false)); in ApplyCeedOperatorGlobalToLocal()
H A Dmat-ceed.c380 Vec X; in MatCreateCeed() local
382 PetscCall(DMGetGlobalVector(dm_x, &X)); in MatCreateCeed()
383 PetscCall(VecGetSize(X, &X_g_size)); in MatCreateCeed()
384 PetscCall(VecGetLocalSize(X, &X_l_size)); in MatCreateCeed()
385 PetscCall(DMRestoreGlobalVector(dm_x, &X)); in MatCreateCeed()
1437 Vec X; in MatCeedContextCreate() local
1439 PetscCall(DMGetLocalVector(dm_x, &X)); in MatCeedContextCreate()
1440 PetscCall(VecGetArrayReadAndMemType(X, &x, &(*ctx)->mem_type)); in MatCeedContextCreate()
1441 PetscCall(VecRestoreArrayReadAndMemType(X, &x)); in MatCeedContextCreate()
1442 PetscCall(DMRestoreLocalVector(dm_x, &X)); in MatCeedContextCreate()
[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);
29 PetscErrorCode ApplyCeedOperatorGlobalToLocal(Vec X, Vec Y_loc, OperatorApplyContext ctx);
/honee/doc/
H A Dtheory.md91 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 …m u$ can be pulled back to the reference element as $\bm u_{\bm X} = \nabla_{\bm x}\bm X \cdot \bm…
229 …y layer element of dimension $(1, \epsilon)$, for which $\nabla_{\bm x} \bm X = \bigl(\begin{small…
231 The ratio $\lVert \bm u \rVert / \lVert \bm u_{\bm X} \rVert$ is a covariant measure of (half) the …
232 …ection of a unit vector $\hat{\bm n}$ is given by $\lVert \bigl(\nabla_{\bm X} \bm x\bigr)^T \hat{…
233 While $\nabla_{\bm X} \bm x$ is readily computable, its inverse $\nabla_{\bm x} \bm X$ is needed di…
240 \mathrm{Pe} = \frac{\lVert \bm u \rVert^2}{\lVert \bm u_{\bm X} \rVert \kappa}.
246 \tau = \frac{\xi(\mathrm{Pe})}{\lVert \bm u_{\bm X} \rVert},
254 …m u \tau \bm u \cdot \nabla q = \nabla_{\bm X} v \cdot (\bm u_{\bm X} \tau \bm u_{\bm X}) \cdot \n…
284 where $\bm g = \nabla_{\bm x} \bm{X}^T \cdot \nabla_{\bm x} \bm{X}$ is the metric tensor and $\Vert…
[all …]
H A Dreferences.bib13 author = {Giraldo, F. X. and Restelli, M. and Läuter, M.},
170 …title = {A new finite element formulation for computational fluid dynamics: {X}. The compressi…
H A Dauxiliary.md26 …\{ \bm v(\bm x) \in H^{1}(\Omega_e^\mathrm{parent}) \,|\, \bm v(\bm x_e(\bm X)) \in P_p(\bm{I}), e…
209 …bm{\Delta} = \bm{g}^{-1/2}$ where $\bm g = \nabla_{\bm x} \bm{X} \cdot \nabla_{\bm x} \bm{X}$ is t…