freESAdaptiveNormalMutation.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: freESAdaptiveNormalMutation.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 _freESAdaptiveNormalMutation_txx
#define _freESAdaptiveNormalMutation_txx

#include "freESAdaptiveNormalMutation.h"

namespace FREE
{
namespace ES
{

template <class TIndividual>
AdaptiveNormalMutation<TIndividual>
::AdaptiveNormalMutation()
{
  m_StrategicMutationVariance = 1.0;
  m_MutateValueVariance = 1.0;
  m_RandomGenerator = VariateGeneratorType::New();
}

template <class TIndividual>
void
AdaptiveNormalMutation<TIndividual>
::RegisterStrategicIndividualParameters(IndividualType& individual) const
{
  for (unsigned long index = 0; index<individual.ObjectiveParameters().size(); index++)
  {
    typename IndividualType::StrategicParameterPointer parameter = IndividualType::StrategicParameterType::New();

    parameter->SetValue(1.0);
    parameter->SetHandling(IndividualType::StrategicParameterType::HDynamic);
    parameter->SetOrigin(this->GetNameOfClass());

    individual.StrategicParameters().push_back(parameter);
  }
};

template <class TIndividual>
void
AdaptiveNormalMutation<TIndividual>
::MutateStrategicIndividualParameters(IndividualType& individual, PopulationType& population) const
{
  if (m_StrategicMutationVariance<=0) throwExceptionMacro("Error. Strategic mutation variance should not be Zero or less. Actual value: "<<m_StrategicMutationVariance);

  //get all strategic values of the mutation
  typename IndividualType::StrategicParametersType::SelectedParametersType parameters =
    individual.StrategicParameters().GetStrategicParameters(this->GetNameOfClass());

  for (unsigned long index = 0; index<parameters.size(); index++)
  {
    typename IndividualType::StrategicParameterValueType variance = (parameters[index])->GetValue();

    //mutate parameter
    variance = variance * (exp(m_RandomGenerator->GetVariate()*m_StrategicMutationVariance));
    
    //set new value;
    (parameters[index])->SetValue(variance);
  }
};

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

template <class TIndividual>
void
AdaptiveNormalMutation<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());

  for (unsigned long index = 0; index<pIndividual->ObjectiveParameters().size(); index++)
  {
    typename IndividualType::StrategicParameterValueType variance = (parameters[index])->GetValue();
    typename IndividualType::ObjectiveParameterType mutatedParameter = (pIndividual->ObjectiveParameters())[index];

    //mutate parameter
    mutatedParameter += (m_RandomGenerator->GetVariate()*variance)/this->m_GeneralObjectiveScales[index];
    
    //set new value;
    (pIndividual->ObjectiveParameters())[index] = mutatedParameter;
  }
};


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

#endif

---