xref: /libCEED/tests/t100-vector.c (revision 4fee36f0a30516a0b5ad51bf7eb3b32d83efd623)
10436c2adSjeremylt /// @file
20436c2adSjeremylt /// Test creation, setting, reading, restoring, and destroying of a vector
30436c2adSjeremylt /// \test Test creation, setting, reading, restoring, and destroying of a vector
40436c2adSjeremylt #include <ceed.h>
50436c2adSjeremylt 
60436c2adSjeremylt int main(int argc, char **argv) {
70436c2adSjeremylt   Ceed       ceed;
80436c2adSjeremylt   CeedVector x;
9*4fee36f0SJeremy L Thompson   CeedInt    len = 10;
10*4fee36f0SJeremy L Thompson   CeedScalar array[len];
110436c2adSjeremylt 
120436c2adSjeremylt   CeedInit(argv[1], &ceed);
130436c2adSjeremylt 
14*4fee36f0SJeremy L Thompson   CeedVectorCreate(ceed, len, &x);
15*4fee36f0SJeremy L Thompson   for (CeedInt i = 0; i < len; i++) array[i] = len + i;
16*4fee36f0SJeremy L Thompson   CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, array);
170436c2adSjeremylt 
18*4fee36f0SJeremy L Thompson   {
19*4fee36f0SJeremy L Thompson     const CeedScalar *read_array;
20*4fee36f0SJeremy L Thompson 
21*4fee36f0SJeremy L Thompson     CeedVectorGetArrayRead(x, CEED_MEM_HOST, &read_array);
22*4fee36f0SJeremy L Thompson     for (CeedInt i = 0; i < len; i++) {
23*4fee36f0SJeremy L Thompson       if (read_array[i] != len + i) printf("Error reading array[%" CeedInt_FMT "] = %f\n", i, (CeedScalar)read_array[i]);
242b730f8bSJeremy L Thompson     }
25*4fee36f0SJeremy L Thompson     CeedVectorRestoreArrayRead(x, &read_array);
26*4fee36f0SJeremy L Thompson   }
270436c2adSjeremylt 
280436c2adSjeremylt   CeedVectorDestroy(&x);
290436c2adSjeremylt   CeedDestroy(&ceed);
300436c2adSjeremylt   return 0;
310436c2adSjeremylt }
32