source: lemon-0.x/src/work/akos/SimAnn.h @ 941:186aa53d2802

#ifndef SIMANN_H
#define SIMANN_H

#ifndef INFTY
#define INFTY 1e24
#endif

#include <iostream>

class SimAnnBase {
public:
  friend class Controller;
  class Controller {
  public:
    SimAnnBase *sab;
    double T;
    void init(SimAnnBase *_sab) { sab = _sab; }
    bool next() {}
    bool accept(double cost) {}
  };
protected:
  double CurrCost;
  double PrevCost;
  double BestCost;

  long Iter;
  long LastImpr;

  //friend class Controller;
  Controller *Ctrl;
public:
  long getIter() { return Iter; }
  long getLastImpr() { return LastImpr; }
  double getBestCost() { return BestCost; }
  virtual void mutate() = 0;
  virtual void revert() = 0;
  virtual void saveAsBest() = 0;
  virtual void init() = 0;

  void run() {
    init();
    while (Ctrl->next()) {
      Iter++;
      std::cout << Iter << std::endl;
      mutate();
      if (Ctrl->accept(PrevCost - CurrCost) && (CurrCost < BestCost)) {
        LastImpr = Iter;
        saveAsBest();
      }
      else {
        revert();
      }
    }
  }
};

template <typename E>
class SimAnn : public SimAnnBase {
private:
  E *CurrEnt;
  E *PrevEnt;
  E *BestEnt;
public:
  SimAnn() {
    CurrCost = PrevCost = BestCost = INFTY;
    Iter = LastImpr = 0;
  }
  void setController(Controller &_Ctrl) { Ctrl = &_Ctrl; }
  void setEnt(E &Ent) {
    CurrEnt = new E(Ent);
    PrevEnt = new E(Ent);
    BestEnt = new E(Ent);
  }
  E getBestEnt() { return *BestEnt; }
  void mutate() {
    *PrevEnt = *CurrEnt;
    PrevCost = CurrCost;
    CurrCost = CurrEnt->mutate();
  }
  void revert() {
    E *tmp = CurrEnt;
    CurrEnt = PrevEnt;
    PrevEnt = tmp;
    CurrCost = PrevCost;
  }
  void saveAsBest() {
    *BestEnt = *CurrEnt;
    BestCost = CurrCost;
  }
  void init() {
    Ctrl->init(this);
    CurrEnt->init();
  }
};

class EntitySkeleton {
public:
  void init() {}
  // returns the new cost
  double mutate() { return 0.0; }
};

#endif