223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

test more dimensions

This commit is contained in:
Martin Diehl 2023-11-26 12:38:35 +01:00
parent 3c62af3fe5
commit 1cf1f9df31
1 changed files with 21 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/test/test_HDF5_utilities.f90
View File

@ -22,21 +22,36 @@ subroutine read_write()
integer(HID_T) :: f integer(HID_T) :: f
real(pREAL), dimension(3) :: d1_in,d1_out real(pREAL), dimension(3) :: d1_in,d1_out
real(pREAL), dimension(3,3) :: d2_in,d2_out real(pREAL), dimension(3,3) :: d2_in,d2_out
real(pREAL), dimension(3,3,3) :: d3_in,d3_out
real(pREAL), dimension(3,3,3,3) :: d4_in,d4_out
real(pREAL), dimension(3,3,3,3,3) :: d5_in,d5_out
call random_number(d1_in) call random_number(d1_in)
call random_number(d2_in) call random_number(d2_in)
call random_number(d3_in)
call random_number(d4_in)
call random_number(d5_in)
f = HDF5_openFile('test.hdf5','w') f = HDF5_openFile('test.hdf5','w')
call HDF5_write(d1_in,f,'d1') call HDF5_write(d1_in,f,'d1')
call HDF5_write(d2_in,f,'d2') call HDF5_write(d2_in,f,'d2')
call HDF5_write(d3_in,f,'d3')
call HDF5_write(d4_in,f,'d4')
call HDF5_write(d5_in,f,'d5')
call HDF5_read(d1_out,f,'d1') call HDF5_read(d1_out,f,'d1')
call HDF5_read(d2_out,f,'d2') call HDF5_read(d2_out,f,'d2')
call HDF5_read(d3_out,f,'d3')
call HDF5_read(d4_out,f,'d4')
call HDF5_read(d5_out,f,'d5')
if (any(d1_in /= d1_out)) error stop 'test_read_write(w)/d1' if (any(d1_in /= d1_out)) error stop 'test_read_write(w)/d1'
if (any(d2_in /= d2_out)) error stop 'test_read_write(w)/d2' if (any(d2_in /= d2_out)) error stop 'test_read_write(w)/d2'
if (any(d3_in /= d3_out)) error stop 'test_read_write(w)/d3'
if (any(d4_in /= d4_out)) error stop 'test_read_write(w)/d4'
if (any(d5_in /= d5_out)) error stop 'test_read_write(w)/d5'
call HDF5_closeFile(f) call HDF5_closeFile(f)
@ -45,9 +60,15 @@ subroutine read_write()
call HDF5_read(d1_out,f,'d1') call HDF5_read(d1_out,f,'d1')
call HDF5_read(d2_out,f,'d2') call HDF5_read(d2_out,f,'d2')
call HDF5_read(d3_out,f,'d3')
call HDF5_read(d4_out,f,'d4')
call HDF5_read(d5_out,f,'d5')
if (any(d1_in /= d1_out)) error stop 'test_read_write(r)/d1' if (any(d1_in /= d1_out)) error stop 'test_read_write(r)/d1'
if (any(d2_in /= d2_out)) error stop 'test_read_write(r)/d2' if (any(d2_in /= d2_out)) error stop 'test_read_write(r)/d2'
if (any(d3_in /= d3_out)) error stop 'test_read_write(r)/d3'
if (any(d4_in /= d4_out)) error stop 'test_read_write(r)/d4'
if (any(d5_in /= d5_out)) error stop 'test_read_write(r)/d5'
call HDF5_closeFile(f) call HDF5_closeFile(f)