On May 16, 2012, Raymond Toy wrote:
-------------------------
>> testing Raymond Toy's bug #3526359
>> fix contained in my gamma-new.lisp.
>>
>> float and bfloat of gamma_incomplete(1/5,-32)
>> good to 15 digits (gcl).
>>
>> compile-file complains about two items.
>>
>> the bfloat evaluation now contains the
>> comment "new stuff", and also should
>>
>
>Oops. I'll fix that.
>
>> have an extra expansion of the result.
>>
>
>Oh, I guess previously you didn't have to call expand to get the expected
>result. I can fix that.
-------------------------------------------------------
another problem, found while playing with some integrals,
is that bfloat is not processing (-1)^(4/5) with domain
equal to complex.
-------------------------------------------------------
(%i1) display2d:false$
(%i2) load("c:/work2/gamma-new.o");
(%o2) "c:/work2/gamma-new.o"
(%i3) domain:complex;
(%o3) complex
(%i4) ee : integrate(exp(x^5),x,1,2);
(%o4) ((-1)^(4/5)*gamma_incomplete(1/5,-32)-(-1)^(4/5)
*gamma_incomplete(1/5,-1))/5
(%i5) block([fpprec:20],expand(bfloat(ee)));
new stuff
(%o5) -5.9556722908249527029b11*(-1)^(4/5)*%i-8.1972796653426019708b11
*(-1)^(4/5)
(%i6) build_info();
(%o6) ?%build_info("5.27.0","2012-05-08 11:27:57","i686-pc-mingw32",
"GNU Common Lisp (GCL)","GCL 2.6.8")
-----------------------------------------------
Ted