xref: /libCEED/tests/t105-vector-f.f90 (revision 0436c2ad4dd95c9ccd89f5346568014600de10a7)
1*0436c2adSjeremylt!-----------------------------------------------------------------------
2*0436c2adSjeremylt      program test
3*0436c2adSjeremylt
4*0436c2adSjeremylt      include 'ceedf.h'
5*0436c2adSjeremylt
6*0436c2adSjeremylt      integer ceed,err
7*0436c2adSjeremylt      integer x,y,n
8*0436c2adSjeremylt      real*8 a(10)
9*0436c2adSjeremylt      real*8 b(10)
10*0436c2adSjeremylt      real*8 c(10)
11*0436c2adSjeremylt      real*8 diff
12*0436c2adSjeremylt      integer*8 aoffset,boffset,coffset
13*0436c2adSjeremylt      character arg*32
14*0436c2adSjeremylt
15*0436c2adSjeremylt      call getarg(1,arg)
16*0436c2adSjeremylt
17*0436c2adSjeremylt      call ceedinit(trim(arg)//char(0),ceed,err)
18*0436c2adSjeremylt
19*0436c2adSjeremylt      n=10
20*0436c2adSjeremylt
21*0436c2adSjeremylt      call ceedvectorcreate(ceed,n,x,err)
22*0436c2adSjeremylt      call ceedvectorcreate(ceed,n,y,err)
23*0436c2adSjeremylt
24*0436c2adSjeremylt      do i=1,10
25*0436c2adSjeremylt        a(i)=9+i
26*0436c2adSjeremylt      enddo
27*0436c2adSjeremylt
28*0436c2adSjeremylt      aoffset=0
29*0436c2adSjeremylt      call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,a,aoffset,err)
30*0436c2adSjeremylt      call ceedvectorgetarrayread(x,ceed_mem_device,b,boffset,err)
31*0436c2adSjeremylt      call ceedvectorsetarray(y,ceed_mem_device,ceed_copy_values,b,boffset,err)
32*0436c2adSjeremylt      call ceedvectorrestorearrayread(x,b,boffset,err)
33*0436c2adSjeremylt      call ceedvectorgetarrayread(y,ceed_mem_host,c,coffset,err)
34*0436c2adSjeremylt      do i=1,10
35*0436c2adSjeremylt        diff = c(i+coffset)-(9+i)
36*0436c2adSjeremylt        if (abs(diff)>1.0D-15) then
37*0436c2adSjeremylt! LCOV_EXCL_START
38*0436c2adSjeremylt          write(*,*) 'Error reading array c(',i,') = ',c(i+coffset),' != ',9+i
39*0436c2adSjeremylt! LCOV_EXCL_STOP
40*0436c2adSjeremylt        endif
41*0436c2adSjeremylt      enddo
42*0436c2adSjeremylt      call ceedvectorrestorearrayread(y,c,coffset,err)
43*0436c2adSjeremylt
44*0436c2adSjeremylt      call ceedvectordestroy(x,err)
45*0436c2adSjeremylt      call ceedvectordestroy(y,err)
46*0436c2adSjeremylt      call ceeddestroy(ceed,err)
47*0436c2adSjeremylt
48*0436c2adSjeremylt      end
49*0436c2adSjeremylt!-----------------------------------------------------------------------
50