xref: /libCEED/tests/t105-vector-f.f90 (revision 874019bc89a44b0eae6b48d78442abcb0851055b)
10436c2adSjeremylt!-----------------------------------------------------------------------
20436c2adSjeremylt      program test
31f9a83abSJed Brown      implicit none
4*ec3da8bcSJed Brown      include 'ceed/fortran.h'
50436c2adSjeremylt
60436c2adSjeremylt      integer ceed,err
71f9a83abSJed Brown      integer x,y,i,n
80436c2adSjeremylt      real*8 a(10)
90436c2adSjeremylt      real*8 b(10)
100436c2adSjeremylt      real*8 c(10)
110436c2adSjeremylt      real*8 diff
120436c2adSjeremylt      integer*8 aoffset,boffset,coffset
130436c2adSjeremylt      character arg*32
140436c2adSjeremylt
150436c2adSjeremylt      call getarg(1,arg)
160436c2adSjeremylt
170436c2adSjeremylt      call ceedinit(trim(arg)//char(0),ceed,err)
180436c2adSjeremylt
190436c2adSjeremylt      n=10
200436c2adSjeremylt
210436c2adSjeremylt      call ceedvectorcreate(ceed,n,x,err)
220436c2adSjeremylt      call ceedvectorcreate(ceed,n,y,err)
230436c2adSjeremylt
240436c2adSjeremylt      do i=1,10
250436c2adSjeremylt        a(i)=9+i
260436c2adSjeremylt      enddo
270436c2adSjeremylt
280436c2adSjeremylt      aoffset=0
290436c2adSjeremylt      call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,a,aoffset,err)
300436c2adSjeremylt      call ceedvectorgetarrayread(x,ceed_mem_device,b,boffset,err)
310436c2adSjeremylt      call ceedvectorsetarray(y,ceed_mem_device,ceed_copy_values,b,boffset,err)
320436c2adSjeremylt      call ceedvectorrestorearrayread(x,b,boffset,err)
330436c2adSjeremylt      call ceedvectorgetarrayread(y,ceed_mem_host,c,coffset,err)
340436c2adSjeremylt      do i=1,10
350436c2adSjeremylt        diff = c(i+coffset)-(9+i)
360436c2adSjeremylt        if (abs(diff)>1.0D-15) then
370436c2adSjeremylt! LCOV_EXCL_START
380436c2adSjeremylt          write(*,*) 'Error reading array c(',i,') = ',c(i+coffset),' != ',9+i
390436c2adSjeremylt! LCOV_EXCL_STOP
400436c2adSjeremylt        endif
410436c2adSjeremylt      enddo
420436c2adSjeremylt      call ceedvectorrestorearrayread(y,c,coffset,err)
430436c2adSjeremylt
440436c2adSjeremylt      call ceedvectordestroy(x,err)
450436c2adSjeremylt      call ceedvectordestroy(y,err)
460436c2adSjeremylt      call ceeddestroy(ceed,err)
470436c2adSjeremylt
480436c2adSjeremylt      end
490436c2adSjeremylt!-----------------------------------------------------------------------
50