Free cookie consent management tool by TermsFeed Policy Generator

source: branches/2994-AutoDiffForIntervals/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis.Symbolic-3.4/DeriveTest.cs @ 17217

Last change on this file since 17217 was 17217, checked in by gkronber, 5 years ago

#2994: more test cases to understand specific problems when optimizing with constraints

File size: 8.0 KB
Line 
1#region License Information
2/* HeuristicLab
3 * Copyright (C) Heuristic and Evolutionary Algorithms Laboratory (HEAL)
4 *
5 * This file is part of HeuristicLab.
6 *
7 * HeuristicLab is free software: you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation, either version 3 of the License, or
10 * (at your option) any later version.
11 *
12 * HeuristicLab is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15 * GNU General Public License for more details.
16 *
17 * You should have received a copy of the GNU General Public License
18 * along with HeuristicLab. If not, see <http://www.gnu.org/licenses/>.
19 */
20#endregion
21
22
23using System;
24using HeuristicLab.Encodings.SymbolicExpressionTreeEncoding;
25using Microsoft.VisualStudio.TestTools.UnitTesting;
26namespace HeuristicLab.Problems.DataAnalysis.Symbolic.Tests {
27
28
29  [TestClass]
30  public class DerivativeTest {
31    [TestMethod]
32    [TestCategory("Problems.DataAnalysis.Symbolic")]
33    [TestProperty("Time", "short")]
34    public void DeriveExpressions() {
35      var formatter = new InfixExpressionFormatter();
36      var parser = new InfixExpressionParser();
37      Assert.AreEqual("0", Derive("3", "x"));
38      Assert.AreEqual("1", Derive("x", "x"));
39      Assert.AreEqual("10", Derive("10*x", "x"));
40      Assert.AreEqual("10", Derive("x*10", "x"));
41      Assert.AreEqual("(2*'x')", Derive("x*x", "x"));
42      Assert.AreEqual("((('x' * 'x') * 2) + ('x' * 'x'))", Derive("x*x*x", "x")); // simplifier does not merge (x*x)*2 + x*x  to 3*x*x
43      Assert.AreEqual("0", Derive("10*x", "y"));
44      Assert.AreEqual("20", Derive("10*x+20*y", "y"));
45      Assert.AreEqual("6", Derive("2*3*x", "x"));
46      Assert.AreEqual("(10*'y')", Derive("10*x*y+20*y", "x"));
47      Assert.AreEqual("(1 / (SQR('x') * (-1)))", Derive("1/x", "x"));
48      Assert.AreEqual("('y' / (SQR('x') * (-1)))", Derive("y/x", "x"));
49      Assert.AreEqual("((((-2*'x') + (-1)) * ('a' + 'b')) / SQR(('x' + ('x' * 'x'))))",
50        Derive("(a+b)/(x+x*x)", "x"));
51      Assert.AreEqual("((((-2*'x') + (-1)) * ('a' + 'b')) / SQR(('x' + SQR('x'))))", Derive("(a+b)/(x+SQR(x))", "x"));
52      Assert.AreEqual("EXP('x')", Derive("exp(x)", "x"));
53      Assert.AreEqual("(EXP((3*'x')) * 3)", Derive("exp(3*x)", "x"));
54      Assert.AreEqual("(1 / 'x')", Derive("log(x)", "x"));
55      Assert.AreEqual("(1 / 'x')", Derive("log(3*x)", "x"));   // 3 * 1/(3*x)
56      Assert.AreEqual("(1 / ('x' + (0.333333333333333*'y')))", Derive("log(3*x+y)", "x"));  // simplifier does not try to keep fractions
57      Assert.AreEqual("(1 / (SQRT(((3*'x') + 'y')) * 0.666666666666667))", Derive("sqrt(3*x+y)", "x"));   // 3 / (2 * sqrt(3*x+y)) = 1 / ((2/3) * sqrt(3*x+y))
58      Assert.AreEqual("(COS((3*'x')) * 3)", Derive("sin(3*x)", "x"));
59      Assert.AreEqual("(SIN((3*'x')) * (-3))", Derive("cos(3*x)", "x"));
60      Assert.AreEqual("(1 / (SQR(COS((3*'x'))) * 0.333333333333333))", Derive("tan(3*x)", "x")); // diff(tan(f(x)), x) = 1.0 / cos²(f(x)), simplifier puts constant factor into the denominator
61
62      Assert.AreEqual("((9*'x') / ABS((3*'x')))", Derive("abs(3*x)", "x"));
63      Assert.AreEqual("(SQR('x') * 3)", Derive("cube(x)", "x"));
64      Assert.AreEqual("(1 / (SQR(CUBEROOT('x')) * 3))", Derive("cuberoot(x)", "x"));
65
66      Assert.AreEqual("0", Derive("(a+b)/(x+SQR(x))", "y")); // df(a,b,x) / dy = 0
67
68
69      Assert.AreEqual("('a' * 'b' * 'c')", Derive("a*b*c*d", "d"));
70      Assert.AreEqual("('a' / ('b' * 'c' * SQR('d') * (-1)))", Derive("a/b/c/d", "d"));
71
72      Assert.AreEqual("('x' * ((SQR(TANH(SQR('x'))) * (-1)) + 1) * 2)", Derive("tanh(sqr(x))", "x")); // (2*'x'*(1 - SQR(TANH(SQR('x'))))
73
74      {
75        // special case: Inv(x) using only one argument to the division symbol
76        // f(x) = 1/x
77        var root = new ProgramRootSymbol().CreateTreeNode();
78        var start = new StartSymbol().CreateTreeNode();
79        var div = new Division().CreateTreeNode();
80        var varNode = (VariableTreeNode)(new Variable().CreateTreeNode());
81        varNode.Weight = 1.0;
82        varNode.VariableName = "x";
83        div.AddSubtree(varNode);
84        start.AddSubtree(div);
85        root.AddSubtree(start);
86        var t = new SymbolicExpressionTree(root);
87        Assert.AreEqual("(1 / (SQR('x') * (-1)))",
88          formatter.Format(DerivativeCalculator.Derive(t, "x")));
89      }
90
91      {
92        // special case: multiplication with only one argument
93        var root = new ProgramRootSymbol().CreateTreeNode();
94        var start = new StartSymbol().CreateTreeNode();
95        var mul = new Multiplication().CreateTreeNode();
96        var varNode = (VariableTreeNode)(new Variable().CreateTreeNode());
97        varNode.Weight = 3.0;
98        varNode.VariableName = "x";
99        mul.AddSubtree(varNode);
100        start.AddSubtree(mul);
101        root.AddSubtree(start);
102        var t = new SymbolicExpressionTree(root);
103        Assert.AreEqual("3",
104          formatter.Format(DerivativeCalculator.Derive(t, "x")));
105      }
106
107      {
108        // division with multiple arguments
109        // div(x, y, z) is interpreted as (x / y) / z
110        var root = new ProgramRootSymbol().CreateTreeNode();
111        var start = new StartSymbol().CreateTreeNode();
112        var div = new Division().CreateTreeNode();
113        var varNode1 = (VariableTreeNode)(new Variable().CreateTreeNode());
114        varNode1.Weight = 3.0;
115        varNode1.VariableName = "x";
116        var varNode2 = (VariableTreeNode)(new Variable().CreateTreeNode());
117        varNode2.Weight = 4.0;
118        varNode2.VariableName = "y";
119        var varNode3 = (VariableTreeNode)(new Variable().CreateTreeNode());
120        varNode3.Weight = 5.0;
121        varNode3.VariableName = "z";
122        div.AddSubtree(varNode1); div.AddSubtree(varNode2); div.AddSubtree(varNode3);
123        start.AddSubtree(div);
124        root.AddSubtree(start);
125        var t = new SymbolicExpressionTree(root);
126
127        Assert.AreEqual("(('y' * 'z' * 60) / (SQR('y') * SQR('z') * 400))", // actually 3 / (4y  5z) but simplifier is not smart enough to cancel numerator and denominator
128                                                                            // 60 y z / y² z² 20² == 6 / y z 40 == 3 / y z 20
129          formatter.Format(DerivativeCalculator.Derive(t, "x")));
130        Assert.AreEqual("(('x' * 'z' * (-60)) / (SQR('y') * SQR('z') * 400))", // actually 3x * -(4 5 z) / (4y 5z)² = -3x / (20 y² z)
131                                                                               // -3 4 5 x z / 4² y² 5² z² = -60 x z / 20² z² y² ==    -60 x z / y² z² 20²
132          formatter.Format(DerivativeCalculator.Derive(t, "y")));
133        Assert.AreEqual("(('x' * 'y' * (-60)) / (SQR('y') * SQR('z') * 400))",
134          formatter.Format(DerivativeCalculator.Derive(t, "z")));
135      }
136
137      // more cases from applications
138      Assert.AreEqual("(2*'R')", Derive("R*R", "R"));
139      Assert.AreEqual("(2*'r')", Derive("sqr(r)", "r"));
140
141      Assert.AreEqual("('R' / (SQR('r') * (-1)))", Derive("R/r", "r"));
142      Assert.AreEqual("(('r' * 'R' * (-2)) / (SQR('r') * SQR('r')))", Derive("R/(r*r)", "r"));
143      Assert.AreEqual("(('r' * 'R' * (-2)) / SQR(SQR('r')))", Derive("R/sqr(r)", "r"));
144
145      Assert.AreEqual("(('r' * ('R' * 'R') * (-2)) / SQR(SQR('r')))", Derive("R*R/sqr(r)", "r"));
146      Assert.AreEqual("(('r' * SQR('R') * (-2)) / SQR(SQR('r')))", Derive("sqr(R)/sqr(r)", "r"));
147      Assert.AreEqual("((('R' * 'R') * 2) / ('r' * SQR('r') * (-1)))", Derive("sqr(R/r)", "r"));
148      Assert.AreEqual("((('R' * 'R') / ('r' * SQR('r') * (-1))) + (('R' * 'R') / ('r' * SQR('r') * (-1))))", Derive("(R/r)*(R/r)", "r"));
149    }
150
151    private string Derive(string expr, string variable) {
152      var parser = new InfixExpressionParser();
153      var formatter = new InfixExpressionFormatter();
154
155      var t = parser.Parse(expr);
156      var tPrime = DerivativeCalculator.Derive(t, variable);
157
158      return formatter.Format(tPrime);
159    }
160  }
161}
Note: See TracBrowser for help on using the repository browser.