Changeset 9776 for branches/HeuristicLab.DataAnalysis.Symbolic.LinearInterpreter/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Interpreter/SymbolicDataAnalysisExpressionTreeLinearInterpreter.cs

Timestamp:

07/26/13 13:15:05 (11 years ago)

Author:

bburlacu

Message:

#2021: Added support for lagged symbols: TimeLag, Derivative, Integral.

File:

• branches/HeuristicLab.DataAnalysis.Symbolic.LinearInterpreter/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Interpreter/SymbolicDataAnalysisExpressionTreeLinearInterpreter.cs (modified) (11 diffs)

branches/HeuristicLab.DataAnalysis.Symbolic.LinearInterpreter/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Interpreter/SymbolicDataAnalysisExpressionTreeLinearInterpreter.cs

@@ -37 +37 @@
     private const string EvaluatedSolutionsParameterName = "EvaluatedSolutions";
 
+    private SymbolicDataAnalysisExpressionTreeInterpreter interpreter;
+
     public override bool CanChangeName {
       get { return false; }
@@ -74 +76 @@
       SymbolicDataAnalysisExpressionTreeLinearInterpreter original, Cloner cloner)
       : base(original, cloner) {
+      interpreter = original.interpreter;
     }
 
@@ -84 +87 @@
       Parameters.Add(new ValueParameter<BoolValue>(CheckExpressionsWithIntervalArithmeticParameterName, "Switch that determines if the interpreter checks the validity of expressions with interval arithmetic before evaluating the expression.", new BoolValue(false)));
       Parameters.Add(new ValueParameter<IntValue>(EvaluatedSolutionsParameterName, "A counter for the total number of solutions the interpreter has evaluated", new IntValue(0)));
+      interpreter = new SymbolicDataAnalysisExpressionTreeInterpreter();
     }
 
@@ -90 +94 @@
       Parameters.Add(new ValueParameter<BoolValue>(CheckExpressionsWithIntervalArithmeticParameterName, "Switch that determines if the interpreter checks the validity of expressions with interval arithmetic before evaluating the expression.", new BoolValue(false)));
       Parameters.Add(new ValueParameter<IntValue>(EvaluatedSolutionsParameterName, "A counter for the total number of solutions the interpreter has evaluated", new IntValue(0)));
+      interpreter = new SymbolicDataAnalysisExpressionTreeInterpreter();
     }
 
@@ -96 +101 @@
       if (!Parameters.ContainsKey(EvaluatedSolutionsParameterName))
         Parameters.Add(new ValueParameter<IntValue>(EvaluatedSolutionsParameterName, "A counter for the total number of solutions the interpreter has evaluated", new IntValue(0)));
+      if (interpreter == null) interpreter = new SymbolicDataAnalysisExpressionTreeInterpreter();
     }
 
@@ -122 +128 @@
     }
 
+    private static LinearInstruction[] GetPrefixSequence(LinearInstruction[] code, int startIndex) {
+      var instr = code[startIndex];
+      var list = new List<LinearInstruction> { instr };
+      for (int i = 0; i != instr.nArguments; ++i) {
+        list.AddRange(GetPrefixSequence(code, instr.childIndex + i));
+      }
+      return list.ToArray();
+    }
+
     private static void PrepareInstructions(LinearInstruction[] code, Dataset dataset) {
-      for (int i = code.Length - 1; i >= 0; --i) {
+      Action<LinearInstruction> setSkip = null;
+      setSkip = instruction => {
+        instruction.skip = true;
+        for (int j = 0; j != instruction.nArguments; ++j) {
+          setSkip(code[instruction.childIndex + j]);
+        }
+      };
+
+      for (int i = 0; i != code.Length; ++i) {
         var instr = code[i];
         #region opcode switch
@@ -130 +153 @@
               var constTreeNode = (ConstantTreeNode)instr.dynamicNode;
               instr.value = constTreeNode.Value;
+              instr.skip = true; // the value is already set so this instruction should be skipped in the evaluation phase
             }
             break;
@@ -147 +171 @@
             }
             break;
+          case OpCodes.TimeLag:
+          case OpCodes.Integral:
+          case OpCodes.Derivative: {
+              for (int j = 0; j != instr.nArguments; ++j) {
+                setSkip(code[instr.childIndex + j]);
+              }
+            }
+            break;
         }
         #endregion
@@ -152 +184 @@
     }
 
-    private static double Evaluate(Dataset dataset, ref int row, LinearInstruction[] code) {
+    private double Evaluate(Dataset dataset, ref int row, LinearInstruction[] code) {
       for (int i = code.Length - 1; i >= 0; --i) {
-        if (code[i].opCode == OpCodes.Constant) continue;
+        if (code[i].skip) continue;
         #region opcode switch
         var instr = code[i];
@@ -304 +336 @@
               else {
                 alglib.sinecosineintegrals(x, out si, out ci);
-                instr.value = si;
+                instr.value = ci;
               }
             }
@@ -431 +463 @@
             }
             break;
+          case OpCodes.TimeLag:
+          case OpCodes.Integral:
+          case OpCodes.Derivative: {
+              if (instr.iArg0 == null) {
+                instr.iArg0 = GetPrefixSequence(code, i); ;
+              }
+              var interpreterState = new InterpreterState((LinearInstruction[])instr.iArg0, 0);
+              instr.value = interpreter.Evaluate(dataset, ref row, interpreterState);
+            }
+            break;
           default:
             var errorText = string.Format("The {0} symbol is not supported by the linear interpreter. To support this symbol, please use another interpreter.", instr.dynamicNode.Symbol.Name);