Commit [8f1085] Maximize Restore History

Work around bugs in COMPLEX-ERF. Bug #1: sqrt(z^2)/z incorrect

when Re(z) = 0 and Im(z) < 0. Bug #2: (EXPT Z 2.0) yields spurious imaginary part
in some Lisps. Bug #3: GAMMA-INCOMPLETE returns different values (conjugates) for
imaginary part of z^2 equal to 0.0 or -0.0, if the Lisp implementation recognizes
signed zero.
Resolution of #1 and #3 is to work in upper half plane. Resolution of #2 is to
replace (EXPT Z 2.0) with (* Z Z).
Incidentally need to increase tolerance for erf(-0.75+%i) test as (* Z Z) yields
slightly different results for some Lisps.

robert_dodier robert_dodier 2012-12-03

changed src
changed src/gamma.lisp
changed tests
changed tests/rtest_gamma.mac
src
Directory.
src/gamma.lisp Diff Switch to side-by-side view
Loading...
tests
Directory.
tests/rtest_gamma.mac Diff Switch to side-by-side view
Loading...