March 2008 - GiNaC-devel

[PATCH 7/8] Implicit conversion from cl_N to numeric considered harmful.
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Finally, mark the numeric(const cl_N&) ctor as explicit. This allows one to mix the code using GiNaC and CLN, i.e. cl_N x, y; // initialize them cl_N z = sin(x) + y*exp(y); symbol a("a"); ex e = exp(a); --- ginac/numeric.h | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) diff --git a/ginac/numeric.h b/ginac/numeric.h index f52ffc0..7363c9e 100644 --- a/ginac/numeric.h +++ b/ginac/numeric.h @@ -184,7 +184,7 @@ public: const numeric denom() const; int int_length() const; // converting routines for interfacing with CLN: - numeric(const cln::cl_N &z); + explicit numeric(const cln::cl_N &z); protected: void print_numeric(const print_context & c, const char *par_open, const char *par_close, const char *imag_sym, const char *mul_sym, unsigned level) const; -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[PATCH 6/8] bernoulli, fibonacci, iquo: explicitly convert return values to numeric.
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Implicit conversion from cl_N to numeric considered harmful, part 6. --- ginac/numeric.cpp | 10 +++++----- 1 files changed, 5 insertions(+), 5 deletions(-) diff --git a/ginac/numeric.cpp b/ginac/numeric.cpp index 704854c..4cbaa7d 100644 --- a/ginac/numeric.cpp +++ b/ginac/numeric.cpp @@ -2161,7 +2161,7 @@ const numeric bernoulli(const numeric &nn) next_r = 4; } if (n<next_r) - return results[n/2-1]; + return numeric(results[n/2-1]); results.reserve(n/2); for (unsigned p=next_r; p<=n; p+=2) { @@ -2186,7 +2186,7 @@ const numeric bernoulli(const numeric &nn) results.push_back(-b/(p+1)); } next_r = n+2; - return results[n/2-1]; + return numeric(results[n/2-1]); } @@ -2243,9 +2243,9 @@ const numeric fibonacci(const numeric &n) if (n.is_even()) // Here we don't use the squaring formula because one multiplication // is cheaper than two squarings. - return u * ((v << 1) - u); + return numeric(u * ((v << 1) - u)); else - return cln::square(u) + cln::square(v); + return numeric(cln::square(u) + cln::square(v)); } @@ -2363,7 +2363,7 @@ const numeric iquo(const numeric &a, const numeric &b, numeric &r) const cln::cl_I_div_t rem_quo = cln::truncate2(cln::the<cln::cl_I>(a.to_cl_N()), cln::the<cln::cl_I>(b.to_cl_N())); r = numeric(rem_quo.remainder); - return rem_quo.quotient; + return numeric(rem_quo.quotient); } else { r = *_num0_p; return *_num0_p; -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[PATCH 5/8] tgamma, lgamma: convert to take cl_N as an argument; provide wrappers for GiNaC::numeric
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Implicit conversion from cl_N to numeric considered harmful, part 5. --- ginac/numeric.cpp | 58 +++++++++++++++++++++++++++++++++++++--------------- 1 files changed, 41 insertions(+), 17 deletions(-) diff --git a/ginac/numeric.cpp b/ginac/numeric.cpp index a06be42..704854c 100644 --- a/ginac/numeric.cpp +++ b/ginac/numeric.cpp @@ -1962,24 +1962,34 @@ lanczos_coeffs::lanczos_coeffs() coeffs[3].swap(coeffs_120); } +static const cln::float_format_t guess_precision(const cln::cl_N& x) +{ + cln::float_format_t prec = cln::default_float_format; + if (!instanceof(realpart(x), cln::cl_RA_ring)) + prec = cln::float_format(cln::the<cln::cl_F>(realpart(x))); + if (!instanceof(imagpart(x), cln::cl_RA_ring)) + prec = cln::float_format(cln::the<cln::cl_F>(imagpart(x))); + return prec; +} /** The Gamma function. * Use the Lanczos approximation. If the coefficients used here are not * sufficiently many or sufficiently accurate, more can be calculated * using the program doc/examples/lanczos.cpp. In that case, be sure to * read the comments in that file. */ -const numeric lgamma(const numeric &x) +const cln::cl_N lgamma_(const cln::cl_N &x) { + cln::float_format_t prec = guess_precision(x); lanczos_coeffs lc; - if (lc.sufficiently_accurate(Digits)) { - cln::cl_N pi_val = cln::pi(cln::default_float_format); - if (x.real() < 0.5) - return log(pi_val) - log(sin(pi_val*x.to_cl_N())) - - lgamma(numeric(1).sub(x)); - cln::cl_N A = lc.calc_lanczos_A(x.to_cl_N()); - cln::cl_N temp = x.to_cl_N() + lc.get_order() - cln::cl_N(1)/2; + if (lc.sufficiently_accurate(prec)) { + cln::cl_N pi_val = cln::pi(prec); + if (realpart(x) < 0.5) + return cln::log(pi_val) - cln::log(sin(pi_val*x)) + - lgamma_(1 - x); + cln::cl_N A = lc.calc_lanczos_A(x); + cln::cl_N temp = x + lc.get_order() - cln::cl_N(1)/2; cln::cl_N result = log(cln::cl_I(2)*pi_val)/2 - + (x.to_cl_N()-cln::cl_N(1)/2)*log(temp) + + (x-cln::cl_N(1)/2)*log(temp) - temp + log(A); return result; @@ -1988,17 +1998,25 @@ const numeric lgamma(const numeric &x) throw dunno(); } -const numeric tgamma(const numeric &x) +const numeric lgamma(const numeric &x) +{ + const cln::cl_N x_ = x.to_cl_N(); + const cln::cl_N result = lgamma_(x_); + return numeric(result); +} + +const cln::cl_N tgamma_(const cln::cl_N &x) { + cln::float_format_t prec = guess_precision(x); lanczos_coeffs lc; - if (lc.sufficiently_accurate(Digits)) { - cln::cl_N pi_val = cln::pi(cln::default_float_format); - if (x.real() < 0.5) - return pi_val/(sin(pi_val*x))/(tgamma(numeric(1).sub(x)).to_cl_N()); - cln::cl_N A = lc.calc_lanczos_A(x.to_cl_N()); - cln::cl_N temp = x.to_cl_N() + lc.get_order() - cln::cl_N(1)/2; + if (lc.sufficiently_accurate(prec)) { + cln::cl_N pi_val = cln::pi(prec); + if (realpart(x) < 0.5) + return pi_val/(cln::sin(pi_val*x))/tgamma_(1 - x); + cln::cl_N A = lc.calc_lanczos_A(x); + cln::cl_N temp = x + lc.get_order() - cln::cl_N(1)/2; cln::cl_N result - = sqrt(cln::cl_I(2)*pi_val) * expt(temp, x.to_cl_N()-cln::cl_N(1)/2) + = sqrt(cln::cl_I(2)*pi_val) * expt(temp, x - cln::cl_N(1)/2) * exp(-temp) * A; return result; } @@ -2006,6 +2024,12 @@ const numeric tgamma(const numeric &x) throw dunno(); } +const numeric tgamma(const numeric &x) +{ + const cln::cl_N x_ = x.to_cl_N(); + const cln::cl_N result = tgamma_(x_); + return numeric(result); +} /** The psi function (aka polygamma function). * This is only a stub! */ -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[PATCH 4/8] Li2, zeta, sqrt, abs, gcd, etc.: explicitly convert return value to numeric.
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Implicit conversion from cl_N to numeric considered harmful, part 4. --- ginac/numeric.cpp | 61 ++++++++++++++++++++++++++++++---------------------- 1 files changed, 35 insertions(+), 26 deletions(-) diff --git a/ginac/numeric.cpp b/ginac/numeric.cpp index 0f1b394..a06be42 100644 --- a/ginac/numeric.cpp +++ b/ginac/numeric.cpp @@ -1611,24 +1611,24 @@ static cln::cl_N Li2_projection(const cln::cl_N &x, return Li2_series(x, prec); } + /** Numeric evaluation of Dilogarithm. The domain is the entire complex plane, * the branch cut lies along the positive real axis, starting at 1 and * continuous with quadrant IV. * * @return arbitrary precision numerical Li2(x). */ -const numeric Li2(const numeric &x) +const cln::cl_N Li2_(const cln::cl_N& value) { - if (x.is_zero()) - return *_num0_p; + if (zerop(value)) + return 0; // what is the desired float format? // first guess: default format cln::float_format_t prec = cln::default_float_format; - const cln::cl_N value = x.to_cl_N(); // second guess: the argument's format - if (!x.real().is_rational()) + if (!instanceof(realpart(value), cln::cl_RA_ring)) prec = cln::float_format(cln::the<cln::cl_F>(cln::realpart(value))); - else if (!x.imag().is_rational()) + else if (!instanceof(imagpart(value), cln::cl_RA_ring)) prec = cln::float_format(cln::the<cln::cl_F>(cln::imagpart(value))); if (value==1) // may cause trouble with log(1-x) @@ -1640,7 +1640,16 @@ const numeric Li2(const numeric &x) - cln::zeta(2, prec) - Li2_projection(cln::recip(value), prec)); else - return Li2_projection(x.to_cl_N(), prec); + return Li2_projection(value, prec); +} + +const numeric Li2(const numeric &x) +{ + const cln::cl_N x_ = x.to_cl_N(); + if (zerop(x_)) + return *_num0_p; + const cln::cl_N result = Li2_(x_); + return numeric(result); } @@ -1656,7 +1665,7 @@ const numeric zeta(const numeric &x) if (x.is_real()) { const int aux = (int)(cln::double_approx(cln::the<cln::cl_R>(x.to_cl_N()))); if (cln::zerop(x.to_cl_N()-aux)) - return cln::zeta(aux); + return numeric(cln::zeta(aux)); } throw dunno(); } @@ -2219,7 +2228,7 @@ const numeric fibonacci(const numeric &n) /** Absolute value. */ const numeric abs(const numeric& x) { - return cln::abs(x.to_cl_N()); + return numeric(cln::abs(x.to_cl_N())); } @@ -2233,8 +2242,8 @@ const numeric abs(const numeric& x) const numeric mod(const numeric &a, const numeric &b) { if (a.is_integer() && b.is_integer()) - return cln::mod(cln::the<cln::cl_I>(a.to_cl_N()), - cln::the<cln::cl_I>(b.to_cl_N())); + return numeric(cln::mod(cln::the<cln::cl_I>(a.to_cl_N()), + cln::the<cln::cl_I>(b.to_cl_N()))); else return *_num0_p; } @@ -2248,8 +2257,8 @@ const numeric smod(const numeric &a, const numeric &b) { if (a.is_integer() && b.is_integer()) { const cln::cl_I b2 = cln::ceiling1(cln::the<cln::cl_I>(b.to_cl_N()) >> 1) - 1; - return cln::mod(cln::the<cln::cl_I>(a.to_cl_N()) + b2, - cln::the<cln::cl_I>(b.to_cl_N())) - b2; + return numeric(cln::mod(cln::the<cln::cl_I>(a.to_cl_N()) + b2, + cln::the<cln::cl_I>(b.to_cl_N())) - b2); } else return *_num0_p; } @@ -2267,8 +2276,8 @@ const numeric irem(const numeric &a, const numeric &b) if (b.is_zero()) throw std::overflow_error("numeric::irem(): division by zero"); if (a.is_integer() && b.is_integer()) - return cln::rem(cln::the<cln::cl_I>(a.to_cl_N()), - cln::the<cln::cl_I>(b.to_cl_N())); + return numeric(cln::rem(cln::the<cln::cl_I>(a.to_cl_N()), + cln::the<cln::cl_I>(b.to_cl_N()))); else return *_num0_p; } @@ -2289,8 +2298,8 @@ const numeric irem(const numeric &a, const numeric &b, numeric &q) if (a.is_integer() && b.is_integer()) { const cln::cl_I_div_t rem_quo = cln::truncate2(cln::the<cln::cl_I>(a.to_cl_N()), cln::the<cln::cl_I>(b.to_cl_N())); - q = rem_quo.quotient; - return rem_quo.remainder; + q = numeric(rem_quo.quotient); + return numeric(rem_quo.remainder); } else { q = *_num0_p; return *_num0_p; @@ -2308,8 +2317,8 @@ const numeric iquo(const numeric &a, const numeric &b) if (b.is_zero()) throw std::overflow_error("numeric::iquo(): division by zero"); if (a.is_integer() && b.is_integer()) - return cln::truncate1(cln::the<cln::cl_I>(a.to_cl_N()), - cln::the<cln::cl_I>(b.to_cl_N())); + return numeric(cln::truncate1(cln::the<cln::cl_I>(a.to_cl_N()), + cln::the<cln::cl_I>(b.to_cl_N()))); else return *_num0_p; } @@ -2329,7 +2338,7 @@ const numeric iquo(const numeric &a, const numeric &b, numeric &r) if (a.is_integer() && b.is_integer()) { const cln::cl_I_div_t rem_quo = cln::truncate2(cln::the<cln::cl_I>(a.to_cl_N()), cln::the<cln::cl_I>(b.to_cl_N())); - r = rem_quo.remainder; + r = numeric(rem_quo.remainder); return rem_quo.quotient; } else { r = *_num0_p; @@ -2345,8 +2354,8 @@ const numeric iquo(const numeric &a, const numeric &b, numeric &r) const numeric gcd(const numeric &a, const numeric &b) { if (a.is_integer() && b.is_integer()) - return cln::gcd(cln::the<cln::cl_I>(a.to_cl_N()), - cln::the<cln::cl_I>(b.to_cl_N())); + return numeric(cln::gcd(cln::the<cln::cl_I>(a.to_cl_N()), + cln::the<cln::cl_I>(b.to_cl_N()))); else return *_num1_p; } @@ -2359,8 +2368,8 @@ const numeric gcd(const numeric &a, const numeric &b) const numeric lcm(const numeric &a, const numeric &b) { if (a.is_integer() && b.is_integer()) - return cln::lcm(cln::the<cln::cl_I>(a.to_cl_N()), - cln::the<cln::cl_I>(b.to_cl_N())); + return numeric(cln::lcm(cln::the<cln::cl_I>(a.to_cl_N()), + cln::the<cln::cl_I>(b.to_cl_N()))); else return a.mul(b); } @@ -2376,7 +2385,7 @@ const numeric lcm(const numeric &a, const numeric &b) * where imag(x)>0. */ const numeric sqrt(const numeric &x) { - return cln::sqrt(x.to_cl_N()); + return numeric(cln::sqrt(x.to_cl_N())); } @@ -2386,7 +2395,7 @@ const numeric isqrt(const numeric &x) if (x.is_integer()) { cln::cl_I root; cln::isqrt(cln::the<cln::cl_I>(x.to_cl_N()), &root); - return root; + return numeric(root); } else return *_num0_p; } -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[PATCH 3/8] elementary functions: explicitly convert return values from cl_N to numeric.
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Implicit conversion from cl_N to numeric considered harmful, part 3. --- ginac/numeric.cpp | 34 +++++++++++++++++----------------- 1 files changed, 17 insertions(+), 17 deletions(-) diff --git a/ginac/numeric.cpp b/ginac/numeric.cpp index ab6e8b8..0f1b394 100644 --- a/ginac/numeric.cpp +++ b/ginac/numeric.cpp @@ -1379,7 +1379,7 @@ const numeric I = numeric(cln::complex(cln::cl_I(0),cln::cl_I(1))); * @return arbitrary precision numerical exp(x). */ const numeric exp(const numeric &x) { - return cln::exp(x.to_cl_N()); + return numeric(cln::exp(x.to_cl_N())); } @@ -1392,7 +1392,7 @@ const numeric log(const numeric &x) { if (x.is_zero()) throw pole_error("log(): logarithmic pole",0); - return cln::log(x.to_cl_N()); + return numeric(cln::log(x.to_cl_N())); } @@ -1401,7 +1401,7 @@ const numeric log(const numeric &x) * @return arbitrary precision numerical sin(x). */ const numeric sin(const numeric &x) { - return cln::sin(x.to_cl_N()); + return numeric(cln::sin(x.to_cl_N())); } @@ -1410,7 +1410,7 @@ const numeric sin(const numeric &x) * @return arbitrary precision numerical cos(x). */ const numeric cos(const numeric &x) { - return cln::cos(x.to_cl_N()); + return numeric(cln::cos(x.to_cl_N())); } @@ -1419,7 +1419,7 @@ const numeric cos(const numeric &x) * @return arbitrary precision numerical tan(x). */ const numeric tan(const numeric &x) { - return cln::tan(x.to_cl_N()); + return numeric(cln::tan(x.to_cl_N())); } @@ -1428,7 +1428,7 @@ const numeric tan(const numeric &x) * @return arbitrary precision numerical asin(x). */ const numeric asin(const numeric &x) { - return cln::asin(x.to_cl_N()); + return numeric(cln::asin(x.to_cl_N())); } @@ -1437,7 +1437,7 @@ const numeric asin(const numeric &x) * @return arbitrary precision numerical acos(x). */ const numeric acos(const numeric &x) { - return cln::acos(x.to_cl_N()); + return numeric(cln::acos(x.to_cl_N())); } @@ -1452,7 +1452,7 @@ const numeric atan(const numeric &x) x.real().is_zero() && abs(x.imag()).is_equal(*_num1_p)) throw pole_error("atan(): logarithmic pole",0); - return cln::atan(x.to_cl_N()); + return numeric(cln::atan(x.to_cl_N())); } @@ -1468,8 +1468,8 @@ const numeric atan(const numeric &y, const numeric &x) if (x.is_zero() && y.is_zero()) return *_num0_p; if (x.is_real() && y.is_real()) - return cln::atan(cln::the<cln::cl_R>(x.to_cl_N()), - cln::the<cln::cl_R>(y.to_cl_N())); + return numeric(cln::atan(cln::the<cln::cl_R>(x.to_cl_N()), + cln::the<cln::cl_R>(y.to_cl_N()))); // Compute -I*log((x+I*y)/sqrt(x^2+y^2)) // == -I*log((x+I*y)/sqrt((x+I*y)*(x-I*y))) @@ -1485,7 +1485,7 @@ const numeric atan(const numeric &y, const numeric &x) // x-I*y==0 => y/x==-I, so this is a pole (we have x!=0). throw pole_error("atan(): logarithmic pole",0); } - return cln::complex(0,-1)*cln::log(aux_p/cln::sqrt(aux_p*aux_m)); + return numeric(cln::complex(0,-1)*cln::log(aux_p/cln::sqrt(aux_p*aux_m))); } @@ -1494,7 +1494,7 @@ const numeric atan(const numeric &y, const numeric &x) * @return arbitrary precision numerical sinh(x). */ const numeric sinh(const numeric &x) { - return cln::sinh(x.to_cl_N()); + return numeric(cln::sinh(x.to_cl_N())); } @@ -1503,7 +1503,7 @@ const numeric sinh(const numeric &x) * @return arbitrary precision numerical cosh(x). */ const numeric cosh(const numeric &x) { - return cln::cosh(x.to_cl_N()); + return numeric(cln::cosh(x.to_cl_N())); } @@ -1512,7 +1512,7 @@ const numeric cosh(const numeric &x) * @return arbitrary precision numerical tanh(x). */ const numeric tanh(const numeric &x) { - return cln::tanh(x.to_cl_N()); + return numeric(cln::tanh(x.to_cl_N())); } @@ -1521,7 +1521,7 @@ const numeric tanh(const numeric &x) * @return arbitrary precision numerical asinh(x). */ const numeric asinh(const numeric &x) { - return cln::asinh(x.to_cl_N()); + return numeric(cln::asinh(x.to_cl_N())); } @@ -1530,7 +1530,7 @@ const numeric asinh(const numeric &x) * @return arbitrary precision numerical acosh(x). */ const numeric acosh(const numeric &x) { - return cln::acosh(x.to_cl_N()); + return numeric(cln::acosh(x.to_cl_N())); } @@ -1539,7 +1539,7 @@ const numeric acosh(const numeric &x) * @return arbitrary precision numerical atanh(x). */ const numeric atanh(const numeric &x) { - return cln::atanh(x.to_cl_N()); + return numeric(cln::atanh(x.to_cl_N())); } -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[PATCH 2/8] inifcns_nstdsums.cpp: Lin_numeric takes cl_N as an argument instead of numeric.
by Alexei Sheplyakov 19 Mar '08

19 Mar '08

Implicit conversion from cl_N to numeric considered harmful, part 2. --- ginac/inifcns_nstdsums.cpp | 39 +++++++++++++++++++++++---------------- 1 files changed, 23 insertions(+), 16 deletions(-) diff --git a/ginac/inifcns_nstdsums.cpp b/ginac/inifcns_nstdsums.cpp index ea42a6e..951caa4 100644 --- a/ginac/inifcns_nstdsums.cpp +++ b/ginac/inifcns_nstdsums.cpp @@ -379,15 +379,14 @@ cln::cl_N Li_projection(int n, const cln::cl_N& x, const cln::float_format_t& pr } } - // helper function for classical polylog Li -numeric Lin_numeric(int n, const numeric& x) +const cln::cl_N Lin_numeric(const int n, const cln::cl_N& x) { if (n == 1) { // just a log - return -cln::log(1-x.to_cl_N()); + return -cln::log(1-x); } - if (x.is_zero()) { + if (zerop(x)) { return 0; } if (x == 1) { @@ -398,12 +397,11 @@ numeric Lin_numeric(int n, const numeric& x) // [Kol] (2.22) return -(1-cln::expt(cln::cl_I(2),1-n)) * cln::zeta(n); } - if (abs(x.real()) < 0.4 && abs(abs(x)-1) < 0.01) { - cln::cl_N x_ = ex_to<numeric>(x).to_cl_N(); - cln::cl_N result = -cln::expt(cln::log(x_), n-1) * cln::log(1-x_) / cln::factorial(n-1); + if (cln::abs(realpart(x)) < 0.4 && cln::abs(cln::abs(x)-1) < 0.01) { + cln::cl_N result = -cln::expt(cln::log(x), n-1) * cln::log(1-x) / cln::factorial(n-1); for (int j=0; j<n-1; j++) { - result = result + (S_num(n-j-1, 1, 1) - S_num(1, n-j-1, 1-x_)) - * cln::expt(cln::log(x_), j) / cln::factorial(j); + result = result + (S_num(n-j-1, 1, 1) - S_num(1, n-j-1, 1-x)) + * cln::expt(cln::log(x), j) / cln::factorial(j); } return result; } @@ -411,11 +409,11 @@ numeric Lin_numeric(int n, const numeric& x) // what is the desired float format? // first guess: default format cln::float_format_t prec = cln::default_float_format; - const cln::cl_N value = x.to_cl_N(); + const cln::cl_N value = x; // second guess: the argument's format - if (!x.real().is_rational()) + if (!instanceof(realpart(x), cln::cl_RA_ring)) prec = cln::float_format(cln::the<cln::cl_F>(cln::realpart(value))); - else if (!x.imag().is_rational()) + else if (!instanceof(imagpart(x), cln::cl_RA_ring)) prec = cln::float_format(cln::the<cln::cl_F>(cln::imagpart(value))); // [Kol] (5.15) @@ -441,7 +439,7 @@ numeric Lin_numeric(int n, const numeric& x) cln::cl_N add; for (int j=0; j<n-1; j++) { add = add + (1+cln::expt(cln::cl_I(-1),n-j)) * (1-cln::expt(cln::cl_I(2),1-n+j)) - * Lin_numeric(n-j,1).to_cl_N() * cln::expt(cln::log(-value),j) / cln::factorial(j); + * Lin_numeric(n-j,1) * cln::expt(cln::log(-value),j) / cln::factorial(j); } result = result - add; return result; @@ -1375,12 +1373,18 @@ static ex Li_evalf(const ex& m_, const ex& x_) // classical polylogs if (m_.info(info_flags::posint)) { if (x_.info(info_flags::numeric)) { - return Lin_numeric(ex_to<numeric>(m_).to_int(), ex_to<numeric>(x_)); + int m__ = ex_to<numeric>(m_).to_int(); + const cln::cl_N x__ = ex_to<numeric>(x_).to_cl_N(); + const cln::cl_N result = Lin_numeric(m__, x__); + return numeric(result); } else { // try to numerically evaluate second argument ex x_val = x_.evalf(); if (x_val.info(info_flags::numeric)) { - return Lin_numeric(ex_to<numeric>(m_).to_int(), ex_to<numeric>(x_val)); + int m__ = ex_to<numeric>(m_).to_int(); + const cln::cl_N x__ = ex_to<numeric>(x_val).to_cl_N(); + const cln::cl_N result = Lin_numeric(m__, x__); + return numeric(result); } } } @@ -1494,7 +1498,10 @@ static ex Li_eval(const ex& m_, const ex& x_) } } if (m_.info(info_flags::posint) && x_.info(info_flags::numeric) && !x_.info(info_flags::crational)) { - return Lin_numeric(ex_to<numeric>(m_).to_int(), ex_to<numeric>(x_)); + int m__ = ex_to<numeric>(m_).to_int(); + const cln::cl_N x__ = ex_to<numeric>(x_).to_cl_N(); + const cln::cl_N result = Lin_numeric(m__, x__); + return numeric(result); } return Li(m_, x_).hold(); -- 1.5.4.2 -- All science is either physics or stamp collecting.

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-31-gd7c5ca7
by varg＠theor.jinr.ru 18 Mar '08

18 Mar '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via d7c5ca75b0fc9e58a8676681e3801869e47ad911 (commit) from 090348bcd9ac2104e1e387430aecd5c8b7289633 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit d7c5ca75b0fc9e58a8676681e3801869e47ad911 Author: Alexei Sheplyakov <varg(a)theor.jinr.ru> Date: Tue Mar 18 09:45:25 2008 +0300 move rotate_left() function away from the public header. ----------------------------------------------------------------------- Summary of changes: ginac/ex.h | 7 ------- ginac/utils.h | 7 +++++++ 2 files changed, 7 insertions(+), 7 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-30-g090348b
by varg＠theor.jinr.ru 18 Mar '08

18 Mar '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 090348bcd9ac2104e1e387430aecd5c8b7289633 (commit) via 433e73f28bef23bd1e3135613f91c283aca97c59 (commit) from 3cb110c5ca2330a0a155660f32d42e3d0ebc0db1 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 090348bcd9ac2104e1e387430aecd5c8b7289633 Author: Alexei Sheplyakov <varg(a)theor.jinr.ru> Date: Tue Mar 18 09:24:56 2008 +0300 INSTALL: libreadline is not mandatory; ginsh builds just fine without it. commit 433e73f28bef23bd1e3135613f91c283aca97c59 Author: Alexei Sheplyakov <varg(a)theor.jinr.ru> Date: Tue Mar 18 09:04:12 2008 +0300 update the compiler info; tell users to avoid GCC 4.3.0 ----------------------------------------------------------------------- Summary of changes: INSTALL | 30 ++++++++++++++++++++++-------- 1 files changed, 22 insertions(+), 8 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, ginac_1-4, updated. release_1-4-0-35-gc4714cd
by varg＠theor.jinr.ru 18 Mar '08

18 Mar '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, ginac_1-4 has been updated via c4714cdfcc0b4817ec2c3dbfc464566eaf829eef (commit) via 489045eabf6520f370a8dce5414831e1bf34f6df (commit) from 0eeb19e927b18a2ddc6d062a42f5809443dee77a (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit c4714cdfcc0b4817ec2c3dbfc464566eaf829eef Author: Alexei Sheplyakov <varg(a)theor.jinr.ru> Date: Tue Mar 18 09:24:56 2008 +0300 INSTALL: libreadline is not mandatory; ginsh builds just fine without it. commit 489045eabf6520f370a8dce5414831e1bf34f6df Author: Alexei Sheplyakov <varg(a)theor.jinr.ru> Date: Tue Mar 18 09:04:12 2008 +0300 update the compiler info; tell users to avoid GCC 4.3.0 ----------------------------------------------------------------------- Summary of changes: INSTALL | 30 ++++++++++++++++++++++-------- 1 files changed, 22 insertions(+), 8 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

ginac-1.4.1 fails to compile under gcc-4.3.0
by Sergei Steshenko 07 Mar '08

07 Mar '08

Hello All, yesterday I built the brand new gcc-4.3.0 and later tried to compile ginac-1.4.1 using it. Compilation failed with these command line and error messages: " g++ -DHAVE_CONFIG_H -I. -I.. -I/maxtor5/sergei/AppsFromScratchWD/install/autoconf-2.61/include -I/maxtor5/sergei/AppsFromScratchWD/install/autogen-5.8.7/include -I/maxtor5/sergei/AppsFromScratchWD/install/bash-3.0/include -I/maxtor5/sergei/AppsFromScratchWD/install/bison-2.3/include -I/maxtor5/sergei/AppsFromScratchWD/install/cln-1.2.0/include -I/maxtor5/sergei/AppsFromScratchWD/install/dejagnu-1.4.4/include -I/maxtor5/sergei/AppsFromScratchWD/install/expect-5.44.1/include -I/maxtor5/sergei/AppsFromScratchWD/install/flex-2.5.35/include -I/maxtor5/sergei/AppsFromScratchWD/install/gcc-4.3.0/include -I/maxtor5/sergei/AppsFromScratchWD/install/gmp-4.2.2/include -I/maxtor5/sergei/AppsFromScratchWD/install/gperf-3.0.2/include -I/maxtor5/sergei/AppsFromScratchWD/install/guile-1.8.4/include -I/maxtor5/sergei/AppsFromScratchWD/install/libtool-1.5.26/include -I/maxtor5/sergei/AppsFromScratchWD/install/m4-1.4.10/include -I/maxtor5/sergei/AppsFromScratchWD/install/make-3.80/include -I/maxtor5/sergei/AppsFromScratchWD/install/mpfr-2.3.1/include -I/maxtor5/sergei/AppsFromScratchWD/install/ncurses-5.6/include -I/maxtor5/sergei/AppsFromScratchWD/install/readline-5.1/include -I/maxtor5/sergei/AppsFromScratchWD/install/tcl-8.4.17/include -I/maxtor5/sergei/AppsFromScratchWD/install/texinfo-4.9/include -I/maxtor5/sergei/AppsFromScratchWD/install/tk-8.4.17/include -I/maxtor5/sergei/AppsFromScratchWD/install/cln-1.2.0/include -march=native -mtune=native -mfpmath=sse,387 -MT numeric.lo -MD -MP -MF .deps/numeric.Tpo -c numeric.cpp -fPIC -DPIC -o .libs/numeric.o numeric.cpp:422: error: explicit template specialization cannot have a storage class numeric.cpp:434: error: explicit template specialization cannot have a storage class make[2]: *** [numeric.lo] Error 1 make[2]: Leaving directory `/maxtor5/sergei/AppsFromScratchWD/build/ginac-1.4.1/ginac' make[1]: *** [all-recursive] Error 1 make[1]: Leaving directory `/maxtor5/sergei/AppsFromScratchWD/build/ginac-1.4.1' make: *** [all] Error 2 ". I have successfully compiled ginac many times before, using various gcc versions up to gcc-4.2.3, so it looks like the newest gcc has stricter C++ rules. Thanks, Sergei. Applications From Scratch: http://appsfromscratch.berlios.de/ ____________________________________________________________________________________ Looking for last minute shopping deals? Find them fast with Yahoo! Search. http://tools.search.yahoo.com/newsearch/category.php?category=shopping

2 1