Maths.h

// -*- C++ -*-
//
// Maths.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_Math_H
#define HERWIG_Math_H
 
#include <cmath>
#include <vector>
#include "ThePEG/Config/Complex.h"
#include "ThePEG/Utilities/Maths.h"
#include "ThePEG/Config/TemplateTools.h"
 
namespace Herwig {
using ThePEG::Complex;
using std::complex;
 
namespace Math {
 
  Complex Li2(Complex);
 
  long double ReLi2(long double);
  
  inline double angleZeroTo2Pi(double angle) {
    double ret = fmod(angle, 2 * M_PI);
    if (ret < 0) ret += 2 * M_PI;
    return ret;
  }
 
  inline double angleMinusPiToPi(double angle) {
    double ret = angleZeroTo2Pi(angle);
    if (ret > M_PI) ret -= 2 * M_PI;
    return ret;
  }
 
  template <typename T>
  inline T median(std::vector<T> v) {
    if (v.empty()) return T();
    sort ( v.begin(), v.end() );
    const size_t N = v.size();
    return (N % 2) ? v.at((N+1)/2 - 1) : v.at(N/2 - 1);
  }
 
}
 
}
 
#endif /* HERWIG_Math_H */