ndistributions.h

Go to the documentation of this file.

//
// ePDF
//
 
#pragma once
 
#include "ePDF/alphaem.h"
#include "ePDF/matrix.h"
 
#include <yaml-cpp/yaml.h>
#include <complex>
#include <vector>
 
namespace ePDF
{
  class NDistributions
  {
  public:
    NDistributions(YAML::Node const& config);
 
    virtual ~NDistributions() {};
 
    virtual std::vector<std::complex<double>> Evolve(std::complex<double> const& N, double const& Q) const;
 
    std::pair<std::complex<double>, Matrix<std::complex<double>>> UMatrix(std::complex<double> const& N, double const& Q) const;
 
    std::pair<std::complex<double>, Matrix<std::complex<double>>> PathOrdering(std::complex<double> const& N, double const& Q) const;
 
    std::pair<std::complex<double>, Matrix<std::complex<double>>> AlphaFixed(std::complex<double> const& N, double const& Q) const;
 
    std::pair<std::complex<double>, Matrix<std::complex<double>>> InitialDistributions(std::complex<double> const& N) const;
 
    std::complex<double> fde(std::complex<double> const& N) const;
    std::complex<double> fdgm(std::complex<double> const& N) const;
 
    Matrix<std::complex<double>> Integrandsg(Matrix<std::complex<double>> const& G0,
                                             Matrix<std::complex<double>> const& G1,
                                             Matrix<std::complex<double>> const& J,
                                             double const& a) const;
    std::complex<double> Integrandns(std::complex<double> const& G0,
                                     std::complex<double> const& G1,
                                     std::complex<double> const& J,
                                     double const& a) const;
    Matrix<std::complex<double>> Minv(Matrix<std::complex<double>> const& M) const;
 
  private:
    AlphaQED    const _aQED;
    int         const _ipt;
    int         const _nl;
    double      const _Qi;
    double      const _ai;
    double      const _bt0;
    double      const _bt1;
    std::string const _method;
    bool        const _expand;
    std::string const _scheme;
  };
}