10436c2adSjeremylt /// @file 20436c2adSjeremylt /// Test setting one vector from array of another vector 30436c2adSjeremylt /// \test Test setting one vector from array of another vector 40436c2adSjeremylt #include <ceed.h> 50436c2adSjeremylt 60436c2adSjeremylt int main(int argc, char **argv) { 70436c2adSjeremylt Ceed ceed; 80436c2adSjeremylt CeedVector X, Y; 90436c2adSjeremylt CeedInt n; 100436c2adSjeremylt CeedScalar a[10]; 110436c2adSjeremylt CeedScalar *x; 120436c2adSjeremylt const CeedScalar *y; 130436c2adSjeremylt 140436c2adSjeremylt CeedInit(argv[1], &ceed); 150436c2adSjeremylt 160436c2adSjeremylt n = 10; 170436c2adSjeremylt CeedVectorCreate(ceed, n, &X); 180436c2adSjeremylt CeedVectorCreate(ceed, n, &Y); 190436c2adSjeremylt 20*2b730f8bSJeremy L Thompson for (CeedInt i = 0; i < n; i++) a[i] = 10 + i; 210436c2adSjeremylt CeedVectorSetArray(X, CEED_MEM_HOST, CEED_USE_POINTER, a); 220436c2adSjeremylt 230436c2adSjeremylt CeedVectorGetArray(X, CEED_MEM_HOST, &x); 240436c2adSjeremylt CeedVectorSetArray(Y, CEED_MEM_HOST, CEED_COPY_VALUES, x); 250436c2adSjeremylt CeedVectorRestoreArray(X, &x); 260436c2adSjeremylt 270436c2adSjeremylt CeedVectorGetArrayRead(Y, CEED_MEM_HOST, &y); 28*2b730f8bSJeremy L Thompson for (CeedInt i = 0; i < n; i++) { 29*2b730f8bSJeremy L Thompson if (y[i] != 10 + i) printf("Error reading array y[%" CeedInt_FMT "] = %f\n", i, (CeedScalar)y[i]); 30*2b730f8bSJeremy L Thompson } 310436c2adSjeremylt CeedVectorRestoreArrayRead(Y, &y); 320436c2adSjeremylt 330436c2adSjeremylt CeedVectorDestroy(&X); 340436c2adSjeremylt CeedVectorDestroy(&Y); 350436c2adSjeremylt CeedDestroy(&ceed); 360436c2adSjeremylt return 0; 370436c2adSjeremylt } 38