- unsigned result = 0;
-
- cout << "examining several historic failures just out of paranoia" << flush;
- clog << "----------several historic failures:" << endl;
-
- result += exam_paranoia1(); cout << '.' << flush;
- result += exam_paranoia2(); cout << '.' << flush;
- result += exam_paranoia3(); cout << '.' << flush;
- result += exam_paranoia4(); cout << '.' << flush;
- result += exam_paranoia5(); cout << '.' << flush;
- result += exam_paranoia6(); cout << '.' << flush;
- result += exam_paranoia7(); cout << '.' << flush;
- result += exam_paranoia8(); cout << '.' << flush;
- result += exam_paranoia9(); cout << '.' << flush;
- result += exam_paranoia10(); cout << '.' << flush;
- result += exam_paranoia11(); cout << '.' << flush;
- result += exam_paranoia12(); cout << '.' << flush;
-
- if (!result) {
- cout << " passed " << endl;
- clog << "(no output)" << endl;
- } else {
- cout << " failed " << endl;
- }
-
- return result;
+ unsigned result = 0;
+ symbol a("a"), b("b"), c("c");
+
+ ex e = (b*a-c*a)/(4-a);
+ ex d = (c*a-b*a)/(a-4);
+
+ try {
+ ex f = e.normal();
+ if (!(f - d).expand().is_zero()) {
+ clog << "normal(" << e << ") returns " << f
+ << " instead of " << d << endl;
+ ++result;
+ }
+ } catch (const exception &err) {
+ clog << "normal(" << e << ") throws " << err.what() << endl;
+ ++result;
+ }
+ return result;
+}
+
+// A bug introduced on July 19, 2001. quo() and rem() would sometimes call
+// vector::reserve() with a negative argument. Fixed on Dec 20, 2001.
+static unsigned exam_paranoia14()
+{
+ unsigned result = 0;
+ symbol x("x");
+
+ ex q = quo(1, pow(x, 3), x);
+ if (!q.is_zero()) {
+ clog << "quo(1,x^3,x) erroneously returned " << q << " instead of 0\n";
+ ++result;
+ }
+
+ return result;
+}
+
+// Under certain conditions, power::expand_add_2() could produce non-canonical
+// numeric expairs. Fixed on Oct 24, 2002.
+static unsigned exam_paranoia15()
+{
+ unsigned result = 0;
+
+ ex q = (pow(pow(2, numeric(1, 2))*2+1, 2)).expand();
+ // this used to produce "1+4*sqrt(2)+4*2" which would never evaluate
+ // to "9+4*sqrt(2)"
+
+ if (!(q-9-4*pow(2, numeric(1, 2))).is_zero()) {
+ clog << "expand((sqrt(2)*2+1)^2) erroneously returned " << q << " instead of 9-4*sqrt(2)\n";
+ ++result;
+ }
+
+ return result;
+}
+
+unsigned exam_paranoia()
+{
+ unsigned result = 0;
+
+ cout << "examining several historic failures just out of paranoia" << flush;
+ clog << "----------several historic failures:" << endl;
+
+ result += exam_paranoia1(); cout << '.' << flush;
+ result += exam_paranoia2(); cout << '.' << flush;
+ result += exam_paranoia3(); cout << '.' << flush;
+ result += exam_paranoia4(); cout << '.' << flush;
+ result += exam_paranoia5(); cout << '.' << flush;
+ result += exam_paranoia6(); cout << '.' << flush;
+ result += exam_paranoia7(); cout << '.' << flush;
+ result += exam_paranoia8(); cout << '.' << flush;
+ result += exam_paranoia9(); cout << '.' << flush;
+ result += exam_paranoia10(); cout << '.' << flush;
+ result += exam_paranoia11(); cout << '.' << flush;
+ result += exam_paranoia12(); cout << '.' << flush;
+ result += exam_paranoia13(); cout << '.' << flush;
+ result += exam_paranoia14(); cout << '.' << flush;
+ result += exam_paranoia15(); cout << '.' << flush;
+
+ if (!result) {
+ cout << " passed " << endl;
+ clog << "(no output)" << endl;
+ } else {
+ cout << " failed " << endl;
+ }
+
+ return result;