#include "cs_config.h"
#include <unistd.h>
#include "util/neo_misc.h"
#include "util/neo_hdf.h"
#include "util/neo_rand.h"
#include "util/neo_files.h"

int main(int argc, char *argv[])
{
  NEOERR *err;
  HDF *hdf;
  int x;
  double tstart = 0;
  double tend = 0;
  char *file;
  int reps = 1000;
  char *s = NULL;

  if (argc > 1)
    file = argv[1];
  else
    file = "test.hdf";

  if (argc > 2)
    reps = atoi(argv[2]);

  err = hdf_init(&hdf);
  if (err != STATUS_OK) 
  {
    nerr_log_error(err);
    return -1;
  }

  tstart = ne_timef();

  for (x = 0; x < reps; x++)
  {
    err = hdf_read_file(hdf, file);
    if (err != STATUS_OK) 
    {
      nerr_log_error(err);
      return -1;
    }
  }
  tend = ne_timef();
  ne_warn("hdf_read_file test finished in %5.3fs, %5.3fs/rep", tend - tstart, (tend-tstart) / reps);

  tstart = ne_timef();

  for (x = 0; x < reps; x++)
  {
    err = ne_load_file(file, &s);
    if (err != STATUS_OK) 
    {
      nerr_log_error(err);
      return -1;
    }
    err = hdf_read_string(hdf, s);
    if (err != STATUS_OK) 
    {
      nerr_log_error(err);
      return -1;
    }
    free(s);
  }
  tend = ne_timef();
  ne_warn("load/hdf_read_string test finished in %5.3fs, %5.3fs/rep", tend - tstart, (tend-tstart) / reps);


  err = ne_load_file(file, &s);
  if (err != STATUS_OK) 
  {
    nerr_log_error(err);
    return -1;
  }
  tstart = ne_timef();
  for (x = 0; x < reps; x++)
  {
    err = hdf_read_string(hdf, s);
    if (err != STATUS_OK) 
    {
      nerr_log_error(err);
      return -1;
    }
  }
  tend = ne_timef();
  free(s);
  ne_warn("hdf_read_string test finished in %5.3fs, %5.3fs/rep", tend - tstart, (tend-tstart) / reps);
  /* hdf_dump(hdf, NULL);  */

  hdf_destroy(&hdf);

  return 0;
}