diff options
Diffstat (limited to 'tools/h5diff/h5diff.c')
| -rw-r--r-- | tools/h5diff/h5diff.c | 986 |
1 files changed, 812 insertions, 174 deletions
diff --git a/tools/h5diff/h5diff.c b/tools/h5diff/h5diff.c index 15f7ac6..4980eeb 100644 --- a/tools/h5diff/h5diff.c +++ b/tools/h5diff/h5diff.c @@ -2,25 +2,67 @@ #include <stdio.h> #include <stdlib.h> +#include <math.h> + +#include <assert.h> + #include "hdf5.h" +#include "h5trav.h" + + +/* + +Example 1 +dset1 dset2 h5diff_test1.h5 h5diff_test2.h5 +Example 2 +dset1 dset2 -v h5diff_test1.h5 h5diff_test2.h5 +Example 3 +-v h5diff_test1.h5 h5diff_test2.h5 +Example 4 +-v -r h5diff_test1.h5 h5diff_test2.h5 +Example 5 +dset1 dset2 -n 2 h5diff_test1.h5 h5diff_test2.h5 +Example 6 +dset3 dset4 -m 0.01 h5diff_test1.h5 h5diff_test2.h5 +Example 7 +dset5 dset6 -p 0.05 h5diff_test1.h5 h5diff_test2.h5 + +*/ #ifndef FALSE -#define FALSE 0 +#define FALSE 0 #endif #ifndef TRUE -#define TRUE 1 +#define TRUE 1 #endif -int do_test_files(); -int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ); -void array_diff( void *buf1, void *buf2, hsize_t tot_cnt, hid_t type_id ); -herr_t get_ndsets( hid_t loc_id, const char *group_name ); -herr_t count_dsets( hid_t loc_id, const char *name, void *op_data); +typedef struct options_t +{ + int d_; + int g_; + int t_; + int r_; /* report only what objects differ */ + int m_; + double m_limit; + int p_; + double p_percent; + int n_; + int n_number_count; + int a_; + int v_; /* verbose */ +} options_t; +int do_test_files(); + +int diff_dataset( hid_t file1_id, hid_t file2_id, const char *obj1_name, + const char *obj2_name, options_t options ); +int array_diff( void *buf1, void *buf2, hsize_t tot_cnt, hid_t type_id, int rank, + hsize_t *dims, options_t options ); +void print_pos( hsize_t curr_pos, hsize_t *acc, hsize_t *pos, int rank ); @@ -46,12 +88,19 @@ herr_t count_dsets( hid_t loc_id, const char *name, void *op_data); void usage( const char *progname ) { #define USAGE "\ - [-d dset] The name of the dataset to compare\n\ + [-m ] Print a basic help message\n\ + [-v ] Print a verbose output message\n\ + [-r ] Print only what objects differ \n\ + [-d ] Compare only datasets\n\ + [-g ] Compare only groups\n\ + [-n count] Print difference up to count number for each variable\n\ + [-m delta] Print difference when it is greater than limit delta\n\ + [-p percent] Print differences which are within a certain percent value\n\ file1 File name of the first HDF5 file\n\ file2 File name of the second HDF5 file\n" fprintf(stderr, - "%s [-d] file1 file2\n%s", + "%s [-d delta] file1 file2\n%s", progname, USAGE); exit(EXIT_FAILURE); @@ -80,15 +129,30 @@ int main(int argc, const char *argv[]) { const char *progname = argv[0]; - int dset_only = FALSE; - int argno; - const char *s = NULL; - const char *file1_name = NULL; - const char *file2_name = NULL; - const char *dset_name = NULL; - hid_t file1_id, file2_id; - herr_t status; - int ndsets; + int dset_only = FALSE; + int have_group = FALSE; + int argno, i, j; + const char *s = NULL; + const char *file1_name = NULL; + const char *file2_name = NULL; + char *obj1_name = NULL; + char *obj2_name = NULL; + hid_t file1_id, file2_id; + herr_t status; + int nobjects1, nobjects2; + info_t *info1=NULL; + info_t *info2=NULL; + char *attr1_name = NULL; + char *attr2_name = NULL; + options_t options = {0,0,0,0,0,0,0,0,0,0,0,0}; + + + /* string compare */ + char *pdest; + int result; + int len; + + do_test_files(); @@ -98,15 +162,15 @@ int main(int argc, const char *argv[]) *------------------------------------------------------------------------- */ - if (argc < 4) { + if (argc < 2) { usage( progname ); exit(EXIT_FAILURE); } - /* parse command line options */ - for (argno=1; argno<argc ; argno++) + /* last 2 items are the file names */ + for (argno=1; argno<argc-2 ; argno++) { - + /* get the single-letter switches */ if ( '-'==argv[argno][0] ) { @@ -117,54 +181,354 @@ int main(int argc, const char *argv[]) case 'h': usage(progname); exit(EXIT_SUCCESS); - case 'd': - dset_only = TRUE; - dset_name = argv[argno+1]; + case 'v': + options.v_ = 1; + break; + case 'd': + options.d_ = 1; + break; + case 'g': + options.g_ = 1; + break; + case 't': + options.t_ = 1; + break; + case 'r': + options.r_ = 1; + break; + case 'a': + options.a_ = 1; + break; + case 'm': + /* if it is not another option */ + if ( '-' !=argv[argno+1][0] ) + { + options.m_ = 1; + options.m_limit = atof(argv[argno+1]); + } + break; + case 'p': + if ( '-' !=argv[argno+1][0] ) + { + options.p_ = 1; + options.p_percent = atof(argv[argno+1]); + } + break; + case 'n': + if ( '-' !=argv[argno+1][0] ) + { + options.n_ = 1; + options.n_number_count = atoi(argv[argno+1]); + } break; } /*switch*/ } /*for*/ } /*if*/ - + + else + + { + + /* 2 last args are the file names, and it is not a -switch parameter */ + if ( argno < argc-2 && '-' !=argv[argno-1][0] ) + { + if ( obj1_name == NULL ) + obj1_name = argv[argno]; + + if ( obj2_name == NULL ) + { + + /* check if we have a second object name */ + if ( '-' !=argv[argno+1][0] ) + /* yes */ + obj2_name = argv[argno+1]; + else + /* no */ + obj2_name = obj1_name; + } + } + + } + }/*for*/ - + /*------------------------------------------------------------------------- * process the files *------------------------------------------------------------------------- */ - file1_name = argv[argno-2]; - file2_name = argv[argno-1]; + file1_name = argv[argc-2]; + file2_name = argv[argc-1]; /* Open the files */ if ((file1_id=H5Fopen(file1_name,H5F_ACC_RDONLY,H5P_DEFAULT))<0 || (file2_id=H5Fopen(file2_name,H5F_ACC_RDONLY,H5P_DEFAULT))<0) exit(EXIT_FAILURE); - /* Get the number of datasets */ - ndsets = get_ndsets( file1_id, "." ); - if ( dset_only ) +/*------------------------------------------------------------------------- + * get the number of objects in the files + *------------------------------------------------------------------------- + */ + + nobjects1 = H5get_object_info( file1_id, NULL ); + nobjects2 = H5get_object_info( file2_id, NULL ); + +/*------------------------------------------------------------------------- + * get the list of objects in the files + *------------------------------------------------------------------------- + */ + + info1 = (info_t*) malloc( nobjects1 * sizeof(info_t)); + info2 = (info_t*) malloc( nobjects2 * sizeof(info_t)); + + H5get_object_info( file1_id, info1 ); + H5get_object_info( file2_id, info2 ); + + if ( options.v_ ) { + printf("File 1: # of entries = %d\n", nobjects1); + for ( i = 0; i < nobjects1; i++) + { + switch ( info1[i].type ) + { + case H5G_GROUP: + printf("%s \t %s\n", info1[i].name, "group" ); + break; + case H5G_DATASET: + printf("%s \t %s\n", info1[i].name, "dataset" ); + break; + case H5G_TYPE: + printf("%s \t %s\n", info1[i].name, "datatype" ); + break; + } + } + + printf("File 2: # of entries = %d\n", nobjects2); + for ( i = 0; i < nobjects2; i++) + { + switch ( info2[i].type ) + { + case H5G_GROUP: + printf("%s \t %s\n", info2[i].name, "group" ); + break; + case H5G_DATASET: + printf("%s \t %s\n", info2[i].name, "dataset" ); + break; + case H5G_TYPE: + printf("%s \t %s\n", info2[i].name, "datatype" ); + break; + } + } + printf("\n"); + } + + - h5diff_dataset(file1_id,file2_id,dset_name); + +/*------------------------------------------------------------------------- + * object name was supplied + *------------------------------------------------------------------------- + */ + + /* object name was supplied, find obj1_name */ + if ( obj1_name ) + { + + for ( i = 0; i < nobjects1; i++) + { + + pdest = strstr( info1[i].name, obj1_name ); + result = pdest - info1[i].name; + len = strlen(obj1_name); + + /* found at position result */ + if( pdest != NULL && + /* check if it is not a substring */ + info1[i].name[result-1] == '/' && + /* check if it is the last or in the middle */ + (info1[i].name[result+len]=='/' || + info1[i].name[result+len]=='\0') ) + { + if ( options.v_ ) + printf( "%s found in file %s\n\n", info1[i].name, file1_name); + + /* go to second file and find obj2_name */ + for ( j = 0; j < nobjects2; j++) + { + + pdest = strstr( info2[j].name, obj2_name ); + result = pdest - info2[j].name; + + len = strlen(obj2_name); + + /* found at position result */ + if( pdest != NULL && + /* check if it is not a substring */ + info2[j].name[result-1] == '/' && + /* check if it is the last or in the middle */ + (info2[j].name[result+len]=='/' || info2[j].name[result+len]=='\0') && + /* objects are the same type */ + info1[i].type == info2[j].type ) + { + + switch ( info1[i].type ) + { + + /*------------------------------------------------------------------------- + * H5G_GROUP + *------------------------------------------------------------------------- + */ + + case H5G_GROUP: + + if ( options.d_==1 || options.t_==1 ) break; + printf( "%s found in file %s, %s found in file %s\n", info1[i].name, file1_name, + info2[j].name, file2_name ); + + + break; + + /*------------------------------------------------------------------------- + * H5G_DATASET + *------------------------------------------------------------------------- + */ + + case H5G_DATASET: + if ( options.g_==1 || options.t_==1 ) break; + if ( options.v_ ) + printf( "%s found in file %s\n\n", info2[j].name, file2_name); + if ( options.r_==1 ) break; + /* compare with the absolute name */ + diff_dataset(file1_id,file2_id,info1[i].name,info2[j].name,options); + printf("\n"); + break; + + /*------------------------------------------------------------------------- + * H5G_TYPE + *------------------------------------------------------------------------- + */ + + case H5G_TYPE: + + break; + + } /* switch */ + + } + + } /* j */ + + } + + } /* i */ } + + + +/*------------------------------------------------------------------------- + * compare all datasets + *------------------------------------------------------------------------- + */ + else + + { + for ( i = 0; i < nobjects1; i++) + { + obj1_name = info1[i].name; + len = strlen(obj1_name); + + for ( j = 0; j < nobjects2; j++) + { + /* find an object in file2 with same name as in file 1 */ + + pdest = strstr( info2[j].name, obj1_name ); + result = pdest - info2[j].name; + + /* found at position result */ + if( pdest != NULL && + /* check if it is not a substring */ + info2[j].name[result] == '/' && + /* check if it is the last or in the middle */ + (info2[j].name[result+len]=='/' || info2[j].name[result+len]=='\0') && + /* objects are the same type */ + info1[i].type == info2[j].type ) + { + + switch ( info1[i].type ) + { + + /*------------------------------------------------------------------------- + * H5G_GROUP + *------------------------------------------------------------------------- + */ + + case H5G_GROUP: + + if ( options.d_==1 || options.t_==1 ) break; + printf( "%s found in file %s, %s found in file %s\n", info1[i].name, file1_name, + info2[j].name, file2_name ); + + break; + + /*------------------------------------------------------------------------- + * H5G_DATASET + *------------------------------------------------------------------------- + */ + + case H5G_DATASET: + + if ( options.g_==1 || options.t_==1 ) break; + if ( options.v_ ) + printf( "%s found in file %s, %s found in file %s\n", info1[i].name, file1_name, + info2[j].name, file2_name ); + if ( options.r_==1 ) break; + /* compare with the absolute name */ + diff_dataset(file1_id,file2_id,info1[i].name,info2[j].name,options); + printf("\n"); + break; + + /*------------------------------------------------------------------------- + * H5G_TYPE + *------------------------------------------------------------------------- + */ + + case H5G_TYPE: + + break; + + } + + } + + } /* j */ + } /* i */ + + } - /* Close */ - status = H5Fclose(file1_id); - status = H5Fclose(file2_id); + /* close */ + status = H5Fclose(file1_id); + status = H5Fclose(file2_id); + if ( info1 ) + free(info1); + if ( info2 ) + free(info2); + return 0; - return 0; } + + + /*------------------------------------------------------------------------- - * Function: h5diff_dataset + * Function: diff_dataset * * Purpose: * @@ -183,40 +547,41 @@ int main(int argc, const char *argv[]) -int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) +int diff_dataset( hid_t file1_id, hid_t file2_id, const char *obj1_name, + const char *obj2_name, options_t options ) { - hid_t dset1_id, dset2_id; - hid_t space1_id, space2_id; - hid_t type1_id, type2_id; - hid_t rank1, rank2; - void *buf1, *buf2; - hsize_t tot_cnt, tot_cnt1, tot_cnt2; - hsize_t dims1[32], dims2[32]; - int i; - herr_t status; + hid_t dset1_id, dset2_id; + hid_t space1_id, space2_id; + hid_t type1_id, type2_id; + hid_t rank1, rank2; + void *buf1, *buf2; + hsize_t tot_cnt, tot_cnt1, tot_cnt2; + hsize_t dims1[32], dims2[32]; + int i; + herr_t status; /*------------------------------------------------------------------------- * open the handles *------------------------------------------------------------------------- */ - /* Open the datasets */ - if ( (dset1_id = H5Dopen(file1_id,dset_name)) < 0 ) + /* Open the datasets */ + if ( (dset1_id = H5Dopen(file1_id,obj1_name)) < 0 ) return -1; - if ( (dset2_id = H5Dopen(file2_id,dset_name)) < 0 ) + if ( (dset2_id = H5Dopen(file2_id,obj2_name)) < 0 ) return -1; - /* Get the datatype */ + /* Get the datatype */ if ( (type1_id = H5Dget_type(dset1_id)) < 0 ) goto out; - /* Get the datatype */ + /* Get the datatype */ if ( (type2_id = H5Dget_type(dset2_id)) < 0 ) goto out; - /* Get the dataspace handle */ + /* Get the dataspace handle */ if ( (space1_id = H5Dget_space(dset1_id)) < 0 ) return -1; @@ -224,7 +589,7 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) if ( (rank1 = H5Sget_simple_extent_ndims(space1_id)) < 0 ) return -1; - /* Get the dataspace handle */ + /* Get the dataspace handle */ if ( (space2_id = H5Dget_space(dset2_id)) < 0 ) return -1; @@ -232,11 +597,11 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) if ( (rank2 = H5Sget_simple_extent_ndims(space2_id)) < 0 ) return -1; - /* Get dimensions */ + /* Get dimensions */ if ( H5Sget_simple_extent_dims(space1_id,dims1,NULL) < 0 ) goto out; - /* Get dimensions */ + /* Get dimensions */ if ( H5Sget_simple_extent_dims(space2_id,dims2,NULL) < 0 ) goto out; @@ -245,13 +610,13 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) * compare *------------------------------------------------------------------------- */ - - printf("Dataset Name: %s .... Comparing\n", dset_name); - - if ( rank1 != rank2 ) - { + + if ( rank1 != rank2 ) + { goto out; - } + } + + printf("Dataset Names: %s and %s .... Comparing\n", obj1_name, obj2_name); tot_cnt1 = 1; for (i = 0; i < rank1; i++) @@ -265,14 +630,14 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) tot_cnt2 *= dims2[i]; } - buf1 = (void *) malloc((unsigned) (tot_cnt1*H5Tget_size(type1_id))); + buf1 = (void *) malloc((unsigned) (tot_cnt1*H5Tget_size(type1_id))); buf2 = (void *) malloc((unsigned) (tot_cnt2*H5Tget_size(type2_id))); if ( H5Dread(dset1_id,type1_id,H5S_ALL,H5S_ALL,H5P_DEFAULT,buf1) < 0 ) goto out; - if ( H5Dread(dset2_id,type2_id,H5S_ALL,H5S_ALL,H5P_DEFAULT,buf2) < 0 ) + if ( H5Dread(dset2_id,type2_id,H5S_ALL,H5S_ALL,H5P_DEFAULT,buf2) < 0 ) goto out; if (tot_cnt1 > tot_cnt2) @@ -280,11 +645,25 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) else tot_cnt = tot_cnt1; - array_diff(buf1,buf2,tot_cnt,type1_id); + if ( array_diff(buf1,buf2,tot_cnt,type1_id,rank1,dims1,options) == 0 ) + printf("No differences found\n" ); - free((char *) buf1); + free((char *) buf1); free((char *) buf2); + +/*------------------------------------------------------------------------- + * attributes + *------------------------------------------------------------------------- + */ + + if ( options.a_==1 ) + { + + + + } + /*------------------------------------------------------------------------- * close @@ -294,15 +673,15 @@ int h5diff_dataset( hid_t file1_id, hid_t file2_id, const char *dset_name ) out: - /* Close */ - status = H5Dclose(dset1_id); - status = H5Dclose(dset2_id); - status = H5Sclose(space1_id); - status = H5Sclose(space2_id); - status = H5Tclose(type1_id); - status = H5Tclose(type2_id); + /* Close */ + status = H5Dclose(dset1_id); + status = H5Dclose(dset2_id); + status = H5Sclose(space1_id); + status = H5Sclose(space2_id); + status = H5Tclose(type1_id); + status = H5Tclose(type2_id); - return 0; + return 0; } @@ -327,23 +706,31 @@ out: */ -void array_diff( void *buf1, void *buf2, hsize_t tot_cnt, hid_t type_id ) +int array_diff( void *buf1, void *buf2, hsize_t tot_cnt, hid_t type_id, int rank, + hsize_t *dims, options_t options ) { - -#if 0 char *i1ptr1, *i1ptr2; short *i2ptr1, *i2ptr2; + int *i4ptr1, *i4ptr2; float *fptr1, *fptr2; double *dptr1, *dptr2; -#endif + int found = 0; - int *i4ptr1, *i4ptr2; - int i; + /* accumulator and matrix position */ + hsize_t acc[32]; + hsize_t pos[32]; + int i; + + H5T_class_t type_class; + size_t type_size; - H5T_class_t type_class; - size_t type_size; + acc[rank-1]=1; + for(i=(rank-2); i>=0; i--) + { + acc[i]=acc[i+1]*dims[i+1]; + } - /* Get the class. */ + /* Get the class. */ type_class = H5Tget_class( type_id ); /* Get the size. */ @@ -355,172 +742,423 @@ void array_diff( void *buf1, void *buf2, hsize_t tot_cnt, hid_t type_id ) case H5T_INTEGER: - switch(type_size) - { - case 4: + switch(type_size) + { - i4ptr1 = (int *) buf1; - i4ptr2 = (int *) buf2; - for ( i = 0; i < tot_cnt; i++) - { - if (*i4ptr1 != *i4ptr2) - { - printf("Index: %d, File1: %d, File2: %d\n", i, *i4ptr1, *i4ptr2); - } - i4ptr1++; i4ptr2++; - } + /*------------------------------------------------------------------------- + * H5T_INTEGER 1 + *------------------------------------------------------------------------- + */ + + case 1: + i1ptr1 = (char *) buf1; + i1ptr2 = (char *) buf2; + + break; + + /*------------------------------------------------------------------------- + * H5T_INTEGER 2 + *------------------------------------------------------------------------- + */ + + case 2: + i2ptr1 = (short *) buf1; + i2ptr2 = (short *) buf2; + + break; + + /*------------------------------------------------------------------------- + * H5T_INTEGER 4 + *------------------------------------------------------------------------- + */ + + case 4: + i4ptr1 = (int *) buf1; + i4ptr2 = (int *) buf2; + for ( i = 0; i < tot_cnt; i++) + { + + if (*i4ptr1 != *i4ptr2) + { + print_pos( i, acc, pos, rank ); + printf(" %d, %d\n", *i4ptr1, *i4ptr2); + found++; + + if ( options.n_ && found>options.n_number_count-1) + return found; + } + i4ptr1++; i4ptr2++; + } + break; + + + } /*switch*/ + + break; /* H5T_INTEGER */ - break; + case H5T_FLOAT: - } /*switch*/ + switch(type_size) + { - - break; /* H5T_INTEGER */ + /*------------------------------------------------------------------------- + * H5T_FLOAT 4 + *------------------------------------------------------------------------- + */ + case 4: + fptr1 = (float *) buf1; + fptr2 = (float *) buf2; + + for ( i = 0; i < tot_cnt; i++) + { + if ( options.n_ && i>options.n_number_count-1) + return found; + + /* delta but not percentage */ + if ( options.m_ && !options.p_ ) + { + if ( fabs(*fptr1 - *fptr2) > options.m_limit ) + { + print_pos( i, acc, pos, rank ); + printf(" %f, %f\n", *fptr1, *fptr2); + found=1; + } + } + + /* percentage but not delta */ + else if ( !options.m_ && options.p_ ) + { + if ( 1 - *fptr1 / *fptr2 > options.p_percent ) + { + print_pos( i, acc, pos, rank ); + printf(" %f, %f\n", *fptr1, *fptr2); + found=1; + } + } + + /* percentage and delta */ + else if ( options.m_ && options.p_ ) + { + if ( 1 - *fptr1 / *fptr2 > options.p_percent && + fabs(*fptr1 - *fptr2) > options.m_limit ) + { + print_pos( i, acc, pos, rank ); + printf(" %f, %f\n", *fptr1, *fptr2); + found=1; + } + } + + else + + if (*fptr1 != *fptr2) + { + print_pos( i, acc, pos, rank ); + printf(" %f, %f\n", *fptr1, *fptr2); + found=1; + } + fptr1++; fptr2++; + } + break; + + /*------------------------------------------------------------------------- + * H5T_FLOAT 8 + *------------------------------------------------------------------------- + */ + + case 8: + dptr1 = (double *) buf1; + dptr2 = (double *) buf2; + - - } /*switch*/ + break; + + } /*switch*/ + + + + break; /* H5T_FLOAT */ + + } /*switch*/ - + return found; } - - /*------------------------------------------------------------------------- - * Function: get_ndsets + * Function: print_pos + * + * Purpose: conver an array index position to matrix notation * - * Purpose: Counts the number of datasets in the group GROUP_NAME + * Return: pos matrix array * * Programmer: Pedro Vicente, pvn@ncsa.uiuc.edu * - * Date: October 10, 2002 + * Date: December 19, 2002 * - * Return: - * Success: The return value of the first operator that - * returns non-zero, or zero if all members were - * processed with no operator returning non-zero. + * Comments: * - * Failure: Negative if something goes wrong within the - * library, or the negative value returned by one - * of the operators. + * Modifications: * *------------------------------------------------------------------------- */ -herr_t get_ndsets( hid_t loc_id, const char *group_name ) + + +void print_pos( hsize_t curr_pos, hsize_t *acc, hsize_t *pos, int rank ) { + int i; - int ndsets = 0; + for ( i = 0; i < rank; i++) + pos[i]=0; + + for ( i = 0; i < rank; i++) + { + pos[i] = curr_pos/acc[i]; + curr_pos -= acc[i]*pos[i]; + } + assert( curr_pos == 0 ); + + printf("[ " ); + for ( i = 0; i < rank; i++) + { + printf("%d ", pos[i]+1 ); + } + printf("] :" ); - if ( H5Giterate( loc_id, group_name, NULL, count_dsets, (void *)&ndsets ) < 0 ) - return -1; - return ndsets; } /*------------------------------------------------------------------------- - * Function: count_dsets - * - * Purpose: operator function used by get_ndsets - * - * Programmer: Pedro Vicente, pvn@ncsa.uiuc.edu - * - * Date: October 10, 2002 - * - * Comments: - * - * Modifications: - * + * do some test files *------------------------------------------------------------------------- */ -static herr_t count_dsets( hid_t loc_id, const char *name, void *op_data) +int do_test_files() { - H5G_stat_t statbuf; + hid_t file1_id, file2_id; + hid_t dataset_id; + hid_t space_id; + hid_t group_id, group2_id; + hsize_t dims [1] = { 7 }; + hsize_t dims2 [2] = { 3,2 }; + int data1[7] = {1,1,1,1,1,1,1}; + int data2[7] = {1,1,1,4,5,6,7}; + float data3[7] = {1,1,3,4,5,6,7}; + float data4[7] = {1,1,3.02f,4.002f,5.00002f,6,7}; + float data5[3][2] = {1,1,3,4,5,6}; + float data6[3][2] = {1,1.1f,3.02f,4.002f,5.00002f,6}; + + /* attribute */ + size_t size_attr = 5; + float attr_data1[5] = {1,2,3,4,5}; + float attr_data2[5] = {1,2.1f,3.01f,4.001f,5.00001f}; + herr_t status; - if (H5Gget_objinfo( loc_id, name, FALSE, &statbuf) < 0 ) - return -1; +/*------------------------------------------------------------------------- + * Create two files + *------------------------------------------------------------------------- + */ + + /* Create a file */ + file1_id = H5Fcreate ("h5diff_test1.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); - if ( statbuf.type == H5G_DATASET ) - (*(int *)op_data)++; + /* Create a file */ + file2_id = H5Fcreate ("h5diff_test2.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); - /* Define a default zero value for return. This will cause the iterator to continue */ - return 0; -} +/*------------------------------------------------------------------------- + * Make dataset "dset1" on file1 + *------------------------------------------------------------------------- + */ + /* Create a data space */ + space_id = H5Screate_simple(1,dims,NULL); + /* Create a dataset "dset1" */ + dataset_id = H5Dcreate(file1_id,"dset1",H5T_NATIVE_INT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data1); + /* Close */ + status = H5Dclose(dataset_id); + status = H5Sclose(space_id); /*------------------------------------------------------------------------- - * do some test files + * Make dataset "dset3" on file1 *------------------------------------------------------------------------- */ -int do_test_files() -{ + /* Create a data space */ + space_id = H5Screate_simple(1,dims,NULL); + + /* Create a dataset "dset3" */ + dataset_id = H5Dcreate(file1_id,"dset3",H5T_NATIVE_FLOAT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_FLOAT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data3); + + /* Close */ + status = H5Dclose(dataset_id); + status = H5Sclose(space_id); + +/*------------------------------------------------------------------------- + * Make group "g1" on file1 + *------------------------------------------------------------------------- + */ + + /* Create a group. */ + group_id = H5Gcreate(file1_id, "g1", 0); + + /* Close */ + status = H5Gclose(group_id); + + +/*------------------------------------------------------------------------- + * Make dataset "dset1" on file2 + *------------------------------------------------------------------------- + */ + + /* Create a data space */ + space_id = H5Screate_simple(1,dims,NULL); + + /* Create a dataset "dset1" */ + dataset_id = H5Dcreate(file2_id,"dset1",H5T_NATIVE_INT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data2); + + /* Close */ + status = H5Dclose(dataset_id); + status = H5Sclose(space_id); + + + /*------------------------------------------------------------------------- + * Make dataset "dset2" on file2 + *------------------------------------------------------------------------- + */ + + /* Create a data space */ + space_id = H5Screate_simple(1,dims,NULL); + + /* Create a dataset "dset1" */ + dataset_id = H5Dcreate(file2_id,"dset2",H5T_NATIVE_INT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data2); + + /* Close */ + status = H5Dclose(dataset_id); + status = H5Sclose(space_id); - hid_t file_id; - hid_t dataset_id; - hid_t space_id; - hsize_t dims1[1] = { 5 }; - int data1[5] = {1,1,1,1,1}; - hsize_t dims2[1] = { 5 }; - int data2[5] = {1,1,1,2,2}; - herr_t status; /*------------------------------------------------------------------------- - * Create one file + * Make dataset "g1/dset1" on file2 *------------------------------------------------------------------------- */ + + /* Create a data space */ + space_id = H5Screate_simple(1,dims,NULL); + + /* Create a group. */ + group_id = H5Gcreate(file2_id, "g1", 0); + + /* Create a dataset "g1/dset1" */ + dataset_id = H5Dcreate(group_id,"dset1",H5T_NATIVE_INT,space_id,H5P_DEFAULT); - /* Create a file */ - file_id = H5Fcreate ("h5diff_test1.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data2); + + /* Close */ + status = H5Dclose(dataset_id); + status = H5Gclose(group_id); + status = H5Sclose(space_id); + +/*------------------------------------------------------------------------- + * Make group "g2/g1" on file2 + *------------------------------------------------------------------------- + */ + + /* Create a group. */ + group_id = H5Gcreate(file2_id, "g2", 0); + group2_id = H5Gcreate(group_id, "g1", 0); + + /* Close */ + status = H5Gclose(group_id); + status = H5Gclose(group2_id); + +/*------------------------------------------------------------------------- + * Make dataset "dset4" on file2 + *------------------------------------------------------------------------- + */ /* Create a data space */ - space_id = H5Screate_simple(1,dims1,NULL); + space_id = H5Screate_simple(1,dims,NULL); /* Create a dataset "dset" */ - dataset_id = H5Dcreate(file_id,"dset",H5T_NATIVE_INT,space_id,H5P_DEFAULT); + dataset_id = H5Dcreate(file2_id,"dset4",H5T_NATIVE_FLOAT,space_id,H5P_DEFAULT); /* Write the data */ - status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data1); + status = H5Dwrite(dataset_id,H5T_NATIVE_FLOAT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data4); - /* Close */ - status = H5Dclose(dataset_id); + /* Close */ + status = H5Dclose(dataset_id); status = H5Sclose(space_id); - status = H5Fclose(file_id); /*------------------------------------------------------------------------- - * Create another file + * Make dataset "dset5" on file1 *------------------------------------------------------------------------- */ - /* Create a file */ - file_id = H5Fcreate ("h5diff_test2.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); - /* Create a data space */ - space_id = H5Screate_simple(1,dims2,NULL); + space_id = H5Screate_simple(2,dims2,NULL); /* Create a dataset "dset" */ - dataset_id = H5Dcreate(file_id,"dset",H5T_NATIVE_INT,space_id,H5P_DEFAULT); + dataset_id = H5Dcreate(file1_id,"dset5",H5T_NATIVE_FLOAT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_FLOAT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data5); - /* Write the data */ - status = H5Dwrite(dataset_id,H5T_NATIVE_INT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data2); + /* Close */ + status = H5Dclose(dataset_id); + status = H5Sclose(space_id); - /* Close */ - status = H5Dclose(dataset_id); +/*------------------------------------------------------------------------- + * Make dataset "dset6" on file2 + *------------------------------------------------------------------------- + */ + + /* Create a data space */ + space_id = H5Screate_simple(2,dims2,NULL); + + /* Create a dataset "dset" */ + dataset_id = H5Dcreate(file2_id,"dset6",H5T_NATIVE_FLOAT,space_id,H5P_DEFAULT); + + /* Write the data */ + status = H5Dwrite(dataset_id,H5T_NATIVE_FLOAT,H5S_ALL,H5S_ALL,H5P_DEFAULT,data6); + + /* Close */ + status = H5Dclose(dataset_id); status = H5Sclose(space_id); - status = H5Fclose(file_id); + +/*------------------------------------------------------------------------- + * Close files + *------------------------------------------------------------------------- + */ + status = H5Fclose(file1_id); + status = H5Fclose(file2_id); - return 0; + return 0; } + + + |