On Feb 21, 7:58 am, Fredrik Johansson <[email protected]>
wrote:
> Hi,
Hi Fredrik,
> Looking around, it seems Sage does not yet implement harmonic numbers
> (except via SymPy)? If anyone is interested, I benchmarked a few
> different algorithms and blogged about it
> here:http://fredrik-j.blogspot.com/2009/02/how-not-to-compute-harmonic-num...
I am not aware of any other way this would be implemented in Sage. One
could call GAP, but as you point out it is rather slow.
> Besides (generalized) harmonic numbers, there is potential for
> asymptotically fast computation of Stirling numbers with large n and
> small k. GAP (which Sage uses) takes 3.5 seconds to compute
> S1(10000,2) on my system while algorithm #5 from my blog post takes
> 0.05 seconds.
Nice.
> I'm curious if anyone here knows of faster algorithms for harmonic
> numbers / Stirling numbers? This list seems like the right place to
> ask :-)
I don't know, but hopefully someone will see the bump and respond.
Paul Zimmermann might be the right guy, but he doesn't read sage-devel
very intensely AFAIK, so I am CCing him.
> Fredrik
Cheers,
Michael
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at http://groups.google.com/group/sage-devel
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
