RGEvolutorNew.cpp

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#include "RGEvolutorNew.h"
 
RGEvolutorNew::RGEvolutorNew(unsigned int dim, schemes scheme, qcd_orders order_qcd_i, qed_orders order_qed_i)
: WilsonTemplateNew<gslpp::matrix<double> >(dim, scheme, order_qcd_i, order_qed_i)
{}
 
 
void RGEvolutorNew::setEvol(unsigned int i, unsigned int  j, double x, qcd_orders order_qcd_i, qed_orders order_qed_i) 
{    
    if (i > size || j > size) {
        std::stringstream out;
        out << i << " " << j;
        throw std::runtime_error("RGEvolutorNew::setEvol(): matrix indices " + out.str() + " out of range"); 
    }
    if (order_qcd_i > order_qcd || order_qed_i > order_qed) {
        std::stringstream out;
        out << order_qcd_i << " and " << order_qed_i;
        throw std::runtime_error("RGEvolutorNew::setEvol(): order " + out.str() +" not implemented "); 
    }
    gslpp::matrix<double> tmp = wilson.getOrd(order_qcd_i, order_qed_i);
    tmp.assign(i, j, x);
    wilson.setOrd(order_qcd_i, order_qed_i, tmp);
}
 
void RGEvolutorNew::setEvol(const gslpp::matrix<double>& m, qcd_orders order_qcd_i, qed_orders order_qed_i)
{
    setWilson(m, order_qcd_i, order_qed_i);
}
 
const Expanded<gslpp::matrix<double> >& RGEvolutorNew::getEvol() const
{
    return wilson;
}
 
double RGEvolutorNew::getM() const
{
    return M;
}
 
void RGEvolutorNew::setScales(double mu, double M)
{
    this->M = M;
    WilsonTemplateNew::setMu(mu);
    setEvol(gslpp::matrix<double>::Id(size), QCD0, QED0);
}
 
void RGEvolutorNew::setM(double M)
{
    setScales(mu, M);
}
 
void RGEvolutorNew::setMu(double mu)
{
    setScales(mu, M);
}
 
const gslpp::matrix<double>& RGEvolutorNew::Evol(qcd_orders order_qcd_i, qed_orders order_qed_i) const
{
    return getWilson(order_qcd_i, order_qed_i);
}