xref: /libCEED/tests/t100-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,n
8*0436c2adSjeremylt      integer*8 aoffset,boffset
9*0436c2adSjeremylt      real*8 a(10)
10*0436c2adSjeremylt      real*8 b(10)
11*0436c2adSjeremylt      real*8 diff
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)=10+i
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 ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
29*0436c2adSjeremylt
30*0436c2adSjeremylt      do i=1,10
31*0436c2adSjeremylt        diff=b(i+boffset)-10-i
32*0436c2adSjeremylt        if (abs(diff)>1.0D-15) then
33*0436c2adSjeremylt! LCOV_EXCL_START
34*0436c2adSjeremylt          write(*,*) 'Error reading array b(',i,')=',b(i+boffset)
35*0436c2adSjeremylt! LCOV_EXCL_STOP
36*0436c2adSjeremylt        endif
37*0436c2adSjeremylt      enddo
38*0436c2adSjeremylt
39*0436c2adSjeremylt      call ceedvectorrestorearrayread(x,b,boffset,err)
40*0436c2adSjeremylt      call ceedvectordestroy(x,err)
41*0436c2adSjeremylt      call ceeddestroy(ceed,err)
42*0436c2adSjeremylt
43*0436c2adSjeremylt      end
44*0436c2adSjeremylt!-----------------------------------------------------------------------
45