/*

(c) Thales group, 2010

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation;
    version 2 of the License.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

Contact: http://eodev.sourceforge.net

Authors:
Lionel Parreaux <lionel.parreaux@gmail.com>

*/

#include <mo>
#include <eo>
#include <es.h>
#include <edo>

#include <edoBounderUniform.h>

#include "neighborhood/moRealNeighbor.h"
#include "neighborhood/moRealNeighborhood.h"

#include "sampling/moStdDevEstimator.h"

#include "coolingSchedule/moTrikiCoolingSchedule.h"


typedef eoReal< eoMinimizingFitness > EOT;
typedef moRealNeighbor< EOT > Neighbor;

double objective_function(const EOT & sol)
{
    double sum = 0;

    for ( size_t i = 0; i < sol.size(); ++i )
    {
        sum += sol[i] * sol[i];
    }

    return sum;
}

int main( int ac, char** av )
{
    eoParser parser( ac, av );

    eoState state;

    eoEvalFuncPtr< EOT, double > eval( objective_function );
    moFullEvalByCopy< Neighbor > neval( eval );
    
    int dimSize = 15;
    
    //-------------------------------------------------------
    // Parameters
    //-------------------------------------------------------

    std::string section( "Temperature initialization paramaters" );

    unsigned int dimension_size = parser.getORcreateParam( (unsigned int)dimSize, "dimension-size", "Dimension size", 'd', section ).value();
    double jump_bound = parser.getORcreateParam( (double)1, "jump-bound", "Bound of jump", '\0', section ).value();
    unsigned int maxiter = parser.getORcreateParam( (unsigned int)10, "max-iter", "Maximum number of iterations", '\0', section ).value();
    
    //-------------------------------------------------------


    //-------------------------------------------------------
    // Instanciate needed classes
    //-------------------------------------------------------

    edoUniform< EOT > distrib( EOT(dimension_size, -1 * jump_bound), EOT(dimension_size, 1 * jump_bound) );
    
    edoBounder< EOT > bounder_search( EOT(dimension_size, -10), EOT(dimension_size, 10) );
    
    edoSamplerUniform< EOT > sampler( bounder_search );
    
    moRealNeighborhood< edoUniform< EOT >, Neighbor > neighborhood( distrib, sampler, bounder_search );
    
    moStdDevEstimator< EOT, Neighbor > init( maxiter, neighborhood, eval );
    
    //-------------------------------------------------------



    
    //-------------------------------------------------------
    // Help + Verbose routines
    //-------------------------------------------------------

    if (parser.userNeedsHelp())
        {
            parser.printHelp(std::cout);
            exit(1);
        }

    make_help(parser);

    //-------------------------------------------------------

    
    EOT solution(dimSize, 5);
    
    std::cout << "init temp: "
          << init( solution )
          << std::endl;
    
    
    

    moTrueContinuator<Neighbor> continuator;
    moCheckpoint<Neighbor> checkpoint(continuator);
    
    
    moFitnessStat<EOT> fitStat;
    checkpoint.add(fitStat);
    eoFileMonitor monitor("triki.out", "");
    //eoGnuplot1DMonitor monitor2("trikignu.out", true);
    moCounterMonitorSaver countMon(100, monitor);
    checkpoint.add(countMon);
    //moCounterMonitorSaver gnuMon (1, monitor2);
    //checkpoint.add(gnuMon);
    monitor.add(fitStat);
    //monitor2.add(fitStat);
    
    
    double stdDevEst = init( solution );
    
    moTrikiCoolingSchedule<EOT> coolingSchedule (
            stdDevEst,
            stdDevEst,
            //50,
            150,//200, //150,
            //100
            250//350 // 250
        );
    //moSA<Neighbor> localSearch(neighborhood, eval, neval, coolingSchedule, checkpoint);
    //moSA<Neighbor> localSearch(neighborhood, eval, neval);
    //moSA<Neighbor> localSearch(neighborhood, eval, coolingSchedule, neval, checkpoint);
    //moSA<Neighbor> localSearch(neighborhood, eval, coolingSchedule);
    moSA<Neighbor> localSearch(neighborhood, eval, coolingSchedule, neval, checkpoint);
    
    std::cout << "#########################################" << std::endl;
    std::cout << "initial solution1: " << solution << std::endl ;
    
    localSearch(solution);
    
    std::cout << "final solution1: " << solution << std::endl ;
    std::cout << "#########################################" << std::endl;
    
    
    return 0;
}