EventPrimitivesToStringConverter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/

// EventPrimitivesToStringConverter.h, (c) ATLAS Detector software
 
#ifndef EVENTPRIMITIVESTOSTRINGCONVERTER_H_
#define EVENTPRIMITIVESTOSTRINGCONVERTER_H_
 
#include "EventPrimitives/EventPrimitives.h"
#include <iomanip>
#include <iostream>
#include <string>
#include <string_view>
#include <cmath>
#include <cassert>
 
namespace Amg {

 
inline double roundWithPrecision(double val, int precision) {
  assert(precision>=0);
  const double scale = std::pow(10.0, precision);
  const double rounded = std::round(val * scale) / scale;
  return rounded == 0.0 ? 0.0 : rounded; //avoids a -0 result.
}
 
inline std::string toString(const MatrixX& matrix, int precision = 4,
                            std::string_view offset = "") {
  std::ostringstream sout;
 
  sout << std::setiosflags(std::ios::fixed) << std::setprecision(precision);
  if (matrix.cols() == 1) {
    sout << "(";
    for (int i = 0; i < matrix.rows(); ++i) {
      double val = roundWithPrecision(matrix(i, 0), precision);
      sout << val;
      if (i != matrix.rows() - 1)
        sout << ", ";
    }
    sout << ")";
  } else {
    for (int i = 0; i < matrix.rows(); ++i) {
      for (int j = 0; j < matrix.cols(); ++j) {
        if (j == 0)
          sout << "(";
        double val = roundWithPrecision(matrix(i, j), precision);
        sout << val;
        if (j == matrix.cols() - 1)
          sout << ")";
        else
          sout << ", ";
      }
      if (i != matrix.rows() - 1) {  // make the end line and the offset in the next line
        sout << std::endl;
        sout << offset;
      }
    }
  }
  return sout.str();
}
 
}  // namespace Amg
 
#endif /* EVENTPRIMITIVESTOSTRINGCONVERTER_H_ */