| /libCEED/tests/ |
| H A D | t207-elemrestriction.c | 50 const CeedScalar *y_array; in main() local 52 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 55 if (y_array[2 * i] != i * mult) { in main() 57 …n restricted array y[%" CeedInt_FMT "] = %f != %f\n", 2 * i, (CeedScalar)y_array[2 * i], i * mult); in main() 60 if (y_array[2 * i + 1] != (10 + i) * mult) { in main() 62 …stricted array y[%" CeedInt_FMT "] = %f != %f\n", 2 * i + 1, (CeedScalar)y_array[2 * i + 1], (10. … in main() 66 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t206-elemrestriction.c | 50 const CeedScalar *y_array; in main() local 52 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 55 …if (y_array[i] != i * mult) printf("Error in restricted array y[%" CeedInt_FMT "] = %f != %f\n", i… in main() 56 if (y_array[i + num_elem + 1] != (10 + i) * mult) { in main() 58 …d array y[%" CeedInt_FMT "] = %f != %f\n", i + num_elem + 1, (CeedScalar)y_array[i + num_elem + 1], in main() 63 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t219-elemrestriction.c | 51 const CeedScalar *y_array; in main() local 53 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 62 if (i == 0 && 10 + k + 1 != -y_array[index]) { in main() 64 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 66 } else if (i == 1 && 10 + k != -y_array[index]) { in main() 68 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 72 if (y_array[index] != x_array[ind[k * elem_size + i]]) { in main() 74 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 81 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t218-elemrestriction.c | 70 const CeedScalar *y_array, *y_unsigned_array; in main() local 72 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 78 if (j == 2 && 10 + 3 * i + j + 1 != -y_array[k]) { in main() 80 … printf("Error in restricted array y[%" CeedInt_FMT "] = %f\n", k, (CeedScalar)y_array[k]); in main() 82 } else if (j == 3 && 10 + 3 * i + j - 1 != -y_array[k]) { in main() 84 … printf("Error in restricted array y[%" CeedInt_FMT "] = %f\n", k, (CeedScalar)y_array[k]); in main() 97 if (10 + 3 * i + j != y_array[k] || 10 + 3 * i + j != y_unsigned_array[k]) { in main() 99 … printf("Error in restricted array y[%" CeedInt_FMT "] = %f\n", k, (CeedScalar)y_array[k]); in main() 105 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t201-elemrestriction.c | 28 const CeedScalar *y_array; in main() local 30 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 36 …if (y_array[i * e_layout[0] + j * e_layout[1] + k * e_layout[2]] != x_array[i * l_layout[0] + j * … in main() 39 (CeedScalar)y_array[i * strides[0] + j * strides[1] + j * strides[2]]); in main() 45 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t205-elemrestriction.c | 39 const CeedScalar *y_array; in main() local 41 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 46 … if (y_array[i * e_layout[0] + j * e_layout[1] + k * e_layout[2]] != x_array[ind[i + k * 2] + j]) { in main() 49 …(CeedScalar)y_array[i * e_layout[0] + j * e_layout[1] + k * e_layout[2]], x_array[ind[i + k * 2] +… in main() 55 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t204-elemrestriction.c | 39 const CeedScalar *y_array; in main() local 41 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 46 …if (y_array[i * e_layout[0] + j * e_layout[1] + k * e_layout[2]] != x_array[ind[i + k * 2] + j * (… in main() 49 …(CeedScalar)y_array[i * e_layout[0] + j * e_layout[1] + k * e_layout[2]], x_array[ind[i + k * 2] +… in main() 55 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t217-elemrestriction.c | 34 const CeedScalar *y_array; in main() local 36 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 40 if (y_array[k] * CeedIntPow(-1, i % 2) != 10 + (k + 1) / 2) { in main() 42 printf("Error in restricted array y[%" CeedInt_FMT "] = %f\n", k, (CeedScalar)y_array[k]); in main() 47 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t202-elemrestriction.c | 36 const CeedScalar *y_array; in main() local 38 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 46 if (y_array[index] != x_array[ind[k * elem_size + i]]) { in main() 48 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 54 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t208-elemrestriction.c | 35 const CeedScalar *y_array; in main() local 37 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 45 if (y_array[index] != x_array[ind[k * elem_size + i]]) { in main() 47 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 53 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t203-elemrestriction.c | 42 const CeedScalar *y_array; in main() local 44 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 52 if (y_array[index] != x_array[ind[k * elem_size + i] + j * (num_elem + 1)]) { in main() 54 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 60 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| H A D | t213-elemrestriction.c | 46 const CeedScalar *y_array; in main() local 48 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 56 if (y_array[index] != x_array[ind[k * elem_size + i] + j * (num_elem + 1)]) in main() 58 … y[%" CeedInt_FMT "][%" CeedInt_FMT "][%" CeedInt_FMT "] = %f\n", i, j, k, (double)y_array[index]); in main() 63 CeedVectorRestoreArrayRead(y, &y_array); in main()
|
| /libCEED/python/tests/ |
| H A D | test-2-elemrestriction.py | 42 with y.array_read() as y_array: 44 assert 10 + (i + 1) // 2 == y_array[i] 68 with y.array_read() as y_array: 70 assert 10 + i == y_array[i] 102 with y.array_read() as y_array: 111 assert y_array[indx] == a[ind[k * elem_size + i]] 150 with y.array_read() as y_array: 159 assert y_array[indx] == a[ind[k * elem_size + i]] 293 with y.array_read() as y_array: 297 assert 10 + (k + 1) // 2 == y_array[k] * pow(-1, i % 2) [all …]
|
| H A D | test-1-vector.py | 110 with y.array_read() as y_array: 112 assert y_array[i] == 10 + i 367 with y.array_read() as y_array: 369 assert y_array[i] == a[i]
|
| /libCEED/backends/sycl-ref/kernels/ |
| H A D | sycl-ref-vector.cpp | 93 extern "C" int CeedDeviceAXPY_Sycl(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, Ceed… in CeedDeviceAXPY_Sycl() argument 99 axpyValueK<<<grid_size, block_size>>>(y_array, alpha, x_array, length); in CeedDeviceAXPY_Sycl() 116 …PointwiseMult_Sycl(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, CeedInt length) { in CeedDevicePointwiseMult_Sycl() argument 122 pointwiseMultValueK<<<grid_size, block_size>>>(w_array, x_array, y_array, length); in CeedDevicePointwiseMult_Sycl()
|
| /libCEED/backends/cuda-ref/kernels/ |
| H A D | cuda-ref-vector.cu | 139 extern "C" int CeedDeviceAXPY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, Ceed… in CeedDeviceAXPY_Cuda() argument 145 axpyValueK<<<grid_size, block_size>>>(y_array, alpha, x_array, length); in CeedDeviceAXPY_Cuda() 164 extern "C" int CeedDeviceAXPBY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedSca… in CeedDeviceAXPBY_Cuda() argument 170 axpbyValueK<<<grid_size, block_size>>>(y_array, alpha, beta, x_array, length); in CeedDeviceAXPBY_Cuda() 186 …PointwiseMult_Cuda(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, CeedSize length)… in CeedDevicePointwiseMult_Cuda() argument 192 pointwiseMultValueK<<<grid_size, block_size>>>(w_array, x_array, y_array, length); in CeedDevicePointwiseMult_Cuda()
|
| /libCEED/backends/hip-ref/kernels/ |
| H A D | hip-ref-vector.hip.cpp | 139 extern "C" int CeedDeviceAXPY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedS… in CeedDeviceAXPY_Hip() argument 145 …hipLaunchKernelGGL(axpyValueK, dim3(grid_size), dim3(block_size), 0, 0, y_array, alpha, x_array, l… in CeedDeviceAXPY_Hip() 164 extern "C" int CeedDeviceAXPBY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedScal… in CeedDeviceAXPBY_Hip() argument 170 …hipLaunchKernelGGL(axpbyValueK, dim3(grid_size), dim3(block_size), 0, 0, y_array, alpha, beta, x_a… in CeedDeviceAXPBY_Hip() 186 …ePointwiseMult_Hip(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, CeedSize length)… in CeedDevicePointwiseMult_Hip() argument 192 …L(pointwiseMultValueK, dim3(grid_size), dim3(block_size), 0, 0, w_array, x_array, y_array, length); in CeedDevicePointwiseMult_Hip()
|
| /libCEED/interface/ |
| H A D | ceed-vector.c | 773 CeedScalar *y_array = NULL; in CeedVectorAXPY() local 816 CeedCall(CeedVectorGetArray(y, CEED_MEM_HOST, &y_array)); in CeedVectorAXPY() 820 assert(y_array); in CeedVectorAXPY() 822 for (CeedSize i = 0; i < length_y; i++) y_array[i] += alpha * x_array[i]; in CeedVectorAXPY() 824 CeedCall(CeedVectorRestoreArray(y, &y_array)); in CeedVectorAXPY() 844 CeedScalar *y_array = NULL; in CeedVectorAXPBY() local 887 CeedCall(CeedVectorGetArray(y, CEED_MEM_HOST, &y_array)); in CeedVectorAXPBY() 891 assert(y_array); in CeedVectorAXPBY() 893 for (CeedSize i = 0; i < length_y; i++) y_array[i] = alpha * x_array[i] + beta * y_array[i]; in CeedVectorAXPBY() 895 CeedCall(CeedVectorRestoreArray(y, &y_array)); in CeedVectorAXPBY() [all …]
|
| /libCEED/examples/python/ |
| H A D | tutorial-2-elemrestriction.ipynb | 82 "with y.array_read() as y_array:\n", 83 " print('y =', y_array)\n" 159 "with y.array_read() as y_array:\n", 160 " print('y =', y_array)\n" 256 "with y.array_read() as y_array:\n", 257 " print('y =', y_array)\n", 326 "with y.array_read() as y_array:\n", 327 " print('y =', y_array)\n",
|
| H A D | tutorial-1-vector.ipynb | 144 "with y.array_read() as y_array:\n", 145 " print(y_array)\n"
|
| /libCEED/backends/cuda-ref/ |
| H A D | ceed-cuda-ref-vector.c | 713 static int CeedHostAXPY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedSize l… in CeedHostAXPY_Cuda() argument 714 for (CeedSize i = 0; i < length; i++) y_array[i] += alpha * x_array[i]; in CeedHostAXPY_Cuda() 721 int CeedDeviceAXPY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedSize length… 760 static int CeedHostAXPBY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedScalar *x… in CeedHostAXPBY_Cuda() argument 761 for (CeedSize i = 0; i < length; i++) y_array[i] = alpha * x_array[i] + beta * y_array[i]; in CeedHostAXPBY_Cuda() 768 int CeedDeviceAXPBY_Cuda(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedScalar *x_arra… 795 static int CeedHostPointwiseMult_Cuda(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array… in CeedHostPointwiseMult_Cuda() argument 796 for (CeedSize i = 0; i < length; i++) w_array[i] = x_array[i] * y_array[i]; in CeedHostPointwiseMult_Cuda() 803 int CeedDevicePointwiseMult_Cuda(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, Cee…
|
| /libCEED/backends/sycl-ref/ |
| H A D | ceed-sycl-vector.sycl.cpp | 580 static int CeedHostAXPY_Sycl(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedSize l… in CeedHostAXPY_Sycl() argument 581 for (CeedSize i = 0; i < length; i++) y_array[i] += alpha * x_array[i]; in CeedHostAXPY_Sycl() 588 static int CeedDeviceAXPY_Sycl(sycl::queue &sycl_queue, CeedScalar *y_array, CeedScalar alpha, Ceed… in CeedDeviceAXPY_Sycl() argument 592 sycl_queue.parallel_for(length, e, [=](sycl::id<1> i) { y_array[i] += alpha * x_array[i]; }); in CeedDeviceAXPY_Sycl() 627 static int CeedHostPointwiseMult_Sycl(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array… in CeedHostPointwiseMult_Sycl() argument 628 for (CeedSize i = 0; i < length; i++) w_array[i] = x_array[i] * y_array[i]; in CeedHostPointwiseMult_Sycl() 635 …queue &sycl_queue, CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, CeedSize length)… in CeedDevicePointwiseMult_Sycl() argument 639 sycl_queue.parallel_for(length, e, [=](sycl::id<1> i) { w_array[i] = x_array[i] * y_array[i]; }); in CeedDevicePointwiseMult_Sycl()
|
| /libCEED/backends/hip-ref/ |
| H A D | ceed-hip-ref-vector.c | 823 static int CeedHostAXPY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedSize le… in CeedHostAXPY_Hip() argument 824 for (CeedSize i = 0; i < length; i++) y_array[i] += alpha * x_array[i]; in CeedHostAXPY_Hip() 831 int CeedDeviceAXPY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar *x_array, CeedSize length); 873 static int CeedHostAXPBY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedScalar *x_… in CeedHostAXPBY_Hip() argument 874 for (CeedSize i = 0; i < length; i++) y_array[i] = alpha * x_array[i] + beta * y_array[i]; in CeedHostAXPBY_Hip() 881 int CeedDeviceAXPBY_Hip(CeedScalar *y_array, CeedScalar alpha, CeedScalar beta, CeedScalar *x_array… 908 static int CeedHostPointwiseMult_Hip(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array,… in CeedHostPointwiseMult_Hip() argument 909 for (CeedSize i = 0; i < length; i++) w_array[i] = x_array[i] * y_array[i]; in CeedHostPointwiseMult_Hip() 916 int CeedDevicePointwiseMult_Hip(CeedScalar *w_array, CeedScalar *x_array, CeedScalar *y_array, Ceed…
|