[GiNaC-devel] Numerical integration in GiNaC

[Richard B. Kreckel]

Hi Chris!

I was playing with the adaptive Simpson function and found some problems
evaluating numerically a (well-behaved) integral.  I've a couple of
questions:

1) Why the decision to *not* couple the precision to the value of Digits?
   Coupling it to Digits would make it available in Ginsh.  Wouldn't it do
   the right thing in almost all cases?

2) Are variations of Simpson's rule really attractive for numerical
   evaluation?  Have you considered/tried something funkier (maybe
   Romberg's method) or is that a particularly smart implementation?

3) Probably an oversight: the last argument to the adaptivesimpson
   function is unused.  Instead, the code always uses the default.

Best wishes
   -richy.
-- 
Richard B. Kreckel
<http://www.ginac.de/~kreckel/>