#region License Information
/* HeuristicLab
* Copyright (C) 2002-2010 Heuristic and Evolutionary Algorithms Laboratory (HEAL)
*
* This file is part of HeuristicLab.
*
* HeuristicLab is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* HeuristicLab is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with HeuristicLab. If not, see .
*/
#endregion
using System;
using System.Linq;
using HeuristicLab.Common;
using HeuristicLab.Core;
using HeuristicLab.Data;
using HeuristicLab.Operators;
using HeuristicLab.Optimization;
using HeuristicLab.Parameters;
using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
using HeuristicLab.Encodings.SymbolicExpressionTreeEncoding.Symbols;
using System.Collections.Generic;
using System.Text;
using System.Diagnostics;
namespace HeuristicLab.Encodings.SymbolicExpressionTreeEncoding.ArchitectureManipulators {
///
/// Creates a new argument within one function-defining branch of a symbolic expression tree.
/// As described in Koza, Bennett, Andre, Keane, Genetic Programming III - Darwinian Invention and Problem Solving, 1999, pp. 106
///
[Item("ArgumentCreater", "Manipulates a symbolic expression by creating a new argument within one function-defining branch.")]
[StorableClass]
public sealed class ArgumentCreater : SymbolicExpressionTreeArchitectureManipulator {
public override sealed void ModifyArchitecture(
IRandom random,
SymbolicExpressionTree symbolicExpressionTree,
ISymbolicExpressionGrammar grammar,
IntValue maxTreeSize, IntValue maxTreeHeight,
IntValue maxFunctionDefiningBranches, IntValue maxFunctionArguments,
out bool success) {
success = CreateNewArgument(random, symbolicExpressionTree, grammar, maxTreeSize.Value, maxTreeHeight.Value, maxFunctionDefiningBranches.Value, maxFunctionArguments.Value);
}
public static bool CreateNewArgument(
IRandom random,
SymbolicExpressionTree symbolicExpressionTree,
ISymbolicExpressionGrammar grammar,
int maxTreeSize, int maxTreeHeight,
int maxFunctionDefiningBranches, int maxFunctionArguments) {
var functionDefiningBranches = symbolicExpressionTree.IterateNodesPrefix().OfType();
if (functionDefiningBranches.Count() == 0)
// no function defining branch found => abort
return false;
// select a random function defining branch
var selectedDefunBranch = functionDefiningBranches.SelectRandom(random);
var definedArguments = (from symbol in selectedDefunBranch.Grammar.Symbols.OfType()
select symbol.ArgumentIndex).Distinct();
if (definedArguments.Count() >= maxFunctionArguments)
// max number of arguments reached => abort
return false;
// select a random cut point in the function defining branch
// the branch at the cut point is to be replaced by a new argument node
var cutPoints = (from node in selectedDefunBranch.IterateNodesPrefix()
where node.SubTrees.Count > 0
from subtree in node.SubTrees
select new { Parent = node, ReplacedChildIndex = node.SubTrees.IndexOf(subtree), ReplacedChild = subtree }).ToList();
if (cutPoints.Count() == 0)
// no cut point found => abort;
return false;
var selectedCutPoint = cutPoints[random.Next(cutPoints.Count)];
var allowedArgumentIndexes = Enumerable.Range(0, maxFunctionArguments);
var newArgumentIndex = allowedArgumentIndexes.Except(definedArguments).First();
// replace the branch at the cut point with an argument node
var newArgNode = MakeArgumentNode(newArgumentIndex);
var replacedBranch = selectedCutPoint.ReplacedChild;
selectedCutPoint.Parent.RemoveSubTree(selectedCutPoint.ReplacedChildIndex);
selectedCutPoint.Parent.InsertSubTree(selectedCutPoint.ReplacedChildIndex, newArgNode);
// find all invocations of the selected ADF and attach a cloned version of the replaced branch (with all argument-nodes expanded)
var invocationNodes = from node in symbolicExpressionTree.IterateNodesPrefix().OfType()
where node.Symbol.FunctionName == selectedDefunBranch.FunctionName
select node;
foreach (var invocationNode in invocationNodes) {
// append a new argument branch after expanding all argument nodes
var clonedBranch = (SymbolicExpressionTreeNode)replacedBranch.Clone();
clonedBranch = ReplaceArgumentsInBranch(clonedBranch, invocationNode.SubTrees);
invocationNode.InsertSubTree(newArgumentIndex, clonedBranch);
}
// increase expected number of arguments of function defining branch
// it's possible that the number of actually referenced arguments was reduced (all references were replaced by a single new argument)
// but the number of expected arguments is increased anyway
selectedDefunBranch.NumberOfArguments++;
selectedDefunBranch.Grammar.AddSymbol(newArgNode.Symbol);
selectedDefunBranch.Grammar.SetMinSubtreeCount(newArgNode.Symbol, 0);
selectedDefunBranch.Grammar.SetMaxSubtreeCount(newArgNode.Symbol, 0);
// allow the argument as child of any other symbol
foreach (var symb in selectedDefunBranch.Grammar.Symbols)
for (int i = 0; i < selectedDefunBranch.Grammar.GetMaxSubtreeCount(symb); i++) {
selectedDefunBranch.Grammar.SetAllowedChild(symb, newArgNode.Symbol, i);
}
foreach (var subtree in symbolicExpressionTree.Root.SubTrees) {
// when the changed function is known in the branch then update the number of arguments
var matchingSymbol = subtree.Grammar.Symbols.OfType().Where(s => s.FunctionName == selectedDefunBranch.FunctionName).SingleOrDefault();
if (matchingSymbol != null) {
subtree.Grammar.SetMinSubtreeCount(matchingSymbol, selectedDefunBranch.NumberOfArguments);
subtree.Grammar.SetMaxSubtreeCount(matchingSymbol, selectedDefunBranch.NumberOfArguments);
foreach (var child in subtree.GetAllowedSymbols(0)) {
for (int i = 0; i < subtree.Grammar.GetMaxSubtreeCount(matchingSymbol); i++) {
subtree.Grammar.SetAllowedChild(matchingSymbol, child, i);
}
}
}
}
return true;
}
private static SymbolicExpressionTreeNode ReplaceArgumentsInBranch(SymbolicExpressionTreeNode branch, IList argumentTrees) {
if (branch is ArgumentTreeNode) {
var argNode = branch as ArgumentTreeNode;
// replace argument nodes by a clone of the original subtree that provided the result for the argument node
return (SymbolicExpressionTreeNode)argumentTrees[argNode.Symbol.ArgumentIndex].Clone();
} else {
// call recursively for all subtree
List subtrees = new List(branch.SubTrees);
while (branch.SubTrees.Count > 0) branch.RemoveSubTree(0);
foreach (var subtree in subtrees) {
branch.AddSubTree(ReplaceArgumentsInBranch(subtree, argumentTrees));
}
return branch;
}
}
private static SymbolicExpressionTreeNode MakeArgumentNode(int argIndex) {
var node = (ArgumentTreeNode)(new Argument(argIndex)).CreateTreeNode();
return node;
}
}
}