WriterCommon.cxx

Go to the documentation of this file.

// this is -*- C++ -*-
/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
#include "HDF5Utils/Writer.h"
 
#include "H5Cpp.h"
 
#include <cassert>
#include <iostream>             // for printing errors in the destructor
 
 
namespace H5Utils {
  namespace internal {
  // packing utility
    H5::CompType packed(H5::CompType in) {
      // TODO: Figure out why a normal copy constructor doesn't work here.
      //       The normal one seems to create shallow copies.
      H5::CompType out(H5Tcopy(in.getId()));
      out.pack();
      return out;
    }
 
    void printDestructorError(const std::string& msg) {
      std::cerr << "ERROR: an exception was thrown in the destructor of an "
        "HDF5 file, the output buffer may be corrupted";
      std::cerr << " (error message: " << msg << ")" << std::endl;
    }
 
    // new functions
    H5::DataSpace getUnlimitedSpace(const std::vector<hsize_t>& extent) {
      std::vector<hsize_t> initial{0};
      initial.insert(initial.end(), extent.begin(), extent.end());
      std::vector<hsize_t> eventual{H5S_UNLIMITED};
      eventual.insert(eventual.end(), extent.begin(), extent.end());
      return H5::DataSpace(eventual.size(), initial.data(), eventual.data());
    }
    std::vector<hsize_t> getStriding(std::vector<hsize_t> extent) {
      // calculate striding
      extent.push_back(1);
      for (size_t iii = extent.size(); iii - 1 != 0; iii--) {
        extent.at(iii-2) = extent.at(iii-2) * extent.at(iii-1);
      }
      return extent;
    }
    void throwIfExists(const std::string& name, const H5::Group& in_group) {
      if (H5Lexists(in_group.getLocId(), name.c_str(), H5P_DEFAULT)) {
        throw std::logic_error("tried to overwrite '" + name + "'");
      }
    }
 
 
 
  }
 
}