[GiNaC-devel] Numerical integration in GiNaC

Richard B. Kreckel kreckel at ginac.de
Mon Sep 12 23:19:50 CEST 2005


On Mon, 12 Sep 2005, Richard B. Kreckel wrote:
> I was playing with the adaptive Simpson function and found some problems
> evaluating numerically a (well-behaved) integral.  I've a couple of
> questions:
[...]

4) The lookup table does not seem to discriminate among different
   precisions.  Presuming the last argument to adaptivesimpson() would
   be used or the precision be coupled to Digits.  Then the lookup map
   would have to hold the precision as well (on the key side, not the
   value side).  Otherwise it becomes impossible to re-run an
   integration with a higher precision goal.

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





More information about the GiNaC-devel mailing list