xref: /libCEED/tests/t105-vector.c (revision 0436c2ad4dd95c9ccd89f5346568014600de10a7)
1*0436c2adSjeremylt /// @file
2*0436c2adSjeremylt /// Test creation, setting, reading, restoring, and destroying of a vector using CEED_MEM_DEVICE
3*0436c2adSjeremylt /// \test Test creation, setting, reading, restoring, and destroying of a vector
4*0436c2adSjeremylt #include <ceed.h>
5*0436c2adSjeremylt 
6*0436c2adSjeremylt int main(int argc, char **argv) {
7*0436c2adSjeremylt   Ceed ceed;
8*0436c2adSjeremylt   CeedVector x;
9*0436c2adSjeremylt   CeedVector y;
10*0436c2adSjeremylt   CeedInt n;
11*0436c2adSjeremylt   CeedScalar a[10];
12*0436c2adSjeremylt   const CeedScalar *b, *c;
13*0436c2adSjeremylt 
14*0436c2adSjeremylt   CeedInit(argv[1], &ceed);
15*0436c2adSjeremylt 
16*0436c2adSjeremylt   n = 10;
17*0436c2adSjeremylt   CeedVectorCreate(ceed, n, &x);
18*0436c2adSjeremylt   CeedVectorCreate(ceed, n, &y);
19*0436c2adSjeremylt   for (CeedInt i=0; i<n; i++)
20*0436c2adSjeremylt     a[i] = 10 + i;
21*0436c2adSjeremylt   CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, a);
22*0436c2adSjeremylt 
23*0436c2adSjeremylt   CeedVectorGetArrayRead(x, CEED_MEM_DEVICE, &b);
24*0436c2adSjeremylt   CeedVectorSetArray(y, CEED_MEM_DEVICE, CEED_COPY_VALUES, (CeedScalar *)b);
25*0436c2adSjeremylt   CeedVectorRestoreArrayRead(x, &b);
26*0436c2adSjeremylt 
27*0436c2adSjeremylt   CeedVectorGetArrayRead(y, CEED_MEM_HOST, &c);
28*0436c2adSjeremylt   for (CeedInt i=0; i<n; i++)
29*0436c2adSjeremylt     if (c[i] != 10+i)
30*0436c2adSjeremylt       // LCOV_EXCL_START
31*0436c2adSjeremylt       printf("Error reading array c[%d] = %f", i, (double)c[i]);
32*0436c2adSjeremylt   // LCOV_EXCL_STOP
33*0436c2adSjeremylt   CeedVectorRestoreArrayRead(y, &c);
34*0436c2adSjeremylt 
35*0436c2adSjeremylt   CeedVectorDestroy(&x);
36*0436c2adSjeremylt   CeedVectorDestroy(&y);
37*0436c2adSjeremylt   CeedDestroy(&ceed);
38*0436c2adSjeremylt   return 0;
39*0436c2adSjeremylt }
40