source: addons/HeuristicLab.Problems.BioBoost/HeuristicLab.Problems.BioBoost/3.3/Evaluators/SinkEvaluator.cs @ 17777

#region License Information
/* HeuristicLab
 * Copyright (C) 2002-2015 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 HeuristicLab.Common;
using HeuristicLab.Core;
using HeuristicLab.Data;
using HeuristicLab.Optimization;
using HeuristicLab.Parameters;
using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
using System.Linq;
using HeuristicLab.BioBoost.Representation;
using HEAL.Attic;

namespace HeuristicLab.BioBoost.Evaluators {
  [StorableType("EB23EA04-FF16-43C0-902D-9DC974F5A246")]
  public class SinkEvaluator : BioBoostEvaluator {

    #region Construction & Cloning
    [StorableConstructor]
    protected SinkEvaluator(StorableConstructorFlag _) : base(_) { }
    protected SinkEvaluator(SinkEvaluator original, Cloner cloner) : base(original, cloner) { }
    public SinkEvaluator() {
      Parameters.Add(new ScopeParameter("Scope", "The scope to search products in."));
    }
    public override IDeepCloneable Clone(Cloner cloner) {
      return new SinkEvaluator(this, cloner);
    }
    #endregion

    private static LayerDescriptor emptyLayerDescriptor = new LayerDescriptor("", "", "");

    public override IOperation Apply() {
      foreach (var product in ProblemData.Products) {
        InvoiceProduct(product.Label, product.Price, emptyLayerDescriptor, LayerDescriptor.DischargeCosts);
        InvoiceProduct(product.Label, product.PriceAtSource, LayerDescriptor.AmountsAtSource, LayerDescriptor.TotalCostsAtSource);
        InvoiceProduct(product.Label, product.PriceAtTarget, LayerDescriptor.AmountsAtTarget, LayerDescriptor.TotalCostsAtTarget);
      }
      return base.Apply();
    }

    private void InvoiceProduct(string productName, double price, LayerDescriptor amountLayerDescriptor, LayerDescriptor costLayerDescriptor) {
      var variables = ExecutionContext.Scope.Variables;
      IVariable amountVariable, costsVariable;
      if (!variables.TryGetValue(amountLayerDescriptor.NameWithPrefix(productName), out amountVariable)) return;
      var amounts = amountVariable.Value as DoubleArray;
      if (amounts == null) return;
      var relativeProductionCosts = new DoubleArray(amounts.Length);
      DoubleArray productionCosts = null;
      if (variables.TryGetValue(costLayerDescriptor.NameWithPrefix(productName), out costsVariable)) {
        productionCosts = costsVariable.Value as DoubleArray;
      }
      var totalCost = amounts.Sum()*-price;
      var length = amounts.Length;
      var costs = new DoubleArray(length);
      for (int i = 0; i < length; i++) {
        costs[i] = amounts[i]*-price;
        if (productionCosts != null && amounts[i] != 0)
          relativeProductionCosts[i] = productionCosts[i]/amounts[i];
      }
      if (totalCost != 0) {
        AddInScope(LayerDescriptor.DischargeCosts.NameWithPrefix(productName), costs, false);
        Costs.Add(new Result(costLayerDescriptor.NameWithPrefix(productName), new DoubleValue(totalCost)));
        if (totalCost > 0)
          AddInScope(LayerDescriptor.LocalAddedValue.Name, costs, false);
      }
    }
  }
}