Fixed bug in multivariate factorization. Fractional numerical content was not

[ginac.git] / check / exam_factor.cpp

diff --git a/check/exam_factor.cpp b/check/exam_factor.cpp
index 8851b7cc2f2c9fedd59d52f91838866f14cf48fe..71a63f24ba68dce756efd210005b93d5b2b6a4c1 100644 (file)
--- a/check/exam_factor.cpp
+++ b/check/exam_factor.cpp
@@ -41,7 +41,7 @@ static unsigned check_factor(const ex& e)
 static unsigned exam_factor1()
 {
        unsigned result = 0;
 static unsigned exam_factor1()
 {
        unsigned result = 0;

-       ex e, d;
+       ex e;

        symbol x("x");
        lst syms;
        syms.append(x);
        symbol x("x");
        lst syms;
        syms.append(x);


@@ -110,7 +110,7 @@ static unsigned exam_factor1()
 static unsigned exam_factor2()
 {
        unsigned result = 0;
 static unsigned exam_factor2()
 {
        unsigned result = 0;

-       ex e, d;
+       ex e;

        symbol x("x"), y("y"), z("z");
        lst syms;
        syms = x, y, z;
        symbol x("x"), y("y"), z("z");
        lst syms;
        syms = x, y, z;


@@ -166,6 +166,23 @@ static unsigned exam_factor2()
        return result;
 }
 
        return result;
 }
 

+static unsigned exam_factor3()
+{
+       unsigned result = 0;
+       ex e;
+       symbol k("k"), n("n");
+       lst syms;
+       syms = k, n;
+       
+       e = ex("1/2*(-3+3*k-n)*(-2+3*k-n)*(-1+3*k-n)", syms);
+       result += check_factor(e);
+
+       e = ex("1/4*(2*k-n)*(-1+2*k-n)", syms);
+       result += check_factor(e);
+
+       return result;
+}
+

 unsigned exam_factor()
 {
        unsigned result = 0;
 unsigned exam_factor()
 {
        unsigned result = 0;


@@ -174,6 +191,7 @@ unsigned exam_factor()
 
        result += exam_factor1(); cout << '.' << flush;
        result += exam_factor2(); cout << '.' << flush;
 
        result += exam_factor1(); cout << '.' << flush;
        result += exam_factor2(); cout << '.' << flush;

+       result += exam_factor3(); cout << '.' << flush;

 
        return result;
 }
 
        return result;
 }