summaryrefslogtreecommitdiffstats
path: root/tools/h5dump.c
diff options
context:
space:
mode:
authorRuey-Hsia Li <rli@ncsa.uiuc.edu>1998-10-29 06:27:31 (GMT)
committerRuey-Hsia Li <rli@ncsa.uiuc.edu>1998-10-29 06:27:31 (GMT)
commit056814423f714f2e841e113d2d2c111b203886e6 (patch)
tree0eaefc9abd586f7bbad733933b23e7c4deb14bdc /tools/h5dump.c
parent54240b9657d3068757e4017670d776f9732e850f (diff)
downloadhdf5-056814423f714f2e841e113d2d2c111b203886e6.zip
hdf5-056814423f714f2e841e113d2d2c111b203886e6.tar.gz
hdf5-056814423f714f2e841e113d2d2c111b203886e6.tar.bz2
[svn-r825] Added features:
* display selected named data type * display named/unamed data type * hard link Also added testing scripts in testh5dump.sh.
Diffstat (limited to 'tools/h5dump.c')
-rw-r--r--tools/h5dump.c758
1 files changed, 648 insertions, 110 deletions
diff --git a/tools/h5dump.c b/tools/h5dump.c
index 9d77d63..d29696d 100644
--- a/tools/h5dump.c
+++ b/tools/h5dump.c
@@ -2,14 +2,33 @@
#include <stdio.h>
#include "H5private.h"
-static int indent = 0;
+int indent = 0;
+int ischar=0;
static int display_data = 1;
static int status = 0;
+static int unamedtype = 0; /* shared data type with no name */
+
+typedef struct shared_obj_t{
+unsigned long objno[2];
+char objname[1024];
+int displayed;
+} shared_obj_t;
+
+typedef struct table_t{
+int size;
+int nobjs;
+shared_obj_t *objs;
+} table_t;
+
+static int prefix_len = 1024;
+static char *prefix;
+static table_t group_table, dset_table, type_table;
static void dump_group (hid_t , const char* );
static void dump_dataset (hid_t, const char*);
static void dump_data (hid_t, int);
static void dump_named_datatype (hid_t , const char *);
+static int search_obj (table_t, unsigned long *);
/*-------------------------------------------------------------------------
* Function: usage
@@ -29,15 +48,17 @@ usage(void)
fprintf(stderr,
"\nUsage of HDF5 Dumper:\n\n \
h5dump [-h] [-bb] [-header] [-a <names>] [-d <names>] [-g <names>]\n \
- [-l <names>] <file>\n\n\
+ [-l <names>] [-t <names>] <file>\n\n\
-h Print information on this command.\n\
- -bb Display the conent of boot block. The default is not to display.\n\
- -header Display header only; that is, no data displayed.\n\
+ -bb Display the conent of the boot block. The default is not to display.\n\
+ -header Display header only; no data is displayed.\n\
-a <names> Display the specified attribute(s).\n\
-d <names> Display the specified dataset(s).\n\
- -g <names> Display all the objects within the specified group(s).\n\
- -l <names> Display the specified link value(s).\n\n");
-
+ -g <names> Display the specified group(s) and all the members.\n\
+ -l <names> Displays the value(s) of the specified soft link(s).\n\
+ -t <names> Display the specified named data type(s).\n\
+ \n\
+ <names> is one or more appropriate object names.\n\n");
}
@@ -55,18 +76,15 @@ h5dump [-h] [-bb] [-header] [-a <names>] [-d <names>] [-g <names>]\n \
*-----------------------------------------------------------------------*/
static void indentation(int x) {
- while (x>0) {
- printf(" ");
- x--;
- }
+ while (x>0) { printf(" "); x--; }
}
/*-------------------------------------------------------------------------
- * Function: datatype_name
+ * Function: print_datatype
*
- * Purpose: Prints the name of data type.
+ * Purpose: print the data type.
*
* Return: void
*
@@ -76,7 +94,13 @@ static void indentation(int x) {
*
*-----------------------------------------------------------------------*/
static void
-datatype_name(hid_t type) {
+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];
+
+H5G_stat_t statbuf;
switch (H5Tget_class(type)) {
@@ -134,7 +158,10 @@ datatype_name(hid_t type) {
printf( "H5T_NATIVE_LLONG");
else if (H5Tequal(type, H5T_NATIVE_ULLONG))
printf( "H5T_NATIVE_ULLONG");
- else printf( "unknown integer");
+ else {
+ printf( "undefined integer");
+ status = 1;
+ }
break;
case H5T_FLOAT:
@@ -152,11 +179,14 @@ datatype_name(hid_t type) {
printf( "H5T_NATIVE_DOUBLE");
else if (H5Tequal(type, H5T_NATIVE_LDOUBLE))
printf( "H5T_NATIVE_LDOUBLE");
- else printf( "unknown float");
+ else {
+ printf( "undefined float");
+ status = 1;
+ }
break;
case H5T_TIME:
- printf( "time: not yet implemented");
+ printf( "H5T_TIME: not yet implemented");
case H5T_STRING:
@@ -164,7 +194,10 @@ datatype_name(hid_t type) {
printf( "H5T_C_S1");
else if (H5Tequal(type,H5T_FORTRAN_S1))
printf( "H5T_FORTRAN_S1");
- else printf( "unknown string");
+ else {
+ printf( "undefined string");
+ status = 1;
+ }
break;
case H5T_BITFIELD:
@@ -185,20 +218,70 @@ datatype_name(hid_t type) {
printf( "H5T_STD_B64BE");
else if (H5Tequal(type, H5T_STD_B64LE))
printf( "H5T_STD_B64LE");
- else printf( "unknown bitfield");
+ else {
+ printf( "undefined bitfield");
+ status = 1;
+ }
break;
case H5T_OPAQUE:
- printf( "opaque: not yet implemented");
+ printf( "H5T_OPAQUE: not yet implemented");
break;
case H5T_COMPOUND:
- printf( "compound: not yet implemented");
+ if (H5Tcommitted(type) > 0) {
+ H5Gget_objinfo(type, ".", TRUE, &statbuf);
+
+ i = search_obj (type_table, statbuf.objno);
+
+ 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);
+ else
+ printf("%s %s \"%s\" %s\n", HARDLINK, BEGIN,type_table.objs[i].objname, END);
+ } else {
+ printf("h5dump error: unknown committed type.\n");
+ status = 1;
+ }
+
+ } else {
+
+ nmembers = H5Tget_nmembers(type);
+
+ for (i = 0; i < nmembers; i++) {
+
+ fname = H5Tget_member_name(type, i);
+
+ mtype = H5Tget_member_type(type, i);
+
+ ndims = H5Tget_member_dims(type, i, dims, perm);
+
+ indentation (indent+col);
+
+ print_datatype(mtype);
+
+ printf (" %s", fname);
+
+ if (ndims != 1 || dims[0] != 1) {
+ for (j = 0; j < ndims; j++)
+ printf("[%d]",dims[j]);
+ }
+
+ printf (";\n");
+
+ free (fname);
+ }
+ }
+
break;
default:
printf( "unknown data type");
+ status = 1;
break;
}
@@ -229,7 +312,7 @@ dump_bb(void) {
* Function: dump_datatype
*
* Purpose: Dump the data type. Data type can be HDF5 predefined
- * atomic data type, named data type or compound data type.
+ * atomic data type or committed/transient data type.
*
* Return: void
*
@@ -241,10 +324,19 @@ dump_bb(void) {
static void
dump_datatype (hid_t type) {
- indentation (indent+col);
- printf ("%s %s \"", DATATYPE, BEGIN);
- datatype_name(type);
- printf ("\" %s\n", END);
+ indent += col;
+ indentation (indent);
+ if (H5Tget_class(type) == H5T_COMPOUND) {
+ printf ("%s %s\n", DATATYPE, BEGIN);
+ print_datatype(type);
+ indentation (indent);
+ printf ("%s\n", END);
+ } else {
+ printf ("%s %s \"", DATATYPE, BEGIN);
+ print_datatype(type);
+ printf ("\" %s\n", END);
+ }
+ indent -= col;
}
@@ -264,8 +356,8 @@ dump_datatype (hid_t type) {
static void
dump_dataspace (hid_t space)
{
- hsize_t size[64];
- hsize_t maxsize[64]; /* check max dims size */
+ hsize_t size[H5DUMP_MAX_NDIMS];
+ hsize_t maxsize[H5DUMP_MAX_NDIMS];
int ndims = H5Sget_simple_extent_dims(space, size, maxsize);
int i;
@@ -275,26 +367,29 @@ dump_dataspace (hid_t space)
if (H5Sis_simple(space)) {
- HDfprintf (stdout, "%s %s ( %Hu",BEGIN, ARRAY, size[0]);
- for (i = 1; i < ndims; i++)
- HDfprintf (stdout, ", %Hu", size[i]);
- printf(" ) ");
+ 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]);
+ for (i = 1; i < ndims; i++)
+ HDfprintf (stdout, ", %Hu", size[i]);
+ printf(" ) ");
- if (maxsize[0]==H5S_UNLIMITED)
- HDfprintf (stdout, "( %s", "H5S_UNLIMITED");
- else
- HDfprintf (stdout, "( %Hu", maxsize[0]);
-
- for (i = 1; i < ndims; i++)
- if (maxsize[i]==H5S_UNLIMITED)
- HDfprintf (stdout, ", %s", "H5S_UNLIMITED");
- else
- HDfprintf (stdout, ", %Hu", maxsize[i]);
-
- printf(" ) %s\n", END);
+ if (maxsize[0]==H5S_UNLIMITED)
+ HDfprintf (stdout, "( %s", "H5S_UNLIMITED");
+ else
+ HDfprintf (stdout, "( %Hu", maxsize[0]);
+
+ for (i = 1; i < ndims; i++)
+ if (maxsize[i]==H5S_UNLIMITED)
+ HDfprintf (stdout, ", %s", "H5S_UNLIMITED");
+ else
+ HDfprintf (stdout, ", %Hu", maxsize[i]);
+
+ printf(" ) %s\n", END);
+ }
} else
-
printf("%s not yet implemented %s\n", BEGIN, END);
}
@@ -372,9 +467,9 @@ hid_t oid, attr_id, type, space;
H5G_stat_t statbuf;
-
j = strlen(name)-1;
obj_name = malloc ((j+2) * sizeof(char));
+
/* find the last / */
while (name[j] != '/' && j >=0) j--;
/* object name */
@@ -472,6 +567,7 @@ H5G_stat_t statbuf;
return SUCCEED;
}
+
/*-------------------------------------------------------------------------
* Function: dump all
*
@@ -489,14 +585,16 @@ H5G_stat_t statbuf;
static herr_t
dump_all (hid_t group, const char *name, void __unused__ *op_data)
{
- hid_t obj;
- char *buf;
- H5G_stat_t statbuf;
-
+hid_t obj;
+char *buf, *tmp;
+H5G_stat_t statbuf;
+int i;
H5Gget_objinfo(group, name, FALSE, &statbuf);
+ tmp = (char *) malloc ((strlen(prefix)+strlen(name)+2) * sizeof(char));
+ strcpy(tmp, prefix);
switch (statbuf.type) {
case H5G_LINK:
@@ -507,7 +605,7 @@ dump_all (hid_t group, const char *name, void __unused__ *op_data)
begin_obj(SOFTLINK, name);
indentation (indent+col);
if (H5Gget_linkval (group, name, statbuf.linklen, buf)>=0)
- printf ("linktarget \"%s\"\n", buf);
+ printf ("LINKTARGET \"%s\"\n", buf);
else {
printf ("h5dump error: unable to get link value.\n");
status = 1;
@@ -520,7 +618,9 @@ dump_all (hid_t group, const char *name, void __unused__ *op_data)
case H5G_GROUP:
if ((obj=H5Gopen (group, name))>=0) {
+ strcat(strcat(prefix,"/"), name);
dump_group (obj, name);
+ strcpy(prefix, tmp);
H5Gclose (obj);
} else {
printf ("h5dump error: unable to dump group %s\n",name);
@@ -532,6 +632,36 @@ dump_all (hid_t group, const char *name, void __unused__ *op_data)
case H5G_DATASET:
if ((obj=H5Dopen (group, name))>=0) {
+
+ /* hard link */
+ H5Gget_objinfo(obj, ".", TRUE, &statbuf);
+ if (statbuf.nlink > 1) {
+ i = search_obj (dset_table, statbuf.objno);
+ if (i < 0) {
+ indentation (indent);
+ begin_obj(DATASET, name);
+ indentation (indent+col);
+ printf("h5dump error: internal error\n");
+ indentation (indent);
+ end_obj();
+ status = 1;
+ goto done;
+ } else if (dset_table.objs[i].displayed) {
+ indentation (indent);
+ begin_obj(DATASET, name);
+ indentation (indent+col);
+ printf("%s %s \"%s\" %s\n",HARDLINK, BEGIN,
+ dset_table.objs[i].objname,END);
+ indentation (indent);
+ end_obj();
+ goto done;
+ } else {
+ dset_table.objs[i].displayed = 1;
+ strcat(tmp,"/");
+ strcat(tmp,name);
+ strcpy(dset_table.objs[i].objname, tmp);
+ }
+ }
dump_dataset (obj, name);
H5Dclose (obj);
} else {
@@ -542,20 +672,24 @@ dump_all (hid_t group, const char *name, void __unused__ *op_data)
break;
case H5G_TYPE:
- dump_named_datatype (obj, name);
- H5Tclose(obj);
+ if ((obj=H5Topen (group, name)) >= 0) {
+ dump_named_datatype (obj, name);
+ H5Tclose(obj);
+ } else {
+ printf ("h5dump error: unable to dump data type %s\n",name);
+ status = 1;
+ }
break;
default:
- printf ("Unknown Object %s\n", name);
+ printf ("h5dump error: unknown object %s\n", name);
status = 1;
return FAIL;
- break;
}
-
-
+done:
+ free(tmp);
return SUCCEED;
}
@@ -574,25 +708,43 @@ dump_all (hid_t group, const char *name, void __unused__ *op_data)
*
*-----------------------------------------------------------------------*/
static void
-dump_named_datatype (hid_t type_id, const char *name) {
-hid_t id;
-/*
+dump_named_datatype (hid_t type, const char *name) {
char *fname ;
hid_t nmembers, mtype;
-int i, ndims, perm[512];
-*/
-
-id = type_id; /*doesn't like warning message */
+int i, j, ndims, perm[H5DUMP_MAX_NDIMS];
+size_t dims[H5DUMP_MAX_NDIMS];
indentation (indent);
begin_obj(DATATYPE, name);
- indentation(indent+col);
- printf("named data type: not yet implemented.\n");
+ nmembers = H5Tget_nmembers(type);
+
+ for (i = 0; i < nmembers; i++) {
+
+ fname = H5Tget_member_name(type, i);
+
+ mtype = H5Tget_member_type(type, i);
+
+ ndims = H5Tget_member_dims(type, i, dims, perm);
+
+ indentation (indent+col);
+
+ print_datatype(mtype);
+
+ printf (" %s", fname);
+
+ if (ndims != 1 || dims[0] != 1) {
+ for (j = 0; j < ndims; j++)
+ printf("[%d]",dims[j]);
+ }
+ printf (";\n");
+
+ free (fname);
+ }
+
indentation (indent);
end_obj();
-
}
@@ -611,25 +763,67 @@ id = type_id; /*doesn't like warning message */
static void
dump_group (hid_t gid, const char *name) {
H5G_stat_t statbuf;
+hid_t dset, type;
+char typename[1024], *tmp;
+int i;
+
+ tmp = (char *) malloc ((strlen(prefix)+strlen(name)+2) * sizeof(char));
+ strcpy(tmp, prefix);
+
indentation (indent);
begin_obj(GROUP, name);
indent += col;
- /* hard link */
+ if (!strcmp(name,"/") && unamedtype) { /* dump unamed type in root group */
+ for (i = 0; i < type_table.nobjs; i++)
+ if (!type_table.objs[i].displayed) {
+ dset = H5Dopen (gid, type_table.objs[i].objname);
+ type = H5Dget_type (dset);
+ sprintf(typename,"#%lu:%lu", type_table.objs[i].objno[0],
+ type_table.objs[i].objno[1]);
+ dump_named_datatype (type, typename);
+ H5Tclose(type);
+ H5Dclose(dset);
+ }
+ }
+
H5Gget_objinfo(gid, ".", TRUE, &statbuf);
if (statbuf.nlink > 1) {
- indentation (indent);
- printf("%s %s (%s %u) (%s %lu %lu) (%s %lu %lu ) %s\n",
- HARDLINK, BEGIN, NLINK, statbuf.nlink, FILENO, statbuf.fileno[0],
- statbuf.fileno[1], OBJNO, statbuf.objno[0], statbuf.objno[1], END);
+
+ i = search_obj (group_table, statbuf.objno);
+ if (i < 0) {
+
+ indentation (indent);
+ printf("h5dump error: internal error\n");
+ status = 1;
+
+ } else if (group_table.objs[i].displayed) {
+
+ indentation (indent);
+ printf("%s %s \"%s\" %s\n",HARDLINK, BEGIN, group_table.objs[i].objname,END);
+
+ } else {
+
+ strcpy(group_table.objs[i].objname, prefix);
+ group_table.objs[i].displayed=1;
+ H5Aiterate (gid, NULL, dump_attr, NULL);
+ H5Giterate (gid, ".", NULL, dump_all, NULL);
+
+ }
+
+ } else {
+
+ H5Aiterate (gid, NULL, dump_attr, NULL);
+ H5Giterate (gid, ".", NULL, dump_all, NULL);
+
}
- H5Aiterate (gid, NULL, dump_attr, NULL);
- H5Giterate (gid, ".", NULL, dump_all, NULL);
indent -= col;
indentation (indent);
end_obj();
+
+ free(tmp);
}
@@ -649,20 +843,10 @@ H5G_stat_t statbuf;
static void
dump_dataset (hid_t did, const char *name) {
hid_t type, space;
-H5G_stat_t statbuf;
indentation (indent);
begin_obj(DATASET, name);
- /* hard link */
- H5Gget_objinfo(did, ".", TRUE, &statbuf);
- if (statbuf.nlink > 1) {
- indentation (indent+col);
- printf("%s %s (%s %u) (%s %lu %lu) (%s %lu %lu ) %s\n",
- HARDLINK, BEGIN, NLINK, statbuf.nlink, FILENO, statbuf.fileno[0],
- statbuf.fileno[1], OBJNO, statbuf.objno[0], statbuf.objno[1], END);
- }
-
type = H5Dget_type (did);
space = H5Dget_space (did);
dump_datatype(type);
@@ -701,10 +885,7 @@ H5G_stat_t statbuf;
printf("DATA{ not yet implemented.}\n");
break;
case H5T_COMPOUND:
- indent += col;
- indentation (indent);
- indent -= col;
- printf("DATA{ not yet implemented.}\n");
+ dump_data(did, DATASET_DATA);
break;
default: break;
}
@@ -715,6 +896,7 @@ H5G_stat_t statbuf;
H5Tclose(type);
H5Sclose(space);
+
indentation (indent);
end_obj();
@@ -722,6 +904,256 @@ H5G_stat_t statbuf;
/*-------------------------------------------------------------------------
+ * Function: init_table
+ *
+ * Purpose: allocate and initialize tables for shared groups, datasets,
+ * and committed types
+ *
+ * Return: void
+ *
+ * Programmer: Ruey-Hsia Li
+ *
+ * Modifications:
+ *
+ *-----------------------------------------------------------------------*/
+static void
+init_table (void){
+int i;
+
+ group_table.size = dset_table.size = type_table.size = 20;
+ group_table.nobjs = dset_table.nobjs = type_table.nobjs = 0;
+
+ group_table.objs = (shared_obj_t*) malloc(group_table.size*sizeof(shared_obj_t));
+ dset_table.objs = (shared_obj_t*) malloc(dset_table.size*sizeof(shared_obj_t));
+ type_table.objs = (shared_obj_t*) malloc(type_table.size*sizeof(shared_obj_t));
+
+ 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;
+ }
+
+ 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;
+ }
+
+ 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;
+ }
+
+ prefix = (char *) malloc(prefix_len * sizeof (char));
+ strcpy(prefix, "");
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: search_obj
+ *
+ * Purpose: search the object specified by objno in the table
+ *
+ * Return: an integer, the location of the object
+ * -1 if object is not found
+ *
+ *
+ * Programmer: Ruey-Hsia Li
+ *
+ * Modifications:
+ *
+ *-----------------------------------------------------------------------*/
+static int
+search_obj (table_t table, unsigned long *objno) {
+int i=0, found=0;
+
+ while (i < table.nobjs && !found)
+ if (table.objs[i].objno[0] == *(objno) &&
+ table.objs[i].objno[1] == *(objno+1) ) found = 1;
+ else i++;
+
+ if (!found) return -1;
+ else return i;
+
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: add_obj
+ *
+ * Purpose: add a shared object to the table
+ * realloc the table if necessary
+ *
+ * Return: void
+ *
+ * Programmer: Ruey-Hsia Li
+ *
+ * Modifications:
+ *
+ *-----------------------------------------------------------------------*/
+static void
+add_obj (table_t *table, unsigned long *objno, char *objname) {
+int i;
+
+ if (table->nobjs == table->size) {
+ table->size *= 2;
+ table->objs = realloc (table->objs, table->size*sizeof(shared_obj_t));
+ for (i = table->nobjs; i < table->size; i++) {
+ table->objs[i].objno[0] = table->objs[i].objno[1] = 0;
+ table->objs[i].displayed = 0;
+ }
+ }
+
+ i = table->nobjs++;
+ table->objs[i].objno[0] = *objno;
+ table->objs[i].objno[1] = *(objno+1);
+ strcpy (table->objs[i].objname, objname);
+
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: dump_tables
+ *
+ * Purpose: display the contents of tables for debugging purposes
+ *
+ * Return: void
+ *
+ * Programmer: Ruey-Hsia Li
+ *
+ * Modifications:
+ *
+ *-----------------------------------------------------------------------*/
+/*
+static void
+dump_tables(void) {
+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],
+ group_table.objs[i].objno[1],
+ group_table.objs[i].objname,
+ group_table.objs[i].displayed);
+
+ 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],
+ dset_table.objs[i].objno[1],
+ dset_table.objs[i].objname,
+ dset_table.objs[i].displayed);
+
+ 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],
+ type_table.objs[i].objno[1],
+ type_table.objs[i].objname,
+ type_table.objs[i].displayed);
+}
+*/
+
+
+/*-------------------------------------------------------------------------
+ * Function: Find_shared_objs
+ *
+ * Purpose: Find shared objects, committed types and store them in tables
+ *
+ * Return: Success: SUCCEED
+ *
+ * Failure: FAIL
+ *
+ * Programmer: Ruey-Hsia Li
+ *
+ * Modifications:
+ *
+ *-----------------------------------------------------------------------*/
+static herr_t
+find_shared_objs(hid_t group, const char *name, void __unused__ *op_data)
+{
+hid_t obj, type;
+H5G_stat_t statbuf;
+char *tmp;
+int i;
+
+ H5Gget_objinfo(group, name, TRUE, &statbuf);
+
+ tmp = (char *) malloc ((strlen(prefix)+strlen(name)+2) * sizeof(char));
+
+ strcpy(tmp, prefix);
+
+ switch (statbuf.type) {
+
+ case H5G_GROUP:
+ if ((obj=H5Gopen (group, name))>=0) {
+
+ if (prefix_len < (int)(strlen(prefix) + strlen(name) + 2)) {
+ prefix_len *= 2;
+ prefix = realloc (prefix, prefix_len * sizeof(char));
+ }
+ strcat(strcat(prefix,"/"), name);
+
+ if (statbuf.nlink > 1) {
+ if (search_obj (group_table, statbuf.objno) < 0) {
+ add_obj (&group_table, statbuf.objno, prefix);
+ H5Giterate (obj, ".", NULL, find_shared_objs, NULL);
+ }
+ } else
+ H5Giterate (obj, ".", NULL, find_shared_objs, NULL);
+
+ strcpy(prefix, tmp);
+ H5Gclose (obj);
+
+ } else
+ status = 1;
+
+ break;
+
+ case H5G_DATASET:
+
+ strcat(tmp,"/");
+ strcat(tmp,name); /* absolute name of the data set */
+ if (statbuf.nlink > 1 &&
+ search_obj (dset_table, statbuf.objno) < 0)
+ add_obj (&dset_table, statbuf.objno, tmp);
+
+ if ((obj=H5Dopen (group, name))>=0) {
+ type = H5Dget_type (obj);
+ if (H5Tcommitted(type) > 0 ) {
+ H5Gget_objinfo(type, ".", TRUE, &statbuf);
+ 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);
+ } else
+ status = 1;
+
+ break;
+
+ case H5G_TYPE:
+ strcat(tmp,"/");
+ strcat(tmp,name); /* absolute name of the type */
+ 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 */
+ } else {
+ strcpy (type_table.objs[i].objname, tmp);
+ type_table.objs[i].displayed = 1; /* named data type */
+ }
+ break;
+
+ default:
+ break;
+ }
+
+ free (tmp);
+
+ return SUCCEED;
+}
+
+
+/*-------------------------------------------------------------------------
* Function: dump_data
*
* Purpose: Dump attribute or dataset data
@@ -742,7 +1174,7 @@ dump_data (hid_t obj_id, int obj_data) {
indent += col;
indentation (indent);
- printf("%s %s", DATA, BEGIN);
+ printf("%s %s\n", DATA, BEGIN);
if (obj_data == DATASET_DATA)
f_type = H5Dget_type(obj_id);
@@ -754,7 +1186,7 @@ dump_data (hid_t obj_id, int obj_data) {
/* Set to all default values and then override */
memset(&info, 0, sizeof info);
info.idx_fmt = " (%s) ";
- info.line_ncols = 80;
+ info.line_ncols = 70 - indent;
/*
* If the dataset is a 1-byte integer type then format it as an ASCI
@@ -764,8 +1196,11 @@ dump_data (hid_t obj_id, int obj_data) {
info.elmt_suf1 = "";
info.elmt_suf2 = "";
info.idx_fmt = " (%s) \"";
+/*
info.line_suf = "\"";
- printf("\"");
+*/
+ info.line_suf = "";
+ ischar = 1;
}
@@ -774,12 +1209,14 @@ dump_data (hid_t obj_id, int obj_data) {
*/
if (h5dump1(stdout, &info, obj_id, -1, obj_data)<0) {
- printf("Unable to print data.");
+ indentation(indent+col);
+ printf("Unable to print data.\n");
}
if (1==size && H5T_INTEGER==H5Tget_class(f_type))
- printf("\"");
+ ischar = 0;
+ indentation(indent);
printf("%s\n", END);
indent -= col;
@@ -789,7 +1226,7 @@ dump_data (hid_t obj_id, int obj_data) {
/*-------------------------------------------------------------------------
- * Function:
+ * Function: main
*
* Purpose: HDF5 dumper
*
@@ -803,19 +1240,20 @@ dump_data (hid_t obj_id, int obj_data) {
*-----------------------------------------------------------------------*/
int
main(int argc, char *argv[]) {
-hid_t fid, gid, dsetid;
+hid_t fid, gid, dsetid, typeid;
hid_t plist=H5P_DEFAULT;
const char *fname = NULL;
-int i, curr_arg, display_bb=0, display_all=1;
+int i, index, curr_arg, display_bb=0, display_all=1;
int nopts=0, *opts;
-char *buf;
+char *buf, name[128], name1[128];
H5G_stat_t statbuf;
void __unused__ *op_data;
- hid_t (*func)(void*);
- void *edata;
- /* Disable error reporting */
- H5Eget_auto (&func, &edata);
- H5Eset_auto (NULL, NULL);
+void *edata;
+hid_t (*func)(void*);
+
+/* Disable error reporting */
+H5Eget_auto (&func, &edata);
+H5Eset_auto (NULL, NULL);
if (argc < 2 ) {
usage();
@@ -832,8 +1270,8 @@ void __unused__ *op_data;
opts[nopts++] = curr_arg;
if (!strcmp(argv[curr_arg],"-h")) {
-
usage();
+ free(opts);
exit(0);
} else if (!strcmp(argv[curr_arg],"-bb"))
@@ -847,11 +1285,13 @@ void __unused__ *op_data;
else if (strcmp(argv[curr_arg],"-a") &&
strcmp(argv[curr_arg],"-d") &&
strcmp(argv[curr_arg],"-g") &&
- strcmp(argv[curr_arg],"-l")) {
+ strcmp(argv[curr_arg],"-l") &&
+ strcmp(argv[curr_arg],"-t")) {
fprintf(stderr, "h5dump error: illegal option %s \n",
argv[curr_arg]);
usage();
+ free(opts);
exit(1);
} else display_all = 0;
}
@@ -861,6 +1301,7 @@ void __unused__ *op_data;
if (opts[0] == 1) { /* argv[1] is an option */
fprintf(stderr, "h5dump error: no <names> or no <file>\n");
usage();
+ free(opts);
exit(1);
}
} else {
@@ -871,6 +1312,7 @@ void __unused__ *op_data;
fprintf(stderr,"h5dump error: no <names> after option %s\n",
argv[opts[i]]);
usage();
+ free(opts);
exit(1);
}
}
@@ -878,6 +1320,7 @@ void __unused__ *op_data;
if (argc - opts[nopts-1] == 1) {
fprintf(stderr,"h5dump error: no <file>\n");
usage();
+ free(opts);
exit(1);
}
if (argc - opts[nopts-1] == 2) {
@@ -885,6 +1328,7 @@ void __unused__ *op_data;
strcmp(argv[opts[i]], "-header") ) {
fprintf (stderr, "h5dump error: no <file> or no <names> after option %s\n", argv[opts[i]]);
usage();
+ free(opts);
exit(1);
}
}
@@ -896,9 +1340,32 @@ void __unused__ *op_data;
if ((fid = H5Fopen (fname, H5F_ACC_RDONLY, plist)) < 0) {
fprintf (stderr, "h5dump error: unable to open file %s \n", fname);
+ free(opts);
exit(1);
}
+ /* allocate and initialize internal data structure */
+ init_table();
+
+ /* find all shared objects */
+ H5Giterate (fid, "/", NULL, find_shared_objs, NULL);
+ strcpy(prefix, "");
+
+ /* assign names to unamed shared data type */
+ for ( i = 0; i < type_table.nobjs; i++)
+ if (type_table.objs[i].displayed == 0) unamedtype = 1;
+
+/*
+ #ifdef H5DUMP_DEBUG
+ dump_tables();
+ #endif
+*/
+
+ if (status) {
+ printf("internal error! \n");
+ goto done;
+ }
+
begin_obj("HDF5", fname);
if (display_bb) dump_bb();
@@ -908,7 +1375,7 @@ void __unused__ *op_data;
if ((gid = H5Gopen (fid, "/")) < 0 ) {
fprintf(stdout, "h5dump error: unable to open root group\n");
status = 1;
- } else
+ } else
dump_group(gid, "/");
if (H5Gclose (gid) < 0) {
@@ -940,12 +1407,32 @@ void __unused__ *op_data;
end_obj();
status = 1;
} else {
- dump_dataset(dsetid, argv[curr_arg]);
+ H5Gget_objinfo(dsetid, ".", TRUE, &statbuf);
+ if (statbuf.nlink > 1) {
+ index = search_obj (dset_table, statbuf.objno);
+ if (index >= 0) {
+ 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);
+ indentation (indent);
+ end_obj();
+ } else {
+ strcpy(dset_table.objs[index].objname, argv[curr_arg]);
+ dset_table.objs[index].displayed = 1;
+ dump_dataset(dsetid, argv[curr_arg]);
+ }
+ } else status = 1;
+ } else
+ dump_dataset(dsetid, argv[curr_arg]);
if (H5Dclose(dsetid)<1) status = 1;
}
}
+
+
} else if (!strcmp(argv[opts[i]],"-g")) {
for (curr_arg = opts[i]+1;
@@ -959,6 +1446,8 @@ void __unused__ *op_data;
end_obj();
status = 1;
} else {
+ H5Gget_objinfo(gid, ".", TRUE, &statbuf);
+ strcpy(prefix, argv[curr_arg]);
dump_group(gid, argv[curr_arg]);
if (H5Gclose (gid) < 0) status = 1;
}
@@ -984,7 +1473,7 @@ void __unused__ *op_data;
begin_obj(SOFTLINK, argv[curr_arg]);
indentation (col);
if (H5Gget_linkval (fid, argv[curr_arg], statbuf.linklen, buf)>=0)
- printf ("linktarget \"%s\"\n", buf);
+ printf ("LINKTARGET \"%s\"\n", buf);
else {
fprintf (stdout, "h5dump error: unable to get link value\n");
status = 1;
@@ -1002,16 +1491,65 @@ void __unused__ *op_data;
}
+ } else if (!strcmp(argv[opts[i]],"-t")) {
+
+
+ for (curr_arg = opts[i]+1;
+ curr_arg < ((i+1)==nopts?(argc-1):opts[i+1]);
+ curr_arg++) {
+
+ if ((typeid=H5Topen (fid, argv[curr_arg])) < 0) {
+
+ /* 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 */
+ 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],
+ type_table.objs[index].objno[1]);
+ if (!strncmp(name, argv[curr_arg], strlen(argv[curr_arg])) ||
+ !strncmp(name1, argv[curr_arg], strlen(argv[curr_arg]))) {
+ break;
+ }
+ }
+ index++;
+ }
+ if (index == type_table.nobjs) { /* unknown type */
+ begin_obj (DATATYPE, argv[curr_arg]);
+ indentation (col);
+ fprintf (stdout, "h5dump error: unable to open %s\n",
+ argv[curr_arg]);
+ end_obj();
+ status = 1;
+ } else {
+ dsetid = H5Dopen (fid, type_table.objs[index].objname) ;
+ typeid = H5Dget_type (dsetid);
+ dump_named_datatype (typeid, argv[curr_arg]);
+ H5Tclose(typeid);
+ H5Dclose(dsetid);
+ }
+
+ } else {
+ dump_named_datatype (typeid, argv[curr_arg]);
+ if (H5Tclose(typeid) < 0) status = 1;
+ }
+ }
}
}
end_obj();
- free(opts);
+done:
- /* Restore error reporting */
H5Eset_auto (func, edata);
- if (H5Fclose (fid) < 0) exit(1);
+ free(opts);
+ if (H5Fclose (fid) < 0) status = 1;
+
+ free (group_table.objs);
+ free (dset_table.objs);
+ free (type_table.objs);
+ free (prefix);
return status;