Switch to use gp instead, and mark it on the todo list to get to the
bottom of this.
Using gp isn't so bad.
On 10/29/06, David Joyner <[EMAIL PROTECTED]> wrote:
>
> 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
> >
> > >
> >
>
> >
>
--
William Stein
Associate Professor of Mathematics
University of Washington
--~--~---------~--~----~------------~-------~--~----~
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/
-~----------~----~----~----~------~----~------~--~---
