operator.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/expression.h"
#include "apfel/distribution.h"
#include "apfel/matrix.h"
 
namespace apfel
{
  class Operator
  {
  public:
    Operator() = delete;
    Operator(Operator const&) = default;
 
    Operator(Grid const& gr, Expression const& expr, double const& eps = 1e-5, bool const& gpd = false);
 
    virtual ~Operator() {}
 
    Distribution operator *= (Distribution const& d) const;         
    Operator& operator *= (Operator const& o);                      
    Operator& operator  = (Operator const& o);                      
    Operator& operator *= (double const& s);                        
    Operator& operator *= (std::function<double(double const&)> f); 
    Operator& operator /= (double const& s);                        
    Operator& operator += (Operator const& o);                      
    Operator& operator -= (Operator const& o);                      
 
    void BuildOperatorDGLAP();
 
    void BuildOperatorGPD();
 
    Distribution Evaluate(double const& x) const;
 
    Grid const& GetGrid() const { return _grid; }
 
    Expression const& GetExpression() const { return _expr; }
 
    double const& GetIntegrationAccuracy() const { return _eps; }
 
    bool const& IsGPD() const { return _gpd; }
 
    std::vector<matrix<double>> GetOperator() const { return _Operator; }
 
    void Print() const { std::cout << *this << std::endl; }
 
  protected:
    Grid                 const& _grid;      
    Expression           const& _expr;      
    double               const  _eps;       
    bool                 const  _gpd;       
    std::vector<matrix<double>> _Operator;  
 
    friend std::ostream& operator << (std::ostream& os, Operator const& sg);
  };
 
  Distribution operator * (Operator lhs, Distribution const& rhs);                
  Operator     operator * (Operator lhs, Operator const& rhs);                    
  Operator     operator * (double const& s, Operator rhs);                        
  Operator     operator * (Operator lhs, double const& s);                        
  Operator     operator * (std::function<double(double const&)> f, Operator rhs); 
  Operator     operator * (Operator lhs, std::function<double(double const&)> f); 
  Operator     operator / (Operator lhs, double const& s);                        
  Operator     operator + (Operator lhs, Operator const& rhs);                    
  Operator     operator - (Operator lhs, Operator const& rhs);                    
 
  std::ostream& operator << (std::ostream& os, Operator const& op);
}