/*! * \brief Validate the GP individual. * \param ioContext Evolutionary context. * \return True if the GP individual is valid, false if not. */ bool GP::Individual::validate(GP::Context& ioContext) { Beagle_StackTraceBeginM(); bool lResult = true; Beagle_LogDetailedM( ioContext.getSystem().getLogger(), std::string("Validating ")+uint2ordinal(ioContext.getIndividualIndex()+1)+ std::string(" individual") ); // Store original values. GP::Tree::Handle lOldTreeHandle = ioContext.getGenotypeHandle(); unsigned int lOldTreeIndex = ioContext.getGenotypeIndex(); // Loop through each of the trees in the individual for (unsigned int i=0; i<size(); i++) { GP::Tree::Handle lTree = (*this)[i]; if (lTree == NULL) { Beagle_LogVerboseM( ioContext.getSystem().getLogger(), std::string("Skipping ")+uint2ordinal(i+1)+std::string(" tree because it's NULL-valued") ); continue; } Beagle_LogVerboseM( ioContext.getSystem().getLogger(), std::string("Validating ")+uint2ordinal(i+1)+std::string(" tree") ); // Store the new values ioContext.setGenotypeHandle(lTree); ioContext.setGenotypeIndex(i); lTree->setContextToNode(0, ioContext); if(!lTree->validateSubTree(0, ioContext)) { Beagle_LogVerboseM( ioContext.getSystem().getLogger(), std::string("Validation of ")+uint2ordinal(i+1)+std::string(" tree failed.") ); lResult = false; break; } } if(lResult) { Beagle_LogVerboseM( ioContext.getSystem().getLogger(), std::string("Individual passed validation testing.") ); } // Restore the original values. ioContext.setGenotypeHandle(lOldTreeHandle); ioContext.setGenotypeIndex(lOldTreeIndex); return lResult; Beagle_StackTraceEndM(); }
/*! * \brief Insert mutate a GP individual. * \param ioIndividual GP individual to mutate. * \param ioContext Context of the evolution. * \return True if the individual is effectively mutated, false if not. */ bool GP::MutationInsertConstrainedOp::mutate(Beagle::Individual& ioIndividual, Beagle::Context& ioContext) { Beagle_StackTraceBeginM(); Beagle_LogDetailedM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Mutating individual with GP::MutationInsertConstrainedOp") ); GP::Individual& lIndividual = castObjectT<GP::Individual&>(ioIndividual); GP::Context& lContext = castObjectT<GP::Context&>(ioContext); const unsigned int lMaxAttempts = mNumberAttempts->getWrappedValue(); const unsigned int lMaxTreeDepth = mMaxTreeDepth->getWrappedValue(); const Factory& lFactory = ioContext.getSystem().getFactory(); GP::Tree::Alloc::Handle lTreeAlloc = castHandleT<GP::Tree::Alloc>(lFactory.getConceptAllocator("Genotype")); // Store original context values const unsigned int lOldGenotypeIndex = lContext.getGenotypeIndex(); const GP::Tree::Handle lOldGenotypeHandle = lContext.getGenotypeHandle(); Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", "Individual before constrained GP insert mutation" ); Beagle_LogObjectDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", ioIndividual ); // Mutation attempts loop for(unsigned int i=0; i<lMaxAttempts; ++i) { // Choose tree and node to mutate const unsigned int lChosenTree = lIndividual.chooseRandomTree(lContext); const unsigned int lChosenNodeIndex = lIndividual.chooseRandomNode(lChosenTree, lContext); GP::Tree::Handle lOriginalTree = lIndividual[lChosenTree]; // Compute depth of tree generated by mutation lContext.setGenotypeIndex(lChosenTree); lContext.setGenotypeHandle(lOriginalTree); lContext.emptyCallStack(); lOriginalTree->setContextToNode(lChosenNodeIndex, lContext); const unsigned int lMutationDepth = lContext.getCallStackSize() + lOriginalTree->getTreeDepth(lChosenNodeIndex); // Check that mutation will not generate a tree deeper than the maximum allowed depth if(lMutationDepth > lMaxTreeDepth) { Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Constrained insert mutation attempt failed as the generated tree will exceed ")+ std::string("maximum allowed tree depth") ); continue; } // Create new tree Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Creating new tree") ); GP::Tree::Handle lNewTree = castHandleT<GP::Tree>(lTreeAlloc->allocate()); lNewTree->setPrimitiveSetIndex(lOriginalTree->getPrimitiveSetIndex()); lNewTree->setNumberArguments(lOriginalTree->getNumberArguments()); // Replace original tree with new tree lIndividual[lChosenTree] = lNewTree; lContext.setGenotypeHandle(lNewTree); // Copy unchanged part of original tree into new tree Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Copying unchanged part of original tree to new tree") ); lNewTree->insert(lNewTree->end(), lOriginalTree->begin(), lOriginalTree->begin()+lChosenNodeIndex); Beagle_AssertM(lNewTree->size() == lChosenNodeIndex); // Generate new branch primitive to insert. GP::PrimitiveSet& lPrimitiveSet = lNewTree->getPrimitiveSet(lContext); Primitive::Handle lBranchInserted = lPrimitiveSet.select(GP::Primitive::eBranch, lContext); if(lBranchInserted==NULL) { Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Constrained insert mutation attempt failed as it seems impossible ")+ std::string("to select a branch primitive in the actual context") ); lIndividual[lChosenTree] = lOriginalTree; lContext.setGenotypeHandle(lOriginalTree); continue; } Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Branch primitive to be inserted by mutation is primitive '")+ lBranchInserted->getName()+std::string("'") ); // Insert new branch lBranchInserted = lBranchInserted->giveReference(GP::Primitive::eBranch, lContext); lNewTree->push_back(Node(lBranchInserted,1)); if(lBranchInserted->validate(lContext) == false) { Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Constrained insert mutation attempt failed as the selected ")+ std::string("branch doesn't match the constraints") ); lIndividual[lChosenTree] = lOriginalTree; lContext.setGenotypeHandle(lOriginalTree); continue; } const unsigned int lNbArgsInsertedBranch = lBranchInserted->getNumberArguments(); Beagle_AssertM(lNbArgsInsertedBranch != 0); const unsigned int lSubtreeArgIndex = lContext.getSystem().getRandomizer().rollInteger(0, lNbArgsInsertedBranch-1); Beagle_AssertM(lSubtreeArgIndex < lNbArgsInsertedBranch); const unsigned int lSubtreeSize = (*lOriginalTree)[lChosenNodeIndex].mSubTreeSize; // Generate inserted node subtrees bool lArgsGenFailed = false; for(unsigned int j=0; j<lNbArgsInsertedBranch; ++j) { if(j == lSubtreeArgIndex) { const unsigned int lSubtreeIndex = lNewTree->size(); lNewTree->insert(lNewTree->end(), lOriginalTree->begin()+lChosenNodeIndex, lOriginalTree->begin()+lChosenNodeIndex+lSubtreeSize); lContext.pushCallStack(lSubtreeIndex); if(lNewTree->validateSubTree(lSubtreeIndex, lContext) == false) { lArgsGenFailed = true; break; } lContext.popCallStack(); (*lNewTree)[lChosenNodeIndex].mSubTreeSize += lSubtreeSize; } else { Primitive::Handle lArgInserted = lPrimitiveSet.select(GP::Primitive::eTerminal, lContext); if(lArgInserted == NULL) { lArgsGenFailed = true; break; } lArgInserted = lArgInserted->giveReference(GP::Primitive::eTerminal, lContext); const unsigned int lSubtreeIndex = lNewTree->size(); lNewTree->push_back(Node(lArgInserted,1)); lContext.pushCallStack(lSubtreeIndex); if(lArgInserted->validate(lContext) == false) { lArgsGenFailed = true; break; } lContext.popCallStack(); ++(*lNewTree)[lChosenNodeIndex].mSubTreeSize; } } if(lArgsGenFailed) { Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Constrained insert mutation attempt failed as it seems impossible ")+ std::string("to select a terminal primitive under the inserted branch in the actual context") ); lIndividual[lChosenTree] = lOriginalTree; lContext.setGenotypeHandle(lOriginalTree); continue; } // Complete new tree with rest of original tree Beagle_AssertM(lOriginalTree->size() >= (lChosenNodeIndex+lSubtreeSize)); lNewTree->insert(lNewTree->end(), lOriginalTree->begin()+lChosenNodeIndex+lSubtreeSize, lOriginalTree->end()); Beagle_AssertM(lNewTree->size() == (lOriginalTree->size()+lNbArgsInsertedBranch)); // Correct subtree size data and terminate mutation process Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("Correcting the 'mSubTreeSize' fields of tree") ); lContext.popCallStack(); for(unsigned int j=0; j<lContext.getCallStackSize(); ++j) { (*lNewTree)[lContext[j]].mSubTreeSize += lNbArgsInsertedBranch; } // Mutation successful, log messages and return std::ostringstream lOSS; lOSS << "Successfully inserted a new node at index " << lChosenNodeIndex; lOSS << " of the " << uint2ordinal(lChosenTree) << " tree of the actual individual"; Beagle_LogTraceM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", lOSS.str() ); Beagle_LogDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", "Individual after constrained GP insert mutation" ); Beagle_LogObjectDebugM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", ioIndividual ); lContext.emptyCallStack(); lContext.setGenotypeIndex(lOldGenotypeIndex); lContext.setGenotypeHandle(lOldGenotypeHandle); return true; } // Insert mutation failed, return without mutating the individual lContext.emptyCallStack(); lContext.setGenotypeIndex(lOldGenotypeIndex); lContext.setGenotypeHandle(lOldGenotypeHandle); Beagle_LogTraceM( ioContext.getSystem().getLogger(), "mutation", "Beagle::GP::MutationInsertConstrainedOp", std::string("All constrained insert mutation attempts failed; ")+ std::string("resuming from mutation without modifying the individual") ); return false; Beagle_StackTraceEndM("bool GP::MutationInsertConstrainedOp::mutate(Beagle::Individual& ioIndividual, Beagle::Context& ioContext)"); }