source: branches/symbreg-factors-2650/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Symbols/FactorVariableTreeNode.cs @ 14259

#region License Information
/* HeuristicLab
 * Copyright (C) 2002-2016 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 <http://www.gnu.org/licenses/>.
 */
#endregion

using System;
using System.Linq;
using HeuristicLab.Common;
using HeuristicLab.Core;
using HeuristicLab.Encodings.SymbolicExpressionTreeEncoding;
using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
using HeuristicLab.Random;
namespace HeuristicLab.Problems.DataAnalysis.Symbolic {
  [StorableClass]
  public class FactorVariableTreeNode : SymbolicExpressionTreeTerminalNode, IVariableTreeNode {
    public new FactorVariable Symbol {
      get { return (FactorVariable)base.Symbol; }
    }
    [Storable]
    private double[] weights;
    public double[] Weights {
      get { return weights; }
      set { weights = value; }
    }
    [Storable]
    private string variableName;
    public string VariableName {
      get { return variableName; }
      set { variableName = value; }
    }

    [StorableConstructor]
    protected FactorVariableTreeNode(bool deserializing) : base(deserializing) { }
    protected FactorVariableTreeNode(FactorVariableTreeNode original, Cloner cloner)
      : base(original, cloner) {
      variableName = original.variableName;
      if (original.weights != null) {
        this.weights = new double[original.Weights.Length];
        Array.Copy(original.Weights, weights, weights.Length);
      }
    }

    public FactorVariableTreeNode(FactorVariable variableSymbol)
      : base(variableSymbol) {
    }

    public override bool HasLocalParameters {
      get { return true; }
    }

    public override void ResetLocalParameters(IRandom random) {
      base.ResetLocalParameters(random);
      variableName = Symbol.VariableNames.SampleRandom(random);
      weights =
        Symbol.GetVariableValues(variableName)
        .Select(_ => NormalDistributedRandom.NextDouble(random, 0, 1)).ToArray();
    }

    public override void ShakeLocalParameters(IRandom random, double shakingFactor) {
      if (random.NextDouble() < 0.2) {
        VariableName = Symbol.VariableNames.SampleRandom(random);
        if (weights.Length != Symbol.GetVariableValues(VariableName).Count()) {
          // if the length of the weight array does not match => re-initialize weights
          weights =
            Symbol.GetVariableValues(variableName)
              .Select(_ => NormalDistributedRandom.NextDouble(random, 0, 1))
              .ToArray();
        }
      } else {
        // mutate only one randomly selected weight
        var idx = random.Next(weights.Length);
        // 50% additive & 50% multiplicative
        if (random.NextDouble() < 0.5) {
          double x = NormalDistributedRandom.NextDouble(random, Symbol.WeightManipulatorMu,
            Symbol.WeightManipulatorSigma);
          weights[idx] = weights[idx] + x * shakingFactor;
        } else {
          double x = NormalDistributedRandom.NextDouble(random, 1.0, Symbol.MultiplicativeWeightManipulatorSigma);
          weights[idx] = weights[idx] * x;
        }
      }
    }

    public override IDeepCloneable Clone(Cloner cloner) {
      return new FactorVariableTreeNode(this, cloner);
    }

    public double GetValue(string cat) {
      // TODO: perf
      var s = Symbol;
      int idx = 0;
      foreach (var val in s.GetVariableValues(VariableName)) {
        if (cat == val) return weights[idx];
        idx++;
      }
      throw new ArgumentOutOfRangeException("Found unknown value " + cat + " for variable " + VariableName);
    }

    public override string ToString() {
      var weightStr = string.Join("; ",
        Symbol.GetVariableValues(VariableName).Select(value => value + ": " + GetValue(value).ToString("E4")));
      return VariableName + " (factor) "
        + "[" + weightStr + "]";
    }
  }
}