Fredrik, > Thanks for the input. Unfortunately, I don't see how Schönhage's > factorial algorithm can be adapted to harmonic numbers or Stirling > numbers, due to the way the partial products need to be summed. > > A slight improvement to the plain binary splitting method is to factor > out 1/2 to eliminate 25% of the terms. But factoring out larger primes > seems to give rapidly diminished returns, unless there is some clever > way to do it that eludes me.
hint: look at http://hal.inria.fr/inria-00177850 Paul --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---