diff options
Diffstat (limited to 'tools/h5toh4')
| -rw-r--r-- | tools/h5toh4/Dependencies | 0 | ||||
| -rw-r--r-- | tools/h5toh4/Makefile.in | 60 | ||||
| -rw-r--r-- | tools/h5toh4/h5toh4.c | 2592 | ||||
| -rw-r--r-- | tools/h5toh4/h5toh4.h | 43 | ||||
| -rwxr-xr-x | tools/h5toh4/testh5toh4.sh | 299 |
5 files changed, 2994 insertions, 0 deletions
diff --git a/tools/h5toh4/Dependencies b/tools/h5toh4/Dependencies new file mode 100644 index 0000000..e69de29 --- /dev/null +++ b/tools/h5toh4/Dependencies diff --git a/tools/h5toh4/Makefile.in b/tools/h5toh4/Makefile.in new file mode 100644 index 0000000..4c529d6 --- /dev/null +++ b/tools/h5toh4/Makefile.in @@ -0,0 +1,60 @@ +## HDF5 Library Makefile(.in) +## +## Copyright (C) 2001 National Center for Supercomputing Applications. +## All rights reserved. +## +## +top_srcdir=@top_srcdir@ +top_builddir=../.. +srcdir=@srcdir@ +SUBDIRS= +@COMMENCE@ + +## Add include directory to the C preprocessor flags, add -lh5tools and +## -lhdf5 to the list of libraries. +## +CPPFLAGS=-I. -I$(srcdir) -I$(top_builddir)/src -I$(top_srcdir)/src \ + -I$(top_srcdir)/tools/lib @CPPFLAGS@ + +## Test programs and scripts. +## +TEST_PROGS= +TEST_SCRIPTS=@TESTH5TOH4@ + +## These are our main targets: library and tools. +## +LIBTOOLS=../lib/libh5tools.la +LIBHDF5=$(top_builddir)/src/libhdf5.la +PUB_PROGS=@H5TOH4@ +PROGS=$(PUB_PROGS) $(TEST_PROGS) + +## Source and object files for the library; do not install +## +LIB_SRC= +LIB_OBJ=$(LIB_SRC:.c=.lo) +PUB_LIB= + +## Source and object files for programs... +## +PROG_SRC=h5toh4.c +PROG_OBJ=$(PROG_SRC:.c=.lo) +PRIVATE_HDR=h5toh4.h + +## Source and object files for the tests +## +TEST_SRC= +TEST_OBJ=$(TEST_SRC:.c=.lo) + +## Programs have to be built before they can be tested! +## +check test _test: $(PROGS) + +## How to build the programs... They all depend on the hdf5 library and +## the tools library compiled in this directory. +## +$(PROGS): $(LIBTOOLS) $(LIBHDF5) + +h5toh4: h5toh4.lo + @$(LT_LINK_EXE) $(CFLAGS) -o $@ h5toh4.lo $(LIBTOOLS) $(LIBHDF5) $(LDFLAGS) $(LIBS) + +@CONCLUDE@ diff --git a/tools/h5toh4/h5toh4.c b/tools/h5toh4/h5toh4.c new file mode 100644 index 0000000..30b381c --- /dev/null +++ b/tools/h5toh4/h5toh4.c @@ -0,0 +1,2592 @@ +/****************************************************************************** + + Description: This file contains routines to translate H5 files to H4 files. + + Author: Paul Harten for the University of Illinois, NCSA + + Creation Date: 04 October 1998 + + History: + + +*****************************************************************************/ + +#include "h5toh4.h" +#include <errno.h> +#include <string.h> +#include <fcntl.h> +#include <h5tools.h> + + +#ifdef WIN32 +typedef unsigned int mode_t; +#endif + +#ifndef S_ISDIR +#define S_ISDIR(mode) (((mode)&0xF000) == S_IFDIR) +#endif +extern void PrintOptions_h5toh4(void); +extern char *BuildFilename(char *h5_filename, char *h4_extension); +extern int test_file(char *filename, int oflag, mode_t mode); +extern int test_dir(char *); +extern int h5toh4(char *, char *); + +extern herr_t convert_group(hid_t, char *, op_data_t *); +extern herr_t convert_dataset(hid_t, char *, op_data_t *); +extern herr_t convert_all(hid_t, char *, op_data_t *); +extern herr_t convert_attr(hid_t, char *, op_data_t *); +extern herr_t convert_shared_dataset(hid_t, int, op_data_t *); +extern herr_t convert_shared_group(hid_t, int, op_data_t *); +extern herr_t convert_dataset_string(hid_t, char *, op_data_t *); +extern int32 h5type_to_h4type(hid_t); +extern hid_t h4type_to_memtype(int32); + +extern void free_table(table_t **temp); +extern void dump_tables(char* name, table_t* table); +extern int get_table_idx(table_t*, unsigned long *); +extern int get_tableflag(table_t*, int); +extern int set_tableflag(table_t*, int); +extern char* get_objectname(table_t*, int); + +typedef herr_t (*H5G_operator_t)(hid_t, const char*, void*); + +static int prefix_len = 1024; +static char *prefix; +static table_t *group_table, *dset_table, *type_table; + +static herr_t h5atomic_type_to_h4type(const hid_t h5type, hid_t* h5memtype, size_t* h5memsize, int32* h4type); + + +/***************************************************************************** + + Routine: main() + + Description: This routine check out arguments sent, makes sure everything is + ok, chooses between the acceptable argument formats and then + calls h5toh4(). + + Input: arg - the number of arguments from the call + 1; + argv - a pointer to an array of strings which are the arguments. + This includes the routine name as the first string. + + Output: function return - error status + +*****************************************************************************/ + +int +main(int argc, char *argv[]) +{ + char *h5_filename=NULL; + char *h4_filename=NULL; + char *h4_extension; + int status = 0; + int status2 = 0; + + + argc--; + argv++; + + if (argc == 0) { + fprintf(stderr,"\nError: Invalid Arguments\n"); + PrintOptions_h5toh4(); + return -1; + } + + /* take care -h (help) option first */ + { int i; + for (i=0; i < argc; i++) + if ( HDstrcmp(argv[i],"-h") == 0 ) { + PrintOptions_h5toh4(); + return 0; + } + } + + /*fargv = argv; */ + + /* Change -m with 1 filename to as 1 filename case */ + if (argc == 2 && HDstrcmp(argv[0],"-m") == 0) { + argv++; + argc--; + } + + switch(argc) { + + case 0: + + PrintOptions_h5toh4(); + break; + + case 1: /* h5toh4 file1 */ + h5_filename = argv[0]; + + if (HDstrcmp(h5_filename,"-m") == 0) { + fprintf(stderr,"\nError: Invalid Arguments\n"); + PrintOptions_h5toh4(); + status = -1; + break; + } +#ifndef WIN32 + if (test_file(h5_filename,O_EXCL,(mode_t)(S_IRUSR|S_IRGRP|S_IROTH)) != 0 ) { /* 292 Decimal - 0444 Octal, a+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } + if (test_dir(h5_filename) != 0 ) { + fprintf(stderr,"%s: Is a directory\n",h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } +#endif + h4_extension = HDstrdup("hdf"); + + h4_filename = BuildFilename(h5_filename,h4_extension); + if (h4_filename == NULL) { + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } +#ifndef WIN32 + if (test_file(h4_filename,O_CREAT|O_EXCL,(mode_t)(S_IRUSR|S_IWUSR|S_IRGRP|S_IWGRP|S_IROTH)) != 0) { /* 436 Decimal - 0664 Octal, ug+rw,o+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } +#endif + status = h5toh4(h5_filename, h4_filename); + if ( status != 0 ) { + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + } + if (h4_filename != NULL) { + HDfree(h4_filename); + } + + break; + + case 2: /* h5toh4 file_in file_out */ + h5_filename = argv[0]; + h4_filename = argv[1]; + + if (HDstrcmp(h4_filename,"-m") == 0) { + fprintf(stderr,"\nError: Invalid Arguments\n"); + PrintOptions_h5toh4(); + status = -1; + break; + } +#ifndef WIN32 + if (test_file(h5_filename,O_EXCL,292) != 0 ) { /* 292 Decimal - 0444 Octal, a+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } + if (test_dir(h5_filename) != 0 ) { + fprintf(stderr,"%s: Is a directory\n",h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } + if (test_file(h4_filename,O_CREAT|O_RDWR,436) != 0) { /* 436 Decimal - 0664 Octal, ug+rw,o+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } + if (test_dir(h4_filename) != 0 ) { + fprintf(stderr,"%s: Is a directory\n",h4_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status = -1; + break; + } +#endif + status = h5toh4(h5_filename, h4_filename); + if ( status != 0 ) { + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + } + break; + + default: /* h5toh4 -m file1 file2 file3 ... */ + if (HDstrcmp(argv[0],"-m") != 0) { + fprintf(stderr,"\nError: Invalid Arguments\n"); + PrintOptions_h5toh4(); + status = -1; + break; + } + argv++; + argc--; + + while ( argc > 0 ) { + h5_filename = argv[0]; + argv++; + argc--; + + if (HDstrcmp(h5_filename,"-m") == 0) { + fprintf(stderr,"\nError: Invalid Arguments\n"); + PrintOptions_h5toh4(); + status2 = -1; + break; + } +#ifndef WIN32 + if (test_file(h5_filename,O_EXCL,292) != 0 ) { /* 292 Decimal - 0444 Octal, a+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status2 = -1; + continue; + } + if (test_dir(h5_filename) != 0 ) { + fprintf(stderr,"%s: Is a directory\n",h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status2 = -1; + continue; + } +#endif + h4_extension = HDstrdup("hdf"); + + h4_filename = BuildFilename(h5_filename,h4_extension); + if (h4_filename == NULL) { + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status2 = -1; + continue; + } +#ifndef WIN32 + if (test_file(h4_filename,O_CREAT|O_EXCL,436) != 0) { /* 436 Decimal - 0664 Octal, ug+rw,o+r */ + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status2 = -1; + continue; + } +#endif + status = h5toh4(h5_filename, h4_filename); + if ( status != 0 ) { + fprintf(stderr,"Error: Problem with %s\n",h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "main", __FILE__, __LINE__); + status2 = status; + } + if (h4_filename != NULL) { + HDfree(h4_filename); + } + + } + status = status2; + break; + } + + return status; + +} + + +/***************************************************************************** + + Routine: PrintOptions_h5toh4() + + Description: This routine prints the acceptable argument formats out to stderr. + + Input: None + + Output: output to stderr + +*****************************************************************************/ + +void PrintOptions_h5toh4(void) +{ + fprintf(stderr,"\nh5toh4 -h (gives this print-out)\n"); + fprintf(stderr,"h5toh4 file.h5 file.hdf\n"); + fprintf(stderr,"h5toh4 file.h5\n"); + fprintf(stderr,"h5toh4 -m file1.h5 file2.h5 ...\n\n"); +} + + +/***************************************************************************** + + Routine: h5toh4() + + Description: This routine translates information from a HDF5 file into a + HDF4 file. + + Input: h5_filename - filename of HDF5 file, + h4_filename - filename of HDF4 file + + Output: function return - error status + +*****************************************************************************/ + +int h5toh4(char *h5_filename, char *h4_filename) +{ + + hid_t fid, gid; + hid_t plist=H5P_DEFAULT; + int status = 0; + int32 hfile_id; + int32 sd_id; + op_data_t op_data; + void *edata; + hid_t (*func)(void*); + + + find_objs_t *info = malloc(sizeof(find_objs_t)); + + + + /* open hdf5 file */ + if ((fid = H5Fopen (h5_filename, H5F_ACC_RDONLY, plist)) <= 0) { + fprintf(stderr,"Error: Unable to open file %s\n",h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + return (fid); + } + + /* open root group */ + if ((gid = H5Gopen (fid, "/")) <= 0 ) { + fprintf(stderr,"Error: Unable to open root group\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + return (gid); + } else { + + /* open hdf4 file */ + if ((hfile_id = Hopen(h4_filename, DFACC_CREATE, 0)) <= 0) { + fprintf(stderr,"Error: Unable to open file %s\n", h4_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + return (hfile_id); + } else { + + if ((sd_id = SDstart(h4_filename, DFACC_RDWR)) <= 0) { + fprintf(stderr,"Error: Unable to open file %s, using SDstart\n", h4_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + return (sd_id); + } else { + + /* Initialize Vgroup interface */ + if ((status = Vstart(hfile_id)) != SUCCEED) { + fprintf(stderr,"Error: Unable to initialize Vgroup interface\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + + /* allocate and initialize internal data structure */ + init_table(&group_table); + init_table(&type_table); + init_table(&dset_table); + init_prefix(&prefix, prefix_len); + + /* init the find_objs_t*/ + info->threshold = 0; + info->prefix_len = prefix_len; + info->prefix = prefix; + info->group_table = group_table; + info->type_table = type_table; + info->dset_table = dset_table; + info->status = status; + + + + /* Disable error reporting */ + H5Eget_auto (&func, &edata); + H5Eset_auto (NULL, NULL); + + /* find all objects one time */ + if ((status = H5Giterate(fid, "/", NULL, (H5G_operator_t)find_objs, (void*)info)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the groups\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + } + + /* Enable error reporting */ + H5Eset_auto (func, edata); + +#ifdef H5DUMP_DEBUG + dump_tables(); +#endif + + if (status != SUCCEED) { + fprintf(stderr,"Error: internal error! \n"); + goto done; + } + + op_data.hfile_id = hfile_id; + op_data.sd_id = sd_id; + op_data.vgroup_id = 0; + + /* start at root group */ + if (( status = convert_group (gid, HDstrdup("/"), &op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_group did not work for %s\n","/"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + + /* Terminate access to Vgroup interface */ + if ((status = Vend(hfile_id)) != SUCCEED) { + fprintf(stderr,"Error: Unable to terminate Vgroup interface\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + } + + /* close hdf4 SDS */ + if ((status = SDend(sd_id)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close file %s, from SDstart\n", h4_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + + } + + /* close hdf4 file */ + if ((status = Hclose(hfile_id)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close file %s\n", h4_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + + } + + if ((status = H5Gclose (gid)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close root group\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + + if ((status = H5Fclose (fid)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close file %s\n", h5_filename); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5toh4", __FILE__, __LINE__); + status = FAIL; + } + +done: + free_table(&group_table); + free_table(&dset_table); + free_table(&type_table); + + return status; + +} + +/*------------------------------------------------------------------------- + * Function: convert_group + * + * Purpose: Dump everything within the specified group + * + * Return: status + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_group (hid_t gid, char *name, op_data_t *op_data) { +H5G_stat_t statbuf; + + int32 hfile_id; + int32 vgroup_id; + int32 obj_idx; + int32 status; + int idx, flag; + + hfile_id = op_data->hfile_id; + + if ((vgroup_id = Vattach(hfile_id, -1, "w")) <= 0 ) { + fprintf(stderr,"Error: Unable to create new vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(vgroup_id); + } + + if ((status = Vsetname(vgroup_id, name)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to set name on new vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + if ((status = Vsetclass(vgroup_id, "HDF5")) != SUCCEED ) { + fprintf(stderr,"Error: Unable to set class on new vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + if (op_data->vgroup_id == 0) { + obj_idx = -1; /* Indexes assigned below start from 0 */ + } else { + if ((obj_idx = Vinsert(op_data->vgroup_id, vgroup_id)) < 0) { + fprintf(stderr,"Error: Index %d of the new vgroup is illegal\n",(int)obj_idx); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(obj_idx); + } + } + op_data->vgroup_id = vgroup_id; + op_data->sds_id = 0; + op_data->vdata_id = 0; + op_data->obj_idx = obj_idx; + + + /* hard link */ + if ((status = H5Gget_objinfo(gid, ".", TRUE, &statbuf)) != SUCCEED ) { + fprintf(stderr,"Error: H5Gget_objinfo() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return (status); + } + + if (HDstrcmp(name,"/") == 0) { /* this is the root group, just iterate */ + + if ((status = H5Aiterate(gid, NULL, (H5A_operator_t)convert_attr, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the attributes\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + if ((status = H5Giterate(gid, ".", NULL, (H5G_operator_t)convert_all, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the groups\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + } else { + + if ((idx = get_table_idx(group_table, statbuf.objno)) < 0 ) { + + fprintf(stderr,"Error: object not found, %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + status = FAIL; + + } else if((flag = get_tableflag(group_table,idx)) < 0 ) { + + fprintf(stderr,"Error: get_tableflag() should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + status = FAIL; + + } else if(flag == TRUE ) { /* this has already been converted, don't convert the attributes again */ + + if ((status = H5Giterate(gid, ".", NULL, (H5G_operator_t)convert_all, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the groups\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + } else { /* flag == FALSE */ + + /* this is now being converted */ + if ((status = set_tableflag(group_table,idx)) < 0 ) { + fprintf(stderr,"Error: set_tableflag should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + if ((status = H5Aiterate(gid, NULL, (H5A_operator_t)convert_attr, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the attributes\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + if ((status = H5Giterate(gid, ".", NULL, (H5G_operator_t)convert_all, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to iterate over all of the groups\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + return(status); + } + + } + + } + + + if ((status = Vdetach(vgroup_id)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to detach the new Vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_group", __FILE__, __LINE__); + status = FAIL; + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_dataset + * + * Purpose: Convert the specified data set + * + * Return: status + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_dataset (hid_t did, char *name, op_data_t *op_data) { +hid_t type, space, class, mem_type, type2; +/* H5G_stat_t statbuf; */ +size_t typesize; +int i, idx; +int32 dim_sizes[32], start[32], edges[32]; +int ndims; +int ndimf; +hsize_t dims[32], maxdims[32]; +hsize_t dimf[4]; +int permf[4]; +int32 hfile_id; +int32 sd_id; +int32 sds_id; +int32 vdata_id; +int32 vgroup_id; +int32 n_values; +int32 status; +int32 h4_type; +int32 recsize; +int32 n_records, num_of_recs, record_pos; +intn nmembers; +char *buffer=NULL; /* read/write buffer*/ +char fieldname_list[4096] = "\0"; +char *fieldname=NULL; +hid_t fieldtype; +int32 order; +off_t offset; +off_t offset_array[512]; +hid_t h4type_array[512], memtype_array[512]; +int32 order_array[512]; + + + /* hard link */ +/* + if ((status = H5Gget_objinfo(did, ".", TRUE, &statbuf)) != SUCCEED ) { + fprintf(stderr,"Error: H5Gget_objinfo() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return (status); + } +*/ + + if ((type = H5Dget_type(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_type() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((space = H5Dget_space(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_space() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((n_values = H5Sget_simple_extent_npoints(space)) <= 0) { + fprintf(stderr, "Error: H5Sget_simple_extent_npoints() returned inappropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((ndims = H5Sget_simple_extent_dims(space,dims,maxdims)) < 0 ) { + fprintf(stderr, "Error: Problems getting ndims, dims, and maxdims of dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = ndims; + return status; + } + + if ((class = H5Tget_class(type)) < 0 ) { + fprintf(stderr,"Error: problem with getting class\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = class; + return status; + } + + switch (class) { + case H5T_INTEGER: + case H5T_FLOAT: + sd_id = op_data->sd_id; +#define NEWWAY +#ifdef NEWWAY + if (FAIL==h5atomic_type_to_h4type(type, &mem_type, &typesize, &h4_type)){ + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +#else + if ((h4_type = h5type_to_h4type(type)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((mem_type = h4type_to_memtype(h4_type)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h4 type to mem type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + if ((typesize = H5Tget_size(mem_type)) <= 0) { + fprintf(stderr, "Error: H5Tget_size() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +#endif + if ((buffer = HDmalloc(n_values*typesize)) == NULL) { + fprintf(stderr, "Error: Problems with HDmalloc of memory space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((status = H5Dread(did, mem_type, space, space, H5P_DEFAULT, buffer)) != SUCCEED) { + fprintf(stderr, "Error: Problems with H5Dread\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + for (i=0;i<ndims;i++) { + if (maxdims[i] == H5S_UNLIMITED) { + if ( i == 0 ) { + dim_sizes[0] = 0; /* this is how HDF4 communicates unlimited dimension */ + } else { + dim_sizes[i] = (int32)dims[i]; + } + } else { + dim_sizes[i] = (int32)maxdims[i]; + } + start[i] = 0; + edges[i] = (int32)dims[i]; + } + if ((sds_id = SDcreate(sd_id, name, h4_type, ndims, dim_sizes)) <= 0 ) { + fprintf(stderr, "Error: Unable to create SDS %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + op_data->sds_id = sds_id; + if ((status = SDwritedata(sds_id, start, NULL, edges, buffer)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to write SDS %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + if ((status = H5Aiterate(did, NULL, (H5A_operator_t)convert_attr, op_data)) < 0 ) { + fprintf(stderr,"Error: iterate over attributes\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + if ((status = SDendaccess(sds_id)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to end access to SDS %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + break; + case H5T_TIME: + fprintf(stderr,"Warning: H5T_TIME not yet implemented.\n"); + break; + case H5T_STRING: + if (ndims==1) { + convert_dataset_string(did, name, op_data); + } else { + fprintf(stderr,"Warning: HDF5 datasets of H5T_STRING type with ndims > 1 are not converted.\n"); + } + break; + case H5T_BITFIELD: + fprintf(stderr,"Warning: H5T_BITFIELD not yet implemented.\n"); + break; + case H5T_OPAQUE: + fprintf(stderr,"Warning: H5T_OPAQUE not yet implemented.\n"); + break; + case H5T_ARRAY: + fprintf(stderr,"Warning: H5T_OPAQUE not yet implemented.\n"); + break; + case H5T_COMPOUND: + if (ndims==1) { + if ((nmembers = H5Tget_nmembers(type)) <= 0 ) { + fprintf(stderr, "Error: Unable to get information about compound datatype %d\n",nmembers); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + + offset = 0; + for (idx=0;idx<nmembers;idx++) { + if ((fieldtype = H5Tget_member_type(type, idx)) < 0 ) { + fprintf(stderr,"Error: H5 datasets of H5T_COMPOUND type with fieldtype %d, idx %d.\n",fieldtype,idx); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + /* Special case for array fields */ + if(H5Tget_class(fieldtype)==H5T_ARRAY) { + hid_t arr_base_type; + + /* Get the number of dimensions */ + if ((ndimf = H5Tget_array_ndims(fieldtype)) < 0 || ndimf > H5S_MAX_RANK ) { + fprintf(stderr, "Error: rank of members of compound type should not be %d\n",ndimf); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + /* Get the dimensions and dimension permutations */ + if (H5Tget_array_dims(fieldtype,dimf,permf) < 0) { + fprintf(stderr, "Error: cannot retrieve dimensions for array\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + /* Get the 'real' field type */ + if ((arr_base_type=H5Tget_super(fieldtype)) < 0) { + fprintf(stderr, "Error: cannot retrieve base type for array\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + if (H5Tclose(fieldtype) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + + /* Assign the array's base type as the field type for future use */ + fieldtype=arr_base_type; + } /* end if */ + else + ndimf=0; +#ifdef NEWWAY + if (FAIL==h5atomic_type_to_h4type(fieldtype, &mem_type, &typesize, &h4_type)){ + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +#else + if ((h4_type = h5type_to_h4type(fieldtype)) < 0 ) { + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((mem_type = h4type_to_memtype(h4_type)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h4 type to mem type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + if ((typesize = H5Tget_size(mem_type)) <= 0) { + fprintf(stderr, "Error: H5Tget_size() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +#endif + order = 1; + if (ndimf > 0) { + order *= dimf[permf[0]]; + if (ndimf > 1) { + order *= dimf[permf[1]]; + if (ndimf > 2) { + order *= dimf[permf[2]]; + if (ndimf > 3) { + order *= dimf[permf[3]]; + } + } + } + } + h4type_array[idx] = h4_type; /* used for VSfdefine */ + memtype_array[idx] = mem_type; /* used during the build of the memory compound type */ + offset_array[idx] = offset; /* used during the build of the memory compound type */ + order_array[idx] = order; /* used for VSfdefine */ + + offset += order * typesize; /* calculating packed size of memory compound type */ + } + hfile_id = op_data->hfile_id; + vgroup_id = op_data->vgroup_id; + if ((vdata_id = VSattach(hfile_id, -1, "w")) <= 0 ) { + fprintf(stderr, "Error: Unable to create vdata %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +/* + * This step is done during the convert_shared_dataset() call instead of here. + * + if ((idx = Vinsert(vgroup_id, vdata_id)) < 0 ) { + fprintf(stderr, "Error: Unable to insert vdata %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +*/ + op_data->vdata_id = vdata_id; + if ((status = VSsetname(vdata_id, name)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set vdata name %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((status = VSsetclass(vdata_id, "HDF5")) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set class on vdata %s\n", name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((type2 = H5Tcreate(H5T_COMPOUND, (size_t)offset)) <= 0 ) { + fprintf(stderr, "Error: unable to execute H5Tcreate()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + for (idx=0;idx<nmembers;idx++) { + /* Special case for array fields */ + if(H5Tget_member_class(type,idx)==H5T_ARRAY) { + hid_t arr_type; + + if ((arr_type = H5Tget_member_type(type, idx)) < 0 ) { + fprintf(stderr,"Error: H5 datasets of H5T_COMPOUND type with fieldtype %d, idx %d.\n",arr_type,idx); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + /* Get the number of dimensions */ + if ((ndimf = H5Tget_array_ndims(arr_type)) < 0 || ndimf > H5S_MAX_RANK ) { + fprintf(stderr, "Error: rank of members of compound type should not be %d\n",ndimf); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + /* Get the dimensions and dimension permutations */ + if (H5Tget_array_dims(arr_type,dimf,permf) < 0) { + fprintf(stderr, "Error: cannot retrieve dimensions for array\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + if (H5Tclose(arr_type) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + } /* end if */ + else + ndimf=0; + + if ((fieldname = H5Tget_member_name(type, idx)) == NULL ) { + fprintf(stderr, "Error: Unable to get fieldname for compound type %d, idx %d\n", type, idx); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((offset = H5Tget_offset(memtype_array[idx])) < 0 || offset >= 128 ) { + fprintf(stderr, "Error: H5Tget_offset() is returning a bad value %d\n",(int)offset); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + if (ndimf == 0 ) { + if ((status = H5Tinsert(type2,fieldname,offset_array[idx]+offset,memtype_array[idx])) != SUCCEED ) { + fprintf(stderr, "Error: Problems inserting field into compound datatype\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + } else { + hid_t arr_type; + + /* Create the array datatype */ + if ((arr_type = H5Tarray_create(memtype_array[idx], ndimf, dimf, permf)) < 0 ) { + fprintf(stderr,"Error: cannot create array datatype\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + + if ((status = H5Tinsert(type2,fieldname,offset_array[idx]+offset,arr_type)) != SUCCEED ) { + fprintf(stderr, "Error: Problems inserting field into compound datatype\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + + if (H5Tclose(arr_type) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + return FAIL; + } + } + if ((status = VSfdefine(vdata_id, fieldname, h4type_array[idx], order_array[idx])) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set field %d\n", idx); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + strcat(fieldname_list,fieldname); + if (idx<nmembers-1) { + strcat(fieldname_list,", "); + } + + if (fieldname != NULL) { + HDfree(fieldname); + } + } + if ((status = VSsetfields(vdata_id, fieldname_list)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set fieldname list %s\n", name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((status = VSsetinterlace(vdata_id, FULL_INTERLACE)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set FULL_INTERLACE mode, status %d\n", (int)status); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((recsize = H5Tget_size(type2)) <= 0 ) { + fprintf(stderr, "Error: Unable to get record size %d\n", (int)recsize); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } +/* + Since the space is rank 1, n_records does not depend on maxdims. +*/ + n_records = n_values; + if ((buffer = HDmalloc(n_records*recsize)) == NULL) { + fprintf(stderr, "Error: Problems with HDmalloc of memory space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((status = H5Dread(did, type2, space, space, H5P_DEFAULT, buffer)) != SUCCEED) { + fprintf(stderr, "Error: Problems with H5Dread\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((status = H5Tclose(type2)) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } + + if ((record_pos = VSseek(vdata_id, 0)) != 0 ) { + fprintf(stderr, "Error: Could not seek the beginning of the Vdata, %d\n", (int)record_pos); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = record_pos; + break; + } + if ((num_of_recs = VSwrite(vdata_id, (void *)buffer, n_records, FULL_INTERLACE)) != n_records ) { + fprintf(stderr, "Error: Only able to write %d of %d records\n", (int)num_of_recs, (int)n_records); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = num_of_recs; + break; + } + + /* there are only vdata attributes, no field attributes */ + if ((status = H5Aiterate(did, NULL, (H5A_operator_t)convert_attr, op_data)) < 0 ) { + fprintf(stderr,"Error: iterate over attributes\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + if ((status = VSdetach(vdata_id)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to detach to vdata %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + break; + } + + } else { + + fprintf(stderr,"Warning: HDF5 datasets of H5T_COMPOUND type with ndims > 1 are not converted.\n"); + + } + break; + default: + fprintf(stderr,"Error: %d class not found\n",class); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + + if ((status = H5Tclose(type)) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + + if ((status = H5Sclose(space)) < 0 ) { + fprintf(stderr,"Error: closing space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + } + + if (buffer != NULL) { + HDfree(buffer); + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_attr + * + * Purpose: dump the attribute + * + * Return: Success: SUCCEED + * + * Failure: FAIL + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ + +herr_t +convert_attr (hid_t attr, char *attr_name, op_data_t *op_data) +{ +hid_t attr_id, type, space, mem_type, class; +size_t typesize; +char *attr_values=NULL; +int32 status; +int32 h4_type; +int32 sds_id; +int32 vdata_id; +int32 vgroup_id; +int32 n_values; +int32 order; +hid_t fxdlenstr; +size_t lenstr; +H5T_cset_t cset; +H5T_str_t strpad; + + sds_id = op_data->sds_id; + vdata_id = op_data->vdata_id; + vgroup_id = op_data->vgroup_id; + + if ((attr_id = H5Aopen_name (attr, attr_name))>= 0) { + + if ((type = H5Aget_type(attr_id)) <= 0) { + fprintf(stderr, "Error: H5Aget_type() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((class = H5Tget_class(type)) < 0 ) { + fprintf(stderr,"Error: problem with getting class\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = class; + return status; + } + + switch (class) { + case H5T_INTEGER: + case H5T_FLOAT: + + + if ((space = H5Aget_space(attr_id)) <= 0) { + fprintf(stderr, "Error: H5Dget_space() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((n_values = H5Sget_simple_extent_npoints(space)) <= 0) { + fprintf(stderr, "Error: H5sget_simple_extent_npoints() didn't return correct value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + +#ifdef NEWWAY + if (FAIL==h5atomic_type_to_h4type(type, &mem_type, &typesize, &h4_type)){ + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + break; + } +#else + if ((h4_type = h5type_to_h4type(type)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h5 type to h4 type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((mem_type = h4type_to_memtype(h4_type)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h4 type to mem type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((typesize = H5Tget_size(mem_type)) <= 0) { + fprintf(stderr, "Error: H5Tget_size() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } +#endif + + if ((attr_values = HDmalloc(n_values*typesize)) == NULL) { + fprintf(stderr, "Error: Problems with HDmalloc of memory space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((status = H5Aread(attr_id, mem_type, attr_values)) != SUCCEED) { + fprintf(stderr, "Error: Problems with H5Aread\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if (sds_id != 0) { + if ((status = SDsetattr(sds_id, attr_name, h4_type, n_values, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } else if (vdata_id != 0) { + if ((status = VSsetattr(vdata_id, -1, attr_name, h4_type, n_values, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } else { + if ((status = Vsetattr(vgroup_id, attr_name, h4_type, n_values, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } + + if ((status = H5Sclose(space)) != SUCCEED ) { + fprintf(stderr, "Error: Problems closing H5Sclose\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + if ((status = H5Aclose(attr_id)) != SUCCEED ) { + fprintf(stderr, "Error: Problems closing H5Aclose\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if (attr_values != NULL) { + HDfree(attr_values); + } + + status = SUCCEED; + break; + + case H5T_TIME: + fprintf(stderr,"Warning: H5T_TIME attribute not yet implemented.\n"); + break; + case H5T_STRING: + + fxdlenstr = type; + h4_type = DFNT_CHAR; + + if ((space = H5Aget_space(attr_id)) <= 0) { + fprintf(stderr, "Error: H5Dget_space() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((n_values = H5Sget_simple_extent_npoints(space)) <= 0) { + fprintf(stderr, "Error: H5sget_simple_extent_npoints() didn't return correct value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((mem_type = H5Tcopy(H5T_C_S1)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h4 type to mem type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = mem_type; + return status; + } + if ((lenstr = H5Tget_size(fxdlenstr)) <= 0 ) { + fprintf(stderr, "Error: size of fixed length string type should not be %d\n",(int)lenstr); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return lenstr; + } + if ((status = H5Tset_size(mem_type,lenstr)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_size()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return status; + } + + strpad = H5Tget_strpad(fxdlenstr); + if ((strpad != H5T_STR_NULLTERM) && (strpad != H5T_STR_NULLPAD) && (strpad != H5T_STR_SPACEPAD)) { + fprintf(stderr, "Error: Invalid string padding value, %d\n",(int)strpad); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return strpad; + } + if ((status = H5Tset_strpad(mem_type, H5T_STR_SPACEPAD)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_strpad()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return status; + } + + if ((cset = H5Tget_cset(fxdlenstr)) != H5T_CSET_ASCII ) { + fprintf(stderr, "Error: cset value != %d\n",(int)H5T_CSET_ASCII); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return cset; + } + if ((status = H5Tset_cset(mem_type,cset)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_cset()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + return status; + } + + order = n_values * lenstr; + if ((attr_values = HDmalloc(order)) == NULL) { + fprintf(stderr, "Error: Problems with HDmalloc of memory space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((status = H5Aread(attr_id, mem_type, attr_values)) != SUCCEED) { + fprintf(stderr, "Error: Problems with H5Aread\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if (sds_id != 0) { + if ((status = SDsetattr(sds_id, attr_name, h4_type, order, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } else if (vdata_id != 0) { + if ((status = VSsetattr(vdata_id, -1, attr_name, h4_type, order, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } else { + if ((status = Vsetattr(vgroup_id, attr_name, h4_type, order, attr_values)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set %s attribute.\n",attr_name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + } + + if ((status = H5Sclose(space)) != SUCCEED ) { + fprintf(stderr, "Error: Problems closing H5Sclose\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + if ((status = H5Aclose(attr_id)) != SUCCEED ) { + fprintf(stderr, "Error: Problems closing H5Aclose\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if (attr_values != NULL) { + HDfree(attr_values); + } + + status = SUCCEED; + break; + case H5T_BITFIELD: + fprintf(stderr,"Warning: H5T_BITFIELD attribute not yet implemented.\n"); + break; + case H5T_OPAQUE: + fprintf(stderr,"Warning: H5T_OPAQUE attribute not yet implemented.\n"); + break; + case H5T_COMPOUND: + fprintf(stderr,"Warning: H5T_COMPOUND attribute not implemented.\n"); + break; + default: + fprintf(stderr,"Error: %d class not found\n",class); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + } + + if ((status = H5Tclose(type)) != SUCCEED ) { + fprintf(stderr, "Error: Problems closing H5Tclose\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_attr", __FILE__, __LINE__); + status = FAIL; + return status; + } + + } else { + + status = FAIL; + + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_all + * + * Purpose: Dump everything in the specified object + * + * Return: Success: SUCCEED + * + * Failure: FAIL + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_all (hid_t group, char *name, op_data_t *op_data) +{ + hid_t obj; + H5G_stat_t statbuf, statbuf2; + int status; + op_data_t op_data_save; + int32 vgroup_id; + int32 sd_id; + int32 sds_id; + int idx, flag; + void *edata; + hid_t (*func)(void*); + + op_data_save = *op_data; + + vgroup_id = op_data->vgroup_id; + sd_id = op_data->sd_id; + sds_id = op_data->sds_id; + + if ((status = H5Gget_objinfo(group, name, FALSE, &statbuf)) != SUCCEED ) { + fprintf(stderr,"Error: H5Gget_objinfo() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (status); + } + statbuf2 = statbuf; + + switch (statbuf.type) { + + case H5G_LINK: /* this is a soft link only */ + + + /* Disable error reporting */ + H5Eget_auto (&func, &edata); + H5Eset_auto (NULL, NULL); + + /* test to see if object exists */ + if ((status = H5Gget_objinfo(group, name, TRUE, NULL)) != SUCCEED ) { + fprintf(stderr,"Warning: the object pointed to by the symbolic link \"%s\" does not exist.\n",name); + } + + /* Enable error reporting */ + H5Eset_auto (func, edata); + + if (status != SUCCEED) { + break; + } + + /* follow link for type */ + if ((status = H5Gget_objinfo(group, name, TRUE, &statbuf)) != SUCCEED ) { + fprintf(stderr,"Error: H5Gget_objinfo() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (status); + } + + if (statbuf.type==H5G_DATASET ) { + + if ((idx = get_table_idx(dset_table, statbuf.objno)) < 0 ) { + + fprintf(stderr,"Error: object not found\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if((flag = get_tableflag(dset_table,idx)) < 0 ) { + + fprintf(stderr,"Error: get_tableflag() should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if(flag == TRUE ) { /* this has already been converted, add as a tag/ref */ + + if ((obj = H5Dopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open H5 dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if ((status = convert_shared_dataset(obj, idx, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to convert to tag/ref\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + } + } + if ((status = H5Dclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close H5 dataset %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } else { /* flag == FALSE */ + + if ((obj = H5Dopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open H5 dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if (( status = convert_dataset (obj, name, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_dataset did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((status = convert_shared_dataset(obj, idx, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to convert to tag/ref\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + } + if(( status = set_tableflag(dset_table,idx)) != SUCCEED ) { + fprintf(stderr,"Error: set_tableflag() did not work for %s\n", name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + break; + } + } + if ((status = H5Dclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close H5 dataset %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } + + } else if (statbuf.type==H5G_GROUP ) { + + if ((idx = get_table_idx(group_table, statbuf.objno)) < 0 ) { + + fprintf(stderr,"Error: object not found\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if((flag = get_tableflag(group_table,idx)) < 0 ) { + + fprintf(stderr,"Error: get_tableflag() should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if(flag == TRUE ) { + + if (( status = convert_shared_group (group, idx, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_group did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } else { /* flag == FALSE */ + + if ((obj = H5Gopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open group\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if (( status = convert_group (obj, name, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_group did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + } + if ((status = H5Gclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close group %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } + + } + + break; + + case H5G_GROUP: + + if ((idx = get_table_idx(group_table, statbuf.objno)) < 0 ) { + + fprintf(stderr,"Error: object not found\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if((flag = get_tableflag(group_table,idx)) < 0 ) { + + fprintf(stderr,"Error: get_tableflag() should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if(flag == TRUE ) { + + if (( status = convert_shared_group (group, idx, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_group did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } else { /* flag == FALSE */ + + if ((obj = H5Gopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open group\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if (( status = convert_group (obj, name, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_group did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + } + if ((status = H5Gclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close group %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } + + break; + + case H5G_DATASET: + + if ((idx = get_table_idx(dset_table, statbuf.objno)) < 0 ) { + + fprintf(stderr,"Error: object not found\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if((flag = get_tableflag(dset_table,idx)) < 0 ) { + + fprintf(stderr,"Error: get_tableflag() should never return < 0\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + + } else if(flag == TRUE ) { /* this has already been converted, add as a tag/ref */ + + if ((obj = H5Dopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open H5 dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if ((status = convert_shared_dataset(obj, idx, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to convert to tag/ref\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + } + } + if ((status = H5Dclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close H5 dataset %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } else { /* flag == FALSE */ + + if ((obj = H5Dopen (group, name)) <= 0 ) { + fprintf(stderr,"Error: Unable to open H5 dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + return (obj); + } else { + if (( status = convert_dataset (obj, name, op_data)) != SUCCEED) { + fprintf(stderr,"Error: convert_dataset did not work for %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + if ((status = convert_shared_dataset(obj, idx, op_data)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to convert to tag/ref\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + } + if(( status = set_tableflag(dset_table,idx)) != SUCCEED ) { + fprintf(stderr,"Error: set_tableflag() did not work for %s\n", name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + break; + } + } + if ((status = H5Dclose(obj)) != SUCCEED) { + fprintf(stderr,"Error: Unable to close H5 dataset %s\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_all", __FILE__, __LINE__); + status = FAIL; + break; + } + + } + + break; + + case H5G_TYPE: + /* object is ignored */ + break; + + default: + fprintf (stderr,"Unknown Object %s\n", name); + status = FAIL; + break; + + } + + *op_data = op_data_save; + + return SUCCEED; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_shared_dataset + * + * Purpose: Handle a shared dataset which has already been converted. + * + * Return: Success: SUCCEED + * + * Failure: FAIL + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_shared_dataset(hid_t did, int idx, op_data_t *op_data) +{ + int status=SUCCEED; + int32 vgroup_id; + char *dataset_name=NULL; + char *dataset_name2=NULL; + int32 hfile_id; + int32 sd_id; + int32 sds_id; + int32 sds_ref; + int32 vdata_ref; + int32 sds_index; + int32 numtagref; + hid_t type, space, class; + hsize_t dims[32], maxdims[32]; + int n_values, ndims; + + vgroup_id = op_data->vgroup_id; + + if ((dataset_name = get_objectname(dset_table, idx)) == NULL ) { + fprintf(stderr,"Error: get_objectname() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (status); + } + + if ((dataset_name2 = strrchr(dataset_name,'/')) == NULL) { /* find last "/" in dataset_name */ + dataset_name2 = dataset_name; /* no "/"s were found */ + } else { + dataset_name2 = dataset_name2 + sizeof(char); /* 1 character past last "/" */ + } + + if ((type = H5Dget_type(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_type() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((space = H5Dget_space(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_space() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((n_values = H5Sget_simple_extent_npoints(space)) <= 0) { + fprintf(stderr, "Error: H5Sget_simple_extent_npoints() returned inappropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((ndims = H5Sget_simple_extent_dims(space,dims,maxdims)) < 0 ) { + fprintf(stderr, "Error: Problems getting ndims, dims, and maxdims of dataset\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = ndims; + return status; + } + + if ((class = H5Tget_class(type)) < 0 ) { + fprintf(stderr,"Error: problem with getting class\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = class; + return status; + } + + switch (class) { + case H5T_INTEGER: + case H5T_FLOAT: + sd_id = op_data->sd_id; + if ((sds_index = SDnametoindex(sd_id, dataset_name2)) < 0 ) { + fprintf (stderr,"Error: Problem with SDnametoindex().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (sds_index); + } + if ((sds_id = SDselect(sd_id, sds_index)) < 0 ) { + fprintf (stderr,"Error: Problem with SDselect().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (sds_id); + } + if ((sds_ref = SDidtoref(sds_id)) < 0 ) { + fprintf (stderr,"Error: Problem with SDidtoref().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (sds_ref); + } + if ((numtagref = Vaddtagref(vgroup_id, DFTAG_NDG, sds_ref)) < 0 ) { + fprintf (stderr,"Error: Problem with Vaddtagref().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (numtagref); + } + break; + case H5T_TIME: + fprintf(stderr,"Error: H5T_TIME not yet implemented.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + break; + case H5T_STRING: + hfile_id = op_data->hfile_id; + if (ndims==1) { + if ((vdata_ref = VSfind(hfile_id,dataset_name2)) <= 0 ) { + fprintf (stderr,"Error: Problem with VSfind().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (vdata_ref); + } + if ((numtagref = Vaddtagref(vgroup_id, DFTAG_VH, vdata_ref)) < 0 ) { + fprintf (stderr,"Error: Problem with Vaddtagref().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (numtagref); + } + } + break; + case H5T_BITFIELD: + fprintf(stderr,"Error: H5T_BITFIELD not yet implemented.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + break; + case H5T_OPAQUE: + fprintf(stderr,"Error: H5T_OPAQUE not yet implemented.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + break; + case H5T_COMPOUND: + hfile_id = op_data->hfile_id; + if (ndims==1) { + if ((vdata_ref = VSfind(hfile_id,dataset_name2)) <= 0 ) { + fprintf (stderr,"Error: Problem with VSfind().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (vdata_ref); + } + if ((numtagref = Vaddtagref(vgroup_id, DFTAG_VH, vdata_ref)) < 0 ) { + fprintf (stderr,"Error: Problem with Vaddtagref().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + return (numtagref); + } + } + break; + default: + fprintf(stderr,"Error: %d class not found\n",class); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = FAIL; + } + + if ((status = H5Tclose(type)) < 0 ) { + fprintf(stderr,"Error: closing type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset", __FILE__, __LINE__); + status = FAIL; + } + + if (dataset_name != NULL) { + HDfree (dataset_name); + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_shared_group + * + * Purpose: Handle a shared group which has already been converted. + * + * Return: status + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_shared_group (hid_t group, int idx, op_data_t *op_data) { + + int32 hfile_id; + int32 vgroup_id; + int32 vgroup_ref; + int32 numtagref; + int32 status = SUCCEED; + hid_t group2; + char *group_name=NULL; + char *group_name2=NULL; + char vgroup_name[VGNAMELENMAX]; + + group2 = group; + hfile_id = op_data->hfile_id; + + if ((group_name = get_objectname(group_table, idx)) == NULL ) { + fprintf(stderr,"Error: get_objectname() did not work\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + status = FAIL; + return (status); + } + + if ((group_name2 = strrchr(group_name,'/')) == NULL) { /* find last "/" in group_name */ + group_name2 = group_name; /* no "/"s were found */ + } else { + group_name2 = group_name2 + sizeof(char); /* 1 character past last "/" */ + } + + if ((status = Vgetname(op_data->vgroup_id,vgroup_name)) < 0 ) { + fprintf (stderr,"Error: Problem with Vfind().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + return (status); + } + + if ((vgroup_ref = Vfind(hfile_id,vgroup_name)) <= 0 ) { + fprintf (stderr,"Error: Problem with Vfind().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + return (vgroup_ref); + } + + if ((vgroup_id = Vattach(hfile_id, vgroup_ref, "w")) <= 0 ) { + fprintf(stderr,"Error: Unable to create new vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + return(vgroup_id); + } + + if ((vgroup_ref = Vfind(hfile_id,group_name2)) <= 0 ) { + fprintf (stderr,"Error: Problem with Vfind().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + return (vgroup_ref); + } + + if ((numtagref = Vaddtagref(vgroup_id, DFTAG_VG, vgroup_ref)) < 0 ) { + fprintf (stderr,"Error: Problem with Vaddtagref().\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + return (numtagref); + } + + if ((status = Vdetach(vgroup_id)) != SUCCEED ) { + fprintf(stderr,"Error: Unable to detach the new Vgroup\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_group", __FILE__, __LINE__); + status = FAIL; + } + + if (group_name != NULL) { + HDfree(group_name); + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: convert_dataset_string + * + * Purpose: convert a dataset with string elements. + * + * Return: status + * + * Programmer: Paul Harten + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +herr_t +convert_dataset_string (hid_t did, char *name, op_data_t *op_data) { + + int32 hfile_id; + int32 vgroup_id; + int32 vdata_id; +hid_t fxdlenstr, space, class, mem_type; +const char* fieldname = {"string"}; +const char* fieldname_list = fieldname; +char *buffer; +int32 status; +int32 h4_type; +int32 recsize, n_records, n_values, num_of_recs, record_pos; +size_t lenstr; +H5T_cset_t cset; +H5T_str_t strpad; + + if ((fxdlenstr = H5Dget_type(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_type() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((class = H5Tget_class(fxdlenstr)) != H5T_STRING ) { + fprintf(stderr,"Error: problem with getting class\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_shared_dataset", __FILE__, __LINE__); + status = class; + return status; + } + + if ((space = H5Dget_space(did)) <= 0) { + fprintf(stderr, "Error: H5Dget_space() didn't return appropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = FAIL; + return status; + } + + if ((n_values = H5Sget_simple_extent_npoints(space)) <= 0) { + fprintf(stderr, "Error: H5Sget_simple_extent_npoints() returned inappropriate value.\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = FAIL; + return status; + } + + h4_type = DFNT_CHAR; + if ((mem_type = H5Tcopy(H5T_C_S1)) == FAIL ) { + fprintf(stderr, "Error: Problems translating h4 type to mem type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = mem_type; + return status; + } + if ((lenstr = H5Tget_size(fxdlenstr)) <= 0 ) { + fprintf(stderr, "Error: size of fixed length string type should not be %d\n",(int)lenstr); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return lenstr; + } + if ((status = H5Tset_size(mem_type,lenstr)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_size()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + strpad = H5Tget_strpad(fxdlenstr); + if ((strpad != H5T_STR_NULLTERM) && (strpad != H5T_STR_NULLPAD) && (strpad != H5T_STR_SPACEPAD)) { + fprintf(stderr, "Error: Invalid string padding value, %d\n",(int)strpad); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return strpad; + } + if ((status = H5Tset_strpad(mem_type, strpad)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_strpad()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + if ((cset = H5Tget_cset(fxdlenstr)) != H5T_CSET_ASCII ) { + fprintf(stderr, "Error: cset value != %d\n",(int)H5T_CSET_ASCII); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return cset; + } + if ((status = H5Tset_cset(mem_type,cset)) != SUCCEED ) { + fprintf(stderr, "Error: Problem with H5Tset_cset()\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + hfile_id = op_data->hfile_id; + vgroup_id = op_data->vgroup_id; + if ((vdata_id = VSattach(hfile_id, -1, "w")) <= 0 ) { + fprintf(stderr, "Error: Unable to create vdata %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = FAIL; + return status; + } + + op_data->vdata_id = vdata_id; + if ((status = VSsetname(vdata_id, name)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set vdata name %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + if ((status = VSsetclass(vdata_id, "HDF5")) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set class on vdata %s\n", name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + if ((status = VSfdefine(vdata_id, fieldname, h4_type, lenstr)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to VSfdefine() field\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + if ((status = VSsetfields(vdata_id, fieldname_list)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set fieldname list %s\n", fieldname_list); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + if ((status = VSsetinterlace(vdata_id, FULL_INTERLACE)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to set FULL_INTERLACE mode, status %d\n", (int)status); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + if ((recsize = H5Tget_size(mem_type)) <= 0 ) { + fprintf(stderr, "Error: Unable to get record size %d\n", (int)recsize); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = recsize; + return status; + } +/* + Since the space is rank 1, n_records does not depend on maxdims. +*/ + n_records = n_values; + if ((buffer = HDmalloc(n_records*recsize)) == NULL) { + fprintf(stderr, "Error: Problems with HDmalloc of memory space\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = FAIL; + return status; + } + if ((status = H5Dread(did, mem_type, space, space, H5P_DEFAULT, buffer)) != SUCCEED) { + fprintf(stderr, "Error: Problems with H5Dread\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + if ((record_pos = VSseek(vdata_id, 0)) != 0 ) { + fprintf(stderr, "Error: Could not seek the beginning of the Vdata, %d\n", (int)record_pos); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = record_pos; + return status; + } + if ((num_of_recs = VSwrite(vdata_id, (void *)buffer, n_records, FULL_INTERLACE)) != n_records ) { + fprintf(stderr, "Error: Only able to write %d of %d records\n", (int)num_of_recs, (int)n_records); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + status = num_of_recs; + return status; + } + + /* there are only vdata attributes, no field attributes */ + if ((status = H5Aiterate(did, NULL, (H5A_operator_t)convert_attr, op_data)) < 0 ) { + fprintf(stderr,"Error: iterate over attributes\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + if ((status = VSdetach(vdata_id)) != SUCCEED ) { + fprintf(stderr, "Error: Unable to detach to vdata %s.\n",name); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "convert_dataset_string", __FILE__, __LINE__); + return status; + } + + if (buffer != NULL) { + HDfree(buffer); + } + + return status; + +} + + +/*------------------------------------------------------------------------- + * Function: h5atomic_type_to_h4type + * + * Purpose: Match an H5 atomic type to an appropriate H4 type. + * Assign an appropirate H5 memory type that matches the H4 type. + * Return the H5 memory type, H4 type and sizeof H4 via pointers. + * + * Return: SUCCEED on suceed, FAIL on failure. + * When fail, pointer values of h5memtype, h5memsize and h4_type + * may have changed and are undefined. + * + * Programmer: Albert Cheng, March 2000 + * + * Modifications: + * + *-----------------------------------------------------------------------*/ +static herr_t h5atomic_type_to_h4type(const hid_t h5type, hid_t* h5memtype, size_t* h5memsize, int32* h4type) +{ + H5T_class_t class; + size_t h5typesize, h4typesize; + H5T_sign_t sign; + hid_t mem_datatype = FAIL; + + if ((class = H5Tget_class(h5type)) < 0 ) { + fprintf(stderr,"Error: problem with getting type class\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + + switch(class){ + case H5T_INTEGER: + if ((h5typesize = H5Tget_size(h5type)) == 0 ) { + fprintf(stderr,"Error: problem with getting type size\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + if ((sign = H5Tget_sign(h5type)) == H5T_SGN_ERROR) { + fprintf(stderr,"Error: problem with getting type sign\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + /* deduce the proper HDF4 integer type to use according to the size of the HDF5 type. */ + /* Current HDF4 types can be 8, 16, 32 bits (1, 2, 4 bytes). */ + switch(h5typesize){ + case 1: + *h4type = (sign == H5T_SGN_2 ? DFNT_INT8 : DFNT_UINT8); + h4typesize = sizeof(int8); + if (h4typesize == H5Tget_size(H5T_NATIVE_CHAR)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SCHAR : H5T_NATIVE_UCHAR); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_SHORT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SHORT : H5T_NATIVE_USHORT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_INT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_INT : H5T_NATIVE_UINT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_LONG)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_LONG : H5T_NATIVE_ULONG); + }else + return(FAIL); + break; + case 2: + *h4type = (sign == H5T_SGN_2 ? DFNT_INT16 : DFNT_UINT16); + h4typesize = sizeof(int16); + if (h4typesize == H5Tget_size(H5T_NATIVE_CHAR)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SCHAR : H5T_NATIVE_UCHAR); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_SHORT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SHORT : H5T_NATIVE_USHORT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_INT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_INT : H5T_NATIVE_UINT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_LONG)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_LONG : H5T_NATIVE_ULONG); + }else + return(FAIL); + break; + case 4: + *h4type = (sign == H5T_SGN_2 ? DFNT_INT32 : DFNT_UINT32); + h4typesize = sizeof(int32); + if (h4typesize == H5Tget_size(H5T_NATIVE_CHAR)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SCHAR : H5T_NATIVE_UCHAR); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_SHORT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_SHORT : H5T_NATIVE_USHORT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_INT)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_INT : H5T_NATIVE_UINT); + }else if (h4typesize == H5Tget_size(H5T_NATIVE_LONG)){ + mem_datatype = (sign == H5T_SGN_2 ? H5T_NATIVE_LONG : H5T_NATIVE_ULONG); + }else + return(FAIL); + break; + default: + fprintf(stderr,"Error: unmatchable integer type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + break; + /* end of case H5T_INTEGER */ + + case H5T_FLOAT: + if ((h5typesize = H5Tget_size(h5type)) == 0 ) { + fprintf(stderr,"Error: problem with getting type size\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + /* deduce the proper HDF4 floating point type to use according to the size of the HDF5 type. */ + /* Current HDF4 types can be 32 or 64 bits (4 or 8 bytes). */ + switch(h5typesize){ + case 4: + *h4type = DFNT_FLOAT32; + h4typesize = sizeof(float32); + if (h4typesize == H5Tget_size(H5T_NATIVE_FLOAT)){ + mem_datatype = H5T_NATIVE_FLOAT; + }else if (h4typesize == H5Tget_size(H5T_NATIVE_DOUBLE)){ + mem_datatype = H5T_NATIVE_DOUBLE; + }else + return(FAIL); + break; + case 8: + *h4type = DFNT_FLOAT64; + h4typesize = sizeof(float64); + if (h4typesize == H5Tget_size(H5T_NATIVE_FLOAT)){ + mem_datatype = H5T_NATIVE_FLOAT; + }else if (h4typesize == H5Tget_size(H5T_NATIVE_DOUBLE)){ + mem_datatype = H5T_NATIVE_DOUBLE; + }else + return(FAIL); + break; + default: + fprintf(stderr,"Error: unmatchable H5 float type\n"); + DEBUG_PRINT("Error detected in %s() [%s line %d]\n", "h5atomic_type_to_h4type", __FILE__, __LINE__); + return FAIL; + } + break; + /* end of case H5T_FLOAT */ + + default: + return FAIL; + } + + *h5memsize = h4typesize; + *h5memtype = mem_datatype; + return(SUCCEED); +} +#ifndef NEWWAY +/***************************************************************************** + + Routine: h5type_to_h4type(h5type) + + Description: Translate h5 datatype into h4 datatype + + Input: h5 datatype + + Output: function return, h4 datatype + +*****************************************************************************/ + +int32 h5type_to_h4type(hid_t h5_datatype) +{ + int32 h4_datatype; + + if (H5Tequal(h5_datatype,H5T_STD_I8BE)) { + h4_datatype = DFNT_INT8; + } else if (H5Tequal(h5_datatype,H5T_STD_I16BE)) { + h4_datatype = DFNT_INT16; + } else if (H5Tequal(h5_datatype,H5T_STD_I32BE)) { + h4_datatype = DFNT_INT32; +/* + * This is not supported by HDF4 + * + } else if (H5Tequal(h5_datatype,H5T_STD_I64BE)) { + h4_datatype = DFNT_INT64; +*/ + } else if (H5Tequal(h5_datatype,H5T_STD_U8BE)) { + h4_datatype = DFNT_UINT8; + } else if (H5Tequal(h5_datatype,H5T_STD_U16BE)) { + h4_datatype = DFNT_UINT16; + } else if (H5Tequal(h5_datatype,H5T_STD_U32BE)) { + h4_datatype = DFNT_UINT32; + } else if (H5Tequal(h5_datatype,H5T_STD_U64BE)) { + h4_datatype = DFNT_UINT64; + } else if (H5Tequal(h5_datatype,H5T_IEEE_F32BE)) { + h4_datatype = DFNT_FLOAT32; + } else if (H5Tequal(h5_datatype,H5T_IEEE_F64BE)) { + h4_datatype = DFNT_FLOAT64; + } else if (H5Tequal(h5_datatype,H5T_STD_I8LE)) { + h4_datatype = DFNT_INT8; + } else if (H5Tequal(h5_datatype,H5T_STD_I16LE)) { + h4_datatype = DFNT_INT16; + } else if (H5Tequal(h5_datatype,H5T_STD_I32LE)) { + h4_datatype = DFNT_INT32; +/* + * This is not supported by HDF4 + * + } else if (H5Tequal(h5_datatype,H5T_STD_I64LE)) { + h4_datatype = DFNT_INT64; +*/ + } else if (H5Tequal(h5_datatype,H5T_STD_U8LE)) { + h4_datatype = DFNT_UINT8; + } else if (H5Tequal(h5_datatype,H5T_STD_U16LE)) { + h4_datatype = DFNT_UINT16; + } else if (H5Tequal(h5_datatype,H5T_STD_U32LE)) { + h4_datatype = DFNT_UINT32; + } else if (H5Tequal(h5_datatype,H5T_STD_U64LE)) { + h4_datatype = DFNT_UINT64; + } else if (H5Tequal(h5_datatype,H5T_IEEE_F32LE)) { + h4_datatype = DFNT_FLOAT32; + } else if (H5Tequal(h5_datatype,H5T_IEEE_F64LE)) { + h4_datatype = DFNT_FLOAT64; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_SCHAR)) { + h4_datatype = DFNT_INT8; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_UCHAR)) { + h4_datatype = DFNT_UINT8; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_SHORT)) { + h4_datatype = DFNT_INT16; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_USHORT)) { + h4_datatype = DFNT_UINT16; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_INT)) { + h4_datatype = DFNT_INT32; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_UINT)) { + h4_datatype = DFNT_UINT32; +/* + * This is not supported by HDF4 + * + } else if (H5Tequal(h5_datatype,H5T_NATIVE_LONG)) { + h4_datatype = DFNT_INT64; +*/ + } else if (H5Tequal(h5_datatype,H5T_NATIVE_ULONG)) { + h4_datatype = DFNT_UINT64; +/* + * This is not supported by HDF4 + * + } else if (H5Tequal(h5_datatype,H5T_NATIVE_LLONG)) { + h4_datatype = DFNT_INT128; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_ULLONG)) { + h4_datatype = DFNT_UINT128; +*/ + } else if (H5Tequal(h5_datatype,H5T_NATIVE_FLOAT)) { + h4_datatype = DFNT_FLOAT32; + } else if (H5Tequal(h5_datatype,H5T_NATIVE_DOUBLE)) { + h4_datatype = DFNT_FLOAT64; +/* + * This is not supported by HDF4 + * + } else if (H5Tequal(h5_datatype,H5T_NATIVE_LDOUBLE)) { + h4_datatype = DFNT_FLOAT128; +*/ + } else { + h4_datatype = FAIL; + } + + return h4_datatype; + +} + + +/***************************************************************************** + + Routine: h4type_to_memtype(h4type) + + Description: Translate h4 datatype into mem datatype + + Input: h4 datatype + + Output: function return, mem datatype + +*****************************************************************************/ + +hid_t h4type_to_memtype(int32 h4_datatype) +{ + hid_t mem_datatype; + + switch (h4_datatype) { + case DFNT_INT8: + case DFNT_NINT8: + case DFNT_LINT8: + mem_datatype = H5T_NATIVE_SCHAR; break; + case DFNT_UINT8: + case DFNT_NUINT8: + case DFNT_LUINT8: + mem_datatype = H5T_NATIVE_UCHAR; break; + case DFNT_INT16: + case DFNT_NINT16: + case DFNT_LINT16: + mem_datatype = H5T_NATIVE_SHORT; break; + case DFNT_UINT16: + case DFNT_NUINT16: + case DFNT_LUINT16: + mem_datatype = H5T_NATIVE_USHORT; break; + case DFNT_INT32: + case DFNT_NINT32: + case DFNT_LINT32: + mem_datatype = H5T_NATIVE_INT; break; + case DFNT_UINT32: + case DFNT_NUINT32: + case DFNT_LUINT32: + mem_datatype = H5T_NATIVE_UINT; break; + case DFNT_INT64: + case DFNT_NINT64: + case DFNT_LINT64: + mem_datatype = H5T_NATIVE_LONG; break; + case DFNT_UINT64: + case DFNT_NUINT64: + case DFNT_LUINT64: + mem_datatype = H5T_NATIVE_ULONG; break; + case DFNT_INT128: + case DFNT_NINT128: + case DFNT_LINT128: + mem_datatype = H5T_NATIVE_LLONG; break; + case DFNT_UINT128: + case DFNT_NUINT128: + case DFNT_LUINT128: + mem_datatype = H5T_NATIVE_ULLONG; break; + case DFNT_FLOAT32: + case DFNT_NFLOAT32: + case DFNT_LFLOAT32: + mem_datatype = H5T_NATIVE_FLOAT; break; + case DFNT_FLOAT64: + case DFNT_NFLOAT64: + case DFNT_LFLOAT64: + mem_datatype = H5T_NATIVE_DOUBLE; break; + case DFNT_FLOAT128: + case DFNT_NFLOAT128: + case DFNT_LFLOAT128: + mem_datatype = H5T_NATIVE_LDOUBLE; break; + default: + mem_datatype = FAIL; + } + + return mem_datatype; + +} +#endif + + +/***************************************************************************** + + Routine: test_file + + Description: Test a file for read/write - ability. + + Input: filename - Unix filename + + Output: function return, global variable - errno + +*****************************************************************************/ + +int test_file(char *filename,int oflag,mode_t mode) +{ + int fid; + + errno = 0; + + fid = open(filename, oflag, mode); + if (fid < 0) { + perror(filename); + } + close(fid); + + return errno; + +} + + +/***************************************************************************** + + Routine: test_dir + + Description: Test pathway to determine if it is a directory + + Input: path - pathname given + + Output: function return TRUE/FALSE + +*****************************************************************************/ + +int test_dir(char *path) +{ + + struct stat buf; + struct stat *buf_ptr; + int idir; + + buf_ptr = &buf; + + idir = stat(path, buf_ptr); + if (idir < 0) { + if (errno == 2) { + return 0; + } else { + perror(path); + } + } + + return S_ISDIR(buf_ptr->st_mode); +} + +/***************************************************************************** + + Routine: BuildFilename() + + Description: Build a filename with new extension + + Input: filename - present filename + ext - extension to root of filename + + Output: (filename:r).ext + +*****************************************************************************/ + +char *BuildFilename(char *filename, char *ext) +{ + /* build outgoing filename */ + + char *filename_out; + char *lastper_ptr, *lastdir_ptr; + int root_len; + + lastper_ptr = strrchr(filename,'.'); + lastdir_ptr = strrchr(filename,'/'); + + if ( lastper_ptr <= lastdir_ptr ) { /* no extension */ + root_len = strlen(filename); + } else { /* existing extension */ + root_len = (int)(lastper_ptr - filename); + } + + filename_out = (char *)HDmalloc(root_len + strlen(ext) + 2); + filename_out = strncpy(filename_out, filename, (size_t)root_len); + filename_out[root_len] = '\0'; + filename_out = strcat(filename_out,"."); + filename_out = strcat(filename_out,ext); + + return filename_out; +} + diff --git a/tools/h5toh4/h5toh4.h b/tools/h5toh4/h5toh4.h new file mode 100644 index 0000000..d268060 --- /dev/null +++ b/tools/h5toh4/h5toh4.h @@ -0,0 +1,43 @@ +/* + * Copyright © 1998 NCSA + * All rights reserved. + * + * Programmer: Paul Harten <pharten@ncsa.uiuc.edu> + * Friday, October 16th, 1998 + * + * Purpose: Convert H5 files to H4 files. + */ + +#ifndef _H5TOH4_H +#define _H5TOH4_H + +#include "hdf.h" +#include "mfhdf.h" +#include "hdf5.h" + +#ifdef H5_HAVE_SYS_STAT_H +# include <sys/stat.h> +#endif + +typedef struct op_data_t { + /* + * information being carried between iterations. + * + */ + + int32 hfile_id; + int32 vgroup_id; + int32 sd_id; + int32 sds_id; + int32 vdata_id; + int32 obj_idx; + +} op_data_t; + +#ifdef H5TOH4_DEBUG +#define DEBUG_PRINT(s1,s2,s3,n1) ( fprintf(stderr,s1,s2,s3,n1) ) +#else +#define DEBUG_PRINT(s1,s2,s3,n1) ( fprintf(stderr," ") ) +#endif + +#endif diff --git a/tools/h5toh4/testh5toh4.sh b/tools/h5toh4/testh5toh4.sh new file mode 100755 index 0000000..986a6a2 --- /dev/null +++ b/tools/h5toh4/testh5toh4.sh @@ -0,0 +1,299 @@ +#! /bin/sh +# +# Copyright (C) 1997-2001 National Center for Supercomputing Applications. +# All rights reserved. +# + +H5TOH4=h5toh4 # a relative name +H5TOH4_BIN=`pwd`/$H5TOH4 # an absolute command path + +cmp='cmp -s' +diff='diff -c' + +RM='rm -f' +SED='sed ' +H4DUMP='hdp' + +# Verify if $H4DUMP is a valid command. +tmpfile=/tmp/testh5toh4.$$ +$H4DUMP -H > $tmpfile +if test -s "$tmpfile"; then + # Find out which version of hdp is being used. Over simplified + # algorithm but will do the job for now. + if ( grep -s 'NCSA HDF Version 4.1 Release [3-9]' $tmpfile > /dev/null ) + then + H4DUMPVER=413 + else + H4DUMPVER=0 + echo " Some tests maybe skipped because your version of $H4DUMP does" + echo " not handle loops in Vgroups correctly. You need version" + echo " 4.1 Release 3 or later. Visit http://hdf.ncsa.uiuc.edu" + echo " or email hdfhelp@ncsa.uiuc.edu for more information." + fi +else + echo " Could not run the '$H4DUMP' command. The test can still proceed" + echo " but it may fail if '$H4DUMP' is needed to verify the output." + echo " You can make sure '$H4DUMP' is among your shell PATH and run" + echo " the test again. You may also visit http://hdf.ncsa.uiuc.edu" + echo " or email hdfhelp@ncsa.uiuc.edu for more information." + H4DUMP=: + H4DUMPVER=0 +fi +$RM $tmpfile + +# The build (current) directory might be different than the source directory. +if test -z "$srcdir"; then + srcdir=. +fi +mkdir ../testfiles >/dev/null 2>&1 + +SRCDIR="$srcdir/../testfiles" +OUTDIR="../testfiles/Results" + +test -d $OUTDIR || mkdir $OUTDIR + +nerrors=0 +verbose=yes + +# Print a line-line message left justified in a field of 70 characters +# beginning with the word "Testing". +TESTING() { + SPACES=" " + echo "Testing $* $SPACES" |cut -c1-70 |tr -d '\012' +} + +# Run a test and print PASS or *FAIL*. If a test fails then increment +# the `nerrors' global variable and (if $verbose is set) display the +# difference between the actual and the expected hdf4 files. The +# expected hdf4 files are in testfiles/Expected directory. +# The actual hdf4 file is not removed if $HDF5_NOCLEANUP is to a non-null +# value. +CONVERT() { + # Run h5toh4 convert. + TESTING $H5TOH4 $@ + + # + # Set up arguments to run the conversion test. + # The converter assumes all hdf5 files has the .h5 suffix as in the form + # of foo.h5. It creates the corresponding hdf4 files with the .hdf suffix + # as in the form of foo.hdf. One exception is that if exactly two file + # names are given, it treats the first argument as an hdf5 file and creates + # the corresponding hdf4 file with the name as the second argument, WITOUT + # any consideration of the suffix. (For this test script, in order to + # match the output hdf4 file with the expected hdf4 file, it expects the + # second file of the two-files tests has the .hdf suffix too.) + # + # If SRCDIR != OUTDIR, need to copy the input hdf5 files from the SRCDIR + # to the OUTDIR and transform the input file pathname because of the suffix + # convention mentioned above. This way, the hdf4 files are always created + # in the OUTDIR directory. + # + + INFILES="" + OUTFILES="" + MULTIRUN="" + + case "$1" in + "-m") # multiple files conversion + MULTIRUN="-m" + shift + for f in $* + do + if test "$SRCDIR" != "$OUTDIR"; then + cp $SRCDIR/$f $OUTDIR/$f + fi + INFILES="$INFILES $f" + OUTFILES="$OUTFILES `basename $f .h5`.hdf" + shift + done + ;; + * ) # Single file conversion + case $# in + 1) if test "$SRCDIR" != "$OUTDIR"; then + cp $SRCDIR/$1 $OUTDIR/$1 + fi + INFILES="$1" + OUTFILES="`basename $1 .h5`.hdf" + ;; + 2) # hdf4 file specified + if test "$SRCDIR" != "$OUTDIR"; then + cp $SRCDIR/$1 $OUTDIR/$1 + fi + INFILES="$1" + OUTFILES="$2" + ;; + *) # Illegal + echo "Illegal arguments" + exit 1 + ;; + esac + ;; + esac + + # run the conversion and remove input files that have been copied over + ( + cd $OUTDIR + $H5TOH4_BIN $MULTIRUN $INFILES 2>/dev/null + if test "$SRCDIR" != "$OUTDIR"; then + $RM $INFILES + fi + ) + + # Verify results + result="passed" + for f in $OUTFILES + do + if $cmp $SRCDIR/Expected/$f $OUTDIR/$f + then + : + else + # Use hdp to dump the files and verify the output. + # Filter out the output of "reference = ..." because + # reference numbers are immaterial in general. + outfile=`basename $f .hdf` + expect_out=$outfile.expect + actual_out=$outfile.actual + + if [ $outfile = "tloop" -a $H4DUMPVER -lt 413 ] + then + echo " -SKIP-" + result="skipped" + touch $expect_out $actual_out # fake them + else + (cd $SRCDIR/Expected + $H4DUMP dumpvg $outfile.hdf + $H4DUMP dumpvd $outfile.hdf + $H4DUMP dumpsds $outfile.hdf ) | + sed -e 's/reference = [0-9]*;//' > $expect_out + (cd $OUTDIR + $H4DUMP dumpvg $outfile.hdf + $H4DUMP dumpvd $outfile.hdf + $H4DUMP dumpsds $outfile.hdf ) | + sed -e 's/reference = [0-9]*;//' > $actual_out + fi + + if [ "passed" = $result -a ! -s $actual_out ] ; then + echo "*FAILED*" + nerrors="`expr $nerrors + 1`" + result=failed + test yes = "$verbose" && + echo " H4DUMP failed to produce valid output" + elif $cmp $expect_out $actual_out; then + : + else + if test "passed" = $result; then + echo "*FAILED*" + nerrors="`expr $nerrors + 1`" + result=failed + fi + test yes = "$verbose" && + echo " Actual result (*.actual) differs from expected result (*.expect)" && + $diff $expect_out $actual_out |sed 's/^/ /' + fi + fi + + # Clean up output file + if test -z "$HDF5_NOCLEANUP"; then + $RM $expect_out $actual_out + $RM $OUTDIR/$f + fi + done + if test "passed" = "$result"; then + echo " PASSED" + fi +} + + + +############################################################################## +############################################################################## +### T H E T E S T S ### +############################################################################## +############################################################################## + +$RM $OUTDIR/*.hdf $OUTDIR/*.tmp + +# +# The HDF4 filenames are created based upon the HDF5 filenames +# without the extension. +# + +# test for converting H5 groups to H4 Vgroups. +CONVERT tgroup.h5 + +# test for converting H5 datasets to H4 SDS's. +CONVERT tdset.h5 + +# test for converting H5 attributes to H4 attributes. +CONVERT tattr.h5 + +# test for converting H5 soft links. +CONVERT tslink.h5 + +# test for converting H5 hard links. +CONVERT thlink.h5 + +# test for converting H5 compound data type to H4 Vdata. +CONVERT tcompound.h5 + +# test for converting all H5 objects at in same file. +CONVERT tall.h5 + +# tests for converting H5 objects with loops. +CONVERT tloop.h5 + +# test for converting extendable H5 datasets to H4 SDS's. +CONVERT tdset2.h5 + +# test for converting extendable H5 datasets with compound data type to H4 Vdata. +CONVERT tcompound2.h5 + +# tests for converting H5 objects from many different pathways. +CONVERT tmany.h5 + +# tests for converting H5 string objects. +CONVERT tstr.h5 + +# tests for converting more H5 string objects. +CONVERT tstr2.h5 + +# +# The test for conversion are the same as above with the only difference +# being that the HDF4 filenames are given explicitly. +# + +$RM $OUTDIR/*.tmp +CONVERT tgroup.h5 tgroup.hdf +CONVERT tdset.h5 tdset.hdf +CONVERT tattr.h5 tattr.hdf +CONVERT tslink.h5 tslink.hdf +CONVERT thlink.h5 thlink.hdf +CONVERT tcompound.h5 tcompound.hdf +CONVERT tall.h5 tall.hdf +CONVERT tloop.h5 tloop.hdf +CONVERT tdset2.h5 tdset2.hdf +CONVERT tcompound2.h5 tcompound2.hdf +CONVERT tmany.h5 tmany.hdf +CONVERT tstr.h5 tstr.hdf +CONVERT tstr2.h5 tstr2.hdf + +# +# Again, the test for conversion are the same as the first set of test. +# Here, multiple conversion are done on HDF5 files at one time. +# + +$RM $OUTDIR/*.hdf $OUTDIR/*.tmp +CONVERT -m tgroup.h5 tdset.h5 tattr.h5 tslink.h5 thlink.h5 +CONVERT -m tcompound.h5 tall.h5 +CONVERT -m tloop.h5 +CONVERT -m tdset2.h5 tcompound2.h5 tmany.h5 +CONVERT -m tstr.h5 tstr2.h5 + +if test $nerrors -eq 0 ; then + echo "All h5toh4 tests passed." +fi + +if test -z "$HDF5_NOCLEANUP"; then + $RM -r $OUTDIR +fi +exit $nerrors |