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
have an extra expansion of the result.
maxima output was compared to wolfram
alpha N[Gamma[1/5,-32]] output which is
-4.098639832671665..x 10^12
-2.977836145416076..x 10^12 i
=======================
(%i7) :lisp (compile-file "c:/work2/gamma-new.lisp")
Compiling c:/work2/gamma-new.lisp.
; (IN-PACKAGE :BIGFLOAT) is being compiled.
;; Warning: The package operation (IN-PACKAGE :BIGFLOAT) was in a bad place.
; (IN-PACKAGE :MAXIMA) is being compiled.
;; Warning: The package operation (IN-PACKAGE :MAXIMA) was in a bad place.
End of Pass 1.
End of Pass 2.
OPTIMIZE levels: Safety=2, Space=3, Speed=3
Finished compiling c:/work2/gamma-new.lisp.
#pc:/work2/gamma-new.o
(%i7) load("c:/work2/gamma-new.o");
(%o7) "c:/work2/gamma-new.o"
(%i8) fpprec;
(%o8) 16
(%i9) float(gamma_incomplete(1/5,-32));
(%o9) -2.9778361454160732E+12*%i-4.0986398326716611E+12
(%i10) bfloat(gamma_incomplete(1/5,-32));
new stuff
(%o10) 4.590843711998803b0-2.533098724238109b12
*(1.175570504584946b0*%i+1.618033988749895b0)
(%i11) expand(%);
(%o11) -2.977836145416077b12*%i-4.098639832671666b12
=====================
Ted Woollett