freESDominantRecombinator.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: freESDominantRecombinator.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 _freESDominantRecombinator_txx
#define _freESDominantRecombinator_txx

#include "freESDominantRecombinator.h"

namespace FREE
{
  namespace ES
  {


    template <class TIndividual>
    DominantRecombinator<TIndividual>
      ::DominantRecombinator()
    {
      m_RandomGenerator = VariateGeneratorType::New();
      m_RandomGenerator->Initialize((int)clock());
    }

    template <class TIndividual>
    typename DominantRecombinator<TIndividual>::IndividualPointer
      DominantRecombinator<TIndividual>
      ::Recombine(const ParentSelectionType& parents) const
    {
      if (!(parents.size()>0)) throwExceptionMacro("Error. At least one parent is needed for recombination.");

      IndividualPointer newIndividual = parents[0]->Clone();

      for (unsigned int index = 0; index<newIndividual->ObjectiveParameters().size(); index++)
      {
        VariateGeneratorType::IntegerType parentIndex = m_RandomGenerator->GetIntegerVariate(parents.size()-1);
        (newIndividual->ObjectiveParameters())[index] = (parents[parentIndex]->ObjectiveParameters())[index];
      }
      
      for (unsigned int index = 0; index<newIndividual->StrategicParameters().size(); index++)
      {
        if ( newIndividual->StrategicParameters().at(index)->GetHandling()!= IndividualType::StrategicParameterType::HDynamic)
        {
          VariateGeneratorType::IntegerType parentIndex = m_RandomGenerator->GetIntegerVariate(parents.size()-1);
          (newIndividual->StrategicParameters())[index] = (parents[parentIndex]->StrategicParameters())[index];
        }
      }
      return newIndividual;
    }

    template <class TIndividual>
    typename DominantRecombinator<TIndividual>::WeightVectorType
    DominantRecombinator<TIndividual>
    ::GetWeights(const unsigned long lParentCount) const
    {
      WeightVectorType weights;
      weights.resize(lParentCount,0.0);
      if (lParentCount) weights[0] = 1.0;
      return weights;
    };

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

#endif

---