projects/oct/oct.git
4 months agoAdd implementation for tan(x) using duplication. master
Raymond Toy [Wed, 27 Nov 2013 17:09:50 +0000] 
Add implementation for tan(x) using duplication.

Includes Lentz' algorithm for evaluating continued fractions and
updates the timing test to include tan(x) using duplication.  Tests
show this is much slower than existing algorithms.

4 months agoUpdate the versions.
Raymond Toy [Tue, 26 Nov 2013 19:54:59 +0000] 
Update the versions.

4 months agoSwitch to using atan-qd/taylor as the default atan implementation.
Raymond Toy [Tue, 26 Nov 2013 17:55:46 +0000] 
Switch to using atan-qd/taylor as the default atan implementation.

qd-fun.lisp:
o Use atan-qd/taylor and atan2-qd/taylor as the atan-qd and atan2-qd
  implementation.
o Fix bug in atan-qd/taylor. We need a special case for x small where
  we don't want to use the table. (Found by our rt-test.  Yay!)

rt-tests.lisp:
o Add comment for oct.atan.4 for why it works.

4 months agoReduce size of the atan partition and node tables by half.
Raymond Toy [Tue, 26 Nov 2013 17:36:02 +0000] 
Reduce size of the atan partition and node tables by half.

qd-const.lisp:
o Update tables and generator to use tables with only half as many
  entries as before. This is because atan-qd/taylor never needs to
  compute directly for x > 1.

qd-fun.lisp:
o Use table length for binary search, so we're independent of other
  constants.
o Remove FIXME.

4 months agoFix atan/taylor bug and add atan2/taylor.
Raymond Toy [Tue, 26 Nov 2013 17:11:45 +0000] 
Fix atan/taylor bug and add atan2/taylor.

qd-fun.lisp:
o Fix issue with atan-qd/taylor for large args. We use the
  relationship atan(x) = pi/2-atan(1/x) for x > 1.
o Add atan2-qd/taylor.  This appears to be much faster than
  atan2-qd/newton.

timing.lisp:
o Add timing test for atan2-qd.

4 months agoImplement a new atan method that is 21 times faster.
Raymond Toy [Tue, 26 Nov 2013 07:10:26 +0000] 
Implement a new atan method that is 21 times faster.

The new method uses a table and the taylor series for atan to compute
atan for real args.

qd-const.lisp:
o Add the two tables used for the atan computation.
o Include clisp code that generates the tables.

qd-fun.lisp:
o Implement atan using a set of tables and a Taylor series.

timing.lisp:
o Add a simple timing program to evaluate how fast each atan method
  is.

4 months agoClean up rem-pi/2.
Raymond Toy [Mon, 25 Nov 2013 04:53:22 +0000] 
Clean up rem-pi/2.

o Remove old version
o Update new version to skip reduction if the arg <= pi/4.

4 months agoUpdate the version.
Raymond Toy [Mon, 25 Nov 2013 04:11:29 +0000] 
Update the version.

4 months agoFix ticket:5
Raymond Toy [Sun, 24 Nov 2013 18:41:27 +0000] 
Fix ticket:5

o Implement better pi reduction and also always do pi reduction.
o Add test for sin(pi) and cos(2^120), to test pi reduction.
o Accuracy for erfc is only 198 bits.

11 months agoAdd continued fractions for erf, erfc, w and Dawson's integral.
Raymond Toy [Fri, 26 Apr 2013 03:11:27 +0000] 
Add continued fractions for erf, erfc, w and Dawson's integral.

 qd-gamma.lisp::
 Add continued fractions

 rt-tests.lisp::
 Add a test for erfc

23 months agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Tue, 17 Apr 2012 23:16:07 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

23 months agoAdd function to compute J(1/2,z) and use it in the tests.
Raymond Toy [Tue, 17 Apr 2012 23:15:37 +0000] 
Add function to compute J(1/2,z) and use it in the tests.

2 years agosum-an was summing one too many terms.
Raymond Toy [Tue, 17 Apr 2012 16:32:27 +0000] 
sum-an was summing one too many terms.

2 years agoAdd iterative versions for some functions.
Raymond Toy [Tue, 17 Apr 2012 03:45:12 +0000] 
Add iterative versions for some functions.

2 years agoMake big-n a defvar so we can change it easily.
Raymond Toy [Tue, 17 Apr 2012 02:30:12 +0000] 
Make big-n a defvar so we can change it easily.

2 years agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Mon, 16 Apr 2012 17:05:21 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

2 years ago * Add comments for {{{integer-bessel-j-exp-arc}}}.
Raymond Toy [Sun, 15 Apr 2012 17:50:16 +0000] 
 * Add comments for {{{integer-bessel-j-exp-arc}}}.
 * Simplify {{{sum-an}}} so we stop the sum when the terms no longer
   contribute to the sum.
 * Change {{{big-n}}}.  This still needs work.

2 years agoAdd more tests for Bessel J. Not all of them pass.
Raymond Toy [Fri, 13 Apr 2012 06:24:03 +0000] 
Add more tests for Bessel J.  Not all of them pass.

2 years agoAdd more tests.
Raymond Toy [Fri, 13 Apr 2012 01:27:08 +0000] 
Add more tests.

2 years agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Thu, 12 Apr 2012 22:32:00 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

2 years agoEnable test of complex args.
Raymond Toy [Thu, 12 Apr 2012 21:37:38 +0000] 
Enable test of complex args.

2 years agoFix for complex args.
Raymond Toy [Thu, 12 Apr 2012 21:37:27 +0000] 
Fix for complex args.

2 years agoUpdate accuracy for bessel-j-1/2.q test.
Raymond Toy [Thu, 12 Apr 2012 16:46:59 +0000] 
Update accuracy for bessel-j-1/2.q test.

2 years agoAdd tests for bessel-j.
Raymond Toy [Thu, 12 Apr 2012 16:30:28 +0000] 
Add tests for bessel-j.

2 years ago * Use the g[k] formula instead of r[2*k+1] because we fail to
Raymond Toy [Thu, 12 Apr 2012 16:30:10 +0000] 
 * Use the g[k] formula instead of r[2*k+1] because we fail to
   converge when v = 0.
 * Clear hash tables in bessel-j.

2 years agoFirst cut at Bessel Y. Not working yet.
Raymond Toy [Thu, 12 Apr 2012 03:32:23 +0000] 
First cut at Bessel Y.  Not working yet.

2 years agoAdd TODO list for bessel-j.
Raymond Toy [Thu, 12 Apr 2012 03:15:47 +0000] 
Add TODO list for bessel-j.

2 years agoCorrect the computation of alpha and beta.
Raymond Toy [Thu, 12 Apr 2012 03:06:38 +0000] 
Correct the computation of alpha and beta.

2 years agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Wed, 11 Apr 2012 20:53:07 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

2 years agoCorrect some comments, remove unused code.
Raymond Toy [Wed, 11 Apr 2012 16:18:31 +0000] 
Correct some comments, remove unused code.

2 years agoUpdate accuracy for oct.elliptic-pi.n2.d and add new test that
Raymond Toy [Wed, 11 Apr 2012 15:47:28 +0000] 
Update accuracy for oct.elliptic-pi.n2.d and add new test that
caused oct.elliptic-pi.n2.d to fail.

2 years agoVALUE-OR-TINY returned a value that was too tiny.
Raymond Toy [Wed, 11 Apr 2012 03:36:11 +0000] 
VALUE-OR-TINY returned a value that was too tiny.

 * qd-gamma.lisp::
   * Return sqrt(least-positive-normalized-double) instead of
     least-positive-normalized-double.
 * rt-tests.lisp::
   * Add test for this case.

2 years agoFix bug in s-bessel-j. Microoptimize integer-bessel-j-exp-arc for the
Raymond Toy [Tue, 10 Apr 2012 16:51:55 +0000] 
Fix bug in s-bessel-j.  Microoptimize integer-bessel-j-exp-arc for the
case where v is an integer.

2 years ago * Implement sum-big-ia.
Raymond Toy [Tue, 10 Apr 2012 07:45:13 +0000] 
 * Implement sum-big-ia.
 * Add series for Bessel J.  (Not working yet.)

2 years agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Mon, 9 Apr 2012 18:14:16 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

2 years agoFix two typos in the names of tests.
Raymond Toy [Mon, 9 Apr 2012 15:55:27 +0000] 
Fix two typos in the names of tests.

2 years agoFix some mistakes in exp-integral-e when changing algorithm to use
Raymond Toy [Mon, 9 Apr 2012 15:52:08 +0000] 
Fix some mistakes in exp-integral-e when changing algorithm to use
series.

2 years agoMerge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct
Raymond Toy [Mon, 9 Apr 2012 15:37:51 +0000] 
Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

2 years agoDefine FLOATP, fix bugs in FLOAT.
Raymond Toy [Sun, 8 Apr 2012 17:14:37 +0000] 
Define FLOATP, fix bugs in FLOAT.

qd-methods.lisp:
* Define FLOATP
* Fix bugs in FLOAT:
  * (FLOAT float nil) is an error
  * (FLOAT float) returns the float
  * (FLOAT rational) returns a single-float.

qd-package.lisp:
o Export FLOATP, shadowing CL:FLOAT.

rt-tests.lisp:
o Add a few tests for FLOAT.

2 years agoOops. The second arg to FLOAT is optional!
Raymond Toy [Sun, 8 Apr 2012 16:59:39 +0000] 
Oops. The second arg to FLOAT is optional!

2 years agoFix typo in %big-a. Just use incomplete-gamma-tail in big-i.
Raymond Toy [Sun, 8 Apr 2012 16:57:12 +0000] 
Fix typo in %big-a.  Just use incomplete-gamma-tail in big-i.

2 years agoT is not a variable. Use a different name.
Raymond Toy [Sun, 8 Apr 2012 16:37:22 +0000] 
T is not a variable.  Use a different name.

2 years agoBuild qd-bessel now. (qd-bessel is a work in progress.)
Raymond Toy [Sun, 8 Apr 2012 16:36:46 +0000] 
Build qd-bessel now.  (qd-bessel is a work in progress.)

2 years agoDefine macro WITH-FLOATING-POINT-CONTAGION.
Raymond Toy [Sun, 8 Apr 2012 16:35:35 +0000] 
Define macro WITH-FLOATING-POINT-CONTAGION.

qd-methods.lisp:
o Define the macro

qd-gamma.lisp:
o Use it.

2 years agoFix some comments.
Raymond Toy [Sun, 8 Apr 2012 16:04:18 +0000] 
Fix some comments.

2 years agoFix a divide by zero error in s-exp-integral-e for v = 1. We need to
Raymond Toy [Sun, 8 Apr 2012 06:28:08 +0000] 
Fix a divide by zero error in s-exp-integral-e for v = 1.  We need to
skip the first term in the series.

2 years agoAdd more parts of the exp-arc algorithm. Needs lots of work, but it
Raymond Toy [Sat, 7 Apr 2012 16:28:16 +0000] 
Add more parts of the exp-arc algorithm.  Needs lots of work, but it
seems that bessel_j(n,z) mostly works.

2 years agoAdd some comments, rename bessel-j-exp-arc to integer-bessel-j-exp-arc.
Raymond Toy [Sat, 7 Apr 2012 02:28:23 +0000] 
Add some comments, rename bessel-j-exp-arc to integer-bessel-j-exp-arc.

2 years agoFirst cut at Bessel functions. Needs lots of work.
Raymond Toy [Sat, 7 Apr 2012 02:21:46 +0000] 
First cut at Bessel functions.  Needs lots of work.

2 years agoqd-gamma depends on qd-complex.
Raymond Toy [Sat, 24 Mar 2012 03:48:40 +0000] 
qd-gamma depends on qd-complex.

2 years agoNeed to take ABS before testing for 1/2.
Raymond Toy [Sat, 24 Mar 2012 02:10:16 +0000] 
Need to take ABS before testing for 1/2.

2 years agoOops. Need to take the realpart before call FLOAT-NAN-P.
Raymond Toy [Sat, 24 Mar 2012 02:05:24 +0000] 
Oops.  Need to take the realpart before call FLOAT-NAN-P.

2 years agoNeed to apply contagion for exp-integral-e.
Raymond Toy [Fri, 23 Mar 2012 20:00:24 +0000] 
Need to apply contagion for exp-integral-e.

2 years agoUse exp-integral-e to evaluate incomplete-gamma-tail for real,
Raymond Toy [Fri, 23 Mar 2012 19:48:54 +0000] 
Use exp-integral-e to evaluate incomplete-gamma-tail for real,
negative values of the parameter.

2 years agoDocument CHECK-ACCURACY
Raymond Toy [Fri, 23 Mar 2012 17:38:31 +0000] 
Document CHECK-ACCURACY

2 years agoBug fix for check-accuracy, and more tests for exp-integral-e.
Raymond Toy [Fri, 23 Mar 2012 17:21:40 +0000] 
Bug fix for check-accuracy, and more tests for exp-integral-e.

* Don't let NaN's fool check-accuracy
* Add tests for exp-integral-e with v = 1.

2 years agoCleanups for psi and bug fix for exp-integral-e.
Raymond Toy [Fri, 23 Mar 2012 17:20:27 +0000] 
Cleanups for psi and bug fix for exp-integral-e.

* Allow (exp-integral-e 1 z) to work.
* psi
  * Handle psi(1) specially.
  * Do a better job with cot(%pi*z) when z is an odd multiple of 1/2
    where cot is 0.
  * Fib bug in computing the number of terms when we try to float a
    complex.  Just float the realpart.

2 years agoOops. In FLOAT, if it's already a float, don't change it if the
Raymond Toy [Fri, 23 Mar 2012 17:06:28 +0000] 
Oops.  In FLOAT, if it's already a float, don't change it if the
second arg is not given.

2 years agoFix bug in FLOAT: second arg is optional! Add FLOAT-NAN-P method.
Raymond Toy [Fri, 23 Mar 2012 16:47:16 +0000] 
Fix bug in FLOAT:  second arg is optional!  Add FLOAT-NAN-P method.

qd-methods.lisp:
* Second arg to {{{FLOAT}}} is optional.
* Add {{{FLOAT-NAN-P}}}.

qd-package.lisp:
* Need to shadow {{{EXT:FLOAT-NAN-P}}} on cmucl.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Fri, 23 Mar 2012 15:38:53 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

2 years agoClean up implemenation of s-exp-integral-e a bit.
Raymond Toy [Fri, 23 Mar 2012 05:58:43 +0000] 
Clean up implemenation of s-exp-integral-e a bit.

2 years agoCorrect the computation of s-exp-integral-e when v is an integer. Add
Raymond Toy [Fri, 23 Mar 2012 04:13:07 +0000] 
Correct the computation of s-exp-integral-e when v is an integer.  Add
tests.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Thu, 22 Mar 2012 16:02:16 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

Conflicts:
qd-gamma.lisp

2 years agoHandle the case of negative v in exp-integral-e(v, z). Add tests.
Raymond Toy [Thu, 22 Mar 2012 05:24:48 +0000] 
Handle the case of negative v in exp-integral-e(v, z).  Add tests.

2 years agoFix bug in psi for -n/2 for n odd which was causing an overflow. Add
Raymond Toy [Thu, 22 Mar 2012 04:53:37 +0000] 
Fix bug in psi for -n/2 for n odd which was causing an overflow.  Add
tests too.

2 years agoAdd a few tests for exp-integral-e. More needed.
Raymond Toy [Thu, 22 Mar 2012 01:58:02 +0000] 
Add a few tests for exp-integral-e.  More needed.

2 years agoImplement psi and fix exp-integral-e for integral values of v. Needs
Raymond Toy [Thu, 22 Mar 2012 01:45:44 +0000] 
Implement psi and fix exp-integral-e for integral values of v.  Needs
some more work.

2 years agoAllow exp-integral-e to work with negative integer orders. First
Raymond Toy [Wed, 21 Mar 2012 23:20:04 +0000] 
Allow exp-integral-e to work with negative integer orders.  First
cut....

2 years agoFirst cut at psi function.
Raymond Toy [Wed, 21 Mar 2012 22:29:07 +0000] 
First cut at psi function.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Wed, 21 Mar 2012 16:37:15 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

2 years agoAdd Euler's constant and fix dependencies.
Raymond Toy [Wed, 21 Mar 2012 16:35:22 +0000] 
Add Euler's constant and fix dependencies.

oct.asd:
o qd-const2 depends on qd-const.

qd-const.lisp:
o Add qd value for Euler's constant

qd-const2.lisp:
o Add qd-real value for Euler's constant.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Tue, 6 Dec 2011 16:50:45 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

2 years agoBetter exp-integral-e computation and fix for incomplete-gamma-tail.
Raymond Toy [Tue, 6 Dec 2011 03:58:44 +0000] 
Better exp-integral-e computation and fix for incomplete-gamma-tail.

For exp-integral-e, use the series for small z and the
incomplete-gamma-tail for near the negative real axis.  Otherwise, use
the continued fraction.

In incomplete-gamma-tail, we were using the continued fraction instead
of the incomplete-gamma function for the region just below the negative
real axis.  We should use the cf except in that region.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Mon, 5 Dec 2011 22:07:49 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

2 years agoOops. Put back QEXPT method for QD-COMPLEX, QD-REAL.
Raymond Toy [Mon, 5 Dec 2011 22:06:39 +0000] 
Oops.  Put back QEXPT method for QD-COMPLEX, QD-REAL.

2 years agoMerge branch 'master' of git://common-lisp.net/projects/oct/oct
Raymond Toy [Mon, 5 Dec 2011 21:44:51 +0000] 
Merge branch 'master' of git://common-lisp.net/projects/oct/oct

2 years agoClean up code:
Raymond Toy [Mon, 5 Dec 2011 20:48:30 +0000] 
Clean up code:

o Get rid of unused var in FLOAT-CONTAGION.
o QEXPT for qd-complex first arg should be defined in
  qd-complex.lisp.

2 years agoGet rid of compiler warning about unused variable in theta-1 and
Raymond Toy [Mon, 5 Dec 2011 20:47:02 +0000] 
Get rid of compiler warning about unused variable in theta-1 and
theta-3.

2 years agoReduce required accuracy so CCL can pass oct.theta2.1.d and
Raymond Toy [Mon, 5 Dec 2011 19:29:30 +0000] 
Reduce required accuracy so CCL can pass oct.theta2.1.d and
oct.elliptic-pi.n1.d

2 years agoUpdate version (and make it compatible with latest asdf2).
Raymond Toy [Mon, 5 Dec 2011 17:45:54 +0000] 
Update version (and make it compatible with latest asdf2).

2 years agoUpdate asdf version (and make it compatible with latest asdf).
Raymond Toy [Mon, 5 Dec 2011 17:43:15 +0000] 
Update asdf version (and make it compatible with latest asdf).

2 years agoRemove some random junk that we left in.
Raymond Toy [Mon, 5 Dec 2011 07:24:06 +0000] 
Remove some random junk that we left in.

2 years agoUse exp-integral-e for (incomplete-gamma-tail 0 z).
Raymond Toy [Mon, 5 Dec 2011 07:22:01 +0000] 
Use exp-integral-e for (incomplete-gamma-tail 0 z).

We already have exp-integral-e function so move expintegral-e
implementation to exp-integral-e.

2 years agoSignal error for gamma of negative integers.
Raymond Toy [Mon, 5 Dec 2011 05:29:17 +0000] 
Signal error for gamma of negative integers.

2 years agoFix bug in CEILING and FCEILING; add tests.
Raymond Toy [Mon, 5 Dec 2011 05:28:27 +0000] 
Fix bug in CEILING and FCEILING; add tests.

2 years agoRevert previous blank line insertions.
Raymond Toy [Fri, 13 May 2011 14:07:54 +0000] 
Revert previous blank line insertions.

The blanks were for testing integration of git and trac to close trac
tickets.  It's working now.

2 years agoFix ticket:1.
Raymond Toy [Fri, 13 May 2011 14:05:29 +0000] 
Fix ticket:1.

Attempt 9.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 13:59:46 +0000] 
Fix ticket #1.

Attempt 8.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 13:57:37 +0000] 
Fix ticket #1.

Attempt 7.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 11:39:22 +0000] 
Fix ticket #1.

Attempt 6.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 02:46:57 +0000] 
Fix ticket #1.

Attempt 5.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 02:39:11 +0000] 
Fix ticket #1.

Attempt 4.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 02:33:34 +0000] 
Fix ticket #1.

Attempt 3.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 01:45:55 +0000] 
Fix ticket #1.

Attempt 2.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

2 years agoFix ticket #1.
Raymond Toy [Fri, 13 May 2011 01:38:53 +0000] 
Fix ticket #1.

This is a white space change to test integration of git and trac to
see if the ticket is automatically closed or updated on commit.

3 years agoMake gamma accurate for integers; add precision test for
Raymond Toy [Tue, 29 Mar 2011 13:58:45 +0000] 
Make gamma accurate for integers; add precision test for
incomplete-gamma-tail.

qd-gamma.lisp:
o For integer values, just compute the gamma value directly by
  multiplication.  This works around the problem that the current
  algorithm is not as accurate as we would like.

rt-test.lisp:
o Reduce required accuracy in gamma-incomplete-tail.3.d.
o Add precision test for gamm incomplete tail near the negative real
  axis.

3 years agoMake INCOMPLETE-GAMMA-TAIL more accurate.
Raymond Toy [Tue, 29 Mar 2011 02:54:32 +0000] 
Make INCOMPLETE-GAMMA-TAIL more accurate.

o Use the new continued fraction for the incomplete-gamma when the
  argument z is close enough to the negative real axis.

3 years agoMore accurate incomplete-gamma function, add debugging to lentz, and
Raymond Toy [Tue, 29 Mar 2011 02:37:32 +0000] 
More accurate incomplete-gamma function, add debugging to lentz, and
some random clean ups.

o Add *DEBUG-CF-EVAL* to enable debugging prints in LENTZ.
o Modify LENTZ to terminate with an error if *MAX-CF-ITERATIONS* is
  reached.
o Modify LENTZ to return the function value, the number of iterations,
  and the number of times a zero value had to be replaced.
o Adjust cf-incomplete-gamma and cf-incomplete-gamma-tail not to
  signal overflow prematurely when calculating z^a*exp(-z).
o Fix doc bug in reference for continued fraction for (original)
  cf-incomplete-gamma.
o Add new version of cf-incomplete-gamma using a different continued
  fraction.  This appears to converge faster and to be more accurate
  than the original, especially for points near the negative real
  axis.

3 years agoFix typo in #q() number in fresnel-s.2q test.
Raymond Toy [Mon, 28 Mar 2011 13:08:37 +0000] 
Fix typo in #q() number in fresnel-s.2q test.

3 years agoUpdate required accuracy for elliptic-pi.n0.q.
Raymond Toy [Fri, 25 Mar 2011 03:01:51 +0000] 
Update required accuracy for elliptic-pi.n0.q.