make_genotype_real.h

//-----------------------------------------------------------------------------
#ifndef EO_make_genotype_h
#define EO_make_genotype_h

#include <iostream>
#include <sstream>
#include <vector>

#include "es/eoReal.h"
#include "es/eoEsChromInit.h"
#include "utils/eoParser.h"
#include "utils/eoRealVectorBounds.h"
#include "utils/eoState.h"


template <class EOT>
eoEsChromInit<EOT> & do_make_genotype(eoParser& _parser, eoState& _state, EOT)
{
    // the fitness type
    typedef typename EOT::Fitness FitT;
    eoEsChromInit<EOT> *init;

    // for eoReal, only thing needed is the size - but might have been created elswhere ...
    eoValueParam<unsigned>& vecSize
        = _parser.getORcreateParam(unsigned(10), "vecSize",
                                   "The number of variables ",
                                   'n', "Genotype Initialization");
    // to build an eoReal Initializer, we need bounds: [-1,1] by default
    eoValueParam<eoRealVectorBounds>& boundsParam
        = _parser.getORcreateParam(eoRealVectorBounds(vecSize.value(), -1, 1),
                                   "initBounds",
                                   "Bounds for initialization (MUST be bounded)",
                                   'B', "Genotype Initialization");
    // now some initial value for sigmas - even if useless?
    // should be used in Normal mutation
    eoValueParam<std::string>& sigmaParam
        = _parser.getORcreateParam(std::string("0.3"), "sigmaInit",
                                   "Initial value for Sigmas (with a '%' -> scaled by the range of each variable)",
                                   's', "Genotype Initialization");
    // check for %
    bool to_scale = false;
    size_t pos =  sigmaParam.value().find('%');
    if(pos < sigmaParam.value().size()) {
        //  found a % - use scaling and get rid of '%'
        to_scale = true;
        sigmaParam.value().resize(pos);
    }
    std::istringstream is(sigmaParam.value());
    double sigma;
    is >> sigma;
    // minimum check
    if(sigma < 0)
        throw std::runtime_error("Negative sigma in make_genotype");
    if(to_scale)
        init = new eoEsChromInit<EOT>(boundsParam.value(), sigma, to_scale);
    else {
        // define parameter
        eoValueParam<std::vector<double> >& vecSigmaParam
            = _parser.getORcreateParam(std::vector<double>(vecSize.value(), sigma), "vecSigmaInit",
                                       "Initial value for Sigmas (only used when initSigma is not scaled)",
                                       'S', "Genotype Initialization");
        init = new eoEsChromInit<EOT>(boundsParam.value(), vecSigmaParam.value());
    }
    // store in state
    _state.storeFunctor(init);
    return *init;
}

#endif // EO_make_genotype_h



// Local Variables:
// coding: iso-8859-1
// mode:C++
// c-file-style: "Stroustrup"
// comment-column: 35
// fill-column: 80
// End: