- Apr 12, 2012
-
-
Raymond Toy authored
-
Raymond Toy authored
-
Raymond Toy authored
-
- Apr 11, 2012
-
-
Raymond Toy authored
caused oct.elliptic-pi.n2.d to fail.
-
Raymond Toy authored
* qd-gamma.lisp:: * Return sqrt(least-positive-normalized-double) instead of least-positive-normalized-double. * rt-tests.lisp:: * Add test for this case.
-
- Apr 09, 2012
-
-
Raymond Toy authored
-
- Apr 08, 2012
-
-
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.
-
- Mar 24, 2012
-
-
Raymond Toy authored
-
- Mar 23, 2012
-
-
Raymond Toy authored
-
Raymond Toy authored
* Don't let NaN's fool check-accuracy * Add tests for exp-integral-e with v = 1.
-
Raymond Toy authored
tests.
-
- Mar 22, 2012
-
-
Raymond Toy authored
-
Raymond Toy authored
tests too.
-
Raymond Toy authored
-
- Dec 05, 2011
-
-
Raymond Toy authored
oct.elliptic-pi.n1.d
-
Raymond Toy authored
-
- Mar 29, 2011
-
-
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.
-
- Mar 28, 2011
-
-
Raymond Toy authored
-
- Mar 25, 2011
-
-
Raymond Toy authored
-
- Mar 24, 2011
-
-
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.
-
- Mar 17, 2011
-
-
Raymond Toy authored
qd-gamma.lisp: o Add implementation for Lentz's algorithm for evaluating continued fractions. o Implement incomplete-gamma and incomplete-gamma-tail using continued fractions. o Implement erf rt-tests.lisp: o Add tests
-
- Mar 16, 2011
-
-
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.
-
- Mar 14, 2011
-
-
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.
-
- Mar 13, 2011
-
-
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.
-
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.
-
- Mar 12, 2011
-
-
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.
-
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.
-
Raymond Toy authored
-
- Mar 11, 2011
-
-
Raymond Toy authored
-
- Mar 08, 2011
-
-
Raymond Toy authored
-
Raymond Toy authored
qd-elliptic.lisp: o Add support for the complete elliptic integral K using Carlson's Rf function. rt-tests.lisp: o Fix indentation for oct.atan.5 o Add tests for elliptic K o Add tests for Jacobi sn, cn, and dn functions.
-
- Mar 07, 2011
-
-
Raymond Toy authored
CHECK-SIGN just returns T if Lisp thinks -0d0 and 0d0 are the same number. This is only true if Lisp doesn't support signed zeroes.
-
- Feb 11, 2011
-
-
Raymond Toy authored
o Fix typo in number in for test oct.atan-branch-pos.2 o Add tests for atanh branch cuts.
-
- Feb 09, 2011
-
-
Raymond Toy authored
-
- Oct 15, 2007
-
-
Raymond Toy authored
OCTI instead of QDI. o Rename OCT package to NET.COMMON-LISP.OCT, with a nickname of OCT o Remove nickname of QD. (Conflicts with other packages dealing with quad-doubles.) o Update all uses of QDI: to OCTI: qd-fun.lisp: o Add REM-PI/2 to do a simpler computation if the arg is small enough. Otherwise, use the accurate but expensive rem operation. o Renamed ACCURATE-SIN-QD to SIN-QD, etc. o Update SIN-QD etc to use REM-PI/2.
-
Raymond Toy authored
o Oops. In INTEGER-DECODE-QD, the signs of the parts were not computed correctly when combining them into the final integer result. rt-tests.lisp: o Add a test for INTEGER-DECODE-QD. o Use OCT as the package, not QD.
-
- Oct 13, 2007
-
-
Raymond Toy authored
-
- Sep 18, 2007
-
-
Raymond Toy authored
o Add optional arg to enable/disable printing of results. Default is on. rt-tests.lisp: o Don't print results.
-
- Aug 27, 2007
-
-
Raymond Toy authored
-