// -*- C++ -*- // // This file is part of LHAPDF // Copyright (C) 2012-2021 The LHAPDF collaboration (see AUTHORS for details) // #include "LHAPDF/KnotArray.h" #include namespace LHAPDF { void KnotArray::initPidLookup(){ _lookup.clear(); if (_pids.size() == 0){ // Is there a better LHAPDF error for that? std::cerr << "Internal error when constructing lookup table; need to fill pids before construction"<< std::endl; throw; } for (int i(-6); i<0; i++) _lookup.push_back(findPidInPids(i, _pids)); _lookup.push_back(findPidInPids(21, _pids)); for (int i(1); i<=6; i++) _lookup.push_back(findPidInPids(i, _pids)); _lookup.push_back(findPidInPids(22, _pids)); } void KnotArray::fillLogKnots() { _logxs.resize(_xs.size()); for (size_t i(0); i<_xs.size(); ++i) _logxs[i] = log(_xs[i]); _logq2s.resize(_q2s.size()); for (size_t i(0); i<_q2s.size(); ++i) _logq2s[i] = log(_q2s[i]); } }