rational function. Just assembling it is boring...
for (int i=1; i<=10; ++i)
s += i*y*pow(t,i)/pow(y + i*t,i);
- if (s.nops()!=10) {
+ s = s.normal();
+
+ if (s.subs(t==0)!=0) {
clog << "something very strange happened" << endl;
return 1;
}
timer rolex;
double time = .0;
- cout << "timing Lewis-Wester test D (sum of rational fcns)" << flush;
- clog << "-------Lewis-Wester test D (sum of rational fcns)" << endl;
+ cout << "timing Lewis-Wester test D (normalized sum of rational fcns)" << flush;
+ clog << "-------Lewis-Wester test D (normalized sum of rational fcns)" << endl;
rolex.start();
// correct for very small times:
} else {
cout << " failed ";
}
- cout << int(100000*(time/count))*0.00001 << 's' << endl;
+ cout << int(1000*(time/count))*0.001 << 's' << endl;
return result;
}
for (int i=1; i<=10; ++i)
s += i*y*pow(t,i)/pow(y + abs(5-i)*t,i);
- if (s.nops()!=10) {
+ s = s.normal();
+
+ if (s.subs(t==0)!=0) {
clog << "something very strange happened" << endl;
return 1;
}
timer rolex;
double time = .0;
- cout << "timing Lewis-Wester test E (sum of rational fcns)" << flush;
- clog << "-------Lewis-Wester test E (sum of rational fcns)" << endl;
+ cout << "timing Lewis-Wester test E (normalized sum of rational fcns)" << flush;
+ clog << "-------Lewis-Wester test E (normalized sum of rational fcns)" << endl;
rolex.start();
// correct for very small times:
} else {
cout << " failed ";
}
- cout << int(100000*(time/count))*0.00001 << 's' << endl;
+ cout << int(1000*(time/count))*0.001 << 's' << endl;
return result;
}
(no output)
-------Lewis-Wester test C (gcd of big integers)
(no output)
--------Lewis-Wester test D (sum of rational fcns)
+-------Lewis-Wester test D (normalized sum of rational fcns)
(no output)
--------Lewis-Wester test E (sum of rational fcns)
+-------Lewis-Wester test E (normalized sum of rational fcns)
(no output)
-------Lewis-Wester test F (gcd of 2-var polys)
(no output)