| /libCEED/python/tests/ |
| H A D | test-2-elemrestriction.py | 24 num_elem = 3
26 x = ceed.Vector(num_elem + 1)
27 a = np.arange(10, 10 + num_elem + 1, dtype=ceed.scalar_type())
30 ind = np.zeros(2 * num_elem, dtype="int32")
31 for i in range(num_elem):
34 r = ceed.ElemRestriction(num_elem, 2, 1, 1, num_elem + 1, ind,
37 y = ceed.Vector(2 * num_elem)
43 for i in range(2 * num_elem):
54 num_elem = 3
56 x = ceed.Vector(2 * num_elem)
[all …]
|
| /libCEED/include/ceed/jit-source/sycl/ |
| H A D | sycl-shared-basis-tensor.h | 22 kernel void Interp(const CeedInt num_elem, global const CeedScalar *restrict d_interp_1d, global co… in Interp() argument
37 …ReadElementStrided1d(BASIS_NUM_COMP, BASIS_P_1D, num_elem, 1, BASIS_NUM_NODES * num_elem, BASIS_NU… in Interp()
39 …WriteElementStrided1d(BASIS_NUM_COMP, BASIS_Q_1D, num_elem, 1, BASIS_NUM_QPTS * num_elem, BASIS_NU… in Interp()
42 …ReadElementStrided2d(BASIS_NUM_COMP, BASIS_P_1D, num_elem, 1, BASIS_NUM_NODES * num_elem, BASIS_NU… in Interp()
44 …WriteElementStrided2d(BASIS_NUM_COMP, BASIS_Q_1D, num_elem, 1, BASIS_NUM_QPTS * num_elem, BASIS_NU… in Interp()
47 …ReadElementStrided3d(BASIS_NUM_COMP, BASIS_P_1D, num_elem, 1, BASIS_NUM_NODES * num_elem, BASIS_NU… in Interp()
49 …WriteElementStrided3d(BASIS_NUM_COMP, BASIS_Q_1D, num_elem, 1, BASIS_NUM_QPTS * num_elem, BASIS_NU… in Interp()
53 kernel void InterpTranspose(const CeedInt num_elem, global const CeedScalar *restrict d_interp_1d, … in InterpTranspose() argument
71 …ReadElementStrided1d(BASIS_NUM_COMP, BASIS_Q_1D, num_elem, 1, BASIS_NUM_QPTS * num_elem, BASIS_NUM… in InterpTranspose()
73 …WriteElementStrided1d(BASIS_NUM_COMP, BASIS_P_1D, num_elem, 1, BASIS_NUM_NODES * num_elem, BASIS_N… in InterpTranspose()
[all …]
|
| /libCEED/include/ceed/jit-source/hip/ |
| H A D | hip-ref-basis-nontensor.h | 21 extern "C" __global__ void Interp(const CeedInt num_elem, const CeedScalar *__restrict__ d_B, const… in Interp() argument
23 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Interp()
24 …IS_Q_COMP_INTERP, BASIS_P, BASIS_Q>(elem, BASIS_P, BASIS_Q, BASIS_P * num_elem, BASIS_Q * num_elem, in Interp()
25 … BASIS_NUM_COMP * BASIS_Q * num_elem, d_B, d_U, d_V); in Interp()
29 extern "C" __global__ void InterpTranspose(const CeedInt num_elem, const CeedScalar *__restrict__ d… in InterpTranspose() argument
31 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpTranspose()
32 …IS_Q_COMP_INTERP, BASIS_P, BASIS_Q>(elem, BASIS_Q, BASIS_P, BASIS_Q * num_elem, BASIS_P * num_elem, in InterpTranspose()
33 … BASIS_NUM_COMP * BASIS_Q * num_elem, d_B, d_U, d_V); in InterpTranspose()
40 extern "C" __global__ void Deriv(const CeedInt num_elem, const CeedScalar *__restrict__ d_B, const … in Deriv() argument
42 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Deriv()
[all …]
|
| H A D | hip-shared-basis-tensor.h | 19 …void Interp(const CeedInt num_elem, const CeedScalar *c_B, const CeedScalar *__restrict__ d_U, Cee… in __launch_bounds__()
38 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in __launch_bounds__()
40 …ReadElementStrided1d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * num_elem, BASIS_P_1D,… in __launch_bounds__()
42 …WriteElementStrided1d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * num_elem, BASIS_Q_1D… in __launch_bounds__()
44 …2d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * BASIS_P_1D * num_elem, BASIS_P_1D * BAS… in __launch_bounds__()
46 …2d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * BASIS_Q_1D * num_elem, BASIS_Q_1D * BAS… in __launch_bounds__()
48 …rided3d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * BASIS_P_1D * BASIS_P_1D * num_elem, in __launch_bounds__()
51 …rided3d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * BASIS_Q_1D * BASIS_Q_1D * num_elem, in __launch_bounds__()
58 …void InterpCollocated(const CeedInt num_elem, const CeedScalar *c_B, const CeedScalar *__restrict_… in __launch_bounds__()
71 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in __launch_bounds__()
[all …]
|
| /libCEED/include/ceed/jit-source/cuda/ |
| H A D | cuda-ref-basis-nontensor.h | 21 extern "C" __global__ void Interp(const CeedInt num_elem, const CeedScalar *__restrict__ d_B, const… in Interp() argument
23 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Interp()
24 …IS_Q_COMP_INTERP, BASIS_P, BASIS_Q>(elem, BASIS_P, BASIS_Q, BASIS_P * num_elem, BASIS_Q * num_elem, in Interp()
25 … BASIS_NUM_COMP * BASIS_Q * num_elem, d_B, d_U, d_V); in Interp()
29 extern "C" __global__ void InterpTranspose(const CeedInt num_elem, const CeedScalar *__restrict__ d… in InterpTranspose() argument
31 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpTranspose()
32 …IS_Q_COMP_INTERP, BASIS_P, BASIS_Q>(elem, BASIS_Q, BASIS_P, BASIS_Q * num_elem, BASIS_P * num_elem, in InterpTranspose()
33 … BASIS_NUM_COMP * BASIS_Q * num_elem, d_B, d_U, d_V); in InterpTranspose()
40 extern "C" __global__ void Deriv(const CeedInt num_elem, const CeedScalar *__restrict__ d_B, const … in Deriv() argument
42 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Deriv()
[all …]
|
| H A D | cuda-shared-basis-tensor.h | 18 extern "C" __global__ void Interp(const CeedInt num_elem, const CeedScalar *c_B, const CeedScalar *… in Interp() argument
37 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in Interp()
39 …ReadElementStrided1d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * num_elem, BASIS_P_1D,… in Interp()
41 …WriteElementStrided1d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * num_elem, BASIS_Q_1D… in Interp()
43 …2d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * BASIS_P_1D * num_elem, BASIS_P_1D * BAS… in Interp()
45 …2d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * BASIS_Q_1D * num_elem, BASIS_Q_1D * BAS… in Interp()
47 …rided3d<BASIS_NUM_COMP, BASIS_P_1D>(data, elem, 1, BASIS_P_1D * BASIS_P_1D * BASIS_P_1D * num_elem, in Interp()
50 …rided3d<BASIS_NUM_COMP, BASIS_Q_1D>(data, elem, 1, BASIS_Q_1D * BASIS_Q_1D * BASIS_Q_1D * num_elem, in Interp()
56 extern "C" __global__ void InterpCollocated(const CeedInt num_elem, const CeedScalar *c_B, const Ce… in InterpCollocated() argument
70 …for (CeedInt elem = blockIdx.x * blockDim.z + threadIdx.z; elem < num_elem; elem += gridDim.x * bl… in InterpCollocated()
[all …]
|
| /libCEED/tests/ |
| H A D | t213-elemrestriction.c | 13 CeedInt num_elem = 8; in main() local
18 CeedInt ind[elem_size * num_elem]; in main()
19 CeedInt *ceed_ind = malloc(sizeof(CeedInt) * elem_size * num_elem); in main()
20 CeedScalar x_array[num_comp * (num_elem + 1)]; in main()
26 CeedVectorCreate(ceed, num_comp * (num_elem + 1), &x); in main()
27 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
28 x_array[i + 0 * (num_elem + 1)] = 10 + i; in main()
29 x_array[i + 1 * (num_elem + 1)] = 20 + i; in main()
30 x_array[i + 2 * (num_elem + 1)] = 30 + i; in main()
35 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t203-elemrestriction.c | 11 CeedInt num_elem = 8; in main() local
16 CeedInt ind[elem_size * num_elem]; in main()
17 CeedScalar x_array[num_comp * (num_elem + 1)]; in main()
23 CeedVectorCreate(ceed, num_comp * (num_elem + 1), &x); in main()
24 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
25 x_array[i + 0 * (num_elem + 1)] = 10 + i; in main()
26 x_array[i + 1 * (num_elem + 1)] = 20 + i; in main()
27 x_array[i + 2 * (num_elem + 1)] = 30 + i; in main()
32 for (CeedInt i = 0; i < num_elem; i++) { in main()
36 …CeedElemRestrictionCreateBlocked(ceed, num_elem, elem_size, blk_size, num_comp, num_elem + 1, num_… in main()
[all …]
|
| H A D | t206-elemrestriction.c | 11 CeedInt num_elem = 5; in main() local
12 CeedInt ind[2 * num_elem]; in main()
20 CeedVectorCreate(ceed, 2 * (num_elem * 2), &x); in main()
21 CeedVectorCreate(ceed, 2 * (num_elem + 1), &y); in main()
24 for (CeedInt i = 0; i < num_elem; i++) { in main()
28 …CeedElemRestrictionCreate(ceed, num_elem, 2, 2, num_elem + 1, 2 * (num_elem + 1), CEED_MEM_HOST, C… in main()
33 CeedScalar x_array[2 * (num_elem * 2)]; in main()
37 for (CeedInt k = 0; k < num_elem; k++) { // Element in main()
53 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
54 mult = i > 0 && i < num_elem ? 2 : 1; in main()
[all …]
|
| H A D | t232-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
18 CeedInt offset = num_elem + 1; in main()
19 CeedInt point_index = num_elem; in main()
21 for (CeedInt i = 0; i < num_elem; i++) { in main()
22 CeedInt num_points_in_elem = (i + 1) % num_elem + 1; in main()
31 ind[num_elem] = offset; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
38 CeedInt point_index = num_elem; in main()
41 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t204-elemrestriction.c | 11 CeedInt num_elem = 3; in main() local
12 CeedInt ind[2 * num_elem]; in main()
14 CeedScalar x_array[2 * (num_elem + 1)]; in main()
20 CeedVectorCreate(ceed, 2 * (num_elem + 1), &x); in main()
21 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
23 x_array[i + num_elem + 1] = 20 + i; in main()
26 CeedVectorCreate(ceed, 2 * (num_elem * 2), &y); in main()
28 for (CeedInt i = 0; i < num_elem; i++) { in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 2, num_elem + 1, 2 * (num_elem + 1), CEED_MEM_HOST, C… in main()
45 for (CeedInt k = 0; k < num_elem; k++) { // Element in main()
[all …]
|
| H A D | t233-elemrestriction.c | 9 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
10 CeedInt ind[(num_elem + 1) + num_points]; in main()
17 CeedInt offset = num_elem + 1; in main()
18 CeedInt point_index = num_elem; in main()
20 for (CeedInt i = 0; i < num_elem; i++) { in main()
21 CeedInt num_points_in_elem = (i + 1) % num_elem + 1; in main()
30 ind[num_elem] = offset; in main()
32 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
36 CeedInt point_index = num_elem; in main()
39 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t234-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
18 CeedInt offset = num_elem + 1; in main()
19 CeedInt point_index = num_elem; in main()
21 for (CeedInt i = 0; i < num_elem; i++) { in main()
22 CeedInt num_points_in_elem = (i + 1) % num_elem + 1; in main()
31 ind[num_elem] = offset; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
42 if (min_points != 1 || max_points != num_elem) { in main()
52 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t593-operator.c | 13 CeedInt num_elem = 3, dim = 1, p = 3, q = 5; in main() local
14 …CeedInt num_nodes_x = num_elem + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem =… in main()
15 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points]; in main()
28 for (CeedInt i = 0; i < num_elem; i++) { in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_nodes_x, CEED_MEM_HOST, CEED_USE_POINTER, i… in main()
37 for (CeedInt i = 0; i < num_elem; i++) { in main()
42 …CeedElemRestrictionCreate(ceed, num_elem, p, 1, 1, num_nodes_u, CEED_MEM_HOST, CEED_USE_POINTER, i… in main()
51 ind_x_points[0] = num_elem + 1; in main()
53 ind_x_points[num_elem + 1 + current_index] = current_index; in main()
56 for (CeedInt e = 1; e < num_elem - 1; e++) { in main()
[all …]
|
| H A D | t231-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
18 CeedInt offset = num_elem + 1; in main()
19 CeedInt point_index = num_elem; in main()
21 for (CeedInt i = 0; i < num_elem; i++) { in main()
22 CeedInt num_points_in_elem = (i + 1) % num_elem + 1; in main()
31 ind[num_elem] = offset; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
38 CeedInt point_index = num_elem; in main()
41 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t217-elemrestriction.c | 10 CeedInt num_elem = 6, elem_size = 2; in main() local
11 CeedInt ind[elem_size * num_elem]; in main()
12 bool orients[elem_size * num_elem]; in main()
13 CeedScalar x_array[num_elem + 1]; in main()
18 CeedVectorCreate(ceed, num_elem + 1, &x); in main()
19 for (CeedInt i = 0; i < num_elem + 1; i++) x_array[i] = 10 + i; in main()
21 CeedVectorCreate(ceed, num_elem * elem_size, &y); in main()
23 for (CeedInt i = 0; i < num_elem; i++) { in main()
29 …CeedElemRestrictionCreateOriented(ceed, num_elem, elem_size, 1, 1, num_elem + 1, CEED_MEM_HOST, CE… in main()
37 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t220-elemrestriction.c | 12 CeedInt num_elem = 6, p = 2, dim = 1; in main() local
13 CeedInt ind[p * num_elem]; in main()
14 bool orient[p * num_elem]; in main()
15 CeedScalar x_array[num_elem + 1]; in main()
20 CeedVectorCreate(ceed, num_elem + 1, &x); in main()
21 for (CeedInt i = 0; i < num_elem + 1; i++) x_array[i] = 10 + i; in main()
23 CeedVectorCreate(ceed, num_elem * 2, &y_oriented); in main()
24 CeedVectorCreate(ceed, num_elem * 2, &y_unsigned); in main()
25 CeedVectorCreate(ceed, num_elem * 2, &y_unsigned_copy); in main()
27 for (CeedInt i = 0; i < num_elem; i++) { in main()
[all …]
|
| H A D | t230-elemrestriction.c | 8 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
9 CeedInt ind[(num_elem + 1) + num_points]; in main()
15 CeedInt offset = num_elem + 1; in main()
16 CeedInt point_index = num_elem; in main()
18 for (CeedInt i = 0; i < num_elem; i++) { in main()
19 CeedInt num_points_in_elem = (i + 1) % num_elem + 1; in main()
28 ind[num_elem] = offset; in main()
30 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
|
| H A D | t594-operator.c | 13 CeedInt num_elem = 3, dim = 1, p = 3, q = 5; in main() local
14 …CeedInt num_nodes_x = num_elem + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem =… in main()
15 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points]; in main()
28 for (CeedInt i = 0; i < num_elem; i++) { in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_nodes_x, CEED_MEM_HOST, CEED_USE_POINTER, i… in main()
37 for (CeedInt i = 0; i < num_elem; i++) { in main()
42 …CeedElemRestrictionCreate(ceed, num_elem, p, 1, 1, num_nodes_u, CEED_MEM_HOST, CEED_USE_POINTER, i… in main()
51 ind_x_points[0] = num_elem + 1; in main()
53 ind_x_points[num_elem + 1 + current_index] = current_index; in main()
56 for (CeedInt e = 1; e < num_elem - 1; e++) { in main()
[all …]
|
| H A D | t200-elemrestriction.c | 10 CeedInt num_elem = 3; in main() local
11 CeedInt ind[2 * num_elem]; in main()
16 CeedVectorCreate(ceed, num_elem + 1, &x); in main()
18 CeedScalar array[num_elem + 1]; in main()
20 for (CeedInt i = 0; i < num_elem + 1; i++) array[i] = 10 + i; in main()
23 CeedVectorCreate(ceed, num_elem * 2, &y); in main()
25 for (CeedInt i = 0; i < num_elem; i++) { in main()
29 …CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_elem + 1, CEED_MEM_HOST, CEED_USE_POINTER, … in main()
35 for (CeedInt i = 0; i < num_elem * 2; i++) { in main()
|
| H A D | t205-elemrestriction.c | 11 CeedInt num_elem = 3; in main() local
12 CeedInt ind[2 * num_elem]; in main()
14 CeedScalar x_array[2 * (num_elem + 1)]; in main()
20 CeedVectorCreate(ceed, 2 * (num_elem + 1), &x); in main()
21 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
26 CeedVectorCreate(ceed, 2 * (num_elem * 2), &y); in main()
28 for (CeedInt i = 0; i < num_elem; i++) { in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 2, 1, 2 * (num_elem + 1), CEED_MEM_HOST, CEED_USE_POI… in main()
45 for (CeedInt k = 0; k < num_elem; k++) { // Element in main()
49 … * e_layout[0] + j * e_layout[1] + k * e_layout[2]], x_array[ind[i + k * 2] + j * (num_elem + 1)]); in main()
|
| H A D | t207-elemrestriction.c | 11 CeedInt num_elem = 5; in main() local
12 CeedInt ind[2 * num_elem]; in main()
18 CeedVectorCreate(ceed, 2 * (num_elem + 1), &y); in main()
22 CeedVectorCreate(ceed, 2 * (num_elem * 2), &x); in main()
24 for (CeedInt i = 0; i < num_elem; i++) { in main()
28 …CeedElemRestrictionCreate(ceed, num_elem, 2, 2, 1, 2 * (num_elem + 1), CEED_MEM_HOST, CEED_USE_POI… in main()
33 CeedScalar x_array[2 * (num_elem * 2)]; in main()
37 for (CeedInt k = 0; k < num_elem; k++) { // Element in main()
53 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
54 mult = i > 0 && i < num_elem ? 2 : 1; in main()
|
| H A D | t209-elemrestriction.c | 10 CeedInt num_elem = 3; in main() local
11 CeedInt ind[4 * num_elem]; in main()
16 CeedVectorCreate(ceed, 3 * num_elem + 1, &mult); in main()
18 for (CeedInt i = 0; i < num_elem; i++) { in main()
24 …CeedElemRestrictionCreate(ceed, num_elem, 4, 1, 1, 3 * num_elem + 1, CEED_MEM_HOST, CEED_USE_POINT… in main()
31 for (CeedInt i = 0; i < 3 * num_elem + 1; i++) { in main()
32 if (mult_array[i] != (1 + (i > 0 && i < 3 * num_elem && (i % 3 == 0) ? 1 : 0))) { in main()
|
| H A D | t218-elemrestriction.c | 10 CeedInt num_elem = 6, elem_size = 4; in main() local
11 CeedInt ind[elem_size * num_elem]; in main()
12 CeedInt8 curl_orients[3 * elem_size * num_elem]; in main()
13 CeedScalar x_array[3 * num_elem + 1]; in main()
18 CeedVectorCreate(ceed, 3 * num_elem + 1, &x); in main()
19 for (CeedInt i = 0; i < 3 * num_elem + 1; i++) x_array[i] = 10 + i; in main()
21 CeedVectorCreate(ceed, num_elem * elem_size, &y); in main()
22 CeedVectorCreate(ceed, num_elem * elem_size, &y_unsigned); in main()
24 for (CeedInt i = 0; i < num_elem; i++) { in main()
62 …CeedElemRestrictionCreateCurlOriented(ceed, num_elem, elem_size, 1, 1, 3 * num_elem + 1, CEED_MEM_… in main()
[all …]
|
| H A D | t202-elemrestriction.c | 11 CeedInt num_elem = 8; in main() local
15 CeedInt ind[elem_size * num_elem]; in main()
16 CeedScalar x_array[num_elem + 1]; in main()
22 CeedVectorCreate(ceed, num_elem + 1, &x); in main()
23 for (CeedInt i = 0; i < num_elem + 1; i++) x_array[i] = 10 + i; in main()
27 for (CeedInt i = 0; i < num_elem; i++) { in main()
31 …CeedElemRestrictionCreateBlocked(ceed, num_elem, elem_size, blk_size, 1, 1, num_elem + 1, CEED_MEM… in main()
42 for (CeedInt k = 0; k < num_elem; k++) { // Element in main()
64 for (CeedInt i = 0; i < num_elem + 1; i++) { in main()
65 if (x_array[i] != (10 + i) * (i > 0 && i < num_elem ? 2.0 : 1.0)) { in main()
|