MathExtensions.h

Go to the documentation of this file.

// Copyright András Vukics 2006–2014. Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE.txt)
// -*- C++ -*-
#ifndef CPPQEDCORE_UTILS_MATHEXTENSIONS_H_INCLUDED
#define CPPQEDCORE_UTILS_MATHEXTENSIONS_H_INCLUDED

#include "MathExtensionsFwd.h"

#include "ComplexExtensions.h"
#include "Exception.h"


namespace mathutils {

struct FactOverflow  : public cpputils::Exception {};
  
extern const double PI    ;
extern const double SQRTPI;
extern const double EULER ;
  
int sign(double);
int fcmp(double, double, double);
  
template<class T> const T sqr(T x) {return x*x;}

double sqr(double x);

double sqrAbs(const dcomp&);

double fact  (unsigned          ) throw(FactOverflow);
double choose(unsigned, unsigned);

template<typename T>
bool parity(T n) { return n & 1;}
// even --- false, odd --- true;


inline int round(double r) {
  return (r>0) ? floor(r+0.5) : ceil(r-0.5);
}

template<typename T>
char minusOneToThePowerOf(T n) {return parity(n) ? -1 : 1;}



dcomp coherentElement(unsigned long n, const dcomp& alpha);

} // mathutils

#endif // CPPQEDCORE_UTILS_MATHEXTENSIONS_H_INCLUDED