Home
last modified time | relevance | path

Searched full:s (Results 1 – 25 of 255) sorted by relevance

1234567891011

/libCEED/tests/
H A Dt541-operator.c90 CeedScalar S[16]; in main() local
117 S[0 * 4 + i] = -v_array[0]; in main()
118 S[1 * 4 + i] = -v_array[p - 1]; in main()
119 S[2 * 4 + i] = -v_array[p * p - p]; in main()
120 S[3 * 4 + i] = -v_array[p * p - 1]; in main()
126 …S_inv[0] = S[5] * S[10] * S[15] - S[5] * S[11] * S[14] - S[9] * S[6] * S[15] + S[9] * S[7] * S[14]… in main()
128 …S_inv[4] = -S[4] * S[10] * S[15] + S[4] * S[11] * S[14] + S[8] * S[6] * S[15] - S[8] * S[7] * S[14… in main()
130 …S_inv[8] = S[4] * S[9] * S[15] - S[4] * S[11] * S[13] - S[8] * S[5] * S[15] + S[8] * S[7] * S[13] … in main()
132 …_inv[12] = -S[4] * S[9] * S[14] + S[4] * S[10] * S[13] + S[8] * S[5] * S[14] - S[8] * S[6] * S[13]… in main()
134 …S_inv[1] = -S[1] * S[10] * S[15] + S[1] * S[11] * S[14] + S[9] * S[2] * S[15] - S[9] * S[3] * S[14… in main()
[all …]
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
H A Dt407-qfunction.c46 …if (strcmp(name, "time")) printf("Incorrect context field description for time: \"%s\" != \"time\"… in main()
50 printf("Incorrect context field type for time: \"%s\" != \"%s\"\n", CeedContextFieldTypes[type], in main()
56 …if (strcmp(name, "count")) printf("Incorrect context field description for count: \"%s\" != \"coun… in main()
60 … printf("Incorrect context field type for count: \"%s\" != \"%s\"\n", CeedContextFieldTypes[type], in main()
66 …if (strcmp(name, "is set")) printf("Incorrect context field description for count: \"%s\" != \"cou… in main()
70 … printf("Incorrect context field type for is_set: \"%s\" != \"%s\"\n", CeedContextFieldTypes[type], in main()
H A Dt005-ceed.c18 …if (strcmp(err_msg, "No error message stored")) printf("Unexpected error message received: \"%s\"\… in main()
30 …_msg, "No error message stored\n")) printf("Unexpected error message received: \"%s\"\n", err_msg); in main()
35 …if (strcmp(err_msg, "No error message stored")) printf("Unexpected error message received: \"%s\"\… in main()
/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 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 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.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 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 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 …]
/libCEED/
H A Dcommon.mk12 ansicolor = $(shell echo $(call last_path,$(1)) | cksum | cut -b1-2 | xargs -IS expr 2 \* S + 17)
13 emacs_out = @printf " %10s %s/%s\n" $(1) $(call rule_path,$(2)) $(call rule_file,$(2))
15 printf " %10s \033[38;5;%d;1m%s\033[m/%s\n" \
18 printf " %10s %s\n" $(1) $(2); fi
H A D.gitlab-ci.yml124 …e --output-file coverage.info --ignore-errors source,mismatch --substitute 's#(t*-f.h)#test/(t*.-f…
125 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F interface;
126 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F gallery;
127 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F backends;
128 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F tests;
129 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F examples;
212 …e --output-file coverage.info --ignore-errors source,mismatch --substitute 's#(t*-f.h)#test/(t*.-f…
213 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F interface;
214 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F gallery;
215 … bash <(curl -s https://codecov.io/bash) -f coverage.info -t ${CODECOV_ACCESS_TOKEN} -F backends;
[all …]
/libCEED/rust/libceed-sys/c-src/
H A Dcommon.mk12 ansicolor = $(shell echo $(call last_path,$(1)) | cksum | cut -b1-2 | xargs -IS expr 2 \* S + 17)
13 emacs_out = @printf " %10s %s/%s\n" $(1) $(call rule_path,$(2)) $(call rule_file,$(2))
15 printf " %10s \033[38;5;%d;1m%s\033[m/%s\n" \
18 printf " %10s %s\n" $(1) $(2); fi
/libCEED/benchmarks/
H A Dpostprocess_plot.py19 draw_iter_lines = 0 # draw the "iter/s" lines?
20 ymin_iter_lines = 3e5 # minimal y value for the "iter/s" lines
21 ymax_iter_lines = 8e8 # maximal y value for the "iter/s" lines
88 print('backend: %s, compute nodes: %i, number of MPI tasks = %i' % (
140 plot(d[:, 0], d[:, 2], 's--', color=colors[i],
143 plot(d[:, 0], d[:, 1], 's--', color=colors[i])
152 plot(y / slope1, y, 'k--', label='%g iter/s' % (slope1 / vdim))
153 plot(y / slope2, y, 'k-', label='%g iter/s' % (slope2 / vdim))
156 title(r'%i node%s $\times$ %i ranks, %s, %s, %s' % (
157 num_nodes, '' if num_nodes == 1 else 's',
[all …]
/libCEED/examples/solids/
H A Dindex.md16 …are stress and strain, respectively, in the small strain regime, while $\bm S$ and $\bm E$ are the…
20 {\overbrace{\bm S(\bm E)}^{\text{Finite Strain Hyperelastic}}}
22 {\overbrace{\bm S = \mathsf C \bm E}^{\text{St. Venant-Kirchoff}}} \\
23 …@V{\text{geometric}}V{\begin{smallmatrix}\bm E \to \bm \epsilon \\ \bm S \to \bm \sigma \end{small…
24 …@V{\begin{smallmatrix}\bm E \to \bm \epsilon \\ \bm S \to \bm \sigma \end{smallmatrix}}V{\text{geo…
245 \bm{P} = \bm{F} \, \bm{S},
248 where $\bm S$ is the *second Piola-Kirchhoff stress* tensor, a symmetric tensor defined entirely in…
249 Different constitutive models can define $\bm S$.
266 The constitutive models considered, appropriate for large deformations, express $\bm S$ as a functi…
294 \bm S = \lambda \log J \bm C^{-1} + \mu (\bm I_3 - \bm C^{-1}).
[all …]
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/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/backends/cuda/
H A Dceed-cuda-compile.cpp48 CeedDebug(ceed, "Running command:\n$ %s", command); in CeedCallSystem_Core()
51 …CeedCheck(output_stream != nullptr, ceed, CEED_ERROR_BACKEND, "Failed to %s\ncommand:\n$ %s", mess… in CeedCallSystem_Core()
59 CeedDebug(ceed, "output:\n%s\n", output.c_str()); in CeedCallSystem_Core()
60 …se(output_stream) == 0, ceed, CEED_ERROR_BACKEND, "Failed to %s\ncommand:\n$ %s\nerror:\n%s", mess… in CeedCallSystem_Core()
162 CeedDebug(ceed, "Source:\n%s\n", code.str().c_str()); in CeedCompileCore_Cuda()
172 CeedDebug(ceed, "Option %d: %s", i, opts[i]); in CeedCompileCore_Cuda()
197 return CeedError(ceed, CEED_ERROR_BACKEND, "%s\n%s", nvrtcGetErrorString(result), log); in CeedCompileCore_Cuda()
201 CeedDebug(ceed, "Error: %s\nCompile log:\n%s\n", nvrtcGetErrorString(result), log); in CeedCompileCore_Cuda()
260 …CeedDebug(ceed, "There are %d source dirs, including %s\n", num_rust_source_dirs, rust_source_dirs… in CeedCompileCore_Cuda()
276 // Compile Rust crate(s) needed in CeedCompileCore_Cuda()
[all …]
/libCEED/examples/python/
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/interface/
H A Dceed-jit-tools.c44 CeedDebug(ceed, "%s\n", source_file_path_only); in CeedCheckFilePath()
54 CeedDebug(ceed, "%s\n", source_file_path_only); in CeedCheckFilePath()
99 …normalized_source_file_path, ceed, CEED_ERROR_MAJOR, "Malformed source path %s", source_file_path); in CeedNormalizePath()
132 CeedDebug(ceed, "%s\n", source_file_path); in CeedLoadSourceToInitializedBuffer()
136 CeedCheck(source_file, ceed, CEED_ERROR_MAJOR, "Couldn't open source file: %s", source_file_path); in CeedLoadSourceToInitializedBuffer()
154 return CeedError(ceed, CEED_ERROR_MAJOR, "Couldn't read source file: %s", source_file_path); in CeedLoadSourceToInitializedBuffer()
256 …CeedDebug256(ceed, CEED_DEBUG_COLOR_SUCCESS, "JiT Including: %s\n", normalized_include_source_path… in CeedLoadSourceToInitializedBuffer()
292 CeedDebug(ceed, "%s\n", source_file_path); in CeedLoadSourceToInitializedBuffer()
294 CeedDebug(ceed, "%s\n", *buffer); in CeedLoadSourceToInitializedBuffer()
391 …ERROR_MAJOR, "Couldn't find relative path including 'ceed/jit-source' for: %s", absolute_file_path… in CeedGetJitRelativePath()
[all …]
/libCEED/python/
H A Dceed_vector.py72 # Set Vector's data array
104 # Get Vector's data array
119 # Setup the pointer's pointer
149 # Get Vector's data array in read-only mode
164 # Setup the pointer's pointer
196 # Get Vector's data array in write-only mode
212 # Setup the pointer's pointer
242 # Restore the Vector's data array
246 # Setup the pointer's pointer
257 # Setup the pointer's pointer

1234567891011