X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Fexam_misc.cpp;h=849dc5adabd666ba1009e68ca18922fe926bc06b;hp=d703a93ada5eca54a7fe1966544d1ae80d2276f3;hb=da2ce481c6a4c5f79e3c2d888bbab45b9a476bd4;hpb=0cf43f3096cbcfc7472ff9c8927c6eb74f2eeb8c diff --git a/check/exam_misc.cpp b/check/exam_misc.cpp index d703a93a..849dc5ad 100644 --- a/check/exam_misc.cpp +++ b/check/exam_misc.cpp @@ -3,7 +3,7 @@ */ /* - * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2002 Johannes Gutenberg University Mainz, Germany * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by @@ -199,6 +199,33 @@ static unsigned exam_operator_semantics(void) return result; } +/* This checks whether subs() works as intended in some special cases. */ +static unsigned exam_subs(void) +{ + unsigned result = 0; + symbol x("x"); + ex e1, e2; + + // This used to fail in GiNaC 1.0.5 because it first substituted + // x+1 -> (x-1)+1 -> x, and then substituted again x -> x-1, giving + // the wrong result + e1 = x+1; + e2 = e1.subs(x == x-1); + if (!e2.is_equal(x)) { + clog << "(x+1).subs(x==x-1) erroneously returned " << e2 << " instead of x" << endl; + ++result; + } + + e1 = sin(1+sin(x)); + e2 = e1.subs(sin(wild()) == cos(wild())); + if (!e2.is_equal(cos(1+cos(x)))) { + clog << "sin(1+sin(x)).subs(sin($1)==cos($1)) erroneously returned " << e2 << " instead of cos(1+cos(x))" << endl; + ++result; + } + + return result; +} + unsigned exam_misc(void) { unsigned result = 0; @@ -211,6 +238,7 @@ unsigned exam_misc(void) result += exam_expand_power(); cout << '.' << flush; result += exam_sqrfree(); cout << '.' << flush; result += exam_operator_semantics(); cout << '.' << flush; + result += exam_subs(); cout << '.' << flush; if (!result) { cout << " passed " << endl;