moILS.h

// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-

// "moILS.h"

// (c) OPAC Team, LIFL, 2003-2007

/* LICENCE TEXT
   
   Contact: paradiseo-help@lists.gforge.inria.fr
*/

#ifndef __moILS_h
#define __moILS_h

#include <eoEvalFunc.h>

#include "moHC.h"
#include "moTS.h"
#include "moSA.h"


template < class M > class moILS:public moAlgo < typename M::EOType >
{

  typedef typename M::EOType EOT;
  
  typedef typename EOT::Fitness Fitness;

public:


  moILS (moAlgo<EOT> &__algo, moSolContinue <EOT> &__continue, moComparator<EOT> &__acceptance_criterion, eoMonOp<EOT> &__perturbation, 
         eoEvalFunc<EOT> &__full_eval):
    algo(__algo), cont(__continue), acceptance_criterion(__acceptance_criterion), perturbation(__perturbation), full_eval(__full_eval)
  {}


  moILS (moMoveInit < M > &__move_init, moNextMove < M > &__next_move, moMoveIncrEval < M > &__incr_eval,
         moMoveSelect < M > &__move_select, moSolContinue <EOT> &__continue, moComparator<EOT> &__acceptance_criterion, 
         eoMonOp<EOT> &__perturbation, eoEvalFunc<EOT> &__full_eval):
    algo(*new moHC<M>(__move_init, __next_move, __incr_eval, __move_select, __full_eval)), cont(__continue), 
    acceptance_criterion(__acceptance_criterion), perturbation(__perturbation), full_eval(__full_eval)
  {}


  moILS (moMoveInit <M> &__move_init, moNextMove <M> &__next_move, moMoveIncrEval <M> &__incr_eval,
         moTabuList <M> &__tabu_list, moAspirCrit <M> &__aspir_crit, moSolContinue <EOT> &__moTS_continue, 
         moSolContinue <EOT> &__continue, moComparator<EOT> &__acceptance_criterion, eoMonOp<EOT> &__perturbation, 
         eoEvalFunc<EOT> &__full_eval):
    algo(*new moTS<M>(__move_init, __next_move, __incr_eval, __tabu_list, __aspir_crit, __moTS_continue, __full_eval)), 
    cont(__continue), acceptance_criterion(__acceptance_criterion), perturbation(__perturbation), full_eval(__full_eval)
  {}


  moILS (moRandMove<M> &__move_rand, moMoveIncrEval <M> &__incr_eval, moSolContinue <EOT> &__moSA_continue, double __init_temp,
         moCoolingSchedule & __cool_sched, moSolContinue <EOT> &__continue, moComparator<EOT> &__acceptance_criterion, 
         eoMonOp<EOT> &__perturbation, eoEvalFunc<EOT> &__full_eval):
    algo(*new moSA<M>(__move_rand, __incr_eval, __moSA_continue, init_temp, __cool_sched, __full_eval)), 
    cont(__continue), acceptance_criterion(__acceptance_criterion), perturbation(__perturbation), full_eval(__full_eval)
  {}




  bool operator()(EOT & __sol)
  {
    EOT __sol_saved=__sol;
            
    cont.init ();

    //some code has been duplicated in order to avoid one perturbation and one evaluation without adding a test in the loop.

    //std::cout << "Before intensification: " << __sol.fitness() << std::endl; 
    algo(__sol);
    //std::cout << "After intensification: " << __sol.fitness() << std::endl; 
    
    if(acceptance_criterion(__sol, __sol_saved))
      {
        std::cout << "Accepted: " << __sol.fitness() << std::endl;
        __sol_saved=__sol;
        
      }
    else
      {
        //std::cout << "Refused" << std::endl;
        __sol=__sol_saved;
      }
    
    while (cont (__sol))
    {
      perturbation(__sol);
      full_eval(__sol);
      
      //std::cout << "Before intensification: " << __sol.fitness() << std::endl; 
      algo(__sol);
      //std::cout << "After intensification: " << __sol.fitness() << std::endl; 
      
      if(acceptance_criterion(__sol, __sol_saved))
        {
          std::cout << "Accepted: " << __sol.fitness() << std::endl;
          __sol_saved=__sol;
        }
      else
        {
          //std::cout << "Refused" << std::endl;
          __sol=__sol_saved;
        }
    }
    
    return true;
  }

private:

  moAlgo<EOT> &algo;

  moSolContinue<EOT> &cont;

  moComparator<EOT> &acceptance_criterion;

  eoMonOp<EOT> &perturbation;
  
  eoEvalFunc<EOT> &full_eval;
};

#endif

---