Alternative bernoulli(const numeric &)

Markus Nullmeier markus.nullmeier at urz.uni-heidelberg.de
Sat Jan 19 00:49:53 CET 2002


> Sure, they are notoriously untractable.  But was there really such a
> limit in the old code?  I was under the impression that I once had it
> compute B_{30000} but I might be wrong...

Oh, sorry about this confusion. I overlooked some points while reading
the 1.0.3 code. I rather should have said that the (long) value of (m*n)
overflows 2^32-1 for B_{32768}. I'm guessing that this should impact the
calculated values.



> I do not think the difference in times will be worth the effort, but I
> haven't tried.  Making it safe such that somebody who wants to see it
> break down would have to let it run for a week or so is more important.

Make-safe patch to follow ... :)

Markus



More information about the GiNaC-devel mailing list