Commit 4a7bc6db authored by vbertone's avatar vbertone
Browse files

Moved initialisation of alpha_s and x-space of APFEL for GPD evolution from...

Moved initialisation of alpha_s and x-space of APFEL for GPD evolution from prepareSubModules into initModule
parent 1765db9b
......@@ -97,17 +97,6 @@ void GPDEvolutionApfel::configure(const ElemUtils::Parameters &parameters) {
void GPDEvolutionApfel::prepareSubModules(const std::map<std::string, BaseObjectData>& subModulesData) {
GPDEvolutionModule::prepareSubModules(subModulesData);
apfel::SetVerbosityLevel(0);
// Setup APFEL++ x-space
std::vector<apfel::SubGrid> vsg;
for (int i = 0; i < (int) m_subgridNodes.size(); i++)
vsg.push_back(apfel::SubGrid{m_subgridNodes[i], m_subgridLowerBounds[i], m_subgridInterDegrees[i]});
m_g = std::unique_ptr<apfel::Grid> (new apfel::Grid{vsg});
// Running coupling
m_as = [=] (double const& mu) -> double{ return getRunningAlphaStrongModule()->compute(mu * mu); };
}
PartonDistribution GPDEvolutionApfel::compute(GPDModule* pGPDModule, const GPDType::Type &type) {
......@@ -277,6 +266,17 @@ GPDEvolutionApfel::GPDEvolutionApfel(const GPDEvolutionApfel &other) :
void GPDEvolutionApfel::initModule() {
GPDEvolutionModule::initModule();
apfel::SetVerbosityLevel(0);
// Setup APFEL++ x-space
std::vector<apfel::SubGrid> vsg;
for (int i = 0; i < (int) m_subgridNodes.size(); i++)
vsg.push_back(apfel::SubGrid{m_subgridNodes[i], m_subgridLowerBounds[i], m_subgridInterDegrees[i]});
m_g = std::unique_ptr<apfel::Grid> (new apfel::Grid{vsg});
// Running coupling
m_as = [=] (double const& mu) -> double{ return getRunningAlphaStrongModule()->compute(mu * mu); };
}
void GPDEvolutionApfel::isModuleWellConfigured() {
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment