1*0436c2adSjeremylt!----------------------------------------------------------------------- 2*0436c2adSjeremylt program test 3*0436c2adSjeremylt 4*0436c2adSjeremylt include 'ceedf.h' 5*0436c2adSjeremylt 6*0436c2adSjeremylt integer ceed,err 7*0436c2adSjeremylt integer x,n 8*0436c2adSjeremylt real*8 a(10) 9*0436c2adSjeremylt real*8 b(10) 10*0436c2adSjeremylt real*8 diff 11*0436c2adSjeremylt integer*8 aoffset,boffset 12*0436c2adSjeremylt character arg*32 13*0436c2adSjeremylt 14*0436c2adSjeremylt call getarg(1,arg) 15*0436c2adSjeremylt 16*0436c2adSjeremylt call ceedinit(trim(arg)//char(0),ceed,err) 17*0436c2adSjeremylt 18*0436c2adSjeremylt n=10 19*0436c2adSjeremylt 20*0436c2adSjeremylt call ceedvectorcreate(ceed,n,x,err) 21*0436c2adSjeremylt 22*0436c2adSjeremylt do i=1,10 23*0436c2adSjeremylt a(i)=0 24*0436c2adSjeremylt enddo 25*0436c2adSjeremylt 26*0436c2adSjeremylt aoffset=0 27*0436c2adSjeremylt call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,a,aoffset,err) 28*0436c2adSjeremylt call ceedvectorgetarray(x,ceed_mem_host,b,boffset,err) 29*0436c2adSjeremylt b(boffset+3) = -3.14 30*0436c2adSjeremylt call ceedvectorrestorearray(x,b,boffset,err) 31*0436c2adSjeremylt diff=a(3)+3.14 32*0436c2adSjeremylt if (abs(diff)>1.0D-15) then 33*0436c2adSjeremylt! LCOV_EXCL_START 34*0436c2adSjeremylt write(*,*) 'Error writing array a(3)=',a(3) 35*0436c2adSjeremylt! LCOV_EXCL_STOP 36*0436c2adSjeremylt endif 37*0436c2adSjeremylt 38*0436c2adSjeremylt call ceedvectordestroy(x,err) 39*0436c2adSjeremylt call ceeddestroy(ceed,err) 40*0436c2adSjeremylt 41*0436c2adSjeremylt end 42*0436c2adSjeremylt!----------------------------------------------------------------------- 43