DoubleAggregator.hpp

Go to the documentation of this file.

#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics.hpp>
#include <boost/accumulators/statistics/sum_kahan.hpp>
 
namespace precice::utils {
 
class DoubleAggregator {
private:
  using Acc = boost::accumulators::accumulator_set<double, boost::accumulators::stats<boost::accumulators::tag::sum_kahan>>;
 
public:
  void operator=(double d)
  {
    acc = Acc{};
    acc(d);
  }
 
  void add(double d)
  {
    acc(d);
  }
 
  void operator+=(double d)
  {
    add(d);
  }
 
  void operator-=(double d)
  {
    add(-d);
  }
 
  double value() const
  {
    return boost::accumulators::sum_kahan(acc);
  }
 
  operator double() const
  {
    return value();
  }
 
private:
  Acc acc;
};
 
[[nodiscard]] inline DoubleAggregator operator+(double lhs, DoubleAggregator rhs)
{
  rhs.add(lhs);
  return rhs;
}
 
[[nodiscard]] inline DoubleAggregator operator+(DoubleAggregator lhs, double rhs)
{
  lhs.add(rhs);
  return lhs;
}
 
[[nodiscard]] inline DoubleAggregator operator-(DoubleAggregator lhs, double rhs)
{
  lhs.add(-rhs);
  return lhs;
}
 
} // namespace precice::utils