- Two new timimgs that are interesting for optimizing.

[ginac.git] / check / exam_numeric.cpp

diff --git a/check/exam_numeric.cpp b/check/exam_numeric.cpp
index d3c1264a3aa1fdde148631f3f9f18b4fcfd19dfa..e37c565d019e0f4faa0a40b46138b1eb0e3be367 100644 (file)
--- a/check/exam_numeric.cpp
+++ b/check/exam_numeric.cpp
@@ -301,6 +301,25 @@ static unsigned exam_numeric4(void)
     return result;
 }
 
+/* This test examines that simplifications of the form 5^(3/2) -> 5*5^(1/2)
+ * are carried out properly. */
+static unsigned exam_numeric5(void)
+{
+    unsigned result = 0;
+    
+    // A variation of one of Ramanujan's wonderful identities must be
+    // verifiable with very primitive means:
+    ex e1 = pow(1 + pow(3,numeric(1,5)) - pow(3,numeric(2,5)),3);
+    ex e2 = expand(e1 - 10 + 5*pow(3,numeric(3,5)));
+    if (!e2.is_zero()) {
+        clog << "expand((1+3^(1/5)-3^(2/5))^3-10+5*3^(3/5)) returned "
+             << e2 << " instead of 0." << endl;
+        ++result;
+    }
+    
+    return result;
+}
+
 unsigned exam_numeric(void)
 {
     unsigned result = 0;
@@ -312,6 +331,7 @@ unsigned exam_numeric(void)
     result += exam_numeric2();  cout << '.' << flush;
     result += exam_numeric3();  cout << '.' << flush;
     result += exam_numeric4();  cout << '.' << flush;
+    result += exam_numeric5();  cout << '.' << flush;
     
     if (!result) {
         cout << " passed " << endl;