Primitives.cpp

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
// Author     : M.Frank
//
//==========================================================================
 
// Framework include files
#include <DD4hep/Primitives.h>
#include <DD4hep/Exceptions.h>
#include <DD4hep/Printout.h>
 
// C/C++ include files
#include <cstring>
 
#if defined(__linux) || defined(__APPLE__) || defined(__powerpc64__)
#include <cxxabi.h>
#ifndef __APPLE__
typedef abi::__class_type_info class_t;
using   abi::__dynamic_cast;
#endif
#endif
 
namespace {
 
#if 0
//-----------------------------------------------------------------------------
// MurmurHash2, 64-bit versions, by Austin Appleby
//
// The same caveats as 32-bit MurmurHash2 apply here - beware of alignment 
// and endian-ness issues if used across multiple platforms.
  inline uint64_t murmur_hash_64 ( const void * key, int len)  {
#define seed 0xFEEDBABE
    typedef unsigned long long int uint64;
 
#if INTPTR_MAX == INT32_MAX
    const unsigned int * data = (const unsigned int *)key;
    const unsigned int m = 0x5bd1e995;
    const int r = 24;
 
    unsigned int h1 = seed ^ len;
    unsigned int h2 = 0;
 
    while(len >= 8)
    {
      unsigned int k1 = *data++;
      k1 *= m; k1 ^= k1 >> r; k1 *= m;
      h1 *= m; h1 ^= k1;
      len -= 4;
 
      unsigned int k2 = *data++;
      k2 *= m; k2 ^= k2 >> r; k2 *= m;
      h2 *= m; h2 ^= k2;
      len -= 4;
    }
 
    if(len >= 4)
    {
      unsigned int k1 = *data++;
      k1 *= m; k1 ^= k1 >> r; k1 *= m;
      h1 *= m; h1 ^= k1;
      len -= 4;
    }
 
    switch(len)
    {
    case 3: h2 ^= ((unsigned char*)data)[2] << 16;
    case 2: h2 ^= ((unsigned char*)data)[1] << 8;
    case 1: h2 ^= ((unsigned char*)data)[0];
      h2 *= m;
    };
 
    h1 ^= h2 >> 18; h1 *= m;
    h2 ^= h1 >> 22; h2 *= m;
    h1 ^= h2 >> 17; h1 *= m;
    h2 ^= h1 >> 19; h2 *= m;
 
    uint64 h = h1;
 
    h = (h << 32) | h2;
#elif INTPTR_MAX == INT64_MAX
    const uint64* data = (const uint64*)key;
    const uint64 m = 0xc6a4a7935bd1e995;
    const int r = 47;
 
    uint64 h = seed ^ (len * m);
 
    const uint64 * end = data + (len/8);
 
    while(data != end)
    {
      uint64 k = *data++;
 
      k *= m; 
      k ^= k >> r; 
      k *= m; 
        
      h ^= k;
      h *= m; 
    }
 
    const unsigned char * data2 = (const unsigned char*)data;
 
    switch(len & 7)
    {
    case 7: h ^= uint64(data2[6]) << 48;
    case 6: h ^= uint64(data2[5]) << 40;
    case 5: h ^= uint64(data2[4]) << 32;
    case 4: h ^= uint64(data2[3]) << 24;
    case 3: h ^= uint64(data2[2]) << 16;
    case 2: h ^= uint64(data2[1]) << 8;
    case 1: h ^= uint64(data2[0]);
      h *= m;
    };
 
    h ^= h >> r;
    h *= m;
    h ^= h >> r;
 
#else
#error "Environment not 32 or 64-bit."
#endif
    return h;
  }
#endif
 
  struct FNV1a_64 {
    static const unsigned long long int hashinit = 14695981039346656037ull;
    static constexpr unsigned long long int doByte(unsigned long long int hash,unsigned char val)
    { return (hash ^ val) * 1099511628211ull; }
  };
 
  static unsigned char crc8_table[] =
    { 0, 94,188,226, 97, 63,221,131,194,156,126, 32,163,253, 31, 65,
      157,195, 33,127,252,162, 64, 30, 95,  1,227,189, 62, 96,130,220,
      35,125,159,193, 66, 28,254,160,225,191, 93,  3,128,222, 60, 98,
      190,224,  2, 92,223,129, 99, 61,124, 34,192,158, 29, 67,161,255,
      70, 24,250,164, 39,121,155,197,132,218, 56,102,229,187, 89,  7,
      219,133,103, 57,186,228,  6, 88, 25, 71,165,251,120, 38,196,154,
      101, 59,217,135,  4, 90,184,230,167,249, 27, 69,198,152,122, 36,
      248,166, 68, 26,153,199, 37,123, 58,100,134,216, 91,  5,231,185,
      140,210, 48,110,237,179, 81, 15, 78, 16,242,172, 47,113,147,205,
      17, 79,173,243,112, 46,204,146,211,141,111, 49,178,236, 14, 80,
      175,241, 19, 77,206,144,114, 44,109, 51,209,143, 12, 82,176,238,
      50,108,142,208, 83, 13,239,177,240,174, 76, 18,145,207, 45,115,
      202,148,118, 40,171,245, 23, 73,  8, 86,180,234,105, 55,213,139,
      87,  9,235,181, 54,104,138,212,149,203, 41,119,244,170, 72, 22,
      233,183, 85, 11,136,214, 52,106, 43,117,151,201, 74, 20,246,168,
      116, 42,200,150, 21, 75,169,247,182,232, 10, 84,215,137,107, 53
    };
  static unsigned char pearson_hash(const void *data, size_t len) {
    const char *s = (const char*)data;
    unsigned char c = 0;
    while (--len) c = crc8_table[c ^ *s++];
    return c;
  }
  static unsigned char pearson_hash(const char *data) {
    const char *s = data;
    unsigned char c = 0;
    while (*s) c = crc8_table[c ^ *s++];
    return c;
  }
}
 
std::string dd4hep::volumeID(VolumeID vid)   {
  char text[32];
  unsigned long long id = (unsigned long long)vid;
  ::snprintf(text,sizeof(text), "%016llx", id);
  return text;
}
 
unsigned long long int dd4hep::detail::hash64(const char* key)   {
  return update_hash64(FNV1a_64::hashinit, key);
}
 
unsigned long long int dd4hep::detail::hash64(const std::string& key)  {
  return update_hash64(FNV1a_64::hashinit, key.c_str(), key.length());
}
 
unsigned long long int dd4hep::detail::hash64(const void* key, std::size_t len)  {
  return update_hash64(FNV1a_64::hashinit, key, len);
}
 
unsigned long long int dd4hep::detail::update_hash64(unsigned long long int hash, const std::string& key)  {
  return update_hash64(hash, key.c_str(), key.length());
}
 
unsigned long long int dd4hep::detail::update_hash64(unsigned long long int hash, const char* key)  {
  const unsigned char* str = (const unsigned char*)key;
  for ( ; *str; ++str) hash = FNV1a_64::doByte(hash, *str);
  return hash;
}
 
unsigned long long int dd4hep::detail::update_hash64(unsigned long long int hash, const void* key, std::size_t len)  {
  const unsigned char* str = (const unsigned char*)key;
  if ( len > 0 )  {
    for ( ; --len; ++str) hash = FNV1a_64::doByte(hash, *str);
  }
  return hash;
}
 
unsigned short dd4hep::detail::hash16(const void* key, std::size_t len)   {
  unsigned short value = (unsigned short)hash32(key, len);
  return value;
}
 
unsigned char dd4hep::detail::hash8(const void* key, std::size_t len)   {
  return pearson_hash(key, len);
}
 
unsigned char dd4hep::detail::hash8(const char* key)   {
  return pearson_hash(key);
}
 
std::string dd4hep::detail::str_replace(const std::string& str,
                    const std::string& pattern,
                    const std::string& replacement)   {
  std::string res = str;
  for(size_t id=res.find(pattern); id != std::string::npos; id = res.find(pattern) )
    res.replace(id, pattern.length(), replacement);
  return res;
}
 
std::string dd4hep::detail::str_replace(const std::string& str,
                    char  pattern,
                    const std::string& replacement)   {
  std::string res = str;
  for(size_t id=res.find(pattern); id != std::string::npos; id = res.find(pattern) )
    res.replace(id, 1, replacement);
  return res;
}
 
std::string dd4hep::detail::str_replace(const std::string& str,
                    char  pattern,
                        char  replacement)   {
  std::string res = str;
  for(size_t id=res.find(pattern); id != std::string::npos; id = res.find(pattern) )
    res.replace(id, 1, 1, replacement);
  return res;
}
 
std::string dd4hep::detail::str_lower(const std::string& str) {
  std::string res = str.c_str();
  std::transform(res.begin(), res.end(), res.begin(), ::tolower);
  return res;
}
 
std::string dd4hep::detail::str_upper(const std::string& str) {
  std::string res = str.c_str();
  std::transform(res.begin(), res.end(), res.begin(), ::toupper);
  return res;
}
 
long int dd4hep::detail::makeTime(int year, int month, int day,
                          int hour, int minutes, int seconds)
{
  struct tm tm_init;
  ::memset(&tm_init,0,sizeof(tm_init));
  tm_init.tm_year  = year > 1900 ? year-1900 : year;
  tm_init.tm_mon   = month;
  tm_init.tm_mday  = day;
  tm_init.tm_hour  = hour;
  tm_init.tm_min   = minutes;
  tm_init.tm_sec   = seconds;
  tm_init.tm_isdst = -1;
  long int ti = ::mktime(&tm_init);
  if ( ti >= 0 ) return ti;
  except("dd4hep","Invalid time data given for conversion to epoch: %d-%d-%d %02d:%02d:%02d",
         year, month, day, hour, minutes, seconds);
  return ti;
}
 
long int dd4hep::detail::makeTime(const std::string& date, const char* fmt)  {
  struct tm tm;
  char* c = ::strptime(date.c_str(),fmt,&tm);
  if ( 0 == c )   {
    except("dd4hep",
           "Invalid time format given for update:%s should be: %s",
           date.c_str(), fmt);
  }
  long ti = ::mktime(&tm);
  if ( ti >= 0 ) return ti;
  except("dd4hep",
         "Invalid time string given for conversion to epoch: %s (fmt='%s')",
         date.c_str(), fmt);
  return ti;
}
 
static const std::string __typeinfoName(const std::type_info& tinfo) {
  const char* class_name = tinfo.name();
  std::string result;
#ifdef WIN32
  std::size_t off = 0;
  if ( ::strncmp(class_name, "class ", 6) == 0 ) {
    // The returned name is prefixed with "class "
    off = 6;
  }
  if ( ::strncmp(class_name, "struct ", 7) == 0 ) {
    // The returned name is prefixed with "struct "
    off = 7;
  }
  if ( off != std::string::npos ) {
    std::string tmp = class_name + off;
    size_t loc = 0;
    while( (loc = tmp.find("class ")) != std::string::npos ) {
      tmp.erase(loc, 6);
    }
    loc = 0;
    while( (loc = tmp.find("struct ")) != std::string::npos ) {
      tmp.erase(loc, 7);
    }
    result = tmp;
  }
  else {
    result = class_name;
  }
  // Change any " *" to "*"
  while ( (off=result.find(" *")) != std::string::npos ) {
    result.replace(off, 2, "*");
  }
  // Change any " &" to "&"
  while ( (off=result.find(" &")) != std::string::npos ) {
    result.replace(off, 2, "&");
  }
#elif defined(sun)
  result = class_name;
#elif !defined(__ICC)
  if (::strlen(class_name) == 1) {
    // See http://www.realitydiluted.com/mirrors/reality.sgi.com/dehnert_engr/cxx/abi.pdf
    // for details
    switch (class_name[0]) {
    case 'v':
      result = "void";
      break;
    case 'w':
      result = "wchar_t";
      break;
    case 'b':
      result = "bool";
      break;
    case 'c':
      result = "char";
      break;
    case 'h':
      result = "unsigned char";
      break;
    case 's':
      result = "short";
      break;
    case 't':
      result = "unsigned short";
      break;
    case 'i':
      result = "int";
      break;
    case 'j':
      result = "unsigned int";
      break;
    case 'l':
      result = "long";
      break;
    case 'm':
      result = "unsigned long";
      break;
    case 'x':
      result = "long long";
      break;
    case 'y':
      result = "unsigned long long";
      break;
    case 'n':
      result = "__int128";
      break;
    case 'o':
      result = "unsigned __int128";
      break;
    case 'f':
      result = "float";
      break;
    case 'd':
      result = "double";
      break;
    case 'e':
      result = "long double";
      break;
    case 'g':
      result = "__float128";
      break;
    case 'z':
      result = "ellipsis";
      break;
    }
  }
  else {
    char buff[16 * 1024];
    std::size_t len = sizeof(buff);
    int status = 0;
    result = __cxxabiv1::__cxa_demangle(class_name, buff, &len, &status);
  }
#else
  result = class_name;
  throw std::runtime_error("CXXABI is missing for ICC!");
#endif
  return result;
}
 
std::string dd4hep::typeName(const std::type_info& typ) {
  return __typeinfoName(typ);
}
 
void dd4hep::invalidHandleError(const std::type_info& type)
{
  throw invalid_handle_exception("Attempt to access invalid object of type "+typeName(type)+" [Invalid Handle]");
}
 
void dd4hep::invalidHandleAssignmentError(const std::type_info& from, 
                                          const std::type_info& to)
{
  std::string msg = "Wrong assingment from ";
  msg += typeName(from);
  msg += " to ";
  msg += typeName(to);
  msg += " not possible!!";
  throw invalid_handle_exception(msg);
}
 
void dd4hep::notImplemented(const std::string& msg)
{
  std::string m = "The requested feature " + msg + " is not implemented!";
  throw std::runtime_error(m);
}
 
void dd4hep::typeinfoCheck(const std::type_info& typ1, const std::type_info& typ2, const std::string& text)
{
  if (typ1 != typ2) {
    throw unrelated_type_error(typ1, typ2, text);
  }
}
 
namespace dd4hep   {
  namespace detail  {
    template<> const char* Primitive<bool>::default_format()           { return "%d"; }
    template<> const char* Primitive<char>::default_format()           { return "%c"; }
    template<> const char* Primitive<unsigned char>::default_format()  { return "%02X"; }
    template<> const char* Primitive<short>::default_format()          { return "%d"; }
    template<> const char* Primitive<unsigned short>::default_format() { return "%04X"; }
    template<> const char* Primitive<int>::default_format()            { return "%d"; }
    template<> const char* Primitive<unsigned int>::default_format()   { return "%08X"; }
    template<> const char* Primitive<long>::default_format()           { return "%ld"; }
    template<> const char* Primitive<unsigned long>::default_format()  { return "%016X"; }
    template<> const char* Primitive<float>::default_format()          { return "%f"; }
    template<> const char* Primitive<double>::default_format()         { return "%g"; }
    template<> const char* Primitive<char*>::default_format()          { return "%s"; }
    template<> const char* Primitive<const char*>::default_format()    { return "%s"; }
    template<> const char* Primitive<std::string>::default_format()    { return "%s"; }
 
    template <typename T> std::string Primitive<T>::toString(T value) {
      char text[1024];
      ::snprintf(text,sizeof(text),default_format(),value);
      return text;
    }
 
    template <> std::string Primitive<const char*>::toString(const char* value) {
      if ( value )  {
        return value;
      }
      throw std::runtime_error("Failed to convert (char*)NULL to std::string!");
    }
    template <> std::string Primitive<char*>::toString(char* value) {
      if ( value )  {
        return value;
      }
      throw std::runtime_error("Failed to convert (char*)NULL to std::string!");
    }
    template <> std::string Primitive<std::string>::toString(std::string value) {
      return value;
    }
 
    template std::string Primitive<bool>::toString(bool value);
    template std::string Primitive<char>::toString(char value);
    template std::string Primitive<unsigned char>::toString(unsigned char value);
    template std::string Primitive<short>::toString(short value);
    template std::string Primitive<unsigned short>::toString(unsigned short value);
    template std::string Primitive<int>::toString(int value);
    template std::string Primitive<unsigned int>::toString(unsigned int value);
    template std::string Primitive<long>::toString(long value);
    template std::string Primitive<unsigned long>::toString(unsigned long value);
    template std::string Primitive<float>::toString(float value);
    template std::string Primitive<double>::toString(double value);
  }
}
 
#ifdef __APPLE__
dd4hep::Cast::Cast(const std::type_info& t, cast_t c) : type(t), cast(c)  {
}
#else
dd4hep::Cast::Cast(const std::type_info& t) : type(t)   {
  abi_class = dynamic_cast<const class_t*>(&type);
  if (!abi_class) {
    throw std::runtime_error("Class type " + typeName(type) + " is not an abi object type!");
  }
}
#endif
 
dd4hep::Cast::~Cast() {
}
 
void* dd4hep::Cast::apply_dynCast(const Cast& to, const void* ptr) const
{
  if (&to == this) {
    return (void*) ptr;
  }
#ifdef __APPLE__
  // First try down cast
  void *r = (*to.cast)(ptr);
  if (r)
    return r;
  // Now try the up-cast
  r = (*cast)(ptr);
  if (r)      return r;
  throw unrelated_type_error(type, to.type, "Failed to apply abi dynamic cast operation!");
#else
  void* r = (void*)ptr;
  if ( to.abi_class )  {
    bool cast_worked = type.__do_upcast((const class_t*)to.abi_class,&r);
    if ( cast_worked ) return r;
    r = (void*)ptr;
    cast_worked = to.type.__do_upcast((const class_t*)abi_class,&r);
    if ( cast_worked ) return r;
#if 0
    const class_t* src_type = (const class_t*)to.abi_class;
    if (src_type) {
      // First try down cast
      void *r = cast_wrap(ptr, src_type, (const class_t*) abi_class, -1);
      if ( r ) return r;
      // Now try the up-cast
      r = cast_wrap(ptr, (const class_t*) abi_class, src_type, -1);
      if (r)      return r;
    }
#endif
    throw unrelated_type_error(type, to.type, "Failed to apply abi dynamic cast operation!");
  }
#endif
  throw unrelated_type_error(type, to.type, "Target type is not an abi class type!");
}
 
void* dd4hep::Cast::apply_upCast(const Cast& to, const void* ptr) const
{
  if (&to == this) {
    return (void*) ptr;
  }
  return apply_dynCast(to, ptr);
}
  
void* dd4hep::Cast::apply_downCast(const Cast& to, const void* ptr) const
{
  if (&to == this) {
    return (void*) ptr;
  }
#ifdef __APPLE__
  void *r = (*to.cast)(ptr);
  if (r) return r;
  throw unrelated_type_error(type, to.type, "Failed to apply abi dynamic cast operation!");
#else
  if ( to.abi_class )  {
    // Since we have to cast a 'to' pointer up to the real pointer
    // no virtual inheritance can be supported!
    void* r = (void*)ptr;
    bool cast_worked = type.__do_upcast((const class_t*)to.abi_class,&r);
    if ( cast_worked ) return r;
#if 0
    void *r = cast_wrap(ptr, src_type, (const class_t*)abi_class, -1);
    if (r) return r;
#endif
    throw unrelated_type_error(type, to.type, "Failed to apply abi dynamic cast operation!");
  }
  throw unrelated_type_error(type, to.type, "Target type is not an abi class type!");
#endif
}