xref: /libCEED/tests/t103-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      integer*8 aoffset,xoffset,yoffset
9*0436c2adSjeremylt      real*8 a(10)
10*0436c2adSjeremylt      real*8 xx(10)
11*0436c2adSjeremylt      real*8 yy(10)
12*0436c2adSjeremylt      real*8 diff
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      call ceedvectorcreate(ceed,n,x,err)
21*0436c2adSjeremylt      call ceedvectorcreate(ceed,n,y,err)
22*0436c2adSjeremylt
23*0436c2adSjeremylt      do i=1,10
24*0436c2adSjeremylt        a(i)=10+i
25*0436c2adSjeremylt      enddo
26*0436c2adSjeremylt
27*0436c2adSjeremylt      aoffset=0
28*0436c2adSjeremylt      call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,a,aoffset,err)
29*0436c2adSjeremylt
30*0436c2adSjeremylt      call ceedvectorgetarray(x,ceed_mem_host,xx,xoffset,err)
31*0436c2adSjeremylt      call ceedvectorsetarray(y,ceed_mem_host,ceed_copy_values,xx,xoffset,err)
32*0436c2adSjeremylt      call ceedvectorrestorearray(x,xx,xoffset,err)
33*0436c2adSjeremylt
34*0436c2adSjeremylt      call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err)
35*0436c2adSjeremylt
36*0436c2adSjeremylt      do i=1,10
37*0436c2adSjeremylt        diff=yy(i+yoffset)-10-i
38*0436c2adSjeremylt        if (abs(diff)>1.0D-15) then
39*0436c2adSjeremylt! LCOV_EXCL_START
40*0436c2adSjeremylt          write(*,*) 'Error reading array y(',i,')=',yy(i+yoffset)
41*0436c2adSjeremylt! LCOV_EXCL_STOP
42*0436c2adSjeremylt        endif
43*0436c2adSjeremylt      enddo
44*0436c2adSjeremylt
45*0436c2adSjeremylt      call ceedvectorrestorearrayread(y,yy,yoffset,err)
46*0436c2adSjeremylt      call ceedvectordestroy(x,err)
47*0436c2adSjeremylt      call ceedvectordestroy(y,err)
48*0436c2adSjeremylt      call ceeddestroy(ceed,err)
49*0436c2adSjeremylt
50*0436c2adSjeremylt      end
51*0436c2adSjeremylt!-----------------------------------------------------------------------
52