Changeset 16738 for branches/2994-AutoDiffForIntervals/Tests/AutoDiffTest.cs

Timestamp:

04/01/19 19:19:14 (5 years ago)

Author:

gkronber

Message:

#2994: another unit test, smaller fixes and checks for invalid parameters

File:

• branches/2994-AutoDiffForIntervals/Tests/AutoDiffTest.cs (modified) (1 diff)

branches/2994-AutoDiffForIntervals/Tests/AutoDiffTest.cs

@@ -189 +189 @@
       }
 
+      {
+
+        // derivatives and intervals for flow psi problem
+        var intervals = new Dictionary<string, Interval>();
+        intervals.Add("x1", new Interval(60.0, 65.0));
+        intervals.Add("x2", new Interval(30.0, 40.0));
+        intervals.Add("x3", new Interval(5.0, 10.0));
+        intervals.Add("x4", new Interval(0.5, 0.8));
+        intervals.Add("x5", new Interval(0.2, 0.5));
+
+        var parser = new InfixExpressionParser();
+        var formatter = new InfixExpressionFormatter();
+
+        var expr = parser.Parse("x1*x2*x5*(1 - sqr(x4/x5)) + x3 * log(x5/x4)");
+
+        var dfdx1 = DerivativeCalculator.Derive(expr, "x1");
+        Assert.AreEqual("('x2' * 'x5' * ((SQR(('x4' / 'x5')) * (-1)) + 1))", formatter.Format(dfdx1));
+        // x2 x5 (1 - sqr(x4/x5))
+
+        var dfdx2 = DerivativeCalculator.Derive(expr, "x2");
+        Assert.AreEqual("('x1' * 'x5' * ((SQR(('x4' / 'x5')) * (-1)) + 1))", formatter.Format(dfdx2));
+        // x1 x5 (1 - sqr(x4/x5))
+
+        var dfdx3 = DerivativeCalculator.Derive(expr, "x3");
+        Assert.AreEqual("LOG(('x5' / 'x4'))", formatter.Format(dfdx3));
+        // log(x5/x4)
+
+        var dfdx4 = DerivativeCalculator.Derive(expr, "x4");
+        Assert.AreEqual("((('x1' * 'x2' * 'x5' * 'x4' * 2) / ('x5' * (-1*'x5'))) + (('x4' * 'x5' * 'x3') / ('x5' * SQR('x4') * (-1))))", formatter.Format(dfdx4));
+        // -2*x1*x2*x5*x4/x5*1/x5 + x3*1/(x5/x4)*x5/sqr(x4)
+
+        var dfdx5 = DerivativeCalculator.Derive(expr, "x5");
+        Assert.AreEqual("((('x4' * 'x3') / ('x5' * 'x4')) + ('x1' * 'x2' * ((SQR(('x4' / 'x5')) * (-1)) + 1)) + (('x1' * 'x2' * 'x5' * ('x4' * 'x4') * 2) / ('x5' * SQR('x5') * 1)))", formatter.Format(dfdx5));
+      }
+
+
     }
   }