xref: /libCEED/tests/t104-vector.c (revision 0436c2ad4dd95c9ccd89f5346568014600de10a7)
1*0436c2adSjeremylt /// @file
2*0436c2adSjeremylt /// Test CeedVectorGetArray to modify array
3*0436c2adSjeremylt /// \test Test CeedVectorGetArray to modify array
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   const CeedInt n = 10;
10*0436c2adSjeremylt   CeedScalar a[n];
11*0436c2adSjeremylt   CeedScalar *b;
12*0436c2adSjeremylt 
13*0436c2adSjeremylt   CeedInit(argv[1], &ceed);
14*0436c2adSjeremylt 
15*0436c2adSjeremylt   CeedVectorCreate(ceed, n, &x);
16*0436c2adSjeremylt   for (CeedInt i=0; i<n; i++)
17*0436c2adSjeremylt     a[i] = 0;
18*0436c2adSjeremylt   CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, a);
19*0436c2adSjeremylt 
20*0436c2adSjeremylt   CeedVectorGetArray(x, CEED_MEM_HOST, &b);
21*0436c2adSjeremylt   b[3] = -3.14;
22*0436c2adSjeremylt   CeedVectorRestoreArray(x, &b);
23*0436c2adSjeremylt 
24*0436c2adSjeremylt   if (a[3] != -3.14)
25*0436c2adSjeremylt     // LCOV_EXCL_START
26*0436c2adSjeremylt     printf("Error writing array a[3] = %f", (double)b[3]);
27*0436c2adSjeremylt   // LCOV_EXCL_STOP
28*0436c2adSjeremylt 
29*0436c2adSjeremylt   CeedVectorDestroy(&x);
30*0436c2adSjeremylt   CeedDestroy(&ceed);
31*0436c2adSjeremylt   return 0;
32*0436c2adSjeremylt }
33