!****h* root/fortran/test/TH5P_F03 ! ! NAME ! tH5P_F03.F90 ! ! FUNCTION ! Test FORTRAN HDF5 H5P APIs which are dependent on FORTRAN 2003 ! features. ! ! COPYRIGHT ! * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ! Copyright by The HDF Group. * ! All rights reserved. * ! * ! This file is part of HDF5. The full HDF5 copyright notice, including * ! terms governing use, modification, and redistribution, is contained in * ! the COPYING file, which can be found at the root of the source code * ! distribution tree, or in https://www.hdfgroup.org/licenses. * ! If you do not have access to either file, you may request a copy from * ! help@hdfgroup.org. * ! * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ! ! USES ! test_genprop_cls_cb1_mod ! ! CONTAINS SUBROUTINES ! test_create, test_genprop_class_callback ! !***** ! ***************************************** ! *** H 5 P T E S T S ! ***************************************** MODULE test_genprop_cls_cb1_mod ! Callback subroutine for test_genprop_class_callback ! and the function H5Pcreate_class_f. USE HDF5 USE ISO_C_BINDING IMPLICIT NONE TYPE, BIND(C) :: cop_cb_struct_ ! Struct for iterations INTEGER :: count INTEGER(HID_T) :: id END TYPE cop_cb_struct_ CONTAINS INTEGER FUNCTION test_genprop_cls_cb1_f(list_id, create_data ) bind(C) IMPLICIT NONE INTEGER(HID_T), INTENT(IN), VALUE :: list_id TYPE(cop_cb_struct_) :: create_data create_data%count = create_data%count + 1 create_data%id = list_id test_genprop_cls_cb1_f = 0 END FUNCTION test_genprop_cls_cb1_f END MODULE test_genprop_cls_cb1_mod MODULE TH5P_F03 USE HDF5 USE TH5_MISC USE TH5_MISC_GEN USE ISO_C_BINDING CONTAINS !------------------------------------------------------------------------- ! * Function: test_create ! * ! * Purpose: Tests H5Pset_fill_value_f and H5Pget_fill_value_f ! * ! * Return: Success: 0 ! * ! * Failure: number of errors ! *------------------------------------------------------------------------- ! SUBROUTINE test_create(total_error) IMPLICIT NONE INTEGER, INTENT(INOUT) :: total_error INTEGER(HID_T) :: fapl INTEGER(hid_t) :: file=-1, space=-1, dcpl=-1, comp_type_id=-1 INTEGER(hid_t) :: dset9=-1 INTEGER(hsize_t), DIMENSION(1:5), PARAMETER :: cur_size = (/2, 8, 8, 4, 2/) INTEGER(hsize_t), DIMENSION(1:5), PARAMETER :: ch_size= (/1, 1, 1, 4, 1/) CHARACTER(LEN=14) :: filename ='test_create.h5' TYPE(comp_datatype), TARGET :: rd_c, fill_ctype INTEGER :: error INTEGER(SIZE_T) :: h5off TYPE(C_PTR) :: f_ptr CHARACTER(LEN=1) :: cfill INTEGER :: ifill REAL :: rfill REAL(KIND=dp) :: dpfill INTEGER :: low, high ! ! * Create a file. ! CALL h5fcreate_f(filename,H5F_ACC_TRUNC_F,file,error) CALL check("h5fcreate_f", error, total_error) CALL h5screate_simple_f(5, cur_size, space, error, cur_size) CALL check("h5screate_simple_f", error, total_error) CALL H5Pcreate_f(H5P_DATASET_CREATE_F, dcpl, error) CALL check("H5Pcreate_f", error, total_error) CALL h5pset_chunk_f(dcpl, 5, ch_size, error) CALL check("h5pset_chunk_f",error, total_error) ! Create a compound datatype CALL h5tcreate_f(H5T_COMPOUND_F, H5_SIZEOF(fill_ctype), comp_type_id, error) CALL check("h5tcreate_f", error, total_error) h5off = H5OFFSETOF(C_LOC(fill_ctype), C_LOC(fill_ctype%a)) CALL h5tinsert_f(comp_type_id, "a", h5off , H5T_NATIVE_REAL, error) CALL check("h5tinsert_f", error, total_error) CALL h5tinsert_f(comp_type_id, "x", H5OFFSETOF(C_LOC(fill_ctype), C_LOC(fill_ctype%x)), H5T_NATIVE_INTEGER, error) CALL check("h5tinsert_f", error, total_error) CALL h5tinsert_f(comp_type_id, "y", H5OFFSETOF(C_LOC(fill_ctype), C_LOC(fill_ctype%y)), H5T_NATIVE_DOUBLE, error) CALL check("h5tinsert_f", error, total_error) CALL h5tinsert_f(comp_type_id, "z", & H5OFFSETOF(C_LOC(fill_ctype), C_LOC(fill_ctype%z)), H5T_NATIVE_CHARACTER, error) CALL check("h5tinsert_f", error, total_error) CALL H5Pset_alloc_time_f(dcpl, H5D_ALLOC_TIME_LATE_F,error) CALL check("H5Pset_alloc_time_f",error, total_error) CALL H5Pset_fill_time_f(dcpl, H5D_FILL_TIME_ALLOC_F, error) CALL check("H5Pset_fill_time_f",error, total_error) ! Compound datatype test f_ptr = C_LOC(fill_ctype) CALL H5Pget_fill_value_f(dcpl, comp_type_id, f_ptr, error) CALL check("H5Pget_fill_value_f",error, total_error) fill_ctype%y = 4444.D0 fill_ctype%z = 'S' fill_ctype%a = 5555. fill_ctype%x = 55 ! Test various fill values CALL H5Pset_fill_value_f(dcpl, H5T_NATIVE_CHARACTER, 'X', error) CALL check("H5Pset_fill_value_f",error, total_error) CALL h5pget_fill_value_f(dcpl, H5T_NATIVE_CHARACTER, cfill, error) CALL check("H5Pget_fill_value_f",error, total_error) IF(cfill.NE.'X')THEN PRINT*,"***ERROR: Returned wrong fill value (character)" total_error = total_error + 1 ENDIF CALL H5Pset_fill_value_f(dcpl, H5T_NATIVE_INTEGER, 9, error) CALL check("H5Pset_fill_value_f",error, total_error) CALL h5pget_fill_value_f(dcpl, H5T_NATIVE_INTEGER, ifill, error) CALL check("H5Pget_fill_value_f",error, total_error) IF(ifill.NE.9)THEN PRINT*,"***ERROR: Returned wrong fill value (integer)" total_error = total_error + 1 ENDIF CALL H5Pset_fill_value_f(dcpl, H5T_NATIVE_DOUBLE, 1.0_dp, error) CALL check("H5Pset_fill_value_f",error, total_error) CALL h5pget_fill_value_f(dcpl, H5T_NATIVE_DOUBLE, dpfill, error) CALL check("H5Pget_fill_value_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong fill value (double)", dpfill, 1.0_dp, total_error) CALL H5Pset_fill_value_f(dcpl, H5T_NATIVE_REAL, 2.0, error) CALL check("H5Pset_fill_value_f",error, total_error) CALL h5pget_fill_value_f(dcpl, H5T_NATIVE_REAL, rfill, error) CALL check("H5Pget_fill_value_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong fill value (real)", rfill, 2.0, total_error) ! For the actual compound type f_ptr = C_LOC(fill_ctype) CALL H5Pset_fill_value_f(dcpl, comp_type_id, f_ptr, error) CALL check("H5Pget_fill_value_f",error, total_error) CALL h5dcreate_f(file,"dset9", comp_type_id, space, dset9, error, dcpl_id=dcpl) CALL check("h5dcreate_f", error, total_error) CALL h5dclose_f(dset9, error) CALL check("h5dclose_f", error, total_error) CALL h5fclose_f(file,error) CALL check("h5fclose_f", error, total_error) ! Open the file and get the dataset fill value from each dataset CALL h5pcreate_f(H5P_FILE_ACCESS_F, fapl, error) CALL check("h5pcreate_f",error, total_error) CALL h5pset_libver_bounds_f(fapl, H5F_LIBVER_V18_F, H5F_LIBVER_V18_F, error) CALL check("h5pset_libver_bounds_f",error, total_error) CALL h5pget_libver_bounds_f(fapl, low, high, error) CALL check("h5pget_libver_bounds_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong low libver_bounds", low, H5F_LIBVER_V18_F, total_error) CALL VERIFY("***ERROR: Returned wrong high libver_bounds", high, H5F_LIBVER_V18_F, total_error) CALL h5pset_libver_bounds_f(fapl, H5F_LIBVER_V18_F, H5F_LIBVER_V110_F, error) CALL check("h5pset_libver_bounds_f",error, total_error) CALL h5pget_libver_bounds_f(fapl, low, high, error) CALL check("h5pget_libver_bounds_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong low libver_bounds", low, H5F_LIBVER_V18_F, total_error) CALL VERIFY("***ERROR: Returned wrong high libver_bounds", high, H5F_LIBVER_V110_F, total_error) CALL h5pset_libver_bounds_f(fapl, H5F_LIBVER_V112_F, H5F_LIBVER_V112_F, error) CALL check("h5pset_libver_bounds_f",error, total_error) CALL h5pget_libver_bounds_f(fapl, low, high, error) CALL check("h5pget_libver_bounds_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong low libver_bounds", low, H5F_LIBVER_V112_F, total_error) CALL VERIFY("***ERROR: Returned wrong high libver_bounds", high, H5F_LIBVER_V112_F, total_error) CALL h5pset_libver_bounds_f(fapl, H5F_LIBVER_V114_F, H5F_LIBVER_V114_F, error) CALL check("h5pset_libver_bounds_f",error, total_error) CALL h5pget_libver_bounds_f(fapl, low, high, error) CALL check("h5pget_libver_bounds_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong low libver_bounds", low, H5F_LIBVER_V114_F, total_error) CALL VERIFY("***ERROR: Returned wrong high libver_bounds", high, H5F_LIBVER_V114_F, total_error) CALL H5Pset_libver_bounds_f(fapl, H5F_LIBVER_LATEST_F, H5F_LIBVER_LATEST_F, error) CALL check("H5Pset_libver_bounds_f",error, total_error) CALL h5pget_libver_bounds_f(fapl, low, high, error) CALL check("h5pget_libver_bounds_f",error, total_error) CALL VERIFY("***ERROR: Returned wrong low libver_bounds", low, H5F_LIBVER_LATEST_F, total_error) CALL VERIFY("***ERROR: Returned wrong high libver_bounds", high, H5F_LIBVER_LATEST_F, total_error) IF(H5F_LIBVER_LATEST_F.GE.H5F_LIBVER_NBOUNDS_F)THEN WRITE(*,'(A,I0,A,I0,A)') & "***ERROR: H5F_LIBVER_LATEST_F (",H5F_LIBVER_LATEST_F,") .GE. H5F_LIBVER_NBOUNDS_F (",H5F_LIBVER_NBOUNDS_F,")" total_error = total_error + 1 ENDIF CALL h5fopen_f (FILENAME, H5F_ACC_RDONLY_F, file, error, fapl) CALL check("h5fopen_f", error, total_error) ! Compound datatype test CALL h5dopen_f(file, "dset9", dset9, error) CALL check("h5dopen_f", error, total_error) CALL H5Pclose_f(dcpl, error) CALL check("H5Pclose_f", error, total_error) CALL H5Dget_create_plist_f(dset9, dcpl, error) CALL check("H5Dget_create_plist_f", error, total_error) f_ptr = C_LOC(rd_c) CALL H5Pget_fill_value_f(dcpl, comp_type_id, f_ptr, error) CALL check("H5Pget_fill_value_f", error, total_error) CALL verify("***ERROR: Returned wrong fill value", rd_c%a, fill_ctype%a, total_error) CALL verify("***ERROR: Returned wrong fill value", rd_c%y, fill_ctype%y, total_error) IF( rd_c%x .NE. fill_ctype%x .OR. & rd_c%z .NE. fill_ctype%z )THEN PRINT*,"***ERROR: Returned wrong fill value" total_error = total_error + 1 ENDIF CALL h5dclose_f(dset9, error) CALL check("h5dclose_f", error, total_error) CALL H5Pclose_f(dcpl, error) CALL check("H5Pclose_f", error, total_error) CALL h5fclose_f(file,error) CALL check("h5fclose_f", error, total_error) END SUBROUTINE test_create SUBROUTINE test_genprop_class_callback(total_error) ! ! ! test_genprop_class_callback(): Test basic generic property list code. ! Tests callbacks for property lists in a generic class. ! ! FORTRAN TESTS: ! Tests function H5Pcreate_class_f with callback. ! ! USE test_genprop_cls_cb1_mod IMPLICIT NONE INTEGER, INTENT(INOUT) :: total_error INTEGER(hid_t) :: cid1, cid2 ! Generic Property class ID INTEGER(hid_t) :: lid1, lid2 ! Generic Property list ID INTEGER(size_t) :: nprops ! Number of properties in class TYPE(cop_cb_struct_), TARGET :: crt_cb_struct, cls_cb_struct INTEGER :: CLASS1_NAME_SIZE = 7 ! length of class string CHARACTER(LEN=7) :: CLASS1_NAME = "Class 1", CLASS1_NAME_BUF TYPE(C_FUNPTR) :: f1, f5 TYPE(C_PTR) :: f2, f6 CHARACTER(LEN=10) :: PROP1_NAME = "Property 1" INTEGER(SIZE_T) :: PROP1_SIZE = 10 CHARACTER(LEN=10) :: PROP2_NAME = "Property 2" INTEGER(SIZE_T) :: PROP2_SIZE = 10 CHARACTER(LEN=10) :: PROP3_NAME = "Property 3" INTEGER(SIZE_T) :: PROP3_SIZE = 10 CHARACTER(LEN=10) :: PROP4_NAME = "Property 4" INTEGER(SIZE_T) :: PROP4_SIZE = 10 INTEGER :: PROP1_DEF_VALUE = 10 INTEGER :: PROP2_DEF_VALUE = 10 INTEGER :: PROP3_DEF_VALUE = 10 INTEGER :: PROP4_DEF_VALUE = 10 INTEGER :: error ! Generic RETURN value LOGICAL :: flag ! for tests f1 = C_FUNLOC(test_genprop_cls_cb1_f) f5 = C_FUNLOC(test_genprop_cls_cb1_f) f2 = C_LOC(crt_cb_struct) f6 = C_LOC(cls_cb_struct) ! Create a new generic class, derived from the root of the class hierarchy CALL h5pcreate_class_f(h5p_ROOT_F, CLASS1_NAME, cid1, error, f1, f2, c_null_funptr, c_null_ptr, f5, f6) CALL check("h5pcreate_class_f", error, total_error) ! Insert first property into class (with no callbacks) CALL h5pregister_f(cid1, PROP1_NAME, PROP1_SIZE, PROP1_DEF_VALUE, error) CALL check("h5pregister_f", error, total_error) ! Insert second property into class (with no callbacks) CALL h5pregister_f(cid1, PROP2_NAME, PROP2_SIZE, PROP2_DEF_VALUE, error) CALL check("h5pregister_f", error, total_error) ! Insert third property into class (with no callbacks) CALL h5pregister_f(cid1, PROP3_NAME, PROP3_SIZE, PROP3_DEF_VALUE, error) CALL check("h5pregister_f", error, total_error) ! Insert fourth property into class (with no callbacks) CALL h5pregister_f(cid1, PROP4_NAME, PROP4_SIZE, PROP4_DEF_VALUE, error) CALL check("h5pregister_f", error, total_error) ! Check the number of properties in class CALL h5pget_nprops_f(cid1, nprops, error) CALL check("h5pget_nprops_f", error, total_error) CALL verify("h5pget_nprops_f", INT(nprops), 4, total_error) ! Initialize class callback structs crt_cb_struct%count = 0 crt_cb_struct%id = -1 cls_cb_struct%count = 0 cls_cb_struct%id = -1 ! Create a property list from the class CALL h5pcreate_f(cid1, lid1, error) CALL check("h5pcreate_f", error, total_error) ! Get the list's class CALL H5Pget_class_f(lid1, cid2, error) CALL check("H5Pget_class_f", error, total_error) ! Check that the list's class is correct CALL H5Pequal_f(cid2, cid1, flag, error) CALL check("H5Pequal_f", error, total_error) CALL verify("H5Pequal_f", flag, .TRUE., total_error) ! Check the class name CALL H5Pget_class_name_f(cid2, CLASS1_NAME_BUF, CLASS1_NAME_SIZE, error) CALL check("H5Pget_class_name_f", error, total_error) CALL verify("H5Pget_class_name_f", CLASS1_NAME_BUF, CLASS1_NAME, error) IF(error.NE.0)THEN WRITE(*,*) 'Class names do not match! name=',CLASS1_NAME_BUF, 'CLASS1_NAME=',CLASS1_NAME total_error = total_error + 1 ENDIF ! Close class CALL h5pclose_class_f(cid2, error) CALL check("h5pclose_class_f", error, total_error) ! Verify that the creation callback occurred CALL verify("h5pcreate_f", crt_cb_struct%count, 1, total_error) CALL verify("h5pcreate_f", crt_cb_struct%id, lid1, total_error) ! Check the number of properties in list CALL h5pget_nprops_f(lid1,nprops, error) CALL check("h5pget_nprops_f", error, total_error) CALL verify("h5pget_nprops_f", INT(nprops), 4, total_error) ! Create another property list from the class CALL h5pcreate_f(cid1, lid2, error) CALL check("h5pcreate_f", error, total_error) ! Verify that the creation callback occurred CALL verify("h5pcreate_f", crt_cb_struct%count, 2, total_error) CALL verify("h5pcreate_f", crt_cb_struct%id, lid2, total_error) ! Check the number of properties in list CALL h5pget_nprops_f(lid2,nprops, error) CALL check("h5pget_nprops_f", error, total_error) CALL verify("h5pget_nprops_f", INT(nprops), 4, total_error) ! Close first list CALL h5pclose_f(lid1, error); CALL check("h5pclose_f", error, total_error) ! Verify that the close callback occurred CALL verify("h5pcreate_f", cls_cb_struct%count, 1, total_error) CALL verify("h5pcreate_f", cls_cb_struct%id, lid1, total_error) ! Close second list CALL h5pclose_f(lid2, error); CALL check("h5pclose_f", error, total_error) ! Verify that the close callback occurred CALL verify("h5pcreate_f", cls_cb_struct%count, 2, total_error) CALL verify("h5pcreate_f", cls_cb_struct%id, lid2, total_error) ! Close class CALL h5pclose_class_f(cid1, error) CALL check("h5pclose_class_f", error, total_error) END SUBROUTINE test_genprop_class_callback !------------------------------------------------------------------------- ! Function: test_h5p_file_image ! ! Purpose: Tests APIs: ! h5pget_file_image_f and h5pset_file_image_f ! ! Return: Success: 0 ! Failure: -1 !------------------------------------------------------------------------- SUBROUTINE test_h5p_file_image(total_error) IMPLICIT NONE INTEGER, INTENT(INOUT) :: total_error INTEGER(hid_t) :: fapl_1 = -1 INTEGER, PARAMETER :: count = 10 INTEGER, DIMENSION(1:count), TARGET :: buffer INTEGER, DIMENSION(1:count), TARGET :: temp INTEGER :: i INTEGER(size_t) :: size INTEGER(size_t) :: temp_size INTEGER :: error ! error return value TYPE(C_PTR) :: f_ptr TYPE(C_PTR), DIMENSION(1:count) :: f_ptr1 TYPE(C_PTR), DIMENSION(1:1) :: f_ptr2 ! Initialize file image buffer DO i = 1, count buffer(i) = i*10 ENDDO ! Create fapl CALL h5pcreate_f(H5P_FILE_ACCESS_F, fapl_1, error) CALL check("h5pcreate_f", error, total_error) ! Test with NULL ptr f_ptr2(1) = C_NULL_PTR temp_size = 1 CALL h5pget_file_image_f(fapl_1, f_ptr2, temp_size, error) CALL check("h5pget_file_image_f", error, total_error) CALL verify("h5pget_file_image_f", INT(temp_size), 0, total_error) ! Set file image f_ptr = C_LOC(buffer(1)) size = H5_SIZEOF(buffer(1))*count CALL h5pset_file_image_f(fapl_1, f_ptr, size, error) CALL check("h5pset_file_image_f", error, total_error) ! Get the same data back DO i = 1, count f_ptr1(i) = C_LOC(temp(i)) ENDDO temp_size = 0 CALL h5pget_file_image_f(fapl_1, f_ptr1, temp_size, error) CALL check("h5pget_file_image_f", error, total_error) ! Check that sizes are the same, and that the buffers are identical but separate CALL verify("h5pget_file_image_f", INT(temp_size), INT(size), total_error) ! Verify the image data is correct DO i = 1, count CALL verify("h5pget_file_image_f", temp(i), buffer(i), total_error) ENDDO END SUBROUTINE test_h5p_file_image !------------------------------------------------------------------------- ! Function: external_test_offset ! ! Purpose: Tests APIs: ! h5pset_external_f (with offsets not equal to zero), h5pget_external_f ! ! Return: Success: 0 ! Failure: -1 !------------------------------------------------------------------------- ! SUBROUTINE external_test_offset(cleanup,total_error) IMPLICIT NONE INTEGER, INTENT(INOUT) :: total_error LOGICAL, INTENT(IN) :: cleanup INTEGER(hid_t) :: fapl=-1 ! file access property list INTEGER(hid_t) :: file=-1 ! file to write to INTEGER(hid_t) :: dcpl=-1 ! dataset creation properties INTEGER(hid_t) :: space=-1 ! data space INTEGER(hid_t) :: dset=-1 ! dataset INTEGER(hid_t) :: grp=-1 ! group to emit diagnostics INTEGER(size_t) :: i, j ! miscellaneous counters INTEGER :: k CHARACTER(LEN=180) :: filename ! file names INTEGER, DIMENSION(1:25) :: part INTEGER, DIMENSION(1:100), TARGET :: whole ! raw data buffers INTEGER(hsize_t), DIMENSION(1:1) :: cur_size ! current data space size INTEGER(hid_t) :: hs_space ! hyperslab data space INTEGER(hsize_t), DIMENSION(1:1) :: hs_start = (/30/) ! hyperslab starting offset INTEGER(hsize_t), DIMENSION(1:1) :: hs_count = (/25/) ! hyperslab size CHARACTER(LEN=1) :: ichr1 ! character conversion holder INTEGER :: error ! error status TYPE(C_PTR) :: f_ptr ! fortran pointer INTEGER(HSIZE_T) :: sizeof_part CHARACTER(LEN=1,KIND=C_CHAR), DIMENSION(1:30) :: temparray temparray(1:30)(1:1) = '0' ! 1 byte character ! Write the data to external files directly DO i = 1, 4 DO j = 1, 25 part(j) = INT((i-1_size_t)*25_size_t+(j-1_size_t)) ENDDO WRITE(ichr1,'(I1.1)') i filename = "extern_"//ichr1//"a.raw" OPEN(10, FILE=filename, ACCESS='STREAM', form='UNFORMATTED') WRITE(10) temparray(1:(i-1)*10) WRITE(10) part CLOSE(10) ENDDO ! ! Create the file and an initial group. CALL h5pcreate_f(H5P_FILE_ACCESS_F, fapl, error) CALL h5fcreate_f('extern_raw.h5', H5F_ACC_TRUNC_F, file, error, access_prp=fapl) CALL check("h5fcreate_f",error,total_error) CALL h5gcreate_f(file, "emit-diagnostics", grp, error) CALL check("h5gcreate_f",error, total_error) ! Create the dataset sizeof_part = INT(H5_SIZEOF(part(1))*25, hsize_t) CALL h5pcreate_f(H5P_DATASET_CREATE_F, dcpl, error) CALL check("h5pcreate_f", error, total_error) CALL h5pset_external_f(dcpl, "extern_1a.raw", INT(0,off_t), sizeof_part, error) CALL check("h5pset_external_f",error,total_error) CALL h5pset_external_f(dcpl, "extern_2a.raw", INT(10,off_t), sizeof_part, error) CALL check("h5pset_external_f",error,total_error) CALL h5pset_external_f(dcpl, "extern_3a.raw", INT(20,off_t), sizeof_part, error) CALL check("h5pset_external_f",error,total_error) CALL h5pset_external_f(dcpl, "extern_4a.raw", INT(30,off_t), sizeof_part, error) CALL check("h5pset_external_f",error,total_error) cur_size(1) = 100 CALL h5screate_simple_f(1, cur_size, space, error) CALL check("h5screate_simple_f", error, total_error) CALL h5dcreate_f(file, "dset1", H5T_NATIVE_INTEGER, space, dset,error,dcpl_id=dcpl) CALL check("h5dcreate_f", error, total_error) ! ! Read the entire dataset and compare with the original whole(:) = 0 f_ptr = C_LOC(whole(1)) CALL h5dread_f(dset, H5T_NATIVE_INTEGER, f_ptr, error, mem_space_id=space, file_space_id=space) CALL check("h5dread_f", error, total_error) DO k = 1, 100 CALL verify("h5dread_f", whole(k), k-1, error) IF(error .NE. 0)THEN WRITE(*,*) "Incorrect value(s) read." total_error = total_error + 1 EXIT ENDIF ENDDO ! ! Read the middle of the dataset CALL h5scopy_f(space, hs_space, error) CALL check("h5scopy_f", error, total_error) CALL h5sselect_hyperslab_f(hs_space, H5S_SELECT_SET_F, hs_start, hs_count, error) CALL check("h5sselect_hyperslab_f", error, total_error) whole(:) = 0 f_ptr = C_LOC(whole(1)) CALL h5dread_f(dset, H5T_NATIVE_INTEGER, f_ptr, error, mem_space_id=hs_space, file_space_id=hs_space) CALL check("h5dread_f", error, total_error) CALL h5sclose_f(hs_space, error) CALL check("h5sclose_f", error, total_error) DO k = INT(hs_start(1))+1, INT(hs_start(1)+hs_count(1)) CALL verify("h5dread_f", whole(k), k-1, error) IF(error .NE. 0)THEN WRITE(*,*) "Incorrect value(s) read." total_error = total_error + 1 EXIT ENDIF ENDDO CALL h5dclose_f(dset, error) CALL check("h5dclose_f", error, total_error) CALL h5pclose_f(dcpl, error) CALL check("h5pclose_f", error, total_error) CALL h5sclose_f(space, error) CALL check("h5sclose_f", error, total_error) CALL h5fclose_f(file, error) CALL check("h5fclose_f", error, total_error) ! cleanup DO i = 1, 4 WRITE(ichr1,'(I1.1)') i filename = "extern_"//ichr1//"a.raw" CALL h5_cleanup_f(filename, H5P_DEFAULT_F, error) CALL check("h5_cleanup_f", error, total_error) ENDDO IF(cleanup) CALL h5_cleanup_f("extern_raw.h5", H5P_DEFAULT_F, error) CALL check("h5_cleanup_f", error, total_error) END SUBROUTINE external_test_offset !------------------------------------------------------------------------- ! NAME ! test_vds ! ! FUNCTION ! Tests VDS API wrappers ! ! RETURNS: ! Success: 0 ! Failure: number of errors !------------------------------------------------------------------------- ! SUBROUTINE test_vds(total_error) IMPLICIT NONE INTEGER, INTENT(INOUT) :: total_error CHARACTER(LEN=6), PARAMETER :: VFILENAME="vds.h5" CHARACTER(LEN=3), PARAMETER :: DATASET="VDS" INTEGER(hsize_t) :: VDSDIM0 INTEGER(hsize_t), PARAMETER :: VDSDIM1 = 10 INTEGER(hsize_t), PARAMETER :: VDSDIM2 = 15 INTEGER(hsize_t) :: DIM0 INTEGER, PARAMETER :: DIM0_1= 4 ! Initial size of the source datasets INTEGER, PARAMETER :: DIM1 = 10 INTEGER, PARAMETER :: DIM2 = 15 INTEGER, PARAMETER :: RANK = 3 INTEGER(hsize_t), PARAMETER :: PLANE_STRIDE = 4 CHARACTER(LEN=4), DIMENSION(1:PLANE_STRIDE) :: SRC_FILE = (/"a.h5","b.h5","c.h5","d.h5"/) CHARACTER(LEN=3), DIMENSION(1:PLANE_STRIDE) :: SRC_DATASET = (/"AAA","BBB","CCC","DDD"/) INTEGER(hid_t) :: vfile, file, src_space, mem_space, vspace, vdset, dset !Handles INTEGER(hid_t) :: dcpl, dapl INTEGER :: error INTEGER(hsize_t), DIMENSION(1:3) :: & vdsdims = (/4_hsize_t*INT(DIM0_1,hsize_t), VDSDIM1, VDSDIM2/), & vdsdims_max, & dims = (/DIM0_1, DIM1, DIM2/), & memdims = (/DIM0_1, DIM1, DIM2/), & extdims = (/0, DIM1, DIM2/), & ! Dimensions of the extended source datasets chunk_dims = (/DIM0_1, DIM1, DIM2/), & dims_max, & vdsdims_out, vdsdims_max_out, & start, & ! Hyperslab parameters stride, & count, & src_count, block INTEGER(hsize_t), DIMENSION(1:2,1:3) :: vdsdims_out_correct INTEGER(hsize_t), DIMENSION(1:3) :: start_out, & !Hyperslab PARAMETER out stride_out, count_out, block_out INTEGER(hsize_t), DIMENSION(1:3,1:PLANE_STRIDE) :: start_correct INTEGER :: i, j INTEGER(size_t) :: i_sz INTEGER :: layout ! Storage layout INTEGER(size_t) :: num_map ! Number of mappings INTEGER(size_t) :: len ! Length of the string also a RETURN value ! Different sized character buffers CHARACTER(len=LEN(SRC_FILE(1))-3) :: SRC_FILE_LEN_TINY CHARACTER(len=LEN(SRC_FILE(1))-1) :: SRC_FILE_LEN_SMALL CHARACTER(len=LEN(SRC_FILE(1))) :: SRC_FILE_LEN_EXACT CHARACTER(len=LEN(SRC_FILE(1))+1) :: SRC_FILE_LEN_LARGE CHARACTER(len=LEN(SRC_FILE(1))+10) :: SRC_FILE_LEN_HUGE CHARACTER(len=LEN(SRC_DATASET(1))) :: SRC_DATASET_LEN_EXACT INTEGER(HID_T) :: space_out INTEGER :: s_type, virtual_view INTEGER :: type1, type2 INTEGER, DIMENSION(DIM0_1*DIM1*DIM2), TARGET :: wdata TYPE(C_PTR) :: f_ptr INTEGER(SIZE_T) :: nsize LOGICAL :: IsRegular INTEGER(HSIZE_T) :: gap_size ! For testing against vdsdims_out_correct(1,1) = DIM0_1*5 vdsdims_out_correct(2,1) = DIM0_1*8 vdsdims_out_correct(1:2,2) = VDSDIM1 vdsdims_out_correct(1:2,3) = VDSDIM2 VDSDIM0 = H5S_UNLIMITED_F DIM0 = H5S_UNLIMITED_F vdsdims_max = (/VDSDIM0, VDSDIM1, VDSDIM2/) dims_max = (/INT(DIM0,hsize_t), INT(DIM1,hsize_t), INT(DIM2,hsize_t)/) ! ! Create source files and datasets. ! DO i = 1, PLANE_STRIDE ! ! Initialize data for i-th source dataset. DO j = 1, DIM0_1*DIM1*DIM2 wdata(j) = i ENDDO ! ! Create the source files and datasets. Write data to each dataset and ! close all resources. CALL h5fcreate_f(SRC_FILE(i), H5F_ACC_TRUNC_F, file, error) CALL check("h5fcreate_f", error, total_error) CALL h5screate_simple_f(RANK, dims, src_space, error, dims_max) CALL check("h5screate_simple_f", error, total_error) CALL h5pcreate_f(H5P_DATASET_CREATE_F, dcpl, error) CALL check("h5pcreate_f", error, total_error) CALL h5pset_chunk_f(dcpl, RANK, chunk_dims, error) CALL check("h5pset_chunk_f",error, total_error) CALL h5dcreate_f(file, SRC_DATASET(i), H5T_NATIVE_INTEGER, src_space, dset, error, dcpl, H5P_DEFAULT_F, H5P_DEFAULT_F) CALL check("h5dcreate_f",error, total_error) f_ptr = C_LOC(wdata(1)) CALL H5Dwrite_f(dset, H5T_NATIVE_INTEGER, f_ptr, error) CALL check("H5Dwrite_f",error, total_error) CALL H5Sclose_f(src_space, error) CALL check("H5Sclose_f",error, total_error) CALL H5Pclose_f(dcpl, error) CALL check("H5Pclose_f",error, total_error) CALL H5Dclose_f(dset, error) CALL check("H5Dclose_f",error, total_error) CALL H5Fclose_f(file, error) CALL check("H5Fclose_f",error, total_error) ENDDO CALL h5fcreate_f(VFILENAME, H5F_ACC_TRUNC_F, vfile, error) CALL check("h5fcreate_f", error, total_error) ! Create VDS dataspace. CALL H5Screate_simple_f(RANK, vdsdims, vspace, error, vdsdims_max) CALL check("H5Screate_simple_f", error, total_error) ! Create dataspaces for the source dataset. CALL H5Screate_simple_f(RANK, dims, src_space, error, dims_max) CALL check("H5Screate_simple_f", error, total_error) ! Create VDS creation property CALL H5Pcreate_f (H5P_DATASET_CREATE_F, dcpl, error) CALL check("H5Pcreate_f", error, total_error) ! Initialize hyperslab values start(1:3) = 0 stride(1:3) = (/PLANE_STRIDE,1_hsize_t,1_hsize_t/) ! we will select every fifth plane in VDS count(1:3) = (/H5S_UNLIMITED_F,1_hsize_t,1_hsize_t/) src_count(1:3) = (/H5S_UNLIMITED_F,1_hsize_t,1_hsize_t/) block(1:3) = (/1, DIM1, DIM2/) ! ! Build the mappings ! start_correct = 0 CALL H5Sselect_hyperslab_f(src_space, H5S_SELECT_SET_F, start, src_count, error, block=block) CALL check("H5Sselect_hyperslab_f", error, total_error) DO i = 1, PLANE_STRIDE start_correct(1,i) = start(1) CALL H5Sselect_hyperslab_f(vspace, H5S_SELECT_SET_F, start, count, error, stride=stride, block=block) CALL check("H5Sselect_hyperslab_f", error, total_error) IF(i.eq.1)THEN ! check src_file and src_dataset with trailing blanks CALL H5Pset_virtual_f (dcpl, vspace, SRC_FILE(i)//" ", SRC_DATASET(i)//" ", src_space, error) ELSE CALL H5Pset_virtual_f (dcpl, vspace, SRC_FILE(i), SRC_DATASET(i), src_space, error) ENDIF CALL check("H5Pset_virtual_f", error, total_error) start(1) = start(1) + 1 ENDDO CALL H5Sselect_none_f(vspace, error) CALL check("H5Sselect_none_f", error, total_error) ! Create a virtual dataset CALL H5Dcreate_f(vfile, DATASET, H5T_NATIVE_INTEGER, vspace, vdset, error, dcpl, H5P_DEFAULT_F, H5P_DEFAULT_F) CALL check("H5Dcreate_f", error, total_error) CALL H5Sclose_f(vspace, error) CALL check("H5Sclose_f", error, total_error) CALL H5Sclose_f(src_space, error) CALL check("H5Sclose_f", error, total_error) CALL H5Pclose_f(dcpl, error) CALL check("H5Pclose_f", error, total_error) ! Let's add data to the source datasets and check new dimensions for VDS ! We will add only one plane to the first source dataset, two planes to the ! second one, three to the third, and four to the forth. DO i = 1, PLANE_STRIDE ! ! Initialize data for i-th source dataset. DO j = 1, i*DIM1*DIM2 wdata(j) = 10*i ENDDO ! ! Open the source files and datasets. Append data to each dataset and ! close all resources. CALL H5Fopen_f (SRC_FILE(i), H5F_ACC_RDWR_F, file, error) CALL check("H5Fopen_f", error, total_error) CALL H5Dopen_f (file, SRC_DATASET(i), dset, error) CALL check("H5Dopen_f", error, total_error) extdims(1) = DIM0_1+i CALL H5Dset_extent_f(dset, extdims, error) CALL check("H5Dset_extent_f", error, total_error) CALL H5Dget_space_f(dset, src_space, error) CALL check("H5Dget_space_f", error, total_error) start(1:3) = (/DIM0_1,0,0/) count(1:3) = 1 block(1:3) = (/i, DIM1, DIM2/) memdims(1) = i CALL H5Screate_simple_f(RANK, memdims, mem_space, error) CALL check("H5Screate_simple_f", error, total_error) CALL H5Sselect_hyperslab_f(src_space, H5S_SELECT_SET_F, start,count, error,block=block) CALL check("H5Sselect_hyperslab_f", error, total_error) f_ptr = C_LOC(wdata(1)) CALL H5Dwrite_f(dset, H5T_NATIVE_INTEGER, f_ptr, error, mem_space, src_space, H5P_DEFAULT_F) CALL check("H5Dwrite_f", error, total_error) CALL H5Sclose_f(src_space, error) CALL check("H5Sclose_f", error, total_error) call H5Dclose_f(dset, error) CALL check("H5Dclose_f", error, total_error) call H5Fclose_f(file, error) CALL check("H5Fclose_f", error, total_error) ENDDO call H5Dclose_f(vdset, error) CALL check("H5Dclose_f", error, total_error) call H5Fclose_f(vfile, error) CALL check("H5Fclose_f", error, total_error) ! ! begin the read section ! ! Open file and dataset using the default properties. CALL H5Fopen_f(VFILENAME, H5F_ACC_RDONLY_F, vfile, error) CALL check("H5Fopen_f", error, total_error) ! ! Open VDS using different access properties to use max or ! min extents depending on the sizes of the underlying datasets CALL H5Pcreate_f(H5P_DATASET_ACCESS_F, dapl, error) CALL check("H5Pcreate_f", error, total_error) DO i = 1, 2 IF(i.NE.1)THEN CALL H5Pset_virtual_view_f(dapl, H5D_VDS_LAST_AVAILABLE_F, error) CALL check("H5Pset_virtual_view_f", error, total_error) ELSE CALL H5Pset_virtual_view_f(dapl, H5D_VDS_FIRST_MISSING_F, error) CALL check("H5Pset_virtual_view_f", error, total_error) ENDIF CALL H5Dopen_f(vfile, DATASET, vdset, error, dapl) CALL check("H5Dopen_f", error, total_error) ! Let's get space of the VDS and its dimension we should get 32(or 20)x10x10 CALL H5Dget_space_f(vdset, vspace, error) CALL check("H5Dget_space_f", error, total_error) CALL H5Sget_simple_extent_dims_f(vspace, vdsdims_out, vdsdims_max_out, error) CALL check("H5Sget_simple_extent_dims_f", error, total_error) ! check VDS dimensions DO j = 1, RANK IF(vdsdims_out(j).NE.vdsdims_out_correct(i,j))THEN total_error = total_error + 1 EXIT ENDIF ENDDO CALL H5Pget_virtual_view_f(dapl, virtual_view, error) CALL check("h5pget_virtual_view_f", error, total_error) IF(i.EQ.1)THEN IF(virtual_view .NE. H5D_VDS_FIRST_MISSING_F)THEN total_error = total_error + 1 ENDIF ELSE IF(virtual_view .NE. H5D_VDS_LAST_AVAILABLE_F)THEN total_error = total_error + 1 ENDIF ENDIF ! Close CALL H5Dclose_f(vdset, error) CALL check("H5Dclose_f", error, total_error) CALL H5Sclose_f(vspace, error) CALL check("H5Sclose_f", error, total_error) ENDDO CALL H5Dopen_f(vfile, DATASET, vdset, error) CALL check("H5Dopen_f", error, total_error) ! ! Get creation property list and mapping properties. ! CALL H5Dget_create_plist_f (vdset, dcpl, error) CALL check("H5Dget_create_plist_f", error, total_error) ! ! Get storage layout. CALL H5Pget_layout_f(dcpl, layout, error) CALL check("H5Pget_layout_f", error, total_error) IF (H5D_VIRTUAL_F .NE. layout) THEN PRINT*,"Wrong layout found" total_error = total_error + 1 ENDIF ! ! Find number of mappings. CALL H5Pget_virtual_count_f(dcpl, num_map, error) CALL check("H5Pget_virtual_count_f", error, total_error) IF(num_map.NE.4_size_t)THEN PRINT*,"Number of mappings is incorrect" total_error = total_error + 1 ENDIF ! ! Get mapping parameters for each mapping. ! DO i_sz = 1, num_map CALL H5Pget_virtual_vspace_f(dcpl, INT(i_sz-1,size_t), vspace, error) CALL check("H5Pget_virtual_vspace_f", error, total_error) CALL h5sget_select_type_f(vspace, s_type, error) CALL check("h5sget_select_type_f", error, total_error) IF(s_type.EQ.H5S_SEL_HYPERSLABS_F)THEN CALL H5Sis_regular_hyperslab_f(vspace, IsRegular, error) CALL check("H5Sis_regular_hyperslab_f", error, total_error) IF(IsRegular)THEN CALL H5Sget_regular_hyperslab_f(vspace, start_out, stride_out, count_out, block_out, error) CALL check("H5Sget_regular_hyperslab_f", error, total_error) DO j = 1, 3 IF(start_out(j).NE.start_correct(j,i_sz) .OR. & stride_out(j).NE.stride(j).OR. & count_out(j).NE.src_count(j))THEN total_error = total_error + 1 EXIT ENDIF ENDDO ENDIF END IF ! Get source file name CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_EXACT, error, nsize) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(nsize.NE.LEN(SRC_FILE_LEN_EXACT))THEN PRINT*,"virtual filename size is incorrect" total_error = total_error + 1 ENDIF ! check passing a buffer that is very small CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_TINY, error) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(SRC_FILE_LEN_TINY.NE.SRC_FILE(i)(1:LEN(SRC_FILE_LEN_TINY)))THEN PRINT*,"virtual filename returned is incorrect" total_error = total_error + 1 ENDIF ! check passing a buffer that small by one CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_SMALL, error) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(SRC_FILE_LEN_SMALL.NE.SRC_FILE(i)(1:LEN(SRC_FILE_LEN_SMALL)))THEN PRINT*,"virtual filename returned is incorrect" total_error = total_error + 1 ENDIF ! check passing a buffer that is exact CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_EXACT, error) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(SRC_FILE_LEN_EXACT.NE.SRC_FILE(i)(1:LEN(SRC_FILE_LEN_EXACT)))THEN PRINT*,"virtual filename returned is incorrect" total_error = total_error + 1 ENDIF ! check passing a buffer that bigger by one CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_LARGE, error) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(SRC_FILE_LEN_LARGE(1:LEN(SRC_FILE_LEN_EXACT)).NE.SRC_FILE(i)(1:LEN(SRC_FILE_LEN_EXACT)).AND. & SRC_FILE_LEN_LARGE(LEN(SRC_FILE_LEN_EXACT):).NE.'')THEN PRINT*,"virtual filename returned is incorrect" total_error = total_error + 1 ENDIF ! check passing a buffer that is very big CALL H5Pget_virtual_filename_f(dcpl, INT(i-1, size_t), SRC_FILE_LEN_HUGE, error) CALL check("H5Pget_virtual_filename_f", error, total_error) IF(SRC_FILE_LEN_HUGE(1:LEN(SRC_FILE_LEN_EXACT)).NE.SRC_FILE(i)(1:LEN(SRC_FILE_LEN_EXACT)).AND. & SRC_FILE_LEN_HUGE(LEN(SRC_FILE_LEN_EXACT):).NE.'')THEN PRINT*,"virtual filename returned is incorrect" total_error = total_error + 1 ENDIF ! Get source dataset name CALL H5Pget_virtual_dsetname_f(dcpl, INT(i-1, size_t), SRC_DATASET_LEN_EXACT, error, nsize) CALL check("H5Pget_virtual_dsetname_f", error, total_error) CALL H5Pget_virtual_dsetname_f(dcpl, INT(i-1, size_t), SRC_DATASET_LEN_EXACT, error) CALL check("H5Pget_virtual_dsetname_f", error, total_error) IF(SRC_DATASET_LEN_EXACT(1:LEN(SRC_DATASET_LEN_EXACT)).NE.SRC_DATASET(i)(1:LEN(SRC_DATASET_LEN_EXACT)).AND. & SRC_DATASET_LEN_EXACT(LEN(SRC_DATASET_LEN_EXACT):).NE.'')THEN PRINT*,"virtual dataset returned is incorrect" total_error = total_error + 1 ENDIF CALL h5pget_virtual_srcspace_f(dcpl, i_sz - 1_size_t, space_out, error) CALL check("H5Pget_virtual_srcspace_f", error, total_error) CALL h5sget_select_type_f(space_out, type1, error) CALL check("H5Sget_select_type_f", error, total_error) CALL h5sget_select_type_f(vspace, type2, error) CALL check("H5Sget_select_type_f", error, total_error) IF(type1.NE.type2)THEN total_error = total_error + 1 ENDIF ENDDO ! ! Close and release resources. ! Clear virtual layout in DCPL CALL h5pset_layout_f(dcpl, H5D_VIRTUAL_F,error) CALL check("H5Pset_layout_f", error, total_error) CALL H5Pclose_f(dcpl, error) CALL check("H5Pclose_f", error, total_error) CALL H5Dclose_f(vdset, error) CALL check("H5Dclose_f", error, total_error) ! Reopen VDS with printf gap set to 1 CALL H5Pset_virtual_printf_gap_f(dapl, 1_hsize_t, error) CALL check("H5Pset_virtual_printf_gap_f", error, total_error) CALL H5Dopen_f(vfile, DATASET, vdset, error, dapl) CALL check("H5Dopen_f", error, total_error) CALL H5Pget_virtual_printf_gap_f(dapl, gap_size, error) CALL check("H5Pget_virtual_printf_gap_f", error, total_error) IF(gap_size.NE.1_hsize_t)THEN PRINT*,"gapsize is incorrect" total_error = total_error + 1 ENDIF CALL H5Dclose_f(vdset, error) CALL check("H5Dclose_f", error, total_error) CALL H5Sclose_f(vspace, error) CALL check("H5Sclose_f", error, total_error) CALL H5Pclose_f(dapl, error) CALL check("H5Pclose_f", error, total_error) CALL H5Fclose_f(vfile, error) CALL check("H5Fclose_f", error, total_error) END SUBROUTINE test_vds END MODULE TH5P_F03