Animated plot and pari stack

[Emmanuel]

Hello,

I would like to reproduce some of the film available at Jeffrey Stopple's webpage with Sage.

See http://www.math.ucsb.edu/~stopple/zeta.html

It seems however (see below) that I need to increase the PARI stacks. How can I do that, knowing that I am not even the administrator of the computer I want to use.

With my best regards,

Emmanuel

MyComputer~/tmp>sage

----------------------------------------------------------------------

| Sage Version 5.4.1, Release Date: 2012-11-15                       |

| Type "notebook()" for the browser-based notebook interface.        |

| Type "help()" for help.                                            |

----------------------------------------------------------------------

sage: i = CDF.0

sage: fz=animate([line([zeta(s+i*t) for t in srange(0,50,0.01)]) for s in srange(0,0.9,0.01)],xmin=-1.5,ymin=-2,xmax=2.5,ymax=2)

  ***   the PARI stack overflows !

  current stack size: 16000000 (15.259 Mbytes)

  [hint] you can increase GP stack with allocatemem()

  ***   Error in the PARI system. End of program.

[Jeroen Demeyer]

On 2013-02-07 23:29, Emmanuel wrote:
> ----------------------------------------------------------------------
> | Sage Version 5.4.1, Release Date: 2012-11-15 |
> | Type "notebook()" for the browser-based notebook interface. |
> | Type "help()" for help. |
> ----------------------------------------------------------------------
> sage: i = CDF.0
> sage: fz=animate([line([zeta(s+i*t) for t in srange(0,50,0.01)]) for s
> in srange(0,0.9,0.01)],xmin=-1.5,ymin=-2,xmax=2.5,ymax=2)
> *** the PARI stack overflows !
> current stack size: 16000000 (15.259 Mbytes)
> [hint] you can increase GP stack with allocatemem()
>
> *** Error in the PARI system. End of program.

Confirmed, this is certainly a bug in Sage:
http://trac.sagemath.org/sage_trac/ticket/14082

[Jeroen Demeyer]

As a workaround, using CC instead of CDF seems to work.

[Emmanuel]

Yes indeed! Thanks

[kcrisman]

On Friday, February 8, 2013 6:00:32 AM UTC-5, Emmanuel wrote:

Yes indeed! Thanks

Please post this animation/movie when you have finished it; I always finish my undergrad NT class with a day on the RH, and having the movie would be great.  (I have some graphics but of course they take a LONG time to compute, like the explicit formula approximation with zeros of zeta up to 150 or 200, so I don't make them into interacts.)

[Jeroen Demeyer]

Fixed (the problem was the Sage <-> PARI interface) and needs review:
http://trac.sagemath.org/sage_trac/ticket/14082