Home
last modified time | relevance | path

Searched refs:s (Results 1 – 25 of 81) sorted by relevance

1234

/libCEED/julia/LibCEED.jl/.style/
H A Dceed_style.jl11 @inline JuliaFormatter.getstyle(s::CeedStyle) = s
16 s::State;
21 t = FST(Binary, cst, nspaces(s))
28 !s.opts.whitespace_typedefs
33 nospace_args = s.opts.whitespace_ops_in_indices ? false : nospace
36 n = pretty(style, cst[1], s, nonest=nonest, nospace=nospace_args)
38 n = pretty(style, cst[1], s)
42 s.opts.whitespace_ops_in_indices &&
46 add_node!(t, paren, s)
47 add_node!(t, n, s, join_lines=true)
[all …]
/libCEED/examples/fluids/qfunctions/
H A Dturb_spanstats.h26 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection() local
28 v[TURB_MEAN_DENSITY][i] = wdetJ * s.U.density; in ChildStatsCollection()
29 v[TURB_MEAN_PRESSURE][i] = wdetJ * s.Y.pressure; in ChildStatsCollection()
30 v[TURB_MEAN_PRESSURE_SQUARED][i] = wdetJ * Square(s.Y.pressure); in ChildStatsCollection()
31 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
32 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
33 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
34 v[TURB_MEAN_DENSITY_TEMPERATURE][i] = wdetJ * s.U.density * s.Y.temperature; in ChildStatsCollection()
35 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
36 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
[all …]
H A Dnewtonian_state.h31 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U() argument
32 U[0] = s.density; in UnpackState_U()
33 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
34 U[4] = s.E_total; in UnpackState_U()
37 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y() argument
38 Y[0] = s.pressure; in UnpackState_Y()
39 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
40 Y[4] = s.temperature; in UnpackState_Y()
43 CEED_QFUNCTION_HELPER void UnpackState_V(StateEntropy s, CeedScalar V[5]) { in UnpackState_V() argument
44 V[0] = s.S_density; in UnpackState_V()
[all …]
H A Dstabilization.h28 CEED_QFUNCTION_HELPER void StabilizationMatrix(NewtonianIdealGasContext gas, State s, CeedScalar Ta… in StabilizationMatrix() argument
37 State ds = StateFromY_fwd(gas, s, dY); in StabilizationMatrix()
38 FluxInviscid_fwd(gas, s, ds, dF); in StabilizationMatrix()
46 CEED_QFUNCTION_HELPER void Stabilization(NewtonianIdealGasContext gas, State s, CeedScalar Tau_d[3]… in Stabilization() argument
54 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
57 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
61 StabilizationMatrix(gas, s, Tau_d, strong_residual, stab); in Stabilization()
71 CEED_QFUNCTION_HELPER void Tau_diagPrim(NewtonianIdealGasContext gas, State s, const CeedScalar dXd… in Tau_diagPrim() argument
81 const CeedScalar rho = s.U.density; in Tau_diagPrim()
94 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
[all …]
H A Dnewtonian.h21 …id InternalDampingLayer(const NewtonianIdealGasContext context, const State s, const CeedScalar si… in InternalDampingLayer() argument
24 State damp_s = StateFromY_fwd(context, s, damp_Y); in InternalDampingLayer()
41 State s = StateFromPrimitive(&context->gas, context->reference); in ICsNewtonianIG() local
42 StateToQ(&context->gas, s, q, state_var); in ICsNewtonianIG()
73 const State s = StateFromQ(context, qi, state_var); in MassFunction_Newtonian() local
87 Tau_diagPrim(context, s, dXdx, context->dt, Tau_d); in MassFunction_Newtonian()
88 Stabilization(context, s, Tau_d, grad_s, U_dot, body_force, stab); in MassFunction_Newtonian()
167 State s = StateFromU(context, U); in RHSFunction_Newtonian() local
170 … StatePhysicalGradientFromReference(Q, i, context, s, STATEVAR_CONSERVATIVE, Grad_q, dXdx, grad_s); in RHSFunction_Newtonian()
176 ViscousEnergyFlux(context, s.Y, grad_s, stress, Fe); in RHSFunction_Newtonian()
[all …]
H A Dbc_slip.h24 State s = StateFromQ(newt_ctx, qi, state_var); in Slip() local
31 const CeedScalar vel_normal = Dot3(s.Y.velocity, norm); in Slip()
32 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal; in Slip()
33 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip()
36 StateConservative flux = RiemannFlux_HLLC(newt_ctx, s, s_reflect, norm); in Slip()
80 State s = StateFromQ(newt_ctx, qi, state_var); in Slip_Jacobian() local
81 State ds = StateFromQ_fwd(newt_ctx, s, dqi, state_var); in Slip_Jacobian()
84 const CeedScalar vel_normal = Dot3(s.Y.velocity, norm); in Slip_Jacobian()
85 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal; in Slip_Jacobian()
86 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip_Jacobian()
[all …]
H A Ddifferential_filter.h45 const State s = StateFromQ(gas, qi, state_var); in DifferentialFilter_RHS() local
47 v0[DIFF_FILTER_PRESSURE][i] = wdetJ * s.Y.pressure; in DifferentialFilter_RHS()
48 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
51 v0[DIFF_FILTER_TEMPERATURE][i] = wdetJ * s.Y.temperature; in DifferentialFilter_RHS()
52 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
53 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
54 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
55 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
[all …]
H A Dblasius.h121 State s = BlasiusSolution(context, x, x0, x_inflow, S_infty.U.density, &t12); in ICsBlasius() local
124 StateToQ(gas, s, q, gas->state_var); in ICsBlasius()
155 State s = BlasiusSolution(context, x, x0, context->x_inflow, rho_0, &t12); in Blasius_Inflow() local
162 s.U.density = s_int.U.density; in Blasius_Inflow()
163 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
165 s.U.E_total = s_int.U.E_total; in Blasius_Inflow()
166 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
170 FluxInviscid(&context->newtonian_ctx, s, Flux_inviscid); in Blasius_Inflow()
211 State s = BlasiusSolution(context, x, x0, 0, rho_0, &t12); in Blasius_Inflow_Jacobian() local
219 CeedScalar dE_kinetic = .5 * drho * Dot3(s.Y.velocity, s.Y.velocity); in Blasius_Inflow_Jacobian()
[all …]
H A Dadvection.h227 …calGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, NewtonianIdealGasContext gas, State s, in StatePhysicalGradientFromReference_ND() argument
237 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
250 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
257 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd… in Tau() argument
262 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
269 MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj); in Tau()
270 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * cont… in Tau()
298 const State s = StateFromU(gas, qi); in IFunction_AdvectionGeneric() local
303 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
318 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
[all …]
H A Dbc_freestream.h35 const State s = StateFromQ(newt_ctx, qi, state_var); in Freestream() local
44 flux = RiemannFlux_HLL(newt_ctx, s, context->S_infty, normal); in Freestream()
47 flux = RiemannFlux_HLLC(newt_ctx, s, context->S_infty, normal); in Freestream()
108 State s = StateFromQ(newt_ctx, qi, state_var); in Freestream_Jacobian() local
109 State ds = StateFromQ_fwd(newt_ctx, s, dqi, state_var); in Freestream_Jacobian()
114 dflux = RiemannFlux_HLL_fwd(newt_ctx, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
117 dflux = RiemannFlux_HLLC_fwd(newt_ctx, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
338 State s = StateFromQ(gas, qi, state_var); in PressureOutflow() local
339 s.Y.pressure = outflow->pressure; in PressureOutflow()
346 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in PressureOutflow()
[all …]
H A Dtaylorgreen.h40 State s = StateFromY(gas, Y); in ICsTaylorGreen() local
41 StateToQ(gas, s, q, gas->state_var); in ICsTaylorGreen()
H A Dvelocity_gradient_projection.h27 const State s = StateFromQ(context, qi, state_var); in VelocityGradientProjectionRHS() local
29 …StatePhysicalGradientFromReference(Q, i, context, s, state_var, (CeedScalar *)Grad_q, dXdx, grad_s… in VelocityGradientProjectionRHS()
H A Ddensitycurrent.h143 State s = Exact_DC(3, 0., x, 5, ctx); in ICsDC() local
145 StateToQ(gas, s, q, gas->state_var); in ICsDC()
/libCEED/tests/
H A Dt412-qfunction-f.f909 integer q,s,i variable
11 parameter(s=3)
12 real*8 uu(q*s)
13 real*8 vv(q*s)
20 call ceedqfunctioncreateidentity(ceed,s,ceed_eval_interp,&
23 do i=0,q*s-1
27 call ceedvectorcreate(ceed,q*s,u,err)
30 call ceedvectorcreate(ceed,q*s,v,err)
41 do i=1,q*s
/libCEED/julia/LibCEED.jl/examples/
H A Dcommon.jl6 s = 0 # find s: num_elem/2 < 2^s <= num_elem
9 s += 1
11 r = s%dim
13 sd = div(s, dim)
/libCEED/examples/solids/src/
H A Dboundary.c59 PetscScalar c = cos(theta), s = sin(theta); in BCClamp() local
61 …u[0] = lx + s * (-kz * y + ky * z) + (1 - c) * (-(ky * ky + kz * kz) * x + kx * ky * y + kx * kz *… in BCClamp()
62 …u[1] = ly + s * (kz * x + -kx * z) + (1 - c) * (kx * ky * x + -(kx * kx + kz * kz) * y + ky * kz *… in BCClamp()
63 …u[2] = lz + s * (-ky * x + kx * y) + (1 - c) * (kx * kz * x + ky * kz * y + -(kx * kx + ky * ky) *… in BCClamp()
/libCEED/examples/rust/mesh/src/
H A Dlib.rs18 let mut s = 0; // find s: num_elem / 2 < 2^s <= num_elem localVariable
21 s += 1;
25 let mut r = s % dim;
28 let mut sd = s / dim;
/libCEED/examples/python/
H A Dex_common.py74 s = 0
77 s += 1
80 r = s % dim
83 sd = s // dim
H A Dtutorial-6-shell.ipynb11 …"While libCEED's focus is on high-order finite/spectral element method implementations, the approa…
48 "$1$s. It computes:\n",
135 "Then move to the standalone libCEED's examples folder"
185 " Approx. # unknowns [-s] : 262144\n",
207 … or 3D domain Ω respectively, by applying the mass operator to a vector of 1s. The command line op…
231 "Given Laplace's equation,\n",
276 " Approx. # unknowns [-s] : 262144\n",
298 …lace's (diffusion) operator. The command line option `-d` specifies the dimensionality of the doma…
/libCEED/examples/ceed/
H A Dex1-volume-f.f9050 integer s, r, d, sd local
52 s = 0
60 s = s + 1
62 r = mod(s, fe_dim)
65 sd = s/fe_dim
/libCEED/examples/solids/
H A DREADME.md40 - [Young's modulus](https://en.wikipedia.org/wiki/Young%27s_modulus), $E > 0$
42 - [Poisson's ratio](https://en.wikipedia.org/wiki/Poisson%27s_ratio), $\nu < 0.5$
52 This solver can use any mesh format that PETSc's `DMPlex` can read (Exodus, Gmsh, Med, etc.).
123 - Poisson's ratio for multigrid smoothers, $\nu < 0.5$
169 …on multigrid coarsening to linear elements, with algebraic multigrid (PETSc's `GAMG`) for the coar…
178 Quantities such as the Young's modulus vary over many orders of magnitude, and thus can lead to poo…
212 * - Young's modulus, $E$
213 - $10^{11} \,\mathrm{Pa} = 10^{11} \,\mathrm{kg}\, \mathrm{m}^{-1}\, \mathrm s^{-2}$
216 …- $5 \cdot 10^4 \,\mathrm{kg}\, \mathrm m^{-2} \, \mathrm s^{-2} \cdot (\text{volume} \, \mathrm m…
/libCEED/doc/
H A DREADME.md16 Sphinx is the tool used for libCEED's User Manual. Sphinx can produce documentation in different ou…
18 To be able to contribute to libCEED's User Manual, Sphinx needs to be [installed](http://www.sphinx…
20 The Sphinx API documentation depends on Doxygen's XML output (via the `breathe` plugin). Build the…
/libCEED/doc/sphinx/source/
H A Dindex.md1 # Welcome to libCEED's User Manual!
/libCEED/doc/sphinx/source/api/
H A DCeedVector.rst6 …ctor` constitutes the main data structure and serves as input/output for the :ref:`CeedOperator`\s.
/libCEED/backends/ref/
H A Dceed-ref-operator.c1425 for (CeedInt s = 0; s < e_vec_size; s++) { in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref() local
1431 array[s] = 1.0; in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1432 if (s > 0) array[s - 1] = 0.0; in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1489 current_value = array[s]; in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1493 array[s] = current_value; in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1505 if (s == e_vec_size - 1) { in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1509 array[s] = 0.0; in CeedOperatorLinearAssembleAddDiagonalAtPoints_Ref()
1632 for (CeedInt s = 0; s < e_vec_size; s++) { in CeedOperatorAssembleSingleAtPoints_Ref() local
1633 const CeedInt comp_in = s / elem_size_active; in CeedOperatorAssembleSingleAtPoints_Ref()
1634 const CeedInt node_in = s % elem_size_active; in CeedOperatorAssembleSingleAtPoints_Ref()
[all …]

1234