/*! * \brief Recombine individuals by weighted mean to generate a new individual. * \param inIndivPool Parents being recombined. * \param ioContext Evolutionary context. * \return Children generated by recombination. */ Individual::Handle SAES::RecombinationWeightedOp::recombine(Individual::Bag& inIndivPool, Context& ioContext) { Beagle_StackTraceBeginM(); // Compute recombination weights. std::vector<double> lWeights(inIndivPool.size()); for(unsigned int i=0; i<lWeights.size(); ++i) { lWeights[i] = std::log(double(lWeights.size()+1)); lWeights[i] -= std::log(double(i+1)); } // Recombine parents to generate new individual. const Factory& lFactory = ioContext.getSystem().getFactory(); Individual::Alloc::Handle lIndivAlloc = castHandleT<Individual::Alloc>(lFactory.getConceptAllocator("Individual")); Genotype::Alloc::Handle lGenotypeAlloc = castHandleT<Genotype::Alloc>(lFactory.getConceptAllocator("Genotype")); Individual::Handle lChildIndiv = castHandleT<Individual>(lIndivAlloc->allocate()); std::vector< std::vector<double> > lCountGenoSum; for(unsigned int i=0; i<inIndivPool.size(); ++i) { const unsigned int lPoolISize = inIndivPool[i]->size(); const unsigned int lChildSize = lChildIndiv->size(); if(lPoolISize > lChildSize) { lCountGenoSum.resize(lPoolISize); lChildIndiv->resize(lPoolISize); for(unsigned int j=lChildSize; j<lPoolISize; ++j) { (*lChildIndiv)[j] = castHandleT<Genotype>(lGenotypeAlloc->allocate()); } } for(unsigned int j=0; j<lPoolISize; ++j) { SAES::PairVector::Handle lChildGenoJ = castHandleT<SAES::PairVector>((*lChildIndiv)[j]); SAES::PairVector::Handle lPoolIGenoJ = castHandleT<SAES::PairVector>((*inIndivPool[i])[j]); const unsigned int lPoolIGenoJSize = lPoolIGenoJ->size(); if(lPoolIGenoJSize > lChildGenoJ->size()) { lChildGenoJ->resize(lPoolIGenoJSize,0.0); lCountGenoSum[j].resize(lPoolIGenoJSize,0); } for(unsigned int k=0; k<lPoolIGenoJSize; ++k) { (*lChildGenoJ)[k].mValue += (lWeights[i] * (*lPoolIGenoJ)[k].mValue); (*lChildGenoJ)[k].mStrategy += (lWeights[i] * (*lPoolIGenoJ)[k].mStrategy); lCountGenoSum[j][k] += lWeights[i]; } } } for(unsigned int i=0; i<lChildIndiv->size(); ++i) { SAES::PairVector::Handle lChildGenoI = castHandleT<SAES::PairVector>((*lChildIndiv)[i]); for(unsigned int j=0; j<lChildGenoI->size(); ++j) { (*lChildGenoI)[j].mValue /= lCountGenoSum[i][j]; (*lChildGenoI)[j].mStrategy /= lCountGenoSum[i][j]; } } Beagle_LogDebugM(ioContext.getSystem().getLogger(), *lChildIndiv); return lChildIndiv; Beagle_StackTraceEndM(); }
/*! * \brief Recombine individuals by averaging to generate a new individual. * \param inIndivPool Parents being recombined. * \param ioContext Evolutionary context. * \return Children generated by recombination. */ Individual::Handle Beagle::GA::RecombinationESVecOp::recombine(Individual::Bag& inIndivPool, Context& ioContext) { Beagle_StackTraceBeginM(); // Recombine parents to generate new individual. const Factory& lFactory = ioContext.getSystem().getFactory(); Individual::Alloc::Handle lIndivAlloc = castHandleT<Individual::Alloc>(lFactory.getConceptAllocator("Individual")); Genotype::Alloc::Handle lGenotypeAlloc = castHandleT<Genotype::Alloc>(lFactory.getConceptAllocator("Genotype")); Individual::Handle lChildIndiv = castHandleT<Individual>(lIndivAlloc->allocate()); std::vector< std::vector<unsigned int> > lCountGenoSum; for(unsigned int i=0; i<inIndivPool.size(); ++i) { const unsigned int lPoolISize = inIndivPool[i]->size(); const unsigned int lChildSize = lChildIndiv->size(); if(lPoolISize > lChildSize) { lCountGenoSum.resize(lPoolISize); lChildIndiv->resize(lPoolISize); for(unsigned int j=lChildSize; j<lPoolISize; ++j) { (*lChildIndiv)[j] = castHandleT<Genotype>(lGenotypeAlloc->allocate()); } } for(unsigned int j=0; j<lPoolISize; ++j) { GA::ESVector::Handle lChildGenoJ = castHandleT<GA::ESVector>((*lChildIndiv)[j]); GA::ESVector::Handle lPoolIGenoJ = castHandleT<GA::ESVector>((*inIndivPool[i])[j]); const unsigned int lPoolIGenoJSize = lPoolIGenoJ->size(); if(lPoolIGenoJSize > lChildGenoJ->size()) { lChildGenoJ->resize(lPoolIGenoJSize,0.0); lCountGenoSum[j].resize(lPoolIGenoJSize,0); } for(unsigned int k=0; k<lPoolIGenoJSize; ++k) { (*lChildGenoJ)[k].mValue += (*lPoolIGenoJ)[k].mValue; (*lChildGenoJ)[k].mStrategy += (*lPoolIGenoJ)[k].mStrategy; ++lCountGenoSum[j][k]; } } } for(unsigned int i=0; i<lChildIndiv->size(); ++i) { GA::ESVector::Handle lChildGenoI = castHandleT<GA::ESVector>((*lChildIndiv)[i]); for(unsigned int j=0; j<lChildGenoI->size(); ++j) { (*lChildGenoI)[j].mValue /= double(lCountGenoSum[i][j]); (*lChildGenoI)[j].mStrategy /= double(lCountGenoSum[i][j]); } } Beagle_LogDebugM( ioContext.getSystem().getLogger(), "crossover", "Beagle::GA::RecombinationESVecOp", "Individual generated by recombination" ); Beagle_LogObjectDebugM( ioContext.getSystem().getLogger(), "crossover", "Beagle::GA::RecombinationESVecOp", *lChildIndiv ); return lChildIndiv; Beagle_StackTraceEndM("Individual::Handle Beagle::GA::RecombinationESVecOp::recombine(Individual::Bag& inIndivPool,Context& ioContext)"); }