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

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

Add function to compute J(1/2,z) and use it in the tests.

sum-an was summing one too many terms.

Add iterative versions for some functions.

Make big-n a defvar so we can change it easily.

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

 * 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.

Add more tests for Bessel J.  Not all of them pass.

Add more tests.

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

Enable test of complex args.

Fix for complex args.

Update accuracy for bessel-j-1/2.q test.

Add tests for bessel-j.

 * 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.

First cut at Bessel Y.  Not working yet.

Add TODO list for bessel-j.

Correct the computation of alpha and beta.

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

Correct some comments, remove unused code.

Update accuracy for oct.elliptic-pi.n2.d and add new test that
caused oct.elliptic-pi.n2.d to fail.

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.

Fix bug in s-bessel-j.  Microoptimize integer-bessel-j-exp-arc for the
case where v is an integer.

 * Implement sum-big-ia.
 * Add series for Bessel J.  (Not working yet.)

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

Fix two typos in the names of tests.

Fix some mistakes in exp-integral-e when changing algorithm to use
series.

Merge branch 'master' of ssh://common-lisp.net/var/git/projects/oct/oct

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.

Oops. The second arg to FLOAT is optional!

Fix typo in %big-a.  Just use incomplete-gamma-tail in big-i.

T is not a variable.  Use a different name.

Build qd-bessel now.  (qd-bessel is a work in progress.)

Define macro WITH-FLOATING-POINT-CONTAGION.

qd-methods.lisp:
o Define the macro

qd-gamma.lisp:
o Use it.

Fix some comments.

Fix a divide by zero error in s-exp-integral-e for v = 1.  We need to
skip the first term in the series.

Add more parts of the exp-arc algorithm.  Needs lots of work, but it
seems that bessel_j(n,z) mostly works.

Add some comments, rename bessel-j-exp-arc to integer-bessel-j-exp-arc.

First cut at Bessel functions.  Needs lots of work.

qd-gamma depends on qd-complex.

Need to take ABS before testing for 1/2.

Oops.  Need to take the realpart before call FLOAT-NAN-P.

Need to apply contagion for exp-integral-e.

Use exp-integral-e to evaluate incomplete-gamma-tail for real,
negative values of the parameter.

Document CHECK-ACCURACY

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.

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.

Oops.  In FLOAT, if it's already a float, don't change it if the
second arg is not given.

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.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

Clean up implemenation of s-exp-integral-e a bit.

Correct the computation of s-exp-integral-e when v is an integer.  Add
tests.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

Conflicts:
qd-gamma.lisp

Handle the case of negative v in exp-integral-e(v, z).  Add tests.

Fix bug in psi for -n/2 for n odd which was causing an overflow.  Add
tests too.

Add a few tests for exp-integral-e.  More needed.

Implement psi and fix exp-integral-e for integral values of v.  Needs
some more work.

Allow exp-integral-e to work with negative integer orders.  First
cut....

First cut at psi function.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

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.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

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.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

Oops.  Put back QEXPT method for QD-COMPLEX, QD-REAL.

Merge branch 'master' of git://common-lisp.net/projects/oct/oct

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.

Get rid of compiler warning about unused variable in theta-1 and
theta-3.

Reduce required accuracy so CCL can pass oct.theta2.1.d and
oct.elliptic-pi.n1.d

Update version (and make it compatible with latest asdf2).

Update asdf version (and make it compatible with latest asdf).

Remove some random junk that we left in.

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.

Signal error for gamma of negative integers.

Fix bug in CEILING and FCEILING; add tests.

Revert previous blank line insertions.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Fix typo in #q() number in fresnel-s.2q test.

Update required accuracy for elliptic-pi.n0.q.

Add DOMAIN-ERROR condition; fix bug in FRESNEL-S-SERIES.

qd-methods.lisp:
o Define DOMAIN-ERROR condition to allow signaling errors for
  incorrect domains.

qd-gamma.lisp:
o Signal domain error in CF-INCOMPLETE-GAMMA-TAIL if necessary.
o Fix bug in FRESNEL-S-SERIES.  We were comparing a real against a
  complex.

Fix incomplete-gamma-tail for negative reals, use series for Fresnel S
for small arg and update tests.

qd-gamma.lisp:
o INCOMPLETE-GAMMA-TAIL was hanging for arguments on the negative real
  axis.  Use INCOMPLETE-GAMMA in this case too.
o Add the series expansion for Fresnel S and use it for evaluating it
  for small arguments.  We were losing accuracy with the existing
  algorithm.

rt-tests.lisp:
o Update thresholds for elliptic-pi-n0.d, elliptic-pi.n2.q,
  theta3.1.d.
o Fix typo in test name.  gamma-incomplete-tail.1.q should have been
  2.q.
o Add tests for gamma-incomplete-tail for arguments on the negative
  real axis.
o Add tests for Fresnel S.

qd-gamma depends on qd-reader.

Move the uses of foo-t before the use of foo-t.

Add rem and mod functions that were left out.

Add series for incomplete-gamma for when the fraction is slow.

o Add series for incomplete gamma function for small a and z.  Needed
  because the continued fraction is slow in this range.
o In INCOMPLETE-GAMMA-TAIL, call INCOMPLETE-GAMMA instead of
  CF-INCOMPLETE-GAMMA just in case a and z are small.
o In INCOMPLETE-GAMMA, use the series for small a and z.
o Simplify evaluation of Si(z) when z is real.

Implement Si and Ci, sin and cos integrals.

Add Fresnel integrals; fix issue in incomplete-gamma for large args.

o INCOMPLETE-GAMMA was returning bad values for large (complex)
  arguments.  Fix this by using incomplete gamma tail function since
  the incomplete gamma function approaches gamma for large arguments.

o Implement Fresnel S and C functions.

Implement exponential integral E.