app-bholger/gfem/sumofpowers.h

Go to the documentation of this file.

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

namespace concepts {


class SumOfPowers {
public:
  typedef std::vector<Cmplx> VC;
  typedef std::pair<Real, Real> PRR;
  typedef std::pair<PRR, int> CacheKey;
  typedef std::map<CacheKey, VC> Cache;

  SumOfPowers() {

  }

  ~SumOfPowers() {
  }

  static SumOfPowers& getInstance() {
    return g_instance;
  }

  Cmplx compute(Cmplx a, int k, int n); 

  Real binom(int a, int b);

  const VC& cached(Cmplx a, int kmax, int nrep);

  static inline Cmplx powi(Cmplx x, int power) {
    return BernoulliPol::powi(x, power);
  }

private:
  static SumOfPowers g_instance;

  Cache cache;
};

} // namespace