diff options
Diffstat (limited to 'tools')
-rw-r--r-- | tools/h5dump.c | 185 | ||||
-rw-r--r-- | tools/h5dump.h | 11 | ||||
-rw-r--r-- | tools/h5dumputil.c | 331 | ||||
-rwxr-xr-x | tools/testh5dump.sh | 37 |
4 files changed, 349 insertions, 215 deletions
diff --git a/tools/h5dump.c b/tools/h5dump.c index b0d19a0..043f70e 100644 --- a/tools/h5dump.c +++ b/tools/h5dump.c @@ -12,6 +12,7 @@ typedef struct shared_obj_t{ unsigned long objno[2]; char objname[1024]; int displayed; +int recorded; } shared_obj_t; typedef struct table_t{ @@ -96,10 +97,11 @@ static void indentation(int x) { static void print_datatype(hid_t type) { char *fname ; -hid_t nmembers, mtype; -int i, j, ndims, perm[H5DUMP_MAX_NDIMS]; -size_t dims[H5DUMP_MAX_NDIMS]; - +hid_t nmembers, mtype, str_type; +int i, j, ndims, perm[H5DUMP_MAX_RANK]; +size_t dims[H5DUMP_MAX_RANK], size; +H5T_str_t str_pad; +H5T_cset_t cset; H5G_stat_t statbuf; switch (H5Tget_class(type)) { @@ -187,17 +189,69 @@ H5G_stat_t statbuf; case H5T_TIME: printf( "H5T_TIME: not yet implemented"); + break; case H5T_STRING: - if (H5Tequal(type,H5T_C_S1)) - printf( "H5T_C_S1"); - else if (H5Tequal(type,H5T_FORTRAN_S1)) - printf( "H5T_FORTRAN_S1"); - else { - printf( "undefined string"); + size = H5Tget_size(type); + str_pad = H5Tget_strpad(type) ; + cset = H5Tget_cset(type); + + indentation (indent+col); + printf("%s %s %d;\n", BEGIN, STRSIZE, size); + + indentation (indent+col); + printf(" %s ", STRPAD); + if (str_pad == H5T_STR_NULLTERM ) + + printf("H5T_STR_NULLTERM;\n"); + + else if (str_pad == H5T_STR_NULLPAD ) + + printf("H5T_STR_NULLPAD;\n"); + + else if (str_pad == H5T_STR_SPACEPAD ) + + printf("H5T_STR_SPACEPAD;\n"); + + else + printf("H5T_STR_ERROR;\n"); + + indentation (indent+col); + printf(" %s ", CSET); + if (cset == H5T_CSET_ASCII) + printf("H5T_CSET_ASCII;\n"); + else + printf("unknown_cset;\n"); + + str_type = H5Tcopy(H5T_C_S1); + H5Tset_cset(str_type, cset ) ; + H5Tset_size(str_type, size); + H5Tset_strpad(str_type, str_pad); + + indentation (indent+col); + printf(" %s ", CTYPE); + if (H5Tequal(type,str_type)) { + printf("H5T_C_S1;\n"); + H5Tclose(str_type); + } else { + + H5Tclose(str_type); + str_type = H5Tcopy(H5T_FORTRAN_S1); + H5Tset_cset(str_type, cset ) ; + H5Tset_size(str_type, size); + H5Tset_strpad(str_type, str_pad ); + if (H5Tequal(type,str_type)) + printf( "H5T_FORTRAN_S1;\n"); + else { + printf("unknown_one_character_type;\n "); status = 1; + } + H5Tclose(str_type); } + indentation (indent+col); + printf("%s\n", END); + break; case H5T_BITFIELD: @@ -237,12 +291,11 @@ H5G_stat_t statbuf; indentation (indent+col); if (i >= 0) { - if (!type_table.objs[i].displayed) /* unamed data type */ - printf("%s %s \"#%lu:%lu\" %s\n", HARDLINK, BEGIN, - type_table.objs[i].objno[0], - type_table.objs[i].objno[1], END); + if (!type_table.objs[i].recorded) /* unamed data type */ + printf("\"#%lu:%lu\"\n", type_table.objs[i].objno[0], + type_table.objs[i].objno[1]); else - printf("%s %s \"%s\" %s\n", HARDLINK, BEGIN,type_table.objs[i].objname, END); + printf("\"%s\"\n", type_table.objs[i].objname); } else { printf("h5dump error: unknown committed type.\n"); status = 1; @@ -264,7 +317,7 @@ H5G_stat_t statbuf; print_datatype(mtype); - printf (" %s", fname); + printf (" \"%s\"", fname); if (ndims != 1 || dims[0] != 1) { for (j = 0; j < ndims; j++) @@ -326,15 +379,16 @@ dump_datatype (hid_t type) { indent += col; indentation (indent); - if (H5Tget_class(type) == H5T_COMPOUND) { + if (H5Tget_class(type) == H5T_COMPOUND || + H5Tget_class(type) == H5T_STRING ) { printf ("%s %s\n", DATATYPE, BEGIN); print_datatype(type); indentation (indent); printf ("%s\n", END); } else { - printf ("%s %s \"", DATATYPE, BEGIN); + printf ("%s %s ", DATATYPE, BEGIN); print_datatype(type); - printf ("\" %s\n", END); + printf (" %s\n", END); } indent -= col; } @@ -356,8 +410,8 @@ dump_datatype (hid_t type) { static void dump_dataspace (hid_t space) { - hsize_t size[H5DUMP_MAX_NDIMS]; - hsize_t maxsize[H5DUMP_MAX_NDIMS]; + hsize_t size[H5DUMP_MAX_RANK]; + hsize_t maxsize[H5DUMP_MAX_RANK]; int ndims = H5Sget_simple_extent_dims(space, size, maxsize); int i; @@ -367,13 +421,13 @@ dump_dataspace (hid_t space) if (H5Sis_simple(space)) { - if (ndims == 0) /* scalar space */ - HDfprintf (stdout, "%s %s ( 0 ) ( 0 ) %s\n",BEGIN, ARRAY, END); - else { - HDfprintf (stdout, "%s %s ( %Hu",BEGIN, ARRAY, size[0]); + if (ndims == 0) /* scalar dataspace */ + HDfprintf (stdout, "%s %s %s\n", BEGIN, SCALAR, END); + else { /* simple dataspace */ + HDfprintf (stdout, "%s %s ( %Hu",BEGIN, SIMPLE, size[0]); for (i = 1; i < ndims; i++) HDfprintf (stdout, ", %Hu", size[i]); - printf(" ) "); + printf(" ) / "); if (maxsize[0]==H5S_UNLIMITED) HDfprintf (stdout, "( %s", "H5S_UNLIMITED"); @@ -466,21 +520,23 @@ char *obj_name, *attr_name; hid_t oid, attr_id, type, space; H5G_stat_t statbuf; - j = strlen(name)-1; - obj_name = malloc ((j+2) * sizeof(char)); + obj_name = malloc ((j+2)* sizeof(char)); /* find the last / */ while (name[j] != '/' && j >=0) j--; + /* object name */ if (j == -1) strcpy(obj_name, "/"); - else strncpy(obj_name, name, j+1); - - attr_name = name+j+1; + else { + strncpy(obj_name, name, j+1); + obj_name[j+1] = '\0'; + } + attr_name = name+j+1; begin_obj (ATTRIBUTE, name); - H5Gget_objinfo(loc_id, obj_name, FALSE, &statbuf); + H5Gget_objinfo(loc_id, obj_name, FALSE , &statbuf); switch (statbuf.type) { case H5G_GROUP: if ((oid = H5Gopen (loc_id, obj_name))<0) { @@ -650,8 +706,7 @@ int i; indentation (indent); begin_obj(DATASET, name); indentation (indent+col); - printf("%s %s \"%s\" %s\n",HARDLINK, BEGIN, - dset_table.objs[i].objname,END); + printf("%s \"%s\"\n", HARDLINK, dset_table.objs[i].objname); indentation (indent); end_obj(); goto done; @@ -711,8 +766,8 @@ static void dump_named_datatype (hid_t type, const char *name) { char *fname ; hid_t nmembers, mtype; -int i, j, ndims, perm[H5DUMP_MAX_NDIMS]; -size_t dims[H5DUMP_MAX_NDIMS]; +int i, j, ndims, perm[H5DUMP_MAX_RANK]; +size_t dims[H5DUMP_MAX_RANK]; indentation (indent); begin_obj(DATATYPE, name); @@ -731,7 +786,7 @@ size_t dims[H5DUMP_MAX_NDIMS]; print_datatype(mtype); - printf (" %s", fname); + printf (" \"%s\"", fname); if (ndims != 1 || dims[0] != 1) { for (j = 0; j < ndims; j++) @@ -777,7 +832,7 @@ int i; if (!strcmp(name,"/") && unamedtype) { /* dump unamed type in root group */ for (i = 0; i < type_table.nobjs; i++) - if (!type_table.objs[i].displayed) { + if (!type_table.objs[i].recorded) { dset = H5Dopen (gid, type_table.objs[i].objname); type = H5Dget_type (dset); sprintf(typename,"#%lu:%lu", type_table.objs[i].objno[0], @@ -801,12 +856,12 @@ int i; } else if (group_table.objs[i].displayed) { indentation (indent); - printf("%s %s \"%s\" %s\n",HARDLINK, BEGIN, group_table.objs[i].objname,END); + printf("%s \"%s\"\n",HARDLINK, group_table.objs[i].objname); } else { strcpy(group_table.objs[i].objname, prefix); - group_table.objs[i].displayed=1; + group_table.objs[i].displayed = 1; H5Aiterate (gid, NULL, dump_attr, NULL); H5Giterate (gid, ".", NULL, dump_all, NULL); @@ -930,20 +985,26 @@ int i; for (i = 0; i < group_table.size; i++) { group_table.objs[i].objno[0] = group_table.objs[i].objno[1] = 0; group_table.objs[i].displayed = 0; + group_table.objs[i].recorded = 0; } for (i = 0; i < dset_table.size; i++) { dset_table.objs[i].objno[0] = dset_table.objs[i].objno[1] = 0; dset_table.objs[i].displayed = 0; + dset_table.objs[i].recorded = 0; } for (i = 0; i < type_table.size; i++) { - dset_table.objs[i].objno[0] = dset_table.objs[i].objno[1] = 0; - dset_table.objs[i].displayed = 0; + type_table.objs[i].objno[0] = type_table.objs[i].objno[1] = 0; + type_table.objs[i].displayed = 0; + type_table.objs[i].recorded = 0; } prefix = (char *) malloc(prefix_len * sizeof (char)); + *prefix = '\0'; +/* strcpy(prefix, ""); +*/ } @@ -999,6 +1060,7 @@ int i; for (i = table->nobjs; i < table->size; i++) { table->objs[i].objno[0] = table->objs[i].objno[1] = 0; table->objs[i].displayed = 0; + table->objs[i].recorded = 0; } } @@ -1029,24 +1091,27 @@ int i; printf("group_table: # of entries = %d\n", group_table.nobjs); for ( i = 0; i < group_table.nobjs; i++) - printf ("%ul %ul %s\n %d", group_table.objs[i].objno[0], + printf ("%ul %ul %s\n %d %d", group_table.objs[i].objno[0], group_table.objs[i].objno[1], group_table.objs[i].objname, - group_table.objs[i].displayed); + group_table.objs[i].displayed, + group_table.objs[i].recorded); printf("\ndset_table: # of entries = %d\n", dset_table.nobjs); for ( i = 0; i < dset_table.nobjs; i++) - printf ("%ul %ul %s %d\n", dset_table.objs[i].objno[0], + printf ("%ul %ul %s %d %d\n", dset_table.objs[i].objno[0], dset_table.objs[i].objno[1], dset_table.objs[i].objname, - dset_table.objs[i].displayed); + dset_table.objs[i].displayed, + dset_table.objs[i].recorded); printf("\ntype_table: # of entries = %d\n", type_table.nobjs); for ( i = 0; i < type_table.nobjs; i++) - printf ("%ul %ul %s %d\n", type_table.objs[i].objno[0], + printf ("%ul %ul %s %d %d\n", type_table.objs[i].objno[0], type_table.objs[i].objno[1], type_table.objs[i].objname, - type_table.objs[i].displayed); + type_table.objs[i].displayed, + type_table.objs[i].recorded); } */ @@ -1118,10 +1183,8 @@ int i; type = H5Dget_type (obj); if (H5Tcommitted(type) > 0 ) { H5Gget_objinfo(type, ".", TRUE, &statbuf); - if (search_obj (type_table, statbuf.objno) < 0) { + if (search_obj (type_table, statbuf.objno) < 0) add_obj (&type_table, statbuf.objno, tmp) ; - type_table.objs[type_table.nobjs-1].displayed = 0; - } } H5Tclose(type); H5Dclose (obj); @@ -1136,10 +1199,10 @@ int i; i = search_obj (type_table, statbuf.objno); if (i < 0) { add_obj (&type_table, statbuf.objno, tmp) ; - type_table.objs[type_table.nobjs-1].displayed = 1; /* named data type */ + type_table.objs[type_table.nobjs-1].recorded = 1; /* named data type */ } else { strcpy (type_table.objs[i].objname, tmp); - type_table.objs[i].displayed = 1; /* named data type */ + type_table.objs[i].recorded = 1; } break; @@ -1192,7 +1255,8 @@ dump_data (hid_t obj_id, int obj_data) { * If the dataset is a 1-byte integer type then format it as an ASCI * character string instead of integers. */ - if (1==size && H5T_INTEGER==H5Tget_class(f_type)) { + if ((1==size && H5T_INTEGER==H5Tget_class(f_type)) || + (H5T_STRING == H5Tget_class(f_type))) { info.elmt_suf1 = ""; info.elmt_suf2 = ""; info.idx_fmt = " (%s) \""; @@ -1211,6 +1275,7 @@ dump_data (hid_t obj_id, int obj_data) { if (h5dump1(stdout, &info, obj_id, -1, obj_data)<0) { indentation(indent+col); printf("Unable to print data.\n"); + status = 1; } if (1==size && H5T_INTEGER==H5Tget_class(f_type)) @@ -1351,9 +1416,9 @@ H5Eset_auto (NULL, NULL); H5Giterate (fid, "/", NULL, find_shared_objs, NULL); strcpy(prefix, ""); - /* assign names to unamed shared data type */ + /* does there exist unamed committed data type */ for ( i = 0; i < type_table.nobjs; i++) - if (type_table.objs[i].displayed == 0) unamedtype = 1; + if (type_table.objs[i].recorded == 0) unamedtype = 1; /* #ifdef H5DUMP_DEBUG @@ -1366,6 +1431,7 @@ H5Eset_auto (NULL, NULL); goto done; } + /* start to dump */ begin_obj("HDF5", fname); if (display_bb) dump_bb(); @@ -1391,7 +1457,6 @@ H5Eset_auto (NULL, NULL); for (curr_arg = opts[i]+1; curr_arg < ((i+1)==nopts?(argc-1):opts[i+1]); curr_arg++) - dump_selected_attr (fid, argv[curr_arg]); } else if (!strcmp(argv[opts[i]],"-d")) { @@ -1414,8 +1479,8 @@ H5Eset_auto (NULL, NULL); if (dset_table.objs[index].displayed) { begin_obj(DATASET, argv[curr_arg]); indentation (indent+col); - printf("%s %s \"%s\" %s\n",HARDLINK, BEGIN, - dset_table.objs[index].objname,END); + printf("%s \"%s\"\n", HARDLINK, + dset_table.objs[index].objname); indentation (indent); end_obj(); } else { @@ -1503,7 +1568,7 @@ H5Eset_auto (NULL, NULL); /* check if argv[curr_arg] is unamed data type */ index = 0; while (index < type_table.nobjs ) { - if (!type_table.objs[index].displayed) { /* unamed data type */ + if (!type_table.objs[index].recorded) { /* unamed data type */ sprintf(name,"#%lu:%lu\n", type_table.objs[index].objno[0], type_table.objs[index].objno[1]); sprintf(name1,"/#%lu:%lu\n", type_table.objs[index].objno[0], diff --git a/tools/h5dump.h b/tools/h5dump.h index 02b3166..0a55508 100644 --- a/tools/h5dump.h +++ b/tools/h5dump.h @@ -10,8 +10,9 @@ #define DATATYPE "DATATYPE" #define DATASPACE "DATASPACE" #define DATA "DATA" -#define ARRAY "ARRAY" -#define OTHER "OTHER" +#define SCALAR "SCALAR" +#define SIMPLE "SIMPLE" +#define COMPLEX "COMPLEX" #define STORAGELAYOUT "STORAGELAYOUT" #define COMPRESSION "COMPRESSION" #define EXTERNAL "EXTERNAL" @@ -20,6 +21,10 @@ #define NLINK "NLINK" #define FILENO "FILENO" #define OBJNO "OBJNO" +#define STRSIZE "STRSIZE" +#define STRPAD "STRPAD" +#define CSET "CSET" +#define CTYPE "CTYPE" #define BEGIN "{" #define END "}" @@ -27,7 +32,7 @@ #define ATTRIBUTE_DATA 0 #define DATASET_DATA 1 -#define H5DUMP_MAX_NDIMS 64 +#define H5DUMP_MAX_RANK H5S_MAX_RANK #define begin_obj(obj,name) printf("%s \"%s\" %s\n", obj, name, BEGIN) #define end_obj() printf("%s\n", END); diff --git a/tools/h5dumputil.c b/tools/h5dumputil.c index d033092..4ae9012 100644 --- a/tools/h5dumputil.c +++ b/tools/h5dumputil.c @@ -36,69 +36,6 @@ extern int indent; extern int ischar; - -/*------------------------------------------------------------------------- - * Function: h5dump_prefix - * - * Purpose: Prints the prefix to show up at the begining of the line. - * - * Return: void - * - * Programmer: Robb Matzke - * Thursday, July 23, 1998 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -/* -static void -h5dump_prefix(char *s, const h5dump_t *info, hsize_t elmtno, int ndims, - hsize_t min_idx[], hsize_t max_idx[]) -{ - hsize_t p_prod[8], p_idx[8]; - hsize_t n, i; - char temp[1024]; - -*/ - /* - * Calculate the number of elements represented by a unit change in a - * certain index position. - */ -/* - for (i=ndims-1, p_prod[ndims-1]=1; i>0; --i) { - p_prod[i-1] = (max_idx[i]-min_idx[i]) * p_prod[i]; - } -*/ - /* - * Calculate the index values from the element number. - */ -/* - for (i=0, n=elmtno; i<(hsize_t)ndims; i++) { - p_idx[i] = n / p_prod[i] + min_idx[i]; - n %= p_prod[i]; - } -*/ - /* - * Print the index values. - */ -/* - *temp = '\0'; - for (i=0; i<(hsize_t)ndims; i++) { - if (i) strcat(temp, OPT(info->idx_sep, ",")); - sprintf(temp+strlen(temp), OPT(info->idx_n_fmt, "%lu"), - (unsigned long)p_idx[i]); - } -*/ - - /* - * Add prefix and suffix to the index. - */ -/* - sprintf(s, OPT(info->idx_fmt, "%s: "), temp); -} -*/ - /*------------------------------------------------------------------------- * Function: h5dump_sprint * @@ -121,11 +58,23 @@ h5dump_sprint(char *s/*out*/, const h5dump_t *info, hid_t type, void *vp) char temp[1024], *name; hid_t memb; int nmembs, j, k, ndims; - + if (H5Tequal(type, H5T_NATIVE_DOUBLE)) { sprintf(temp, "%g", *((double*)vp)); } else if (H5Tequal(type, H5T_NATIVE_FLOAT)) { sprintf(temp, "%g", *((float*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_SHORT)) { + sprintf(temp, "%d", *((short*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_USHORT)) { + sprintf(temp, "%u", *((unsigned short*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_INT)) { + sprintf(temp, "%d", *((int*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_UINT)) { + sprintf(temp, "%u", *((unsigned*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_LONG)) { + sprintf(temp, "%ld", *((long*)vp)); + } else if (H5Tequal(type, H5T_NATIVE_ULONG)) { + sprintf(temp, "%lu", *((unsigned long*)vp)); } else if (H5Tequal(type, H5T_NATIVE_SCHAR) || H5Tequal(type, H5T_NATIVE_UCHAR)) { switch (*((char*)vp)) { @@ -155,19 +104,7 @@ h5dump_sprint(char *s/*out*/, const h5dump_t *info, hid_t type, void *vp) else sprintf(temp, "\\%03o", *((unsigned char*)vp)); break; } - } else if (H5Tequal(type, H5T_NATIVE_SHORT)) { - sprintf(temp, "%d", *((short*)vp)); - } else if (H5Tequal(type, H5T_NATIVE_USHORT)) { - sprintf(temp, "%u", *((unsigned short*)vp)); - } else if (H5Tequal(type, H5T_NATIVE_INT)) { - sprintf(temp, "%d", *((int*)vp)); - } else if (H5Tequal(type, H5T_NATIVE_UINT)) { - sprintf(temp, "%u", *((unsigned*)vp)); - } else if (H5Tequal(type, H5T_NATIVE_LONG)) { - sprintf(temp, "%ld", *((long*)vp)); - } else if (H5Tequal(type, H5T_NATIVE_ULONG)) { - sprintf(temp, "%lu", *((unsigned long*)vp)); - } else if (H5T_COMPOUND==H5Tget_class(type)) { + } else if (H5T_COMPOUND==H5Tget_class(type)) { nmembs = H5Tget_nmembers(type); strcpy(temp, OPT(info->cmpd_pre, "{")); for (j=0; j<nmembs; j++) { @@ -196,6 +133,10 @@ h5dump_sprint(char *s/*out*/, const h5dump_t *info, hid_t type, void *vp) H5Tclose(memb); } strcat(temp, OPT(info->cmpd_suf, "}")); + } else if (H5T_STRING==H5Tget_class(type)) { + for (i = 0; i < H5Tget_size(type); i++) + h5dump_sprint(s+i, info, H5T_NATIVE_SCHAR, (char*)vp+i); + goto done; } else { strcpy(temp, "0x"); n = H5Tget_size(type); @@ -205,8 +146,12 @@ h5dump_sprint(char *s/*out*/, const h5dump_t *info, hid_t type, void *vp) } sprintf(s, OPT(info->elmt_fmt, "%s"), temp); +done: + ; + } + /*------------------------------------------------------------------------- * Function: h5dump_simple @@ -226,14 +171,14 @@ h5dump_sprint(char *s/*out*/, const h5dump_t *info, hid_t type, void *vp) *------------------------------------------------------------------------- */ static int -h5dump_simple(FILE *stream, const h5dump_t *info, hid_t dset, hid_t p_type, int data_flag) +h5dump_simple(FILE *stream, const h5dump_t *info, hid_t oid, hid_t p_type, int obj_type) { hid_t f_space; /*file data space */ int ndims; /*dimensionality */ hsize_t elmtno, i; /*counters */ - int carry; /*counter carry value */ + /*int carry;*/ /*counter carry value */ hssize_t zero[8]; /*vector of zeros */ - int need_prefix=1; /*indices need printing */ + /* int need_prefix=1;*/ /*indices need printing */ /* Print info */ hsize_t p_min_idx[8]; /*min selected index */ @@ -241,9 +186,10 @@ h5dump_simple(FILE *stream, const h5dump_t *info, hid_t dset, hid_t p_type, int size_t p_type_nbytes; /*size of memory type */ hsize_t p_nelmts; /*total selected elmts */ char p_buf[256]; /*output string */ - size_t p_column=0; /*output column */ +/* size_t p_column=0;*/ /*output column */ size_t p_ncolumns=80; /*default num columns */ - char p_prefix[1024]; /*line prefix string */ +/* + char p_prefix[1024]; */ /*line prefix string */ /* Stripmine info */ hsize_t sm_size[8]; /*stripmine size */ @@ -256,19 +202,31 @@ h5dump_simple(FILE *stream, const h5dump_t *info, hid_t dset, hid_t p_type, int hssize_t hs_offset[8]; /*starting offset */ hsize_t hs_size[8]; /*size this pass */ hsize_t hs_nelmts; /*elements in request */ -int j, print_once=1; +int j, row_size, ntokens; +char out_buf[80], prefix[80]; + +/* + prefix[0] = NULL; + out_buf[0] = NULL; +*/ + prefix[0] = '\0'; + out_buf[0] = '\0'; + if ((indent+col) > 80) indent = 0; + for (j=0;j<indent+col;j++) strcat (prefix, " "); + + if (obj_type == DATASET_DATA) + f_space = H5Dget_space(oid); + else + f_space = H5Aget_space(oid); + /* * Check that everything looks okay. The dimensionality must not be too * great and the dimensionality of the items selected for printing must * match the dimensionality of the dataset. */ - if (data_flag == DATASET_DATA) - f_space = H5Dget_space(dset); - else - f_space = H5Aget_space(dset); - ndims = H5Sget_simple_extent_ndims(f_space); + if ((size_t)ndims>NELMTS(sm_size)) return -1; /* Assume entire data space to be printed */ @@ -300,6 +258,40 @@ int j, print_once=1; /* The stripmine loop */ memset(hs_offset, 0, sizeof hs_offset); memset(zero, 0, sizeof zero); + + if (ndims == 0) { /* scalar space */ + + if (obj_type == DATASET_DATA) { + if (H5Dread(oid, p_type, sm_space, f_space, H5P_DEFAULT, sm_buf)<0) + return -1; + } else { + if (H5Aread(oid, p_type, sm_buf) < 0) return -1; + } + + if (ischar) { + + for (j = 0; j < (int)p_type_nbytes; j++) { + h5dump_sprint(p_buf, info, H5T_NATIVE_SCHAR, sm_buf+j); + if ((strlen(out_buf)+strlen(p_buf)+5) > (80-strlen(prefix))) { + printf("%s\"%s\" //\n", prefix, out_buf); + strcpy(out_buf, p_buf); + if ( j+1 == (int)p_type_nbytes ) + printf("%s\"%s\"\n", prefix, out_buf); + } else { + strcat(out_buf, p_buf); + if ( j+1 == (int)p_type_nbytes ) + printf("%s\"%s\"\n", prefix, out_buf); + } + } + } else { + h5dump_sprint(p_buf, info, p_type, sm_buf); + printf("%s%s\n", prefix, p_buf); + } + + return 0 ; + } + +row_size = 0; for (elmtno=0; elmtno<p_nelmts; elmtno+=hs_nelmts) { /* Calculate the hyperslab size */ @@ -314,25 +306,81 @@ int j, print_once=1; H5Sselect_hyperslab(sm_space, H5S_SELECT_SET, zero, NULL, &hs_nelmts, NULL); - if (data_flag == DATASET_DATA) { - if (H5Dread(dset, p_type, sm_space, f_space, H5P_DEFAULT, sm_buf)<0) { + if (obj_type == DATASET_DATA) { + if (H5Dread(oid, p_type, sm_space, f_space, H5P_DEFAULT, sm_buf)<0) { return -1; } } else { - if (H5Aread(dset, p_type, sm_buf) < 0) + if (H5Aread(oid, p_type, sm_buf) < 0) return -1; } - /* Print the data */ + + if (ischar) + ntokens = 5; + else + ntokens = 1; + for (i=0; i<hs_nelmts; i++) { - /* Render the element */ + row_size++; + h5dump_sprint(p_buf, info, p_type, sm_buf+i*p_type_nbytes); + + if (H5T_COMPOUND==H5Tget_class(p_type)) + printf("%s%s\n", prefix, p_buf); + + else { + if ((strlen(out_buf)+strlen(p_buf)+ntokens) > (80-strlen(prefix))) { + printf("%s", prefix); + if (ischar) printf("\""); + printf("%s", out_buf); + if (ischar) printf("\" //"); + printf("\n"); + + strcpy(out_buf, p_buf); + if (row_size == (int)p_max_idx[ndims-1]) { + printf("%s", prefix); + if (ischar) printf("\""); + printf("%s", out_buf); + *out_buf = '\0'; + if (ischar) printf("\""); + if ((i+1) != hs_nelmts) { + if (ischar) printf(" //"); + else printf(","); + } + printf("\n"); + row_size = 0; + } else + if (!ischar) strcat(out_buf, ", "); + } else { + strcat(out_buf, p_buf); + if (row_size == (int)p_max_idx[ndims-1]) { + printf("%s", prefix); + if (ischar) printf("\""); + printf("%s", out_buf); + if (ischar) printf("\""); + if ((i+1) != hs_nelmts) { + if (ischar) printf(" //"); + else printf(","); + } + printf("\n"); + *out_buf = '\0'; + row_size = 0; + } else + if (!ischar) strcat(out_buf, ", "); + + } + } + + } +/* + for (i=0; i<hs_nelmts; i++) { + Render the element h5dump_sprint(p_buf, info, p_type, sm_buf+i*p_type_nbytes); if (elmtno+i+1<p_nelmts) { strcat(p_buf, OPT(info->elmt_suf1, ",")); } - /* Print the prefix */ if ((p_column + strlen(p_buf) + strlen(OPT(info->elmt_suf2, " ")) + @@ -340,10 +388,6 @@ int j, print_once=1; need_prefix = 1; } if (need_prefix) { -/* - h5dump_prefix(p_prefix, info, elmtno+i, ndims, - p_min_idx, p_max_idx); -*/ if (p_column) { fputs(OPT(info->line_suf, ""), stream); @@ -356,9 +400,7 @@ int j, print_once=1; putc('"', stream); print_once=0; } -/* fputs(p_prefix, stream); -*/ p_column = strlen(p_prefix); need_prefix = 0; } else { @@ -370,8 +412,10 @@ int j, print_once=1; p_column += strlen(p_buf); } +*/ /* Calculate the next hyperslab offset */ +/* for (i=ndims, carry=1; i>0 && carry; --i) { hs_offset[i-1] += hs_size[i-1]; if (hs_offset[i-1]==(hssize_t)p_max_idx[i-1]) { @@ -380,19 +424,18 @@ int j, print_once=1; carry = 0; } } +*/ } - if (p_column) { /* + if (p_column) { fputs(OPT(info->line_suf, ""), stream); -*/ if (ischar) putc('"',stream); putc('\n', stream); -/* fputs(OPT(info->line_sep, ""), stream); -*/ } +*/ H5Sclose(sm_space); H5Sclose(f_space); return 0; @@ -422,8 +465,10 @@ h5dump_fixtype(hid_t f_type) hid_t m_type=-1, f_memb; hid_t *memb=NULL; char **name=NULL; - int nmembs=0, i, j, *ndims=NULL; + int nmembs=0, i, j, *ndims=NULL, perm[4]; size_t size, offset, *dims=NULL, nelmts; + H5T_str_t strpad; + size = H5Tget_size(f_type); switch (H5Tget_class(f_type)) { @@ -487,7 +532,7 @@ h5dump_fixtype(hid_t f_type) if (memb[i]<0) goto done; /* Get the member dimensions */ - ndims[i] = H5Tget_member_dims(f_type, i, dims+i*4, NULL); + ndims[i] = H5Tget_member_dims(f_type, i, dims+i*4, perm); assert(ndims[i]>=0 && ndims[i]<=4); for (j=0, nelmts=1; j<ndims[i]; j++) nelmts *= dims[i*4+j]; @@ -506,15 +551,33 @@ h5dump_fixtype(hid_t f_type) m_type = H5Tcreate(H5T_COMPOUND, size); for (i=0, offset=0; i<nmembs; i++) { H5Tinsert_array(m_type, name[i], offset, ndims[i], dims+i*4, - NULL, memb[i]); + perm, memb[i]); for (j=0, nelmts=1; j<ndims[i]; j++) nelmts *= dims[i*4+j]; offset = ALIGN(offset, H5Tget_size(memb[i])) + nelmts * H5Tget_size(memb[i]); } +/* free name */ break; case H5T_TIME: case H5T_STRING: + + m_type = H5Tcopy(H5T_C_S1); + H5Tset_size(m_type, size); + strpad = H5Tget_strpad(f_type) ; + H5Tset_strpad(m_type, strpad); + + if (H5Tequal(m_type,f_type) < 0) { + H5Tclose(m_type); + m_type = H5Tcopy(H5T_FORTRAN_S1); + H5Tset_size(m_type, size); + H5Tset_strpad(m_type, strpad); + if (H5Tequal(m_type,f_type) < 0) + m_type = -1; + } + + break; + case H5T_BITFIELD: case H5T_OPAQUE: /* @@ -545,7 +608,7 @@ h5dump_fixtype(hid_t f_type) /*------------------------------------------------------------------------- - * Function: h5dump + * Function: h5dump1 * * Purpose: Print some values from a dataset DSET to the file STREAM * after converting all types to P_TYPE (which should be a @@ -564,17 +627,14 @@ h5dump_fixtype(hid_t f_type) *------------------------------------------------------------------------- */ int -h5dump1(FILE *stream, const h5dump_t *info, hid_t dset, hid_t _p_type, int data_flag) +h5dump1 (FILE *stream, const h5dump_t *info, hid_t oid, hid_t _p_type, int obj_type) { hid_t f_space; hid_t p_type = _p_type; hid_t f_type; - int status; + int status = -1; h5dump_t info_dflt; - char p_buf[256], sm_buf[256]; /*tmp for scala */ - int j; - /* Use default values */ if (!stream) stream = stdout; if (!info) { @@ -582,42 +642,35 @@ h5dump1(FILE *stream, const h5dump_t *info, hid_t dset, hid_t _p_type, int data_ info = &info_dflt; } - if (p_type<0) { - if (data_flag == DATASET_DATA) - f_type = H5Dget_type(dset); + if (p_type < 0) { + + if (obj_type == DATASET_DATA) + f_type = H5Dget_type(oid); else - f_type = H5Aget_type(dset); + f_type = H5Aget_type(oid); + + if (f_type < 0) return status; p_type = h5dump_fixtype(f_type); + H5Tclose(f_type); - if (p_type<0) return -1; + + if (p_type < 0) return status; } /* Check the data space */ - if (data_flag == DATASET_DATA) - f_space = H5Dget_space(dset); + if (obj_type == DATASET_DATA) + f_space = H5Dget_space(oid); else - f_space = H5Aget_space(dset); + f_space = H5Aget_space(oid); + + if (f_space < 0) return status; - if (H5Sis_simple(f_space)<=0) return -1; - - - if ( H5Sget_simple_extent_ndims(f_space) == 0){ - if (data_flag == ATTRIBUTE_DATA) { - if (H5Aread(dset, p_type, sm_buf) < 0) - return -1; - } else return -1; - h5dump_sprint(p_buf, info, p_type, sm_buf); - for (j=0;j<indent+col;j++) putc(' ', stream); - fputs(p_buf, stream); - putc('\n', stream); - H5Sclose(f_space); - } else { + if (H5Sis_simple(f_space) >= 0) + status = h5dump_simple(stream, info, oid, p_type, obj_type); + H5Sclose(f_space); - /* Print the data */ - status = h5dump_simple(stream, info, dset, p_type, data_flag); - } - if (p_type!=_p_type) H5Tclose(p_type); + if (p_type != _p_type) H5Tclose(p_type); return status; } diff --git a/tools/testh5dump.sh b/tools/testh5dump.sh index cb68db7..5645d60 100755 --- a/tools/testh5dump.sh +++ b/tools/testh5dump.sh @@ -69,35 +69,46 @@ DUMP() ############################################################################## ############################################################################## +# test for displaying groups DUMP tgroup-1.ddl tgroup.h5 -DUMP tgroup-2.ddl -g / tgroup.h5 -DUMP tgroup-3.ddl -g /g2 /y tgroup.h5 +# test for displaying the selected groups +DUMP tgroup-2.ddl -g /g2 / /y tgroup.h5 +# test for displaying simple space datasets DUMP tdset-1.ddl tdset.h5 -DUMP tdset-2.ddl -d dset1 /dset2 tdset.h5 -DUMP tdset-3.ddl -d /dset1 -header tdset.h5 -DUMP tdset-4.ddl -d dset3 tdset.h5 +# test for displaying selected datasets +DUMP tdset-2.ddl -header -d dset1 /dset2 dset3 tdset.h5 +# test for displaying attributes DUMP tattr-1.ddl tattr.h5 -DUMP tattr-2.ddl -a attr1 attr3 tattr.h5 -DUMP tattr-3.ddl -header -a attr2 tattr.h5 -DUMP tattr-4.ddl -a attr4 tattr.h5 +# test for displaying the selected attributes of string type and scalar space +DUMP tattr-2.ddl -a attr1 attr4 attr5 tattr.h5 +# test for header and error messages +DUMP tattr-3.ddl -header -a attr2 attr tattr.h5 +# test for displaying soft links DUMP tslink-1.ddl tslink.h5 +# test for displaying the selected link DUMP tslink-2.ddl -l slink2 tslink.h5 +# tests for hard links DUMP thlink-1.ddl thlink.h5 -DUMP thlink-2.ddl -d /g1/link2 /dset /g1/link1/link3 thlink.h5 -DUMP thlink-3.ddl -d /dset /g1/link1/link3 /g1/link2 thlink.h5 +DUMP thlink-2.ddl -d /g1/dset2 /dset1 /g1/g1.1/dset3 thlink.h5 +DUMP thlink-3.ddl -d /g1/g1.1/dset3 /g1/dset2 /dset1 thlink.h5 DUMP thlink-4.ddl -g /g1 thlink.h5 -DUMP thlink-5.ddl -d /dset -g /g2 -d /g1/link2 thlink.h5 +DUMP thlink-5.ddl -d /dset1 -g /g2 -d /g1/dset2 thlink.h5 +# tests for compound data types DUMP tcomp-1.ddl tcompound.h5 +# test for named data types DUMP tcomp-2.ddl -t /type1 /type2 /group1/type3 tcompound.h5 -DUMP tcomp-3.ddl -d /group2/dset5 -g /group1 tcompound.h5 -DUMP tcomp-4.ddl -t /#3432:0 -g /group2 tcompound.h5 +# test for unamed type +DUMP tcomp-3.ddl -t /#5992:0 -g /group2 tcompound.h5 +# test for options DUMP tall-1.ddl tall.h5 DUMP tall-2.ddl -header -g /g1/g1.1 -a attr2 tall.h5 DUMP tall-3.ddl -d /g2/dset2.1 -l /g1/g1.2/g1.2.1/slink tall.h5 +# test for loop detection +DUMP tloop-1.ddl tloop.h5 |