Karim BELABAS on Thu, 3 Jul 2003 08:37:50 +0200 (MEST)


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: using zbrent in C mode


On Wed, 2 Jul 2003, Igor Schein wrote:
> On Wed, Jul 02, 2003 at 11:36:46PM +0200, Karim BELABAS wrote:
> > On Wed, 2 Jul 2003, Bill Allombert wrote:
> > > On Wed, Jul 02, 2003 at 08:58:55PM +0200, Olivier Ramare wrote:
> > >>   I'll be able to program in pari-C one day ...
> > >> Just now, can anyone tell me how to program that
> > >>
> > >> x=solve(t=a,b,sin(t))
> > >
> > > Short answer: You cannot.
> > >
> > > Long answer: You can try to use lisexpr():
> >
> > Intermediate answer: it was a 2 minutes job to convert the existing
> > zbrent into
> >
> > GEN
> > zbrent(void *dat, GEN (*f)(GEN, void*), GEN a, GEN b, long prec)
>
> Here're some instant regression cases:
>
> ? solve(x=-1,0,sin(x)-x);
>   ***   forbidden assignment t_SMALL --> t_REAL.
> ? solve(x=0,1,deriv(x));
>   ***   forbidden  ***   unknown type 52.

Your 'Def' database is not up to date.

You must rebuild 'Def' first [ cf Bill's TODO postings from a few days ago ]:

  cd desc && make

I forgot this side effect of the new function description system.

    Karim.
-- 
Karim Belabas                     Tel: (+33) (0)1 69 15 57 48
Dép. de Mathématiques, Bât. 425   Fax: (+33) (0)1 69 15 60 19
Université Paris-Sud              http://www.math.u-psud.fr/~belabas/
F-91405 Orsay (France)            http://www.parigp-home.de/  [PARI/GP]