On Jul 27, 2:16 am, Bill Hart <[EMAIL PROTECTED]> wrote: > Alternatively you could just use the implementation here: > > http://www.ark.in-berlin.de/part.c
In fact, that was my attempt of a C implementation of this Pari script: * Psi(n, q) = local(a, b, c); a=sqrt(2/3)*Pi/q; b=n-1/24; c=sqrt(b); * (sqrt(q)/(2*sqrt(2)*b*Pi))*(a*cosh(a*c)-(sinh(a*c)/c)) * L(n, q) = if(q==1,1,sum(h=1,q-1,if(gcd(h,q)>1,0,cos((g(h,q)-2*h*n)*Pi/q)))) * g(h, q) = if(q<3,0,sum(k=1,q-1,k*(frac(h*k/q)-1/2))) * part(n) = round(sum(q=1,5 + 0.24*sqrt(n),L(n,q)*Psi(n,q))) which is also of my design. The C version then served as base for the implementation of numbpart() within Pari proper. So, they all use the same algorithm (modulo some changes by Karim Belabas etc). Just to get the history straight, Ralf Stephan --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com 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://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~----------~----~----~----~------~----~------~--~---
