1. 23 Dec, 2020 1 commit
  2. 29 Nov, 2020 1 commit
    • Raymond Toy's avatar
      Fix #6: Support sbcl · e01e26ab
      Raymond Toy authored
      Update code to work with sbcl:
      - Add implementation for `float-infinity-p` and friends.
      - Typos where `-` and `/` were used instead of `sub-qd` and `div-qd`.
      - Change a `defconstant +digits+` to use `defparameter`
  3. 27 Nov, 2020 1 commit
  4. 26 Nov, 2013 4 commits
    • Raymond Toy's avatar
      Switch to using atan-qd/taylor as the default atan implementation. · 0a607f3a
      Raymond Toy authored
      o Use atan-qd/taylor and atan2-qd/taylor as the atan-qd and atan2-qd
      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!)
      o Add comment for oct.atan.4 for why it works.
    • Raymond Toy's avatar
      Reduce size of the atan partition and node tables by half. · 793166a0
      Raymond Toy authored
      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.
      o Use table length for binary search, so we're independent of other
      o Remove FIXME.
    • Raymond Toy's avatar
      Fix atan/taylor bug and add atan2/taylor. · 4d09f5a3
      Raymond Toy authored
      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
      o Add timing test for atan2-qd.
    • Raymond Toy's avatar
      Implement a new atan method that is 21 times faster. · b9a21f1f
      Raymond Toy authored
      The new method uses a table and the taylor series for atan to compute
      atan for real args.
      o Add the two tables used for the atan computation.
      o Include clisp code that generates the tables.
      o Implement atan using a set of tables and a Taylor series.
      o Add a simple timing program to evaluate how fast each atan method
  5. 25 Nov, 2013 1 commit
    • Raymond Toy's avatar
      Clean up rem-pi/2. · 3c136f7c
      Raymond Toy authored
      o Remove old version
      o Update new version to skip reduction if the arg <= pi/4.
  6. 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.
  7. 10 Feb, 2011 1 commit
  8. 16 Jul, 2008 1 commit
  9. 28 Nov, 2007 1 commit
    • Raymond Toy's avatar
      Adjust code so that CMUCL can use arrays to store quad-doubles instead · 3e39d143
      Raymond Toy authored
      of using a (complex double-double-float).
      With these changes, CMUCL uses arrays and (rt:do-tests) passes
      o Push :oct-array onto *FEATURES* to use arrays.  This is the default
        if not building on CMUCL.
      o Fix two erroneous uses of zerop on a quad-double in sinh-qd and
      o Fix two erroneous uses of + on %quad-double; they should have used
        ADD-QD instead.
      o Change conditionalization to allow arrays for CMUCL.
      o Update compiler macros appropriately.
      o Adjust optional target arg appropriately for oct-array feature.
      o Clean up IGNORE declarations.
      o Add some more declarations for the target to make CMUCL happier.
  10. 07 Nov, 2007 1 commit
  11. 18 Oct, 2007 1 commit
  12. 17 Oct, 2007 1 commit
  13. 16 Oct, 2007 3 commits
    • Raymond Toy's avatar
      qd-fun.lisp: · d58b6166
      Raymond Toy authored
      o Remove extra rem-pi/2-int.
      o Fix typo in float-infinity-p for Allegro.
      o Ignore var in float-infinity-p and friends.
    • Raymond Toy's avatar
      o Add default implementation of float-infinity-p, float-nan-p, · 4968eb46
      Raymond Toy authored
        float-trapping-nan-p.  These return NIL by default, unless the Lisp
        implementation has a suitable version.
      o Remove CMU conditionalization for float-infinity-p, float-nan-p,
    • Raymond Toy's avatar
      qd-package.lisp: · 0aaa8918
      Raymond Toy authored
      o Don't :USE CMUCL's EXTENSIONS package anymore.  Import just the
        symbols we need.
      o No need to use package qualifiers.
      o Need package qualifier for MAYBE-INLINE.
      o Need package qualifier for *INLINE-EXPANSION-LIMIT*
  14. 15 Oct, 2007 3 commits
  15. 13 Oct, 2007 1 commit
  16. 10 Oct, 2007 1 commit
  17. 12 Sep, 2007 1 commit
    • Raymond Toy's avatar
      qd-fun.lisp: · 70c699b1
      Raymond Toy authored
      o Remove unused var R1 in EXP-QD/REDUCE.
      o TAN-QD was calling ZEROP instead of ZEROP-QD.
      o Comment out extra copy of ASINH-QD.
      o Ignore unused var in QD-PRINT-EXPONENT and QD-READER.
      o Remove extra version of DIV-QD.
  18. 25 Aug, 2007 1 commit
  19. 14 Jun, 2007 1 commit
  20. 13 Jun, 2007 5 commits
  21. 12 Jun, 2007 3 commits
    • toy's avatar
      Change exit condition from < to <= since very small values of the arg · f4870c1a
      toy authored
      would set the threshold to 0 and we'd never exit the loop.  (Should we
      change the threshold to be something else?  if x is less than about
      1d-33, we don't need any more terms of the taylor series.)
    • toy's avatar
      o Fix debugging prints to use qdi::qd-format. · cd998025
      toy authored
      o Oops.  Desired iteration was chosen incorrectly if xx and yy had
        different signs.
    • toy's avatar
      qd-fun.lisp: · 8aadd62b
      toy authored
      o SQRT-QD has round-off problems for very large and very small
        numbers.  Solve this by scaling the number by a power of 4 before
        computing the square root.
      o Moved LOGB-FINITE and the core of QD-CSSQS to this file from
        qd-methods.lisp.  Rename the core routine as HYPOT-QD.
      o Use it in ATAN-QD/NEWTON to get rid of unnecessary overflow.
      o Moved LOGB-FINITE and the core of QD-CSSQS from this file to
        qd-fun.lisp.  Use HYPOT-QD in QD-CSSQS.
      o Export HYPOT-QD.
  22. 02 Jun, 2007 3 commits
  23. 01 Jun, 2007 3 commits