mpi_hyperslab: check that filters are available first

Author: scivision

src/interface.f90

@@ -13,9 +13,9 @@ module h5mpi
 H5P_DEFAULT_F, H5P_FILE_ACCESS_F, H5P_DATASET_CREATE_F, H5P_DATASET_XFER_F, &
 H5S_ALL_F, H5S_SELECT_SET_F, &
 H5D_CHUNKED_F, H5D_CONTIGUOUS_F, H5D_COMPACT_F, &
-h5dcreate_f, h5dopen_f, h5dclose_f, h5dget_space_f, &
+h5dcreate_f, h5dopen_f, h5dclose_f, h5dget_space_f, h5dget_create_plist_f, &
 h5fopen_f, h5fclose_f, h5fcreate_f, h5fget_filesize_f, h5fflush_f, &
-h5pcreate_f, h5pclose_f, h5pset_chunk_f, h5pset_dxpl_mpio_f, h5pset_fapl_mpio_f, &
+h5pcreate_f, h5pclose_f, h5pset_chunk_f, h5pset_dxpl_mpio_f, h5pset_fapl_mpio_f, h5pall_filters_avail_f, &
 h5sselect_hyperslab_f, h5screate_simple_f, h5sclose_f, &
 h5get_libversion_f, &
 h5open_f, h5close_f
@@ -456,18 +456,36 @@ subroutine mpi_hyperslab(mem_dims, dset_dims, dset_id, filespace, memspace, dnam
 integer(HSIZE_T), dimension(size(istart)), intent(in) :: iend
 
 integer(HSIZE_T), dimension(size(mem_dims)) :: c_mem_dims, i0
+integer(HID_T) :: dcpl
 integer :: ierr
 
+logical :: filters_OK
+
+!> check that all necessary filters to access dataset are available on the system.
+call h5dget_create_plist_f(dset_id, dcpl, ierr)
+if (ierr/=0) error stop "h5fortran:mpi_hyperslab:h5dget_create_plist: " // dname
+
+call h5pall_filters_avail_f(dcpl, filters_OK, ierr)
+if (ierr/=0) error stop "h5fortran:mpi_hyperslab:h5pall_filters_avail: " // dname
+if (.not. filters_OK) then
+  error stop "h5fortran: filter(s) missing necessary for dataset " // dname // " in parallel with MPI. This is " // &
+    "typically caused by missing DEFLATE compression with HDF5-MPI."
+endif
+
+call h5pclose_f(dcpl, ierr)
+if(ierr/=0) error stop "h5fortran:mpi_hyperslab:h5pclose: " // dname
 
 if(filespace == H5S_ALL_F) then
   !> create dataspace
   call h5screate_simple_f(rank=size(dset_dims), dims=dset_dims, space_id=filespace, hdferr=ierr)
-  if (ierr/=0) error stop "h5screate_simple:filespace " // dname
+  if (ierr/=0) error stop "h5fortran:mpi_hyperslab:h5screate_simple:filespace " // dname
 endif
 
 !> Select hyperslab in the file.
 call h5dget_space_f(dset_id, filespace, ierr)
-if (ierr/=0) error stop "h5dget_space: " // dname
+if (ierr/=0) error stop "h5fortran:mpi_hyperslab:h5dget_space: " // dname
+
+
 
 ! blk(1) = 1
 ! blk(2:) = dset_dims(2:)
@@ -493,11 +511,11 @@ subroutine mpi_hyperslab(mem_dims, dset_dims, dset_id, filespace, memspace, dnam
 ! stride=1, &  !< for now we don't stride data
 ! block=blk  !< would this help performance?
 
-if (ierr/=0) error stop "h5sselect_hyperslab: " // dname
+if (ierr/=0) error stop "g5fortran:mpi_hyperslab:h5sselect_hyperslab: " // dname
 
 !> create memory dataspace
 call h5screate_simple_f(rank=size(c_mem_dims), dims=c_mem_dims, space_id=memspace, hdferr=ierr)
-if (ierr/=0) error stop "h5fortran:h5screate_simple:memspace " // dname
+if (ierr/=0) error stop "h5fortran:mpi_hyperslab:h5screate_simple:memspace " // dname
 
 end subroutine mpi_hyperslab
 

src/read/read.f90

@@ -1,7 +1,7 @@
 submodule (h5mpi) hdf5_read
 
 use hdf5, only : h5tget_native_type_f, h5tget_class_f, h5dget_type_f, h5tget_size_f, h5tclose_f, h5topen_f, &
-h5dread_f, h5dget_create_plist_f, &
+h5dread_f, &
 h5pget_nfilters_f, h5pget_filter_f, h5pget_layout_f, h5pget_chunk_f, &
 H5T_DIR_ASCEND_F, &
 H5Z_FILTER_DEFLATE_F

src/tests/unit/test_deflate_read.f90

@@ -64,11 +64,7 @@ subroutine test_read_deflate(fn, N, mpi_id, Nmpi)
 
 call h5f%open(fn, action='r', mpi=.true.)
 call h5f%read('/A', A, istart=i0, iend=i1)
-if (h5f%parallel_compression) then
-  call h5f%read('/noMPI', A, istart=i0, iend=i1)
-else
-  print *, "test_read_deflate: no parallel compression: " // fn
-endif
+call h5f%read('/noMPI', A, istart=i0, iend=i1)  !< compressed with MPI
 call h5f%close()
 
 end subroutine test_read_deflate