freESAdaptiveScaleMutation.txx

Go to the documentation of this file.

/*=========================================================================

  Program:   F.R.E.E. - flexible registration evaluation engine
  Version:   v.1.0.0
  Date:      $Date: 2006/09/01 12:00:00 $
  Module:    $RCSfile: freESAdaptiveScaleMutation.txx,v $
  Language:  C++



  Copyright (c) 2007 Ralf o Floca (Department of Medical Informatics,
  Institute for Medical Biometry and Informatics, University of Heidelberg,
  Germany). All rights reserved.
  See FREECopyright.txt or http://www.mi.med.uni-hd.de/free/copyright.htm
  for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef _freESAdaptiveScaleMutation_txx
#define _freESAdaptiveScaleMutation_txx

#include "freESAdaptiveScaleMutation.h"

namespace FREE
{
namespace ES
{

template <class TIndividual>
AdaptiveScaleMutation<TIndividual>
::AdaptiveScaleMutation()
{
}

template <class TIndividual>
double
AdaptiveScaleMutation<TIndividual>
::MutateValue(const double& value)
{
    double mutatedValue = value;
    mutatedValue += ((this->m_RandomGenerator->GetVariate())*(this->m_MutateValueVariance))/this->m_MutateValueScale;
    return mutatedValue;
};

template <class TIndividual>
void
AdaptiveScaleMutation<TIndividual>
::MutateObjectivParameters(IndividualType* pIndividual, PopulationType* pPopulation) const
{
  if (!pIndividual) throwExceptionMacro("Error. Passed individual is Null.");

  //get the strategic parameters
  typename IndividualType::StrategicParametersType::SelectedParametersType parameters =
    pIndividual->StrategicParameters().GetStrategicParameters(this->GetNameOfClass());

  if (parameters.size()!=pIndividual->ObjectiveParameters().size()) throwExceptionMacro("Error. Number of objective parameters and related strategic parameters are not equal. Ensure that strategic parameters are properly registered. Objective count:" << pIndividual->ObjectiveParameters().size() <<"; strategic count: "<< parameters.size());
  if (parameters.size()!=this->m_GeneralObjectiveScales.size()) throwExceptionMacro("Error. Number of objective parameters and global scales are not equal. Ensure that scales have been set properly. Objective count:" << pIndividual->ObjectiveParameters().size() <<"; objective scales count: "<< this->m_GeneralObjectiveScales.size());

  //mutate all parameters
  for (unsigned long index = 0; index<pIndividual->ObjectiveParameters().size(); index++)
  {
    const typename IndividualType::StrategicParameterValueType variance = (parameters[index])->GetValue();
    const double variate = this->m_RandomGenerator->GetVariate();
    typename IndividualType::ObjectiveParameterType mutatedParameter = (pIndividual->ObjectiveParameters())[index];

    double scaling = 1.0;

    //calculate the scaling
    if (variate<0.0)
    {
      scaling = 1/(1-(variate * variance * (this->m_GeneralObjectiveScales[index])));
    }
    else
    {
      scaling = 1+(variate * variance * (this->m_GeneralObjectiveScales[index]));
    }

    //mutate parameter
    mutatedParameter *= scaling;
    
    //set new value;
    (pIndividual->ObjectiveParameters())[index] = mutatedParameter;
  }
};


} // end of namespace ES
} // end of namespace FREE

#endif

---