freESAdaptiveCovarianceMutation.txx

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: freESAdaptiveCovarianceMutation.txx,v $

  Date:      $Date: 2005/03/07 03:39:58 $
  Version:   $Revision: 1.32 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/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 _freESAdaptiveCovarianceMutation_txx
#define _freESAdaptiveCovarianceMutation_txx

#include "freESAdaptiveCovarianceMutation.h"

namespace FREE
{
namespace ES
{

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

template <class TIndividual>
void
AdaptiveCovarianceMutation<TIndividual>
::RegisterStrategicIndividualParameters(IndividualType& individual) const
{
  for (unsigned long index = 0; index<individual.ObjectiveParameters().size(); index++)
  {
    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
AdaptiveCovarianceMutation<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
  IndividualType::StrategicParametersType::SelectedParametersType parameters =
    individual.StrategicParameters().GetStrategicParameters(this->GetNameOfClass());

  for (unsigned long index = 0; index<parameters.size(); index++)
  {
    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
AdaptiveCovarianceMutation<TIndividual>
::MutateValue(const double& value)
{
    double mutatedValue = value;
    mutatedValue += (m_RandomGenerator->GetVariate()*m_MutateValueVariance)/m_MutateValueScale;
    return mutatedValue;
};

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

  //get the strategic parameters
  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()!=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: "<< m_GeneralObjectiveScales.size());

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

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


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

#endif

---