app-bholger/gfem/bernoulliPol.h

Go to the documentation of this file.

#pragma once
#include "basics/typedefs.hh"
#include <map>
#include <vector>


namespace concepts {


class BernoulliPol {
private:
  // you probably don't want to create a copy of this class
  BernoulliPol& operator=(BernoulliPol& other);
  BernoulliPol(BernoulliPol& other);
public:
  //typedef std::vector<Cmplx> VC;
  //typedef std::pair<Real, Real> PRR;
  typedef std::pair<int, Real> CacheKey;
  typedef std::map<CacheKey, Real> Cache;

  BernoulliPol() {

  }

  ~BernoulliPol() {
  }

  static BernoulliPol& getInstance() {
    return g_instance;
  }

  Real eval(int k, Real x); 

  Real binom(int a, int b);

  template<typename Real>
  static inline Real powi(Real x, int power) {
#if 0
    Real res(1);
    for(int i = 0; i < power; ++i) {
      res *= x;
    }
    return res;
#else
    return pow(x, power);
#endif
  }

private:
  static BernoulliPol g_instance;

  Cache cache;
};

} // namespace