Fix of a bug in sqrfree_parfrac() related to Yun factorisation.

[ginac.git] / check / exam_paranoia.cpp

diff --git a/check/exam_paranoia.cpp b/check/exam_paranoia.cpp
index d96cda59c035b0a37443b9719b370bd56cd75f0f..52466b059b055fe2884636de1aef5a063cfa72c8 100644 (file)
--- a/check/exam_paranoia.cpp
+++ b/check/exam_paranoia.cpp
@@ -636,6 +636,21 @@ unsigned exam_paranoia25()
        return 0;
 }
 
+// Bug in partial fraction expansion
+unsigned exam_paranoia26()
+{
+       symbol x("x");
+       ex ex1=pow(x,4)/(x-1)/4;
+       ex ex2=sqrfree_parfrac(ex1,x);
+       ex e = (ex1-ex2).normal();
+
+       if (! e.is_zero()) {
+               clog << "partial fraction expansion of " << ex1 << " produces error.\n";
+               return 1;
+       }
+       return 0;
+}
+
 unsigned exam_paranoia()
 {
        unsigned result = 0;
@@ -668,6 +683,7 @@ unsigned exam_paranoia()
        result += exam_paranoia23();  cout << '.' << flush;
        result += exam_paranoia24();  cout << '.' << flush;
        result += exam_paranoia25();  cout << '.' << flush;
+       result += exam_paranoia26();  cout << '.' << flush;
        
        return result;
 }