question about writing an array of compound...

• To: netcdf-hdf@xxxxxxxxxxxxxxxx
• Subject: question about writing an array of compound...
• From: Ed Hartnett <ed@xxxxxxxxxxxxxxxx>
• Date: Fri, 25 Mar 2005 15:17:28 -0700

Guys,

What am I doing wrong here? The program below creates a dataset the
way I expect, but the numbers that it writes are not what I expect to
get. Here's the h5dump of the output:

HDF5 "tst_h_compounds.h5" {
GROUP "/" {
   GROUP "group1" {
      DATATYPE "compound_type_1" H5T_COMPOUND {
         H5T_STD_I32LE "f1";
         H5T_STD_I32LE "f2";
      }
      DATASET "v1" {
         DATATYPE  "/group1/compound_type_1"
         DATASPACE  SIMPLE { ( 3 ) / ( 3 ) }
         DATA {
         (0): {
               0,
               0
            },
         (1): {
               24641422,
               -1073743792
            },
         (2): {
               1455352,
               0
            }
         }
      }
   }
}
}


#include "tests.h"

#define FILE_NAME "tst_h_compounds.h5"
#define DIM1_LEN 3
#define GRP_NAME "group1"
#define VAR1_NAME "v1"
#define FIELD1_NAME "f1"
#define FIELD2_NAME "f2"
#define COMPOUND_NAME "compound_type_1"

int
main()
{
   hid_t fileid, grpid, spaceid, typeid, datasetid;
   int bool_out[DIM1_LEN] = {0, 1, 0}, bool_in[DIM1_LEN];
   hsize_t dims[1];
   struct s1 {
         int i1, i2;
   } data[DIM1_LEN];
   int i;

   printf("*** Checking HDF5 compound types...");
   
   for (i; i<DIM1_LEN; i++)
   {
      data[i].i1 = 99;
      data[i].i2 = -99;
   }

   /* Open file and create group. */
   if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT, 
                           H5P_DEFAULT)) < 0) ERR;
   if ((grpid = H5Gcreate(fileid, GRP_NAME, 0)) < 0) ERR;

   /* Create a simple compound type. */
   if ((typeid = H5Tcreate(H5T_COMPOUND, 8))) ERR;
   if (H5Tinsert(typeid, FIELD1_NAME, 0, H5T_NATIVE_INT) < 0) ERR;
   if (H5Tinsert(typeid, FIELD2_NAME, 4, H5T_NATIVE_INT) < 0) ERR;
   if (H5Tcommit(grpid, COMPOUND_NAME, typeid)) ERR;

   /* Create a space. */
   dims[0] = DIM1_LEN;
   if ((spaceid = H5Screate_simple(1, dims, dims)) < 0) ERR;

   /* Create a dataset of this compound type. */
   if ((datasetid = H5Dcreate(grpid, VAR1_NAME, typeid, 
                              spaceid, H5P_DEFAULT)) < 0) ERR;

   if (H5Dwrite(datasetid, typeid, H5S_ALL, H5S_ALL, 
                H5P_DEFAULT, data) < 0) ERR;

   if (H5Dclose(datasetid) < 0 ||
       H5Tclose(typeid) < 0 ||
       H5Sclose(spaceid) < 0 ||
       H5Gclose(grpid) < 0 ||
       H5Fclose(fileid) < 0) ERR;

   SUMMARIZE_ERR;

   return 0;
}

-- 
Ed Hartnett  -- ed@xxxxxxxxxxxxxxxx

• Follow-Ups:
  • Re: question about writing an array of compound...
    • From: Quincey Koziol