10436c2adSjeremylt /// @file 20436c2adSjeremylt /// Test creation, setting, reading, restoring, and destroying of a vector using CEED_MEM_DEVICE 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; 90436c2adSjeremylt CeedVector y; 100436c2adSjeremylt CeedInt n; 110436c2adSjeremylt CeedScalar a[10]; 120436c2adSjeremylt const CeedScalar *b, *c; 130436c2adSjeremylt 140436c2adSjeremylt CeedInit(argv[1], &ceed); 150436c2adSjeremylt 160436c2adSjeremylt n = 10; 170436c2adSjeremylt CeedVectorCreate(ceed, n, &x); 180436c2adSjeremylt CeedVectorCreate(ceed, n, &y); 19*2b730f8bSJeremy L Thompson for (CeedInt i = 0; i < n; i++) a[i] = 10 + i; 200436c2adSjeremylt CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, a); 210436c2adSjeremylt 220436c2adSjeremylt CeedVectorGetArrayRead(x, CEED_MEM_DEVICE, &b); 230436c2adSjeremylt CeedVectorSetArray(y, CEED_MEM_DEVICE, CEED_COPY_VALUES, (CeedScalar *)b); 240436c2adSjeremylt CeedVectorRestoreArrayRead(x, &b); 250436c2adSjeremylt 260436c2adSjeremylt CeedVectorGetArrayRead(y, CEED_MEM_HOST, &c); 27*2b730f8bSJeremy L Thompson for (CeedInt i = 0; i < n; i++) { 28*2b730f8bSJeremy L Thompson if (c[i] != 10 + i) printf("Error reading array c[%" CeedInt_FMT "] = %f\n", i, (CeedScalar)c[i]); 29*2b730f8bSJeremy L Thompson } 300436c2adSjeremylt CeedVectorRestoreArrayRead(y, &c); 310436c2adSjeremylt 320436c2adSjeremylt CeedVectorDestroy(&x); 330436c2adSjeremylt CeedVectorDestroy(&y); 340436c2adSjeremylt CeedDestroy(&ceed); 350436c2adSjeremylt return 0; 360436c2adSjeremylt } 37