10436c2adSjeremylt /// @file 20436c2adSjeremylt /// Test CeedVectorGetArray to modify array 30436c2adSjeremylt /// \test Test CeedVectorGetArray to modify array 40436c2adSjeremylt #include <ceed.h> 50436c2adSjeremylt 60436c2adSjeremylt int main(int argc, char **argv) { 70436c2adSjeremylt Ceed ceed; 80436c2adSjeremylt CeedVector x; 9*4fee36f0SJeremy L Thompson const 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] = 0; 16*4fee36f0SJeremy L Thompson CeedVectorSetArray(x, CEED_MEM_HOST, CEED_USE_POINTER, array); 170436c2adSjeremylt 18*4fee36f0SJeremy L Thompson { 19*4fee36f0SJeremy L Thompson CeedScalar *writable_array; 200436c2adSjeremylt 21*4fee36f0SJeremy L Thompson CeedVectorGetArray(x, CEED_MEM_HOST, &writable_array); 22*4fee36f0SJeremy L Thompson writable_array[3] = -3.14; 23*4fee36f0SJeremy L Thompson CeedVectorRestoreArray(x, &writable_array); 24*4fee36f0SJeremy L Thompson } 25*4fee36f0SJeremy L Thompson 26*4fee36f0SJeremy L Thompson if (array[3] != (CeedScalar)(-3.14)) printf("Error writing array[3] = %f\n", (CeedScalar)array[3]); 270436c2adSjeremylt 280436c2adSjeremylt CeedVectorDestroy(&x); 290436c2adSjeremylt CeedDestroy(&ceed); 300436c2adSjeremylt return 0; 310436c2adSjeremylt } 32