coloring_mcmc.hpp

Go to the documentation of this file.

#ifndef _COLORING_MCMC_HPP
#define _COLORING_MCMC_HPP

#include "properties.hpp"
#include "random.hpp"
#include "coloring.hpp"
#include "arak.hpp"
#include "mcmc.hpp"

namespace Arak {

  class ColoringMove {

  public:

    virtual void execute(Coloring& c) = 0;
  
    virtual void undo(Coloring& c) = 0;
  
    virtual ~ColoringMove() {}
  
  };

  class NullColoringMove : public ColoringMove {

  public:

    NullColoringMove() { }
    void execute(Coloring& c) { }
    void undo(Coloring& c) { }
    ~NullColoringMove() {}

    static const NullColoringMove instance;

  };

  class ColoringProposal {

  public:

    virtual ~ColoringProposal() {}

    virtual void sample(const Coloring& state,
      Arak::Util::Random& random,
      bool reversible = true) = 0;

    virtual ColoringMove& move() = 0;

    virtual double ll(const Coloring& c) = 0;
  
    virtual double rll(const Coloring& c) = 0;

    virtual void result(bool accepted) = 0;
  };

  struct ArakMarkovChainTraits {
    typedef ArakProcess DistributionType;
    typedef Coloring StateType;
    typedef ColoringMove MoveType;
    typedef ColoringProposal ProposalType;
  };

  typedef MarkovChain<ArakMarkovChainTraits> ArakMarkovChain;

  typedef AnnealedMarkovChain<ArakMarkovChainTraits> AnnealedArakMarkovChain;

  typedef StochasticHillClimber<ArakMarkovChainTraits> HillClimbingArakMarkovChain;
}

#endif

---