1*0436c2adSjeremylt /// @file 2*0436c2adSjeremylt /// Test setting one vector from array of another vector 3*0436c2adSjeremylt /// \test Test setting one vector from array of another vector 4*0436c2adSjeremylt #include <ceed.h> 5*0436c2adSjeremylt 6*0436c2adSjeremylt int main(int argc, char **argv) { 7*0436c2adSjeremylt Ceed ceed; 8*0436c2adSjeremylt CeedVector X, Y; 9*0436c2adSjeremylt CeedInt n; 10*0436c2adSjeremylt CeedScalar a[10]; 11*0436c2adSjeremylt CeedScalar *x; 12*0436c2adSjeremylt const CeedScalar *y; 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 20*0436c2adSjeremylt for (CeedInt i=0; i<n; i++) 21*0436c2adSjeremylt a[i] = 10 + i; 22*0436c2adSjeremylt CeedVectorSetArray(X, CEED_MEM_HOST, CEED_USE_POINTER, a); 23*0436c2adSjeremylt 24*0436c2adSjeremylt CeedVectorGetArray(X, CEED_MEM_HOST, &x); 25*0436c2adSjeremylt CeedVectorSetArray(Y, CEED_MEM_HOST, CEED_COPY_VALUES, x); 26*0436c2adSjeremylt CeedVectorRestoreArray(X, &x); 27*0436c2adSjeremylt 28*0436c2adSjeremylt CeedVectorGetArrayRead(Y, CEED_MEM_HOST, &y); 29*0436c2adSjeremylt for (CeedInt i=0; i<n; i++) 30*0436c2adSjeremylt if (y[i] != 10+i) 31*0436c2adSjeremylt // LCOV_EXCL_START 32*0436c2adSjeremylt printf("Error reading array y[%d] = %f",i,(double)y[i]); 33*0436c2adSjeremylt // LCOV_EXCL_STOP 34*0436c2adSjeremylt CeedVectorRestoreArrayRead(Y, &y); 35*0436c2adSjeremylt 36*0436c2adSjeremylt CeedVectorDestroy(&X); 37*0436c2adSjeremylt CeedVectorDestroy(&Y); 38*0436c2adSjeremylt CeedDestroy(&ceed); 39*0436c2adSjeremylt return 0; 40*0436c2adSjeremylt } 41