| /libCEED/python/tests/ |
| H A D | test-1-vector.py | 23 def check_values(ceed, x, value): argument
24 with x.array_read() as b:
37 x = ceed.Vector(n)
40 x.set_array(a, cmode=libceed.USE_POINTER)
42 with x.array_read() as b:
54 x = ceed.Vector(n)
57 x.set_array(a, cmode=libceed.USE_POINTER)
59 with x.array() as b:
63 x.set_value(3.0)
64 check_values(ceed, x, 3.0)
[all …]
|
| /libCEED/examples/python/ |
| H A D | tutorial-2-elemrestriction.ipynb | 55 "# x -- x -- x -- x\n",
59 "# x -- x | x -- x | x -- x\n",
66 "x = ceed.Vector(num_elem+1)\n",
68 "x.set_array(a, cmode=libceed.USE_POINTER)\n",
80 "r.apply(x, y)\n",
101 "# x -- o -- o -- x -- o -- o -- x -- o -- o -- x\n",
137 "# x -- x -- x\n",
141 "# x -- x | x -- x | x -- x\n",
146 "x = ceed.Vector(2*num_elem)\n",
148 "x.set_array(a, cmode=libceed.USE_POINTER)\n",
[all …]
|
| /libCEED/tests/ |
| H A D | t101-vector.c | 11 static int CheckValues(Ceed ceed, CeedVector x, CeedScalar value) { in CheckValues() argument
15 CeedVectorGetLength(x, &len); in CheckValues()
16 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &read_array); in CheckValues()
20 CeedVectorRestoreArrayRead(x, &read_array); in CheckValues()
26 CeedVector x; in main() local
35 CeedVectorCreate(ceed, len, &x); in main()
37 CeedVectorSetArray(x, CEED_MEM_HOST, CEED_COPY_VALUES, array); in main()
45 CeedVectorSyncArray(x, type); in main()
51 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &read_array); in main()
55 CeedVectorRestoreArrayRead(x, &read_array); in main()
[all …]
|
| H A D | t109-vector.c | 10 CeedVector x; in main() local
16 CeedVectorCreate(ceed, len, &x); in main()
19 CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, array); in main()
25 CeedVectorTakeArray(x, CEED_MEM_HOST, &array); in main()
33 CeedVectorGetArrayWrite(x, CEED_MEM_HOST, &writable_array); in main()
35 CeedVectorRestoreArray(x, &writable_array); in main()
39 CeedVectorDestroy(&x); in main()
42 CeedVectorCreate(ceed, 0, &x); in main()
43 CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, NULL); in main()
47 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &read_array); in main()
[all …]
|
| H A D | t110-vector.c | 10 CeedVector x; in main() local
16 CeedVectorCreate(ceed, len, &x); in main()
17 CeedVectorSetValue(x, 0.0); in main()
20 CeedVectorGetArray(x, CEED_MEM_HOST, &a); in main()
21 CeedVectorGetArray(x, CEED_MEM_HOST, &b); in main()
24 CeedVectorRestoreArray(x, &a); in main()
25 CeedVectorRestoreArray(x, &b); in main()
27 CeedVectorDestroy(&x); in main()
|
| H A D | t102-vector.c | 10 CeedVector x; in main() local
15 CeedVectorCreate(ceed, len, &x); in main()
16 CeedVectorSetValue(x, 0.0); in main()
22 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &a); in main()
23 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &b); in main()
25 CeedVectorRestoreArrayRead(x, &a); in main()
26 CeedVectorRestoreArrayRead(x, &b); in main()
29 CeedVectorDestroy(&x); in main()
|
| H A D | t113-vector.c | 10 CeedVector x; in main() local
17 CeedVectorCreate(ceed, len, &x); in main()
18 CeedVectorSetValue(x, 0.0); in main()
19 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &a); in main()
22 CeedVectorGetArray(x, CEED_MEM_HOST, &b); in main()
25 CeedVectorRestoreArrayRead(x, &a); in main()
26 CeedVectorRestoreArray(x, &b); in main()
28 CeedVectorDestroy(&x); in main()
|
| H A D | t115-vector.c | 10 CeedVector x; in main() local
17 CeedVectorCreate(ceed, len, &x); in main()
18 CeedVectorSetValue(x, 0.0); in main()
21 CeedVectorGetArray(x, CEED_MEM_HOST, &a); in main()
22 CeedVectorGetArrayRead(x, CEED_MEM_HOST, &b); in main()
25 CeedVectorRestoreArray(x, &a); in main()
26 CeedVectorRestoreArrayRead(x, &b); in main()
28 CeedVectorDestroy(&x); in main()
|
| H A D | t102-vector-f.f90 | 7 integer x,n variable
19 call ceedvectorcreate(ceed,n,x,err)
20 call ceedvectorsetvalue(x,0.0,err)
22 call ceedvectorgetarrayread(x,ceed_mem_host,a,aoffset,err)
23 call ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
25 call ceedvectorrestorearrayread(x,a,aoffset,err)
26 call ceedvectorrestorearrayread(x,b,boffset,err)
28 call ceedvectordestroy(x,err)
|
| H A D | t130-vector.c | 16 CeedVector x, y; in main() local
18 CeedGetWorkVector(ceed, 20, &x); in main()
20 CeedVector x_copy = x; in main()
22 CeedRestoreWorkVector(ceed, &x); in main()
30 CeedVector x, y; in main() local
32 CeedGetWorkVector(ceed, 20, &x); in main()
34 CeedVector x_copy = x; in main()
36 CeedRestoreWorkVector(ceed, &x); in main()
|
| /libCEED/examples/solids/qfunctions/ |
| H A D | manufactured-force.h | 43 CeedScalar x = coords[i + 0 * Q], y = coords[i + 1 * Q], z = coords[i + 2 * Q]; in SetupMMSForce() local
48 …force[i + 0 * Q] = (-(E * (cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 4.0 - cos(z * 4.0) * sin(y… in SetupMMSForce()
50 …(E * (cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (4.5) + sin(x * 2.0) * sin(z * 4.0) * exp(y * 3… in SetupMMSForce()
52 … (E * nu * cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 8.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
53 … (E * nu * sin(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
54 …(E * cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (nu - 1.0) * 4.0) / ((nu * 2.0 - 1.0) * (nu + 1.… in SetupMMSForce()
58 …force[i + 1 * Q] = (-(E * (cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 3.0 - cos(x * 2.0) * sin(z… in SetupMMSForce()
60 …(E * (cos(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 8.0 + sin(x * 2.0) * sin(y * 3.0) * exp(z * 4.0… in SetupMMSForce()
62 … (E * nu * cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
63 … (E * nu * sin(x * 2.0) * sin(y * 3.0) * exp(z * 4.0) * 12.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
[all …]
|
| /libCEED/backends/sycl-ref/kernels/ |
| H A D | sycl-ref-vector.cpp | 16 int index = threadIdx.x + blockDim.x * blockIdx.x; in setValueK()
39 int index = threadIdx.x + blockDim.x * blockIdx.x; in rcpValueK()
61 __global__ static void scaleValueK(CeedScalar *__restrict__ x, CeedScalar alpha, CeedInt size) { in scaleValueK() argument
62 int index = threadIdx.x + blockDim.x * blockIdx.x; in scaleValueK()
65 x[index] *= alpha; in scaleValueK()
84 …xpyValueK(CeedScalar *__restrict__ y, CeedScalar alpha, CeedScalar *__restrict__ x, CeedInt size) { in axpyValueK() argument
85 int index = threadIdx.x + blockDim.x * blockIdx.x; in axpyValueK()
87 y[index] += alpha * x[index]; in axpyValueK()
106 __global__ static void pointwiseMultValueK(CeedScalar *__restrict__ w, CeedScalar *x, CeedScalar *_… in pointwiseMultValueK() argument
107 int index = threadIdx.x + blockDim.x * blockIdx.x; in pointwiseMultValueK()
[all …]
|
| /libCEED/backends/cuda-ref/kernels/ |
| H A D | cuda-ref-vector.cu | 15 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in copyStridedK()
39 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in setValueK()
61 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in setValueStridedK()
85 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in rcpValueK()
108 __global__ static void scaleValueK(CeedScalar *__restrict__ x, CeedScalar alpha, CeedSize size) { in scaleValueK() argument
109 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in scaleValueK()
111 if (index < size) x[index] *= alpha; in scaleValueK()
130 …pyValueK(CeedScalar *__restrict__ y, CeedScalar alpha, CeedScalar *__restrict__ x, CeedSize size) { in axpyValueK() argument
131 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in axpyValueK()
133 if (index < size) y[index] += alpha * x[index]; in axpyValueK()
[all …]
|
| /libCEED/backends/hip-ref/kernels/ |
| H A D | hip-ref-vector.hip.cpp | 15 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in copyStridedK()
39 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in setValueK()
61 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in setValueStridedK()
85 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in rcpValueK()
108 __global__ static void scaleValueK(CeedScalar *__restrict__ x, CeedScalar alpha, CeedSize size) { in scaleValueK() argument
109 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in scaleValueK()
111 if (index < size) x[index] *= alpha; in scaleValueK()
130 …pyValueK(CeedScalar *__restrict__ y, CeedScalar alpha, CeedScalar *__restrict__ x, CeedSize size) { in axpyValueK() argument
131 const CeedSize index = threadIdx.x + (CeedSize)blockDim.x * blockIdx.x; in axpyValueK()
133 if (index < size) y[index] += alpha * x[index]; in axpyValueK()
[all …]
|
| /libCEED/include/ceed/jit-source/cuda/ |
| H A D | cuda-shared-basis-nontensor.h | 22 data.t_id_x = threadIdx.x; in Interp()
25 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in Interp()
37 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Interp()
49 data.t_id_x = threadIdx.x; in InterpTranspose()
52 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTranspose()
64 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpTranspose()
76 data.t_id_x = threadIdx.x; in InterpTransposeAdd()
79 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTransposeAdd()
91 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpTransposeAdd()
105 data.t_id_x = threadIdx.x; in Grad()
[all …]
|
| H A D | cuda-ref-restriction-offset.h | 16 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OffsetNoTranspose()
32 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OffsetTranspose()
47 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in OffsetTranspose()
|
| H A D | cuda-ref-restriction-oriented.h | 17 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OrientedNoTranspose()
35 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OrientedTranspose()
53 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in OrientedTranspose()
|
| H A D | cuda-ref-restriction-strided.h | 16 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in StridedNoTranspose()
31 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in StridedTranspose()
|
| H A D | cuda-shared-basis-tensor-at-points.h | 28 data.t_id_x = threadIdx.x; in InterpAtPoints()
31 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpAtPoints()
45 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpAtPoints()
60 …const CeedInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.… in InterpAtPoints()
62 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in InterpAtPoints()
84 data.t_id_x = threadIdx.x; in InterpTransposeAtPoints()
87 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTransposeAtPoints()
101 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpTransposeAtPoints()
117 …const CeedInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.… in InterpTransposeAtPoints()
119 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in InterpTransposeAtPoints()
[all …]
|
| /libCEED/include/ceed/jit-source/hip/ |
| H A D | hip-shared-basis-nontensor.h | 23 data.t_id_x = threadIdx.x; in __launch_bounds__()
26 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__()
38 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in __launch_bounds__()
50 data.t_id_x = threadIdx.x; in __launch_bounds__()
53 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__()
65 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in __launch_bounds__()
77 data.t_id_x = threadIdx.x; in __launch_bounds__()
80 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__()
92 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in __launch_bounds__()
107 data.t_id_x = threadIdx.x; in __launch_bounds__()
[all …]
|
| H A D | hip-ref-restriction-offset.h | 16 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OffsetNoTranspose()
32 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OffsetTranspose()
47 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in OffsetTranspose()
|
| H A D | hip-ref-restriction-oriented.h | 17 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OrientedNoTranspose()
35 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in OrientedTranspose()
53 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in OrientedTranspose()
|
| H A D | hip-ref-restriction-strided.h | 16 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in StridedNoTranspose()
31 …r (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; no… in StridedTranspose()
|
| /libCEED/examples/petsc/ |
| H A D | index.md | 10 …ith the exception that the physical coordinates for this problem are $\bm{x}=(x,y,z)\in \mathbb{R}…
23 …e, denoted by $\bar{\bm{x}}=(\bar{x},\bar{y},\bar{z})$, and physical coordinates on the discrete s…
26 …c{\partial \bm{x}}{\partial \bm{X}}_{(2\times2)} = \frac{\partial {\bm{x}}}{\partial \bar{\bm{x}}}…
32 …|col_1\left(\frac{\partial \bar{\bm{x}}}{\partial \bm{X}}\right)\right\| \left\|col_2 \left(\frac{…
35 … ${\partial\bar{\bm{x}}}/{\partial \bm{X}}_{(3\times2)}$ is provided by the library, while ${\part…
38 …x}}}{\partial \bm{X}}\right) / \left\| col_1\left(\frac{\partial\bar{\bm{x}}}{\partial \bm{X}}\rig…
60 …c}{\bm{x}}=(\overset{\circ}{x},\overset{\circ}{y},\overset{\circ}{z})$, and physical coordinates o…
63 …t{\circ}{\bm{x}}}{\partial \bm{X}}_{(3\times2)} = \frac{\partial \overset{\circ}{\bm{x}}}{\partial…
69 …frac{\partial \overset{\circ}{\bm{x}}}{\partial \bm{X}}\right) \times col_2 \left(\frac{\partial \…
77 …d by $\bm x(\bm X)$, are mapped to their corresponding radial projections on the circle, which hav…
[all …]
|
| /libCEED/examples/fluids/qfunctions/ |
| H A D | taylorgreen.h | 30 CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsTaylorGreen() local
32 ScaleN(x, 2 * M_PI / context->lx, 3); in ICsTaylorGreen()
34 …ference.pressure + (density0 * Square(V0) / 16) * (cos(2 * x[0]) + cos(2 * x[1])) * (cos(2 * x[2] … in ICsTaylorGreen()
35 Y[1] = V0 * sin(x[0]) * cos(x[1]) * cos(x[2]); in ICsTaylorGreen()
36 Y[2] = -V0 * cos(x[0]) * sin(x[1]) * cos(x[2]); in ICsTaylorGreen()
|