TransitionKernel.cpp

Go to the documentation of this file.

#include "MUQ/SamplingAlgorithms/TransitionKernel.h"
 
namespace pt = boost::property_tree;
using namespace muq::Modeling;
using namespace muq::SamplingAlgorithms;
 
TransitionKernel::TransitionKernel(pt::ptree                 const& pt,
                                   std::shared_ptr<AbstractSamplingProblem> problem) : blockInd(pt.get("BlockIndex",0)),
                                                                                       problem(problem),
                                                                                       reeval(pt.get<bool>("ReevaluateAcceptedDensity", false)) {}
 
 
std::shared_ptr<TransitionKernel::TransitionKernelMap> TransitionKernel::GetTransitionKernelMap() {
  // define a static map from type to constructor
  static std::shared_ptr<TransitionKernelMap> map;
 
  if( !map ) { // if the map has not yet been created ...
    // ... create the map
    map = std::make_shared<TransitionKernelMap>();
  }
 
  return map;
}
 
std::shared_ptr<TransitionKernel> TransitionKernel::Construct(pt::ptree const& pt, std::shared_ptr<AbstractSamplingProblem> problem) {
 
  // get the name of the kernel
  const std::string& kernelName = pt.get<std::string>("Method");
 
  // construct it from the map
  auto kernelMap = GetTransitionKernelMap();
  auto iter = kernelMap->find(kernelName);
  if(iter == kernelMap->end()){
    std::cerr << "ERROR: Could not find Transition Kernel \"" << kernelName << "\", available types are:\n";
 
    for(auto it=kernelMap->begin(); it!=kernelMap->end(); ++it)
      std::cerr << "  " << it->first << std::endl;
    std::cerr << std::endl;
 
    assert(iter != kernelMap->end());
 
  }
 
  return iter->second(pt, problem);
}
 
#if MUQ_HAS_PARCER
void TransitionKernel::SetCommunicator(std::shared_ptr<parcer::Communicator> newcomm) {
  comm = newcomm;
}
#endif