Home
last modified time | relevance | path

Searched refs:u (Results 1 – 25 of 289) sorted by relevance

12345678910>>...12

/libCEED/examples/fluids/qfunctions/
H A Deulervortex.h87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
98 u[0] = mean_velocity[0] - C * y0; in Exact_Euler()
99 u[1] = mean_velocity[1] + C * x0; in Exact_Euler()
103 q[1] = rho * u[0]; in Exact_Euler()
104 q[2] = rho * u[1]; in Exact_Euler()
105 q[3] = rho * u[2]; in Exact_Euler()
106 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
114 q[1] = rho * u[0]; in Exact_Euler()
115 q[2] = rho * u[1]; in Exact_Euler()
116 q[3] = rho * u[2]; in Exact_Euler()
[all …]
H A Dshocktube.h90 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube() local
103 q[1] = rho * u[0]; in Exact_ShockTube()
104 q[2] = rho * u[1]; in Exact_ShockTube()
105 q[3] = rho * u[2]; in Exact_ShockTube()
106 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
114 …cobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar… in ConvectiveFluxJacobian_Euler()
116 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
119 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
122 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
123 … = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.) * u[i] * u[k]; in ConvectiveFluxJacobian_Euler()
[all …]
H A Dutils.h31 CEED_QFUNCTION_HELPER void ScaleN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in ScaleN() argument
32 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] *= alpha; in ScaleN()
36 CEED_QFUNCTION_HELPER void SetValueN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in SetValueN() argument
37 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] = alpha; in SetValueN()
47 CEED_QFUNCTION_HELPER CeedScalar DotN(const CeedScalar *u, const CeedScalar *v, const CeedInt N) { in DotN() argument
49 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i]; in DotN()
54 …ION_HELPER CeedScalar Dot3(const CeedScalar *u, const CeedScalar *v) { return u[0] * v[0] + u[1] *… in Dot3() argument
57 CEED_QFUNCTION_HELPER void Cross3(const CeedScalar u[3], const CeedScalar v[3], CeedScalar w[3]) { in Cross3()
58 w[0] = (u[1] * v[2]) - (u[2] * v[1]); in Cross3()
59 w[1] = (u[2] * v[0]) - (u[0] * v[2]); in Cross3()
[all …]
/libCEED/examples/petsc/qfunctions/bps/
H A Dcommon.h15 const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; in Error() local
18 error[i] = (u[i] - target[i]) * (u[i] - target[i]) * q_data[i]; in Error()
25 const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; in Error3() local
28 …error[i + 0 * Q] = (u[i + 0 * Q] - target[i + 0 * Q]) * (u[i + 0 * Q] - target[i + 0 * Q]) * q_dat… in Error3()
29 …error[i + 1 * Q] = (u[i + 1 * Q] - target[i + 1 * Q]) * (u[i + 1 * Q] - target[i + 1 * Q]) * q_dat… in Error3()
30 …error[i + 2 * Q] = (u[i + 2 * Q] - target[i + 2 * Q]) * (u[i + 2 * Q] - target[i + 2 * Q]) * q_dat… in Error3()
/libCEED/backends/magma/tuning/
H A Dtuning.cpp52 CeedVector u, v; in main() local
67 CeedVectorCreate(ceed, P * N, &u); in main()
71 CeedVectorSetValue(u, 1.0); in main()
73 ierr = CeedBasisApply(basis, N, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
77 CeedBasisApply(basis, N, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
83 CeedVectorSetValue(u, 1.0); in main()
85 ierr = CeedBasisApply(basis, N, CEED_TRANSPOSE, CEED_EVAL_INTERP, v, u); in main()
89 CeedBasisApply(basis, N, CEED_TRANSPOSE, CEED_EVAL_INTERP, v, u); in main()
94 CeedVectorDestroy(&u); in main()
100 CeedVectorCreate(ceed, P * N, &u); in main()
[all …]
/libCEED/examples/
H A Dbps.md14 …lem, posed as a weak form on a Hilbert space $V^p \subset H^1$, i.e., find $u \in V^p$ such that f…
17 \langle v,u \rangle = \langle v,f \rangle ,
20 …v,u\rangle$ and $\langle v,f\rangle$ express the continuous bilinear and linear forms, respectivel…
23 \begin{aligned} \langle v,u \rangle &:= \int_{\Omega} \, v \, u \, dV ,\\ \langle v,f \rangle &:= \…
30 u(\bm x) &= \sum_{j=1}^n u_j \, \phi_j(\bm x) ,\\
35 The coefficients $\{u_j\}$ and $\{v_i\}$ are the nodal values of $u$ and $v$, respectively.
39 \langle v,u \rangle = \bm v^T M \bm u , \qquad \langle v,f\rangle = \bm v^T \bm b \,.
54 …n BP3-BP6 is defined via the following variational formulation, i.e., find $u \in V^p$ such that f…
57 a(v,u) = \langle v,f \rangle , \,
60 where now $a (v,u)$ expresses the continuous bilinear form defined on $V^p$ for sufficiently regula…
[all …]
/libCEED/tests/
H A Dt332-basis.c12 CeedVector u, v; in main() local
33 CeedVectorCreate(ceed, p, &u); in main()
34 CeedVectorSetValue(u, 1.0); in main()
38 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_DIV, u, v); in main()
51 CeedVectorSetValue(u, 0.0); in main()
53 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_DIV, v, u); in main()
58 CeedVectorGetArrayRead(u, CEED_MEM_HOST, &u_array); in main()
62 CeedVectorRestoreArrayRead(u, &u_array); in main()
66 CeedVectorDestroy(&u); in main()
H A Dt331-basis.c12 CeedVector u, v; in main() local
33 CeedVectorCreate(ceed, p, &u); in main()
34 CeedVectorSetValue(u, 1.0); in main()
38 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
51 CeedVectorSetValue(u, 0.0); in main()
53 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_INTERP, v, u); in main()
58 CeedVectorGetArrayRead(u, CEED_MEM_HOST, &u_array); in main()
64 CeedVectorRestoreArrayRead(u, &u_array); in main()
68 CeedVectorDestroy(&u); in main()
H A Dt342-basis.c11 CeedVector u, v; in main() local
46 CeedVectorCreate(ceed, p, &u); in main()
47 CeedVectorSetValue(u, 1.0); in main()
51 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_CURL, u, v); in main()
64 CeedVectorSetValue(u, 0.0); in main()
66 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_CURL, v, u); in main()
71 CeedVectorGetArrayRead(u, CEED_MEM_HOST, &u_array); in main()
75 CeedVectorRestoreArrayRead(u, &u_array); in main()
79 CeedVectorDestroy(&u); in main()
H A Dt341-basis.c11 CeedVector u, v; in main() local
46 CeedVectorCreate(ceed, p, &u); in main()
47 CeedVectorSetValue(u, 1.0); in main()
51 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
64 CeedVectorSetValue(u, 0.0); in main()
66 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_INTERP, v, u); in main()
71 CeedVectorGetArrayRead(u, CEED_MEM_HOST, &u_array); in main()
75 CeedVectorRestoreArrayRead(u, &u_array); in main()
79 CeedVectorDestroy(&u); in main()
H A Dt363-basis.c12 CeedVector u, v; in main() local
21 CeedVectorCreate(ceed, q * dim, &u); in main()
22 CeedVectorSetValue(u, 1); in main()
29 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_GRAD, u, v); in main()
30 CeedBasisApplyAdd(basis, 1, CEED_TRANSPOSE, CEED_EVAL_GRAD, u, v); in main()
49 CeedVectorDestroy(&u); in main()
H A Dt310-basis.c11 CeedVector u, v; in main() local
16 CeedVectorCreate(ceed, len, &u); in main()
23 CeedVectorSetArray(u, CEED_MEM_HOST, CEED_COPY_VALUES, u_array); in main()
28 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
41 CeedVectorDestroy(&u); in main()
H A Dt411-qfunction.c11 CeedVector u, v; in main() local
18 CeedVectorCreate(ceed, q, &u); in main()
20 CeedVectorSetArray(u, CEED_MEM_HOST, CEED_USE_POINTER, u_array); in main()
26 in[0] = u; in main()
42 CeedVectorDestroy(&u); in main()
H A Dt409-qfunction.c13 CeedVector u, v; in main() local
22 CeedVectorCreate(ceed, q, &u); in main()
23 CeedVectorSetValue(u, 1.0); in main()
36 in[0] = u; in main()
79 in[0] = u; in main()
85 CeedVectorDestroy(&u); in main()
H A Dt412-qfunction.c11 CeedVector u, v; in main() local
18 CeedVectorCreate(ceed, q * size, &u); in main()
20 CeedVectorSetArray(u, CEED_MEM_HOST, CEED_USE_POINTER, u_array); in main()
26 in[0] = u; in main()
42 CeedVectorDestroy(&u); in main()
H A Dt325-basis.c12 CeedVector u, v; in main() local
24 CeedVectorCreate(ceed, q * dim * num_comp, &u); in main()
28 CeedVectorGetArrayWrite(u, CEED_MEM_HOST, &u_array); in main()
34 CeedVectorRestoreArray(u, &u_array); in main()
39 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_GRAD, u, v); in main()
64 CeedVectorDestroy(&u); in main()
H A Dt303-basis.c12 CeedVector u, v; in main() local
17 CeedVectorCreate(ceed, len, &u); in main()
23 CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, u, v); in main()
27 CeedVectorDestroy(&u); in main()
H A Dt509-operator.c14 CeedVector u, v; in main() local
20 CeedVectorCreate(ceed, num_nodes, &u); in main()
35 CeedVectorSetValue(u, 3.0); in main()
36 CeedOperatorApply(op_identity, u, v, CEED_REQUEST_IMMEDIATE); in main()
49 CeedVectorDestroy(&u); in main()
H A Dt324-basis.c12 CeedVector u, v; in main() local
21 CeedVectorCreate(ceed, q * dim, &u); in main()
22 CeedVectorSetValue(u, 1); in main()
29 CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_GRAD, u, v); in main()
48 CeedVectorDestroy(&u); in main()
/libCEED/include/ceed/jit-source/hip/
H A Dhip-ref-restriction-curl-oriented.h16 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedNoTranspose() argument
29 value += loc_node > 0 ? u[ind_dl + comp * RSTR_COMP_STRIDE] * curl_orient_dl : 0.0; in CurlOrientedNoTranspose()
30 value += u[ind_d + comp * RSTR_COMP_STRIDE] * curl_orient_d; in CurlOrientedNoTranspose()
31 …value += loc_node < (RSTR_ELEM_SIZE - 1) ? u[ind_du + comp * RSTR_COMP_STRIDE] * curl_orient_du : … in CurlOrientedNoTranspose()
41 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedUnsignedNoTranspose() argument
54 value += loc_node > 0 ? u[ind_dl + comp * RSTR_COMP_STRIDE] * curl_orient_dl : 0.0; in CurlOrientedUnsignedNoTranspose()
55 value += u[ind_d + comp * RSTR_COMP_STRIDE] * curl_orient_d; in CurlOrientedUnsignedNoTranspose()
56 …value += loc_node < (RSTR_ELEM_SIZE - 1) ? u[ind_du + comp * RSTR_COMP_STRIDE] * curl_orient_du : … in CurlOrientedUnsignedNoTranspose()
67 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedTranspose() argument
78 …value += loc_node > 0 ? u[loc_node - 1 + comp * RSTR_ELEM_SIZE * RSTR_NUM_ELEM + elem * RSTR_ELEM_… in CurlOrientedTranspose()
[all …]
/libCEED/include/ceed/jit-source/cuda/
H A Dcuda-ref-restriction-curl-oriented.h16 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedNoTranspose() argument
29 value += loc_node > 0 ? u[ind_dl + comp * RSTR_COMP_STRIDE] * curl_orient_dl : 0.0; in CurlOrientedNoTranspose()
30 value += u[ind_d + comp * RSTR_COMP_STRIDE] * curl_orient_d; in CurlOrientedNoTranspose()
31 …value += loc_node < (RSTR_ELEM_SIZE - 1) ? u[ind_du + comp * RSTR_COMP_STRIDE] * curl_orient_du : … in CurlOrientedNoTranspose()
41 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedUnsignedNoTranspose() argument
54 value += loc_node > 0 ? u[ind_dl + comp * RSTR_COMP_STRIDE] * curl_orient_dl : 0.0; in CurlOrientedUnsignedNoTranspose()
55 value += u[ind_d + comp * RSTR_COMP_STRIDE] * curl_orient_d; in CurlOrientedUnsignedNoTranspose()
56 …value += loc_node < (RSTR_ELEM_SIZE - 1) ? u[ind_du + comp * RSTR_COMP_STRIDE] * curl_orient_du : … in CurlOrientedUnsignedNoTranspose()
67 … const CeedScalar *__restrict__ u, CeedScalar *__restrict__ v) { in CurlOrientedTranspose() argument
78 …value += loc_node > 0 ? u[loc_node - 1 + comp * RSTR_ELEM_SIZE * RSTR_NUM_ELEM + elem * RSTR_ELEM_… in CurlOrientedTranspose()
[all …]
/libCEED/examples/solids/src/
H A Dboundary.c23 …tscReal load_increment, const PetscReal coords[], PetscInt num_comp_u, PetscScalar *u, void *ctx) { in BCMMS() argument
30 u[0] = exp(2 * x) * sin(3 * y) * cos(4 * z) / 1e8 * load_increment; in BCMMS()
31 u[1] = exp(3 * y) * sin(4 * z) * cos(2 * x) / 1e8 * load_increment; in BCMMS()
32 u[2] = exp(4 * z) * sin(2 * x) * cos(3 * y) / 1e8 * load_increment; in BCMMS()
42 …tscReal load_increment, const PetscReal coords[], PetscInt num_comp_u, PetscScalar *u, void *ctx) { in BCClamp() argument
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/backends/avx/
H A Dceed-avx-tensor.c44 … const CeedScalar *restrict u, CeedScalar *restrict v, const CeedInt JJ, const CeedInt CC) { in CeedTensorContract_Avx_Blocked() argument
64 fmadd(vv[jj][cc], tqv, loadu(&u[(a * B + b) * C + c + cc * 4])); in CeedTensorContract_Avx_Blocked()
88 fmadd(vv[jj][cc], tqv, loadu(&u[(a * B + b) * C + c + cc * 4])); in CeedTensorContract_Avx_Blocked()
106 … const CeedScalar *restrict u, CeedScalar *restrict v, const CeedInt JJ, const CeedInt CC) { in CeedTensorContract_Avx_Remainder() argument
127 if (C - c == 1) tqu = set(0.0, 0.0, 0.0, u[(a * B + b) * C + c + 0]); in CeedTensorContract_Avx_Remainder()
128 … else if (C - c == 2) tqu = set(0.0, 0.0, u[(a * B + b) * C + c + 1], u[(a * B + b) * C + c + 0]); in CeedTensorContract_Avx_Remainder()
129 …else if (C - c == 3) tqu = set(0.0, u[(a * B + b) * C + c + 2], u[(a * B + b) * C + c + 1], u[(a *… in CeedTensorContract_Avx_Remainder()
130 else tqu = loadu(&u[(a * B + b) * C + c]); in CeedTensorContract_Avx_Remainder()
143 … for (CeedInt c = (C / CC) * CC; c < C; c++) v[(a * J + j) * C + c] += tq * u[(a * B + b) * C + c]; in CeedTensorContract_Avx_Remainder()
154 … CeedTransposeMode t_mode, const CeedInt add, const CeedScalar *restrict u, CeedScalar *restrict v, in CeedTensorContract_Avx_Single() argument
[all …]
/libCEED/julia/LibCEED.jl/examples/
H A Dex1-volume.jl18 u = coords[i]
20 u = 1.0 + u
22 coords[i] = u*cos(v)
23 coords[i+num_nodes] = u*sin(v)
109 (u, :in, EVAL_INTERP),
113 v .= qdata*u
125 (:u, sol_rstr, sol_basis, CeedVectorActive()),
135 u = CeedVector(ceed, sol_size)
137 # Initialize 'u' with ones.
138 u[] = 1.0
[all …]
H A Dex3-volume.jl18 u = coords[i]
20 u = 1.0 + u
22 coords[i] = u*cos(v)
23 coords[i+num_nodes] = u*sin(v)
126 (u, :in, EVAL_INTERP),
132 # Apply mass: v = qdata[1] * u
133 v .= qdata[1].*u
162 (:u, sol_rstr, sol_basis, CeedVectorActive()),
174 u = CeedVector(ceed, sol_size)
176 # Initialize 'u' with ones.
[all …]

12345678910>>...12