find_root with bfloats



 On 9/26/10 9:39 PM, Robert Dodier wrote:
> On 9/26/10, Raymond Toy <toy.raymond at gmail.com> wrote:
>
>> As the author of the bigfloat package, I'm biased.  But, as you can see,
>> I barely made any changes to find-root-subr, so it does make it easier
>> to convert existing double-float routines to handle bigfloats.  The
>> difficult parts are usually in figuring out when to convert to and from
>> the bigfloat representation to maxima's internal bfloat representation.
>>
>> But here is my quick attempt at bf_find_root.  I'm pretty sure there are
>> bugs here, but the couple of examples I tried work:
> I don't really have any preference. So if you have a working version,
> by all means please go ahead and commit it.
I hadn't thought about that, but I guess I can add that to intpol.lisp. 
However, on review of the code, there are some issues with
coerce-bfloat-fun.  I'll have to fix those.  And it would be best if we
didn't have to duplicate coerce-float-fun and friends.  I'll have to
think of some way to have them work with both float and bfloat numbers.

Ray