Program Listing for File debug.hpp

Return to documentation for file (src/debug.hpp)

/* This file is part of brille.

Copyright © 2019,2020 Greg Tucker <greg.tucker@stfc.ac.uk>

brille is free software: you can redistribute it and/or modify it under the
terms of the GNU Affero General Public License as published by the Free
Software Foundation, either version 3 of the License, or (at your option)
any later version.

brille is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.
See the GNU Affero General Public License for more details.

You should have received a copy of the GNU Affero General Public License
along with brille. If not, see <https://www.gnu.org/licenses/>.            */
#ifndef BRILLE_DEBUG_HPP_
#define BRILLE_DEBUG_HPP_
#include <stdio.h>
#include <sstream>
#include <iomanip>
#include <iostream>
#include <string>
#include <array>
#include <vector>
#include <complex>
#include <chrono>
// #define VERBOSE_DEBUG
// #define DEBUG // comment-out for no debugging output
namespace brille {

int terminal_width(void);
int terminal_height(void);

template <typename T> struct is_container {
  enum { value = false };
};
template <typename T> struct is_container<std::vector<T>> {
  enum { value = true };
};
template <typename T, size_t N> struct is_container<std::array<T, N>> {
  enum { value = true };
};

template<bool C, typename T> using enable_if_t = typename std::enable_if<C,T>::type;
template<typename T> static enable_if_t< std::is_integral<T>::value && std::is_unsigned<T>::value, T> local_abs(T x) { return x; }
template<typename T> static enable_if_t<!std::is_integral<T>::value ||   std::is_signed<T>::value, T> local_abs(T x) { return std::abs(x); }

template<typename T, typename=typename std::enable_if<!is_container<T>::value>::type>
const std::string my_to_string(const T x, const size_t width=0){
  std::ostringstream streamobj;
  size_t w{width};
  if (!std::is_integral<T>::value){
    streamobj << std::fixed;
    streamobj << std::setprecision(4);
    if (w>4) w -= 5u; // account for the decimal mark and four places
  }
  // char may or may not be signed, depending on the system
  if (std::is_base_of<char,T>::value || (std::is_integral<T>::value && std::is_unsigned<T>::value) ){
    if (w) streamobj << std::setw(w);
    streamobj << x;
  } else {
    if (w) streamobj << std::setw(w-1); // -1 to account for the sign
    streamobj << (x<0 ? "-" : " ") << local_abs(x);
  }
  return streamobj.str();
}
template<typename T, typename=typename std::enable_if<!is_container<T>::value>::type>
const std::string my_to_string(const std::complex<T> x, const size_t width=0){
  T r = std::real(x), i=std::imag(x);
  std::ostringstream streamobj;
  size_t w{width};
  if (!std::is_integral<T>::value){
    streamobj << std::fixed;
    streamobj << std::setprecision(4);
    if (w>9) w -= 10u; // account for the decimal mark and four places
  }
  if (!std::is_integral<T>::value || std::is_signed<T>::value){
    if (w>3) streamobj << std::setw(w-3); // -3 for -±i
    streamobj << (r<0 ? "-" : " ") << local_abs(r);
    streamobj << (std::signbit(i) ? "-i" : "+i") << local_abs(i);
  } else {
    if (w>2) streamobj << std::setw(w-2); // -2 for +i
    streamobj << r << "+i" << i;
  }
  return streamobj.str();
}
template<typename T, template<class> class C,
        typename=typename std::enable_if<!is_container<T>::value>::type,
        typename=typename std::enable_if<is_container<C<T>>::value>::type>
const std::string my_to_string(const std::vector<C<T>>& v, const size_t){
  std::string s;
  for (C<T> x: v) s += my_to_string(x) + "\n";
  return s;
}
template<typename T, typename=typename std::enable_if<is_container<T>::value>::type>
const std::string my_to_string(const T & a, const size_t w=0){
  std::string s;
  for (auto x: a) s += my_to_string(x, w);
  return s;
}

template<typename T> static enable_if_t<!is_container<T>::value, size_t> max_element_length(const T& v){
  return my_to_string(v).size();
}
template<typename T> static enable_if_t<is_container<T>::value, size_t> max_element_length(const T& v){
  size_t l=0;
  for (auto x: v){
    size_t t = max_element_length(x);
    if (t > l) l = t;
  }
  return l;
}

template<typename T>
std::string time_to_string(std::chrono::time_point<T> time) {
  using namespace std;
  using namespace std::chrono;

  time_t curr_time = T::to_time_t(time);
  char sRep[100];
  strftime(sRep, sizeof(sRep), "%Y-%m-%d %H:%M:%S", localtime(&curr_time));

  typename T::duration since_epoch = time.time_since_epoch();
  seconds s = duration_cast<seconds>(since_epoch);
  since_epoch -= s;
  milliseconds milli = duration_cast<milliseconds>(since_epoch);

  stringstream buffer;
  buffer << '[' << sRep << ':' << setfill('0') << setw(3) << milli.count() << "] ";
  return buffer.str();
}

class DebugPrinter{
  std::string last_function; // replace this with the stack?
  bool _silenced;
  bool emit_datetime;
  size_t _before;
public:
  DebugPrinter(const std::string& s)
  : last_function(s), _silenced(false), emit_datetime(false), _before(0)
  {
    #if defined(PROFILING)
    emit_datetime = true;
    #endif
  };
  bool silenced() const {return _silenced;}
  bool silenced(bool slc) {_silenced = slc; return _silenced;}
  bool silence() {_silenced = true; return _silenced;}
  bool unsilence() {_silenced = false; return !_silenced;}
  bool datetime() const {return emit_datetime;}
  bool datetime(bool edt) {emit_datetime = edt; return emit_datetime;}
  template<typename... L> void print(const std::string& fnc, L... l){
    if (!_silenced){
      size_t torem{0};
      if (last_function.compare(fnc)){
        last_function = fnc;
        std::cout << std::endl << fnc << std::endl;
      }
      if (emit_datetime){
        // auto tp = std::chrono::system_clock::now();
        // auto tt = std::chrono::system_clock::to_time_t(tp);
        // std::tm tm = *std::localtime(&tt);
        // std::stringstream buffer;
        // buffer << std::put_time(&tm, "%FT%T ");
        // std::string tstr = buffer.str();
        std::string tstr = time_to_string(std::chrono::system_clock::now());
        this->inner_print(tstr); // print without extra lead in spaces
        _before += (torem = tstr.size()); // add to lead in spacing
      }
      this->inner_print(l...);
      _before -= torem;
    }
  }
  template<typename... L> void println(const std::string& fnc, L... l){
    if (!_silenced){
      this->print(fnc, l...);
      std::cout << std::endl;
    }
  }
private:
  std::string lead_in() const {
    std::stringstream buffer;
    for (size_t i=0; i<_before; ++i) buffer << " ";
    return buffer.str();
  }
  template<typename T, typename... L>
  enable_if_t<!is_container<T>::value, void> inner_print(const T& x, L... l){
    std::cout << x;
    this->inner_print(l...);
  }
  template<typename T, typename... L>
  enable_if_t<!is_container<T>::value, void> inner_print(const std::vector<T>& x, L... args){
    size_t l = max_element_length(x);
    // size_t n = std::sqrt(x.size());
    int w = terminal_width();
    if (static_cast<int>(_before) < w) w -= static_cast<int>(_before);
    if (l) w /= static_cast<int>(l)+1;
    int count = 0;
    std::string s;
    // // if (x.size() == n*n){
    // if (3u == n){
    //   for (auto y: x){
    //     s += " " + my_to_string(y, l);
    //     ++count;
    //     if (!(count % w)||!(count % n)){
    //       s += "\n" + this->lead_in();
    //       count = 0;
    //     }
    //   }
    // } else {
      for (auto y: x){
        s += " " + my_to_string(y, l);
        if (!(++count % w)) s += "\n" + this->lead_in();
      }
    // }
    this->inner_print(s, args...);
  }
  template<typename T, size_t N, typename... L>
  enable_if_t<!is_container<T>::value, void> inner_print(const std::array<T,N>& x, L... args){
    size_t l= max_element_length(x);
    std::string s;
    // if (N==9){
    //   for (int a=0; a<3; ++a){
    //     for (int b=0; b<3; ++b) s += " " + my_to_string(x[a*3+b], l);
    //     s += "\n" + this->lead_in();
    //   }
    // } else {
      size_t w = static_cast<size_t>(terminal_width());
      if (_before < w) w -= _before;
      if (l) w /= l+1;
      size_t count = 0;
      for (size_t i=0; i<N; ++i){
        s += " " + my_to_string(x[i], l);
        if (!(++count % w)) s += "\n" + this->lead_in();
      }
    // }
    this->inner_print(s, args...);
  }
  template<typename T, typename... L>
  void inner_print(const std::vector<std::vector<T>>& vv, L... args){
    size_t l = max_element_length(vv);
    size_t w = static_cast<size_t>(terminal_width());
    if (_before < w) w -= _before;
    size_t num;
    std::string s;
    if (l) w /= l+1;
    for (auto v: vv){
      num = 0;
      for (auto x: v){
        s += my_to_string(x, l);
        if (!(++num %w)) s += "\n" + this->lead_in();
      }
      s += "\n" + this->lead_in();
    }
    this->inner_print(s, args...);
  }
  template<typename T, size_t N, typename... L>
  enable_if_t<!is_container<T>::value, void> inner_print(const std::vector<std::array<T,N>>& x, L... args){
    size_t l = max_element_length(x);
    size_t w = static_cast<size_t>(terminal_width());
    if (_before < w) w -= _before;
    size_t num;
    std::string s;
    // if (N==9){
    //   num = (l) ? w/(3*l+4) : w/3;
    //   for (size_t i=0; i<x.size(); i+=num){
    //     for (int a=0; a<3; ++a){
    //       for (size_t j=0; j<num && (i+j)<x.size(); ++j){
    //         for (int b=0; b<3; ++b) s += my_to_string(x[i+j][a*3+b], l);
    //         s += " ";
    //       }
    //       s += "\n" + this->lead_in();
    //     }
    //     s += "\n" + this->lead_in();
    //   }
    // } else {
      if (l) w /= l+1;
      for (size_t i=0; i<x.size(); num=0, ++i){
        for (auto y: x[i]){
          s += my_to_string(y, l);
          if (!(++num % w)) s += "\n" + this->lead_in();
        }
        s += "\n" + this->lead_in();
      }
    // }
    this->inner_print(s, args...);
  }
  void inner_print(void){};
};

static DebugPrinter printer("");

template<typename TimeT = std::chrono::milliseconds>
class Stopwatch{
  typedef std::chrono::high_resolution_clock ClockT;
private:
    std::chrono::time_point<ClockT> _start, _end, _split;
    size_t presses;
public:
    Stopwatch(): presses(0u){
      tic();
    }
    void tic(){
      presses = 0u;
      _start = _end = ClockT::now();
    }
    double toc(){
      _end = ClockT::now();
      ++presses;
      return elapsed();
    }
    double elapsed() const {
      auto delta = std::chrono::duration_cast<TimeT>(_end - _start);
      return static_cast<double>(delta.count());
    }
    double average() const {
      return elapsed()/static_cast<double>(presses);
    }
    double jitter() const {
      return std::sqrt(elapsed())/static_cast<double>(presses);
    }
    double split(){
      auto new_split = ClockT::now();
      auto delta = std::chrono::duration_cast<TimeT>(new_split - _split);
      _split = new_split;
      ++presses;
      return static_cast<double>(delta.count());
    }
};


} //end namespace brille

// replace GNU __PRETTY_FUNCTION__ by __FUNCSIG__ on Windows (MSVC)
#ifdef _MSC_VER
  #define __PRETTY_FUNCTION__ __FUNCSIG__
#endif
#ifdef VERBOSE_DEBUG
  #define _MY_PRETTY_FUNC_ __PRETTY_FUNCTION__
#else
  #define _MY_PRETTY_FUNC_ ""
#endif

#define info_update(...) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
#define info_update_if(tf, ...) if (tf) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)

#if defined(VERBOSE_DEBUG) || defined(DEBUG)
  #define debug_exec(...) __VA_ARGS__
  #define debug_update(...) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
  #define debug_update_if(tf, ...) if (tf) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
#else
  #define debug_update_if(...)
  #define debug_update(...)
  #define debug_exec(...)
#endif
#ifdef VERBOSE_DEBUG
  #define verbose_update(...) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
  #define verbose_update_if(tf, ...) if (tf) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
#else
  #define verbose_update(...)
  #define verbose_update_if(...)
#endif

#if defined(PROFILING)
  #define profile_update(...) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
  #define profile_update_if(tf, ...) if (tf) brille::printer.println(_MY_PRETTY_FUNC_, __VA_ARGS__)
#else
  #define profile_update(...)
  #define profile_update_if(...)
#endif

#endif //_DEBUG_H_