1. 26 Nov, 2013 1 commit
    • Raymond Toy's avatar
      Switch to using atan-qd/taylor as the default atan implementation. · 0a607f3a
      Raymond Toy authored
      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.
      0a607f3a
  2. 24 Nov, 2013 1 commit
    • Raymond Toy's avatar
      Fix ticket:5 · 6127a79e
      Raymond Toy authored
      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.
      6127a79e
  3. 26 Apr, 2013 1 commit
  4. 17 Apr, 2012 1 commit
  5. 13 Apr, 2012 2 commits
  6. 12 Apr, 2012 3 commits
  7. 11 Apr, 2012 2 commits
  8. 09 Apr, 2012 1 commit
  9. 08 Apr, 2012 1 commit
    • Raymond Toy's avatar
      Define FLOATP, fix bugs in FLOAT. · 8ec0d200
      Raymond Toy authored
      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.
      8ec0d200
  10. 24 Mar, 2012 1 commit
  11. 23 Mar, 2012 3 commits
  12. 22 Mar, 2012 3 commits
  13. 05 Dec, 2011 2 commits
  14. 29 Mar, 2011 1 commit
    • Raymond Toy's avatar
      Make gamma accurate for integers; add precision test for · af9fc6fc
      Raymond Toy authored
      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.
      af9fc6fc
  15. 28 Mar, 2011 1 commit
  16. 25 Mar, 2011 1 commit
  17. 24 Mar, 2011 1 commit
    • Raymond Toy's avatar
      Fix incomplete-gamma-tail for negative reals, use series for Fresnel S · bed5763d
      Raymond Toy authored
      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.
      bed5763d
  18. 17 Mar, 2011 1 commit
  19. 16 Mar, 2011 1 commit
    • Raymond Toy's avatar
      Add gamma and log-gamma functions; work in progress. · a40062c6
      Raymond Toy authored
      oct.asd:
      o Add qd-gamma.lisp.  The implementations need some work.  The
        accuracy is less than desired because gamma(2.0) /= 1.  It's close
        but not quite right.
      
      rt-tests.lisp:
      o Basic tests of the gamma function.  Accuracy is not as good as we
        would ike.
      
      qd-gamma.lisp:
      o New file for implementation of gamma function.
      a40062c6
  20. 14 Mar, 2011 1 commit
    • Raymond Toy's avatar
      Fix typo; update required accuracy. · d6439bd9
      Raymond Toy authored
      o The names of the elliptic functions changed and we forgot to update
        the tests to use the new names.
      o Reduce required accuracy of some tests.
      d6439bd9
  21. 13 Mar, 2011 2 commits
    • Raymond Toy's avatar
      Add Elliptic theta functions and tests. · 8ade177a
      Raymond Toy authored
      oct.asd:
      o Add qd-theta.
      
      qd-theta.lisp:
      o New file with Elliptic theta functions and elliptic nome function.
      
      rt-tests.lisp:
      o Tests for theta functions.
      o Relax accuracy requirements for some of the tests os that they can
        pass.
      8ade177a
    • Raymond Toy's avatar
      Fix typo; update accuracy requirements. · f4a60f8c
      Raymond Toy authored
      qd-elliptic.lisp:
      o Fix missing paren.
      
      rt-tests.lisp:
      o Reduce accuracy requirements so the tests pass for the random
        arguments.
      f4a60f8c
  22. 12 Mar, 2011 3 commits
    • Raymond Toy's avatar
      Fix possible bug in elliptic-pi; add comments. · 5bd5df93
      Raymond Toy authored
      qd-elliptic.lisp:
      o Add some comments
      o Fix a possible bug if n is a complex number or a negative number.
      
      rt-tests.lisp:
      o Remove one broken test.
      o Fix the other tests for elliptic-pi and adjust required precision down
        a bit so the tests can pass.
      5bd5df93
    • Raymond Toy's avatar
      Add prelimary support for integrals of the 3rd kind. · 390a7483
      Raymond Toy authored
      qd-elliptic.lisp:
      o Clean up for unused variable in ELLIPTIC-K
      o Add Carlson's Rj functions
      o Implement elliptic-pi using Carlson's method.
      
      rt-tests.lisp:
      o Add many tests for elliptic-pi.  Some tests pass, and some fail.  The
        failing tests are not enabled because I don't know if the failure is
        because the test itself is wrong or if the integral is wrong.
      390a7483
    • Raymond Toy's avatar
      Add tests for contagion support. · 147fa2c7
      Raymond Toy authored
      147fa2c7
  23. 11 Mar, 2011 1 commit
  24. 08 Mar, 2011 2 commits
  25. 07 Mar, 2011 1 commit
  26. 11 Feb, 2011 1 commit
  27. 09 Feb, 2011 1 commit