added timing test of HDF5 many attributes creation

nc_test4/Makefile.am

@@ -32,7 +32,7 @@ endif # LARGE_FILE_TESTS
 
 if BUILD_BENCHMARKS
 check_PROGRAMS += tst_create_files bm_file tst_chunks tst_ar4	\
-tst_ar4_3d tst_ar4_4d bm_many_objs bm_many_atts
+tst_ar4_3d tst_ar4_4d bm_many_objs tst_h_many_atts bm_many_atts
 TESTS += tst_ar4_3d tst_create_files run_bm_test1.sh run_bm_elena.sh	\
 run_bm_test2.sh run_tst_chunks.sh run_bm_ar4.sh
 

nc_test4/bm_many_atts.c

@@ -1,3 +1,4 @@
+
 /*
 Copyright 2010, UCAR/Unidata
 See COPYRIGHT file for copying and redistribution conditions.

nc_test4/tst_h_many_atts.c

@@ -0,0 +1,104 @@
+/* This is part of the netCDF package.
+   Copyright 2005 University Corporation for Atmospheric Research/Unidata
+   See COPYRIGHT file for conditions of use.
+
+   Test HDF5 file code. These are not intended to be exhaustive tests,
+   but they use HDF5 the same way that netCDF-4 does, so if these
+   tests don't work, than netCDF-4 won't work either.
+
+   $Id$
+*/
+#include <config.h>
+#include <nc_tests.h>
+#include "netcdf.h"
+#include <nc4internal.h>
+#include <H5DSpublic.h>
+#include <time.h>
+#include <sys/time.h> /* Extra high precision time info. */
+
+#define FILE_NAME "tst_h_many_atts.h5"
+#define GRP_NAME "group1"
+
+/* Subtract the `struct timeval' values X and Y, storing the result in
+   RESULT.  Return 1 if the difference is negative, otherwise 0.  This
+   function from the GNU documentation. */
+int
+timeval_subtract (result, x, y)
+   struct timeval *result, *x, *y;
+{
+   /* Perform the carry for the later subtraction by updating Y. */
+   if (x->tv_usec < y->tv_usec) {
+      int nsec = (y->tv_usec - x->tv_usec) / 1000000 + 1;
+      y->tv_usec -= 1000000 * nsec;
+      y->tv_sec += nsec;
+   }
+   if (x->tv_usec - y->tv_usec > 1000000) {
+      int nsec = (x->tv_usec - y->tv_usec) / 1000000;
+      y->tv_usec += 1000000 * nsec;
+      y->tv_sec -= nsec;
+   }
+
+   /* Compute the time remaining to wait.
+      `tv_usec' is certainly positive. */
+   result->tv_sec = x->tv_sec - y->tv_sec;
+   result->tv_usec = x->tv_usec - y->tv_usec;
+
+   /* Return 1 if result is negative. */
+   return x->tv_sec < y->tv_sec;
+}
+
+int
+main()
+{
+   printf("\n*** Checking many attributes in HDF5 file.\n");
+   printf("*** Checking some more simple atts...");
+   {
+#define NUM_ATTS 10000
+      hid_t fcpl_id, hdfid, grpid;
+      hid_t spaceid, attid1;
+      int one = 1;
+      hsize_t dims[1] = {1};
+      int i;
+      char name[NC_MAX_NAME];
+      struct timeval start_time, end_time, diff_time;
+      double sec;
+
+      /* Create a HDF5 file. */
+      if ((fcpl_id = H5Pcreate(H5P_FILE_CREATE)) < 0) ERR;
+      if (H5Pset_link_creation_order(fcpl_id, (H5P_CRT_ORDER_TRACKED |
+					       H5P_CRT_ORDER_INDEXED)) < 0) ERR;
+      if (H5Pset_attr_creation_order(fcpl_id, (H5P_CRT_ORDER_TRACKED |
+					       H5P_CRT_ORDER_INDEXED)) < 0) ERR;
+      if ((hdfid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, fcpl_id, H5P_DEFAULT)) < 0) ERR;
+      if (H5Pclose(fcpl_id) < 0) ERR;
+
+      /* Open the root group. */
+      if ((grpid = H5Gopen2(hdfid, "/", H5P_DEFAULT)) < 0) ERR;
+
+      /* Write an attribute. */
+      if ((spaceid = H5Screate_simple(1, dims, NULL)) < 0) ERR;
+      for (i = 0; i < NUM_ATTS; i++)
+      {
+	 if (gettimeofday(&start_time, NULL)) ERR;
+	 sprintf(name, "att_%d", i);
+	 if ((attid1 = H5Acreate2(grpid, name, H5T_NATIVE_INT, spaceid, 
+				  H5P_DEFAULT, H5P_DEFAULT)) < 0) ERR;
+	 if (H5Awrite(attid1, H5T_NATIVE_INT, &one) < 0) ERR;
+/*	 if (H5Aclose(attid1) < 0) ERR;*/
+	 if(i % 1000 == 0) 
+	 {		/* only print every 100th attribute name */
+	    if (gettimeofday(&end_time, NULL)) ERR;
+	    if (timeval_subtract(&diff_time, &end_time, &start_time)) ERR;
+	    sec = diff_time.tv_sec + 1.0e-6 * diff_time.tv_usec;
+	    printf("%i\t%.3g sec\n", i, sec);
+	 }
+      }
+
+      /* Close everything. */
+      if (H5Sclose(spaceid) < 0) ERR;
+      if (H5Gclose(grpid) < 0) ERR;
+      if (H5Fclose(hdfid) < 0) ERR;
+   }
+   SUMMARIZE_ERR;
+   FINAL_RESULTS;
+}