graeme-winter
101 lines
3.5 KiB
Fortran
101 lines
3.5 KiB
Fortran
! Example test program for existing external library
|
|
! This should be saved in a file called test_generic_host.f90
|
|
! Kay Diederichs 4/2017
|
|
!
|
|
! compile with
|
|
! ifort -qopenmp generic_data_plugin.f90 test_generic_host.f90 -o test_generic_host
|
|
! or
|
|
! gfortran -O -fopenmp -ldl generic_data_plugin.f90 test_generic.f90 -o test_generic_host
|
|
! run with
|
|
! ./test_generic_host < test.in
|
|
! To test the dectris-neggia library, one could use this test.in:
|
|
!/usr/local/lib64/dectris-neggia.so
|
|
!/scratch/data/Eiger_16M_Nov2015/2015_11_10/insu6_1_??????.h5
|
|
!1 900
|
|
!
|
|
! The OMP_NUM_THREADS environment variable may be used for benchmarks!
|
|
|
|
|
|
PROGRAM test_generic_host
|
|
USE generic_data_plugin, ONLY: library, firstqm, lastqm, nx, ny, is_open, &
|
|
generic_open, generic_get_header, generic_get_data, generic_close
|
|
IMPLICIT NONE
|
|
INTEGER :: ier,nxny,ilow,ihigh,nbyte,info_array(1024), &
|
|
number_of_frames,len,numfrm
|
|
INTEGER, ALLOCATABLE :: iframe(:)
|
|
REAL :: qx,qy,avgcounts
|
|
CHARACTER(len=:), ALLOCATABLE :: master_file
|
|
CHARACTER(len=512) :: ACTNAM
|
|
|
|
! what should be done?
|
|
WRITE(*,*)'enter parameter of LIB= keyword:'
|
|
READ(*,'(a)') actnam
|
|
library=TRIM(actnam)
|
|
WRITE(*,*)'enter parameter of NAME_TEMPLATE_OF_DATA_FRAMES= keyword:'
|
|
READ(*,'(a)') actnam
|
|
WRITE(*,*)'enter parameters of the DATA_RANGE= keyword:'
|
|
READ(*,*) ilow,ihigh
|
|
|
|
! set some more module variables
|
|
firstqm=INDEX(actnam,'?') ! qm means question mark
|
|
lastqm =INDEX(actnam,'?',BACK=.TRUE.)
|
|
len =LEN_TRIM(actnam)
|
|
IF (actnam(len-2:len)=='.h5')THEN
|
|
master_file=actnam(:len-9)//'master.h5'
|
|
PRINT*,'master_file=',TRIM(master_file)
|
|
ELSE
|
|
master_file=TRIM(actnam)
|
|
ENDIF
|
|
info_array(1) = 1 ! 1=XDS (generic_open may check this)
|
|
info_array(2) = 123456789 ! better: e.g. 20160510; generic_open may check this
|
|
|
|
! initialize
|
|
CALL generic_open(library, master_file,info_array, ier)
|
|
IF (ier<0) THEN
|
|
WRITE(*,*)'error from generic_open, ier=',ier
|
|
STOP
|
|
END IF
|
|
is_open=.TRUE.
|
|
|
|
! get header and report
|
|
CALL generic_get_header(nx,ny,nbyte,qx,qy,number_of_frames,info_array,ier)
|
|
IF (ier<0) THEN
|
|
WRITE(*,*)'error from generic_get_header, ier=',ier
|
|
STOP
|
|
END IF
|
|
WRITE(*,'(a,3i6,2f10.6,i6)')'nx,ny,nbyte,qx,qy,number_of_frames=', &
|
|
nx,ny,nbyte,qx,qy,number_of_frames
|
|
WRITE(*,'(a,4i4,i12)')'INFO(1:5)=vendor/major version/minor version/patch/timestamp=', &
|
|
info_array(1:5)
|
|
IF (info_array(1)==0) THEN
|
|
WRITE(*,*) 'generic_getfrm: data are not vendor-specific',info_array(1) ! 1=Dectris
|
|
ELSE IF (info_array(1)==1) THEN
|
|
WRITE(*,*) 'generic_getfrm: data are from Dectris'
|
|
END IF
|
|
nxny=nx*ny
|
|
avgcounts=0.
|
|
|
|
! read the data (possibly in parallel)
|
|
!$omp parallel default(shared) private(numfrm,iframe,info_array,ier)
|
|
ALLOCATE(iframe(nxny))
|
|
!$omp do reduction(+:avgcounts)
|
|
DO numfrm=ilow,ihigh
|
|
CALL generic_get_data(numfrm, nx, ny, iframe, info_array, ier)
|
|
IF (ier<0) THEN
|
|
WRITE(*,*)'error from generic_get_data, numfrm, ier=',numfrm,ier
|
|
STOP
|
|
END IF
|
|
avgcounts=avgcounts + SUM(iframe)/REAL(nxny) ! do something with data
|
|
END DO
|
|
!$omp end parallel
|
|
WRITE(*,*)'average counts:',avgcounts/(ihigh-ilow+1)
|
|
|
|
! close
|
|
CALL generic_close(ier)
|
|
IF (ier<0) THEN
|
|
WRITE(*,*)'error from generic_close, ier=',ier
|
|
STOP
|
|
END IF
|
|
|
|
END PROGRAM test_generic_host
|