xref: /libCEED/tests/t103-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      integer*8 aoffset,xoffset,yoffset
90436c2adSjeremylt      real*8 a(10)
100436c2adSjeremylt      real*8 xx(10)
110436c2adSjeremylt      real*8 yy(10)
120436c2adSjeremylt      real*8 diff
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      call ceedvectorcreate(ceed,n,x,err)
210436c2adSjeremylt      call ceedvectorcreate(ceed,n,y,err)
220436c2adSjeremylt
230436c2adSjeremylt      do i=1,10
240436c2adSjeremylt        a(i)=10+i
250436c2adSjeremylt      enddo
260436c2adSjeremylt
270436c2adSjeremylt      aoffset=0
280436c2adSjeremylt      call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,a,aoffset,err)
290436c2adSjeremylt
300436c2adSjeremylt      call ceedvectorgetarray(x,ceed_mem_host,xx,xoffset,err)
310436c2adSjeremylt      call ceedvectorsetarray(y,ceed_mem_host,ceed_copy_values,xx,xoffset,err)
320436c2adSjeremylt      call ceedvectorrestorearray(x,xx,xoffset,err)
330436c2adSjeremylt
340436c2adSjeremylt      call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err)
350436c2adSjeremylt
360436c2adSjeremylt      do i=1,10
370436c2adSjeremylt        diff=yy(i+yoffset)-10-i
380436c2adSjeremylt        if (abs(diff)>1.0D-15) then
390436c2adSjeremylt! LCOV_EXCL_START
400436c2adSjeremylt          write(*,*) 'Error reading array y(',i,')=',yy(i+yoffset)
410436c2adSjeremylt! LCOV_EXCL_STOP
420436c2adSjeremylt        endif
430436c2adSjeremylt      enddo
440436c2adSjeremylt
450436c2adSjeremylt      call ceedvectorrestorearrayread(y,yy,yoffset,err)
460436c2adSjeremylt      call ceedvectordestroy(x,err)
470436c2adSjeremylt      call ceedvectordestroy(y,err)
480436c2adSjeremylt      call ceeddestroy(ceed,err)
490436c2adSjeremylt
500436c2adSjeremylt      end
510436c2adSjeremylt!-----------------------------------------------------------------------
52