qcdloop.h

//
// QCDLoop 2016
//
// Authors: Stefano Carrazza: stefano.carrazza@cern.ch
//          Keith Ellis: keith.ellis@durham.ac.uk
//          Giulia Zanderighi: giulia.zanderighi@cern.ch

#pragma once

#include <vector>
#include <memory>
#include "types.h"
#include "tadpole.h"
#include "bubble.h"
#include "triangle.h"
#include "box.h"
#include "maths.h"
#include "tools.h"
#include "timer.h"
using std::vector;

namespace ql
{
  template<typename TOutput = complex, typename TMass = double, typename TScale = double>
  class QCDLoop
  {
  public:
    QCDLoop();  
    ~QCDLoop(); 

    void integral(vector<TOutput> &res,
                  TScale const& mu2,
                  vector<TMass> const& m,
                  vector<TScale> const& p = {}) const;

    void setCacheSize(int const& size);

  private:
    TadPole<TOutput,TMass,TScale>*  _tp;
    Bubble<TOutput,TMass,TScale>*   _bb;
    Triangle<TOutput,TMass,TScale>* _tr;
    Box<TOutput,TMass,TScale>*      _tb;
  };
}