* functions. */
/*
- * GiNaC Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2001 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
if (errorflag) {
// we don't count each of those errors
clog << "sin(n*Pi) with integer n does not always return exact 0"
- << endl;
+ << endl;
++result;
}
errorflag = false;
for (int n=-10; n<=10; ++n) {
if (!sin((n+numeric(1,2))*Pi).eval().info(info_flags::integer) ||
- !(sin((n+numeric(1,2))*Pi).eval() == numeric(1) ||
- sin((n+numeric(1,2))*Pi).eval() == numeric(-1)))
+ !(sin((n+numeric(1,2))*Pi).eval() == numeric(1) ||
+ sin((n+numeric(1,2))*Pi).eval() == numeric(-1)))
errorflag = true;
}
if (errorflag) {
clog << "sin((n+1/2)*Pi) with integer n does not always return exact {+|-}1"
- << endl;
+ << endl;
++result;
}
argument = n*Pi/60;
if (abs(sin(evalf(argument))-evalf(sin(argument)))>epsilon) {
clog << "sin(" << argument << ") returns "
- << sin(argument) << endl;
+ << sin(argument) << endl;
errorflag = true;
}
}
errorflag = false;
for (int n=-10; n<=10; ++n) {
if (cos((n+numeric(1,2))*Pi).eval() != numeric(0) ||
- !cos((n+numeric(1,2))*Pi).eval().info(info_flags::integer))
+ !cos((n+numeric(1,2))*Pi).eval().info(info_flags::integer))
errorflag = true;
}
if (errorflag) {
clog << "cos((n+1/2)*Pi) with integer n does not always return exact 0"
- << endl;
+ << endl;
++result;
}
errorflag = false;
for (int n=-10; n<=10; ++n) {
if (!cos(n*Pi).eval().info(info_flags::integer) ||
- !(cos(n*Pi).eval() == numeric(1) ||
- cos(n*Pi).eval() == numeric(-1)))
+ !(cos(n*Pi).eval() == numeric(1) ||
+ cos(n*Pi).eval() == numeric(-1)))
errorflag = true;
}
if (errorflag) {
clog << "cos(n*Pi) with integer n does not always return exact {+|-}1"
- << endl;
+ << endl;
++result;
}
argument = n*Pi/60;
if (abs(cos(evalf(argument))-evalf(cos(argument)))>epsilon) {
clog << "cos(" << argument << ") returns "
- << cos(argument) << endl;
+ << cos(argument) << endl;
errorflag = true;
}
}
argument = n*Pi/60;
if (abs(tan(evalf(argument))-evalf(tan(argument)))>epsilon) {
clog << "tan(" << argument << ") returns "
- << tan(argument) << endl;
+ << tan(argument) << endl;
errorflag = true;
}
}
numeric epsilon(double(1e-16));
for (int n=0; n<200; ++n) {
argument = numeric(20.0*rand()/(RAND_MAX+1.0)-10.0)
- + numeric(20.0*rand()/(RAND_MAX+1.0)-10.0)*I;
+ + numeric(20.0*rand()/(RAND_MAX+1.0)-10.0)*I;
if (abs(Li2(pow(argument,2))-2*Li2(argument)-2*Li2(-argument)) > epsilon) {
cout << "Li2(z) at z==" << argument
- << " failed to satisfy Li2(z^2)==2*(Li2(z)+Li2(-z))" << endl;
+ << " failed to satisfy Li2(z^2)==2*(Li2(z)+Li2(-z))" << endl;
errorflag = true;
}
}