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; 8*4fee36f0SJeremy L Thompson CeedVector x, y; 9*4fee36f0SJeremy L Thompson CeedInt len = 10; 100436c2adSjeremylt 110436c2adSjeremylt CeedInit(argv[1], &ceed); 120436c2adSjeremylt 13*4fee36f0SJeremy L Thompson CeedVectorCreate(ceed, len, &x); 14*4fee36f0SJeremy L Thompson CeedVectorCreate(ceed, len, &y); 150436c2adSjeremylt 16*4fee36f0SJeremy L Thompson { 17*4fee36f0SJeremy L Thompson CeedScalar array[len]; 180436c2adSjeremylt 19*4fee36f0SJeremy L Thompson for (CeedInt i = 0; i < len; i++) array[i] = len + i; 20*4fee36f0SJeremy L Thompson CeedVectorSetArray(x, CEED_MEM_HOST, CEED_COPY_VALUES, array); 212b730f8bSJeremy L Thompson } 22*4fee36f0SJeremy L Thompson { 23*4fee36f0SJeremy L Thompson const CeedScalar *read_array; 240436c2adSjeremylt 25*4fee36f0SJeremy L Thompson CeedVectorGetArrayRead(x, CEED_MEM_HOST, &read_array); 26*4fee36f0SJeremy L Thompson CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, (CeedScalar *)read_array); 27*4fee36f0SJeremy L Thompson CeedVectorRestoreArrayRead(x, &read_array); 28*4fee36f0SJeremy L Thompson } 29*4fee36f0SJeremy L Thompson { 30*4fee36f0SJeremy L Thompson const CeedScalar *read_array; 31*4fee36f0SJeremy L Thompson 32*4fee36f0SJeremy L Thompson CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); 33*4fee36f0SJeremy L Thompson for (CeedInt i = 0; i < len; i++) { 34*4fee36f0SJeremy L Thompson if (read_array[i] != len + i) printf("Error reading array[%" CeedInt_FMT "] = %f\n", i, (CeedScalar)read_array[i]); 35*4fee36f0SJeremy L Thompson } 36*4fee36f0SJeremy L Thompson CeedVectorRestoreArrayRead(y, &read_array); 37*4fee36f0SJeremy L Thompson } 38*4fee36f0SJeremy L Thompson 39*4fee36f0SJeremy L Thompson CeedVectorDestroy(&x); 40*4fee36f0SJeremy L Thompson CeedVectorDestroy(&y); 410436c2adSjeremylt CeedDestroy(&ceed); 420436c2adSjeremylt return 0; 430436c2adSjeremylt } 44