On 10/29/06, William Stein <[EMAIL PROTECTED]> wrote:
>
> On Sun, 29 Oct 2006 10:36:15 -0700, David Joyner <[EMAIL PROTECTED]>
> wrote:
> > This is where it hangs. If I just try to plot the function directly
> > then it immediately
> > fails witha PariError:
> >
> > sage: P = plot(hypergeometric_U(1.0,1.0,x),0.1,0.9)
> > ---------------------------------------------------------------------------
> > <class 'gen.PariError'> Traceback (most recent call
> > last)
> >
> > /mnt/hd200/sagefiles/sage-1.4.1/<ipython console> in <module>()
> >
> > /home/wdj/sagefiles/sage-1.4.1/local/lib/python2.5/site-packages/sage/functions/special.py
> > in hypergeometric_U(alpha, beta, x, prec)
> > 551 from sage.libs.pari.all import pari
> > 552 R,a = _setup(prec)
> > --> 553 b = R(pari(alpha).hyperu(beta,x))
> > 554 pari.set_real_precision(a)
> > 555 return b
> >
> > /mnt/hd200/sagefiles/sage-1.4.1/gen.pyx in gen._pari_trap()
> >
> > <class 'gen.PariError'>: incorrect type (20)
> >
> >
> > Do you know what this means?
>
> It means you're trying to evaluate the function at a point for which PARI
> doesn't know how to do the evaluation.
It seems to me that Pari does know how to do the calculation,
since
sage: f = lambda x: gp.eval("hyperu(1,1,%s)"%x)
sage: P = plot(f,0,1)
sage: show(P)
work just fine. If I was to guess, it seems that the pari class is
barfing for some reason I can't figure out.
>
> William
>
> >
>
--~--~---------~--~----~------------~-------~--~----~
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-support
URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/
-~----------~----~----~----~------~----~------~--~---
