static unsigned inifcns_consist_trans()
{
using GiNaC::asin; using GiNaC::acos;
+ using GiNaC::asinh; using GiNaC::acosh; using GiNaC::atanh;
unsigned result = 0;
symbol x("x");
* exists in closed form and check if it's ok. */
static unsigned inifcns_consist_gamma()
{
+ using GiNaC::tgamma;
unsigned result = 0;
ex e;
static unsigned inifcns_consist_psi()
{
using GiNaC::log;
+ using GiNaC::tgamma;
unsigned result = 0;
symbol x;
int digitsbuf = Digits;
// precision of data
Digits = 22;
- ex prec = 5 * pow(10, -(int)Digits);
+ ex prec = 5 * pow(10, -(ex)Digits);
unsigned result = 0;
static unsigned inifcns_test_HLi()
{
+ using GiNaC::log;
int digitsbuf = Digits;
Digits = 17;
- ex prec = 5 * pow(10, -(int)Digits);
+ ex prec = 5 * pow(10, -(ex)Digits);
numeric almostone("0.999999999999999999");
unsigned result = 0;
for (lst::const_iterator it = res.begin(); it != res.end(); it++) {
Digits = 17;
- ex prec = 5 * pow(10, -(int)Digits);
+ ex prec = 5 * pow(10, -(ex)Digits);
ex diff = abs((*it).evalf());
if (diff > prec) {
clog << *it << " seems to be wrong: " << diff << endl;
}
cout << "." << flush;
Digits = 40;
- prec = 5 * pow(10, -(int)Digits);
+ prec = 5 * pow(10, -(ex)Digits);
diff = abs((*it).evalf());
if (diff > prec) {
clog << *it << " seems to be wrong: " << diff << endl;
{
int digitsbuf = Digits;
Digits = 17;
- ex prec = 5 * pow(10, -(int)Digits);
+ ex prec = 5 * pow(10, -(ex)Digits);
numeric almostone("0.99999999999999999999");
unsigned result = 0;
static unsigned inifcns_test_legacy()
{
Digits = 17;
- ex prec = 5 * pow(10, -(int)Digits);
+ ex prec = 5 * pow(10, -(ex)Digits);
unsigned result = 0;
}
Digits = 100;
- prec = 5 * pow(10, -(int)Digits);
+ prec = 5 * pow(10, -(ex)Digits);
ex x0 = 1.;
x1 = exp(Pi*I/3).evalf();
x2 = exp(2*Pi*I/3).evalf();
// Series exponentiation
static unsigned exam_series4()
{
+ using GiNaC::tgamma;
unsigned result = 0;
ex e, d;
// Series expansion of tgamma(-1)
static unsigned exam_series6()
{
+ using GiNaC::tgamma;
ex e = tgamma(2*x);
ex d = pow(x+1,-1)*numeric(1,4) +
pow(x+1,0)*(numeric(3,4) -
static unsigned exam_series12()
{
using GiNaC::log;
+ using GiNaC::atanh;
unsigned result = 0;
ex e, d;
/* F_ab(a, i, b, j, "x") is a common pattern in all vertex evaluators. */
static ex F_ab(int a, int i, int b, int j, const symbol &x)
{
+ using GiNaC::tgamma;
if ((i==0 && a<=0) || (j==0 && b<=0))
return 0;
else
unsigned result = 0;
symbol x;
- ex myseries = series(tgamma(x),x==0,order);
+ ex myseries = series(GiNaC::tgamma(x),x==0,order);
// compute the last coefficient numerically:
ex last_coeff = myseries.coeff(x,order-1).evalf();
// compute a bound for that coefficient using a variation of the leading
// term in Stirling's formula:
ex bound = exp(-.57721566490153286*(order-1))/(order-1);
- if (abs((last_coeff-pow(-1,order))/bound) > 1) {
+ if (abs((last_coeff-pow(-1,ex(order)))/bound) > 1) {
clog << "The " << order-1
<< "th order coefficient in the power series expansion of tgamma(0) was erroneously found to be "
<< last_coeff << ", violating a simple estimate." << endl;
static ex f_ginac_function(const exprseq &es, int serial)
{
- return function(serial, es).eval(1);
+ return GiNaC::function(serial, es).eval(1);
}
// All registered GiNaC functions