Commit a5476401 authored by Liam Healy's avatar Liam Healy

Added tests from GSL for more random distributions

Added chi-squared (which was not loaded at all), exponential,
exponential-power, fdist, gamma-randist, tdist.  Introduced new macros
assert-posinf and assert-neginf to check for expected infinite
results.  Improved the automatic test converter and noted its known
failings.
parent e9017ff3
;; Definition of GSLL system
;; Liam Healy
;; Time-stamp: <2010-05-22 10:24:43EDT gsll-tests.asd>
;; Time-stamp: <2010-05-22 18:08:01EDT gsll-tests.asd>
;;
;; Copyright 2007, 2008, 2009 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
......@@ -51,6 +51,7 @@
(:file "cauchy")
(:file "cdot")
(:file "chebyshev")
(:file "chi-squared")
(:file "cholesky")
(:file "clausen")
(:file "column")
......
;; Mathematical functions
;; Liam Healy, Wed Mar 8 2006 - 22:09
;; Time-stamp: <2010-04-18 00:22:18EDT mathematical.lisp>
;; Time-stamp: <2010-05-22 16:48:47EDT mathematical.lisp>
;;
;; Copyright 2006, 2007, 2008, 2009 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
......@@ -49,6 +49,9 @@
sb-ext:double-float-positive-infinity)))
(export '(+nan+ +positive-infinity+ +negative-infinity+))
(map-name '+positive-infinity+ "gsl_posinf")
(map-name '+negative-infinity+ "gsl_neginf")
(map-name '+nan+ "gsl_nan")
(defconstant +nan+
(ignore-errors
......
;; Additional methods for lisp-unit
;; Liam Healy 2009-04-15 23:23:30EDT augment.lisp
;; Time-stamp: <2010-05-22 10:54:02EDT augment.lisp>
;; Time-stamp: <2010-05-22 17:52:17EDT augment.lisp>
;;
;; Copyright 2009 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
......@@ -42,3 +42,13 @@
(lisp-unit::assert-numerical-equal
,expected-value
,form)))
(defmacro assert-posinf (form)
`(lisp-unit::assert-true
(let ((val ,form))
(and (infinityp val) (plusp val)))))
(defmacro assert-neginf (form)
`(lisp-unit::assert-true
(let ((val ,form))
(and (infinityp val) (minusp val)))))
;; Convert the GSL tests
;; Liam Healy 2010-05-22 13:03:53EDT convert.lisp
;; Time-stamp: <2010-05-22 16:43:04EDT convert.lisp>
;; Time-stamp: <2010-05-22 18:26:00EDT convert.lisp>
;;; This file is not normally loaded; it is only used to convert the
;;; GSL tests in C to CL tests. It requires cl-ppcre, lisp-util, and iterate.
;;; (convert-tests-in-file "/home/liam/mathematics/gsl/cdf/test.c")
;;; Things it won't do yet:
;;; Will not convert fixnum to double.
;;; Will not recognize infinities and substitute the appropriate assert-posinf, assert-neginf test.
;;; Will not parse/convert C math specified as an argument, e.g. 100 * TEST_TOL6.
(in-package :gsl)
;;; (convert-gsl-test "(gsl_cdf_tdist_P, (0.001, 1.0), 5.00318309780080559e-1, TEST_TOL6)")
......@@ -17,7 +22,8 @@
(cons 'assert-to-tolerance
(let ((ppcre-convert
(read-from-string
(remove
(substitute
#\space
#\, ; Get rid of commas
(cl-ppcre:regex-replace ; Replace function name
"\\((\\w*)\\W*\\((.*\\))"
......
......@@ -34,16 +34,49 @@
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 1.8000000000000016d0)
(MULTIPLE-VALUE-LIST (BETA-PDF 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.19000000000000017d0)
(MULTIPLE-VALUE-LIST (BETA-P 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.8099999999999998d0)
(MULTIPLE-VALUE-LIST (BETA-Q 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.09999999999999992d0)
(MULTIPLE-VALUE-LIST (BETA-PINV 0.19d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.09999999999999988d0)
(MULTIPLE-VALUE-LIST (BETA-QINV 0.81d0 1.0d0 2.0d0))))
(ASSERT-TO-TOLERANCE (BETA-P 0.0d0 1.2d0 1.3d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (BETA-P 1.d-100 1.2d0 1.3d0) 1.344349446564896d-120 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.001d0 1.2d0 1.3d0) 3.376300425045358d-4 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.01d0 1.2d0 1.3d0) 0.005343172640389295d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.1d0 1.2d0 1.3d0) 0.08339978283067484d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.325d0 1.2d0 1.3d0) 0.3286986541805839d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.5d0 1.2d0 1.3d0) 0.5297814294512991d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.9d0 1.2d0 1.3d0) 0.9385293972244306d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.99d0 1.2d0 1.3d0) 0.9968864383412543d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 0.999d0 1.2d0 1.3d0) 0.9998437928330677d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-P 1.0d0 1.2d0 1.3d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.0d0 1.2d0 1.3d0) 1.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (BETA-Q 1.d-100 1.2d0 1.3d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.001d0 1.2d0 1.3d0) 0.9996623699574955d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.01d0 1.2d0 1.3d0) 0.9946568273596107d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.1d0 1.2d0 1.3d0) 0.9166002171693252d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.325d0 1.2d0 1.3d0) 0.671301345819416d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.5d0 1.2d0 1.3d0) 0.4702185705487009d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.9d0 1.2d0 1.3d0) 0.061470602775569344d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.99d0 1.2d0 1.3d0) 0.0031135616587456196d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 0.999d0 1.2d0 1.3d0) 1.5620716693236575d-4 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-Q 1.0d0 1.2d0 1.3d0) 0.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.0d0 1.2d0 1.3d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (BETA-PINV 1.344349446564896d-120 1.2d0 1.3d0) 1.d-100 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 3.376300425045358d-4 1.2d0 1.3d0) 0.001d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.005343172640389295d0 1.2d0 1.3d0) 0.01d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.08339978283067484d0 1.2d0 1.3d0) 0.1d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.3286986541805839d0 1.2d0 1.3d0) 0.325d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.5297814294512991d0 1.2d0 1.3d0) 0.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.9385293972244306d0 1.2d0 1.3d0) 0.9d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.9968864383412543d0 1.2d0 1.3d0) 0.99d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.9998437928330677d0 1.2d0 1.3d0) 0.999d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 1.0d0 1.2d0 1.3d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 1.0d0 1.2d0 1.3d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (BETA-QINV 1.0d0 1.2d0 1.3d0) 0.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.9996623699574955d0 1.2d0 1.3d0) 0.001d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.9946568273596107d0 1.2d0 1.3d0) 0.01d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.9166002171693252d0 1.2d0 1.3d0) 0.1d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.671301345819416d0 1.2d0 1.3d0) 0.325d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.4702185705487009d0 1.2d0 1.3d0) 0.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.061470602775569344d0 1.2d0 1.3d0) 0.9d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.0031135616587456196d0 1.2d0 1.3d0) 0.99d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 1.5620716693236575d-4 1.2d0 1.3d0) 0.999d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-QINV 0.0d0 1.2d0 1.3d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (BETA-PINV 0.025d0 2133.0d0 7868.0d0) 0.20530562929915866d0 +TEST-TOL6+))
......@@ -34,18 +34,58 @@
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.43939128946772227d0)
(MULTIPLE-VALUE-LIST (CHI-SQUARED-PDF 0.5d0 1.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.5204998778130462d0)
(MULTIPLE-VALUE-LIST (CHI-SQUARED-P 0.5d0 1.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.4795001221869538d0)
(MULTIPLE-VALUE-LIST (CHI-SQUARED-Q 0.5d0 1.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.5000000000000003d0)
(MULTIPLE-VALUE-LIST
(CHI-SQUARED-PINV 0.5204998778130463d0 1.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.5000000000000003d0)
(MULTIPLE-VALUE-LIST
(CHI-SQUARED-QINV 0.4795001221869537d0 1.0d0))))
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 0.0d0 13.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 1.d-100 13.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 0.001d0 13.0d0) 1.86631102655846d-25 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 0.01d0 13.0d0) 5.878822485045298d-19 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 0.1d0 13.0d0) 1.7879698335855542d-12 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 0.325d0 13.0d0) 3.446113137799052d-9 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 1.0d0 13.0d0) 3.8347347351359515d-6 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 1.5d0 13.0d0) 4.317183892010419d-5 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 2.0d0 13.0d0) 2.2625008465604717d-4 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 10.0d0 13.0d0) 0.30606563201925113d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 20.0d0 13.0d0) 0.9047897439219085d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 100.0d0 13.0d0) 0.9999999999999983d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 1000.0d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-P 10000.0d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 0.0d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 1.d-100 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 0.001d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 0.01d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 0.1d0 13.0d0) 0.999999999998212d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 0.325d0 13.0d0) 0.9999999965538868d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 1.0d0 13.0d0) 0.9999961652652649d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 1.5d0 13.0d0) 0.9999568281610799d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 2.0d0 13.0d0) 0.999773749915344d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 10.0d0 13.0d0) 0.6939343679807489d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 20.0d0 13.0d0) 0.09521025607809151d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 100.0d0 13.0d0) 1.6590260807085882d-15 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 1000.0d0 13.0d0) 1.7485119154486024d-205 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-Q 10000.0d0 13.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.0d0 13.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 1.86631102655846d-25 13.0d0) 0.001d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 5.878822485045298d-19 13.0d0) 0.01d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 1.7879698335855542d-12 13.0d0) 0.1d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 3.446113137799052d-9 13.0d0) 0.325d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 3.8347347351359515d-6 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 4.317183892010419d-5 13.0d0) 1.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 2.2625008465604717d-4 13.0d0) 2.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.30606563201925113d0 13.0d0) 10.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.9047897439219085d0 13.0d0) 20.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.9999999999999983d0 13.0d0) 100.0d0 0.01d0)
(ASSERT-posinf (CHI-SQUARED-PINV 1.0d0 13.0d0))
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.1932381452061236d0 1.5d0) 0.21198009293179954d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 4.83d-8 19.19d0) 1.6322801868602668d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.05d0 1263131.0d0) 1260517.7711333886d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-PINV 0.05d0 2526262.0d0) 2522565.864973351d0 +TEST-TOL6+)
(ASSERT-posinf (CHI-SQUARED-QINV 0.0d0 13.0d0))
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 1.6590260807085882d-15 13.0d0) 100.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.09521025607809151d0 13.0d0) 20.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.6939343679807489d0 13.0d0) 10.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.999773749915344d0 13.0d0) 2.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.9999568281610799d0 13.0d0) 1.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.9999961652652649d0 13.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.9999999965538868d0 13.0d0) 0.325d0 1.d-6)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 0.999999999998212d0 13.0d0) 0.1d0 1.d-5)
(ASSERT-TO-TOLERANCE (CHI-SQUARED-QINV 1.0d0 13.0d0) 0.0d0 +DBL-EPSILON+))
......@@ -35,12 +35,16 @@
(LIST 0.5641895835477557d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-POWER-PDF 0.0d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.9213503964748571d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-POWER-P 1.0d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.07864960352514248d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-POWER-Q 1.0d0 1.0d0 2.0d0))))
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P -1000.0d0 0.7d0 1.8d0) 0.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P -0.1d0 0.7d0 1.8d0) 0.42053490828675155d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P -1.d-32 0.7d0 1.8d0) 0.5d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P 0.1d0 0.7d0 1.8d0) 0.5794650917132484d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P 1.d-32 0.7d0 1.8d0) 0.5d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-P 1000.0d0 0.7d0 1.8d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q -1000.0d0 0.7d0 1.8d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q -0.1d0 0.7d0 1.8d0) 0.5794650917132484d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q -1.d-32 0.7d0 1.8d0) 0.5d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q 0.1d0 0.7d0 1.8d0) 0.42053490828675155d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q 1.d-32 0.7d0 1.8d0) 0.5d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-POWER-Q 1000.0d0 0.7d0 1.8d0) 0.0d0 +TEST-TOL6+))
......@@ -19,36 +19,29 @@
(in-package :gsl)
(LISP-UNIT:DEFINE-TEST EXPONENTIAL
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST 0.0025828444588394794d0 18.145581427987647d0
12.636598054339759d0 0.5424387252062355d0
14.624994234158105d0 7.236607929535993d0
0.4345362449683603d0 2.95303920904529d0
6.161052939065796d0 3.011686333539114d0
2.7451079819355364d0))
(MULTIPLE-VALUE-LIST
(LET ((RNG (MAKE-RANDOM-NUMBER-GENERATOR +MT19937+ 0)))
(LOOP FOR I FROM 0 TO 10 COLLECT
(sample rng :exponential :mu 10.0d0)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.1d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-PDF
0.0d0 10.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.3934693402873666d0)
(MULTIPLE-VALUE-LIST (EXPONENTIAL-P 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.6065306597126334d0)
(MULTIPLE-VALUE-LIST (EXPONENTIAL-Q 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 1.0d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-PINV
0.3934693402873666d0
2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 1.0d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-QINV
0.6065306597126334d0
2.0d0))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST 0.0025828444588394794d0 18.145581427987647d0
12.636598054339759d0 0.5424387252062355d0
14.624994234158105d0 7.236607929535993d0
0.4345362449683603d0 2.95303920904529d0
6.161052939065796d0 3.011686333539114d0
2.7451079819355364d0))
(MULTIPLE-VALUE-LIST
(LET ((RNG (MAKE-RANDOM-NUMBER-GENERATOR +MT19937+ 0)))
(LOOP FOR I FROM 0 TO 10 COLLECT
(sample rng :exponential :mu 10.0d0)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.1d0)
(MULTIPLE-VALUE-LIST
(EXPONENTIAL-PDF
0.0d0 10.0d0)))
(ASSERT-TO-TOLERANCE (EXPONENTIAL-P 0.1d0 0.7d0) 0.13312210024981838d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-P 1.d-32 0.7d0) 1.4285714285714287d-32 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-P 1000.0d0 0.7d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-Q 0.1d0 0.7d0) 0.8668778997501816d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-Q 1.d-32 0.7d0) 1.0d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-Q 1000.0d0 0.7d0) 0.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-PINV 0.13d0 0.7d0) 0.09748344713345536d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-PINV 1.42d-32 0.7d0) 9.94d-33 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-QINV 0.86d0 0.7d0) 0.10557602281420854d0 +TEST-TOL0+)
(ASSERT-TO-TOLERANCE (EXPONENTIAL-QINV 0.99999d0 0.7d0) 7.000035000233335d-6 +TEST-TOL6+))
This diff is collapsed.
......@@ -46,17 +46,111 @@
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.475614712250357d0)
(MULTIPLE-VALUE-LIST (GAMMA-PDF 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.04877057549928599d0)
(MULTIPLE-VALUE-LIST (GAMMA-P 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.951229424500714d0)
(MULTIPLE-VALUE-LIST (GAMMA-Q 0.1d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.10000000000000006d0)
(MULTIPLE-VALUE-LIST
(GAMMA-PINV 0.048770575499286005d0 1.0d0 2.0d0)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.10000000000000003d0)
(MULTIPLE-VALUE-LIST
(GAMMA-QINV 0.951229424500714d0 1.0d0 2.0d0))))
(ASSERT-TO-TOLERANCE (GAMMA-P 0.0d0 1.0d0 1.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.d-100 1.0d0 1.0d0) 1.d-100 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.001d0 1.0d0 1.0d0) 9.995001666250082d-4 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.01d0 1.0d0 1.0d0) 0.009950166250831947d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.1d0 1.0d0 1.0d0) 0.09516258196404043d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.325d0 1.0d0 1.0d0) 0.2774726463579278d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.0d0 1.0d0 1.0d0) 0.6321205588285577d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.5d0 1.0d0 1.0d0) 0.7768698398515702d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 2.0d0 1.0d0 1.0d0) 0.8646647167633873d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10.0d0 1.0d0 1.0d0) 0.9999546000702375d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 20.0d0 1.0d0 1.0d0) 0.9999999979388464d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 100.0d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1000.0d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10000.0d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.0d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.d-100 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.001d0 1.0d0 1.0d0) 0.999000499833375d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.01d0 1.0d0 1.0d0) 0.9900498337491681d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.1d0 1.0d0 1.0d0) 0.9048374180359595d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.325d0 1.0d0 1.0d0) 0.7225273536420722d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.0d0 1.0d0 1.0d0) 0.36787944117144233d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.5d0 1.0d0 1.0d0) 0.22313016014842982d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 2.0d0 1.0d0 1.0d0) 0.1353352832366127d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10.0d0 1.0d0 1.0d0) 4.5399929762484854d-5 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 20.0d0 1.0d0 1.0d0) 2.061153622438558d-9 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 100.0d0 1.0d0 1.0d0) 3.720075976020836d-44 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1000.0d0 1.0d0 1.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10000.0d0 1.0d0 1.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.0d0 1.0d0 10.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.d-100 1.0d0 10.0d0) 1.d-101 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.001d0 1.0d0 10.0d0) 9.99950001666625d-5 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.01d0 1.0d0 10.0d0) 9.995001666250082d-4 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.1d0 1.0d0 10.0d0) 0.009950166250831947d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.325d0 1.0d0 10.0d0) 0.031977550168693955d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.0d0 1.0d0 10.0d0) 0.09516258196404043d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.5d0 1.0d0 10.0d0) 0.1392920235749422d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 2.0d0 1.0d0 10.0d0) 0.18126924692201815d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10.0d0 1.0d0 10.0d0) 0.6321205588285577d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 20.0d0 1.0d0 10.0d0) 0.8646647167633873d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 100.0d0 1.0d0 10.0d0) 0.9999546000702375d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1000.0d0 1.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10000.0d0 1.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.0d0 1.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.d-100 1.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.001d0 1.0d0 10.0d0) 0.9999000049998333d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.01d0 1.0d0 10.0d0) 0.999000499833375d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.1d0 1.0d0 10.0d0) 0.9900498337491681d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.325d0 1.0d0 10.0d0) 0.968022449831306d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.0d0 1.0d0 10.0d0) 0.9048374180359595d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.5d0 1.0d0 10.0d0) 0.8607079764250578d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 2.0d0 1.0d0 10.0d0) 0.8187307530779818d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10.0d0 1.0d0 10.0d0) 0.36787944117144233d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 20.0d0 1.0d0 10.0d0) 0.1353352832366127d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 100.0d0 1.0d0 10.0d0) 4.5399929762484854d-5 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1000.0d0 1.0d0 10.0d0) 3.720075976020836d-44 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10000.0d0 1.0d0 10.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.0d0 17.0d0 10.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.d-100 17.0d0 10.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.001d0 17.0d0 10.0d0) 2.8111917404042285d-83 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.01d0 17.0d0 10.0d0) 2.808803246519859d-66 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.1d0 17.0d0 10.0d0) 2.7850299808749213d-49 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 0.325d0 17.0d0 10.0d0) 1.3728365324512585d-40 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.0d0 17.0d0 10.0d0) 2.5581193229254426d-32 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1.5d0 17.0d0 10.0d0) 2.4042044117542236d-29 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 2.0d0 17.0d0 10.0d0) 3.0509292621789857d-27 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10.0d0 17.0d0 10.0d0) 1.094920130378183d-15 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 20.0d0 17.0d0 10.0d0) 5.606050961731617d-11 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 100.0d0 17.0d0 10.0d0) 0.027041609784801128d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 1000.0d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-P 10000.0d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.0d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.d-100 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.001d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.01d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.1d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 0.325d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.0d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1.5d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 2.0d0 17.0d0 10.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10.0d0 17.0d0 10.0d0) 0.9999999999999989d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 20.0d0 17.0d0 10.0d0) 0.9999999999439395d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 100.0d0 17.0d0 10.0d0) 0.9729583902151989d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 1000.0d0 17.0d0 10.0d0) 2.112009516339486d-25 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-Q 10000.0d0 17.0d0 10.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.0d0 1.0d0 1.0d0) 0.0d0 +DBL-EPSILON+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 1.d-100 1.0d0 1.0d0) 1.d-100 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 9.995001666250082d-4 1.0d0 1.0d0) 0.001d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.009950166250831947d0 1.0d0 1.0d0) 0.01d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.09516258196404043d0 1.0d0 1.0d0) 0.1d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.2774726463579278d0 1.0d0 1.0d0) 0.325d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.6321205588285577d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.7768698398515702d0 1.0d0 1.0d0) 1.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.8646647167633873d0 1.0d0 1.0d0) 2.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.9999546000702375d0 1.0d0 1.0d0) 10.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.9999999979388464d0 1.0d0 1.0d0) 20.0d0 (* 100 +TEST-TOL6+))
(ASSERT-posinf (GAMMA-PINV 1.0d0 1.0d0 1.0d0))
(ASSERT-TO-TOLERANCE (GAMMA-PINV 0.1d0 11.887411491530846D0 1.0D0) 7.73788447848618d0 +TEST-TOL1+)
(ASSERT-posinf (GAMMA-QINV 0.0d0 1.0d0 1.0d0))
(ASSERT-TO-TOLERANCE (GAMMA-QINV 2.061153622438558d-9 1.0d0 1.0d0) 20.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 4.5399929762484854d-5 1.0d0 1.0d0) 10.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.1353352832366127d0 1.0d0 1.0d0) 2.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.22313016014842982d0 1.0d0 1.0d0) 1.5d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.36787944117144233d0 1.0d0 1.0d0) 1.0d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.7225273536420722d0 1.0d0 1.0d0) 0.325d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.9048374180359595d0 1.0d0 1.0d0) 0.1d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.9900498337491681d0 1.0d0 1.0d0) 0.01d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 0.999000499833375d0 1.0d0 1.0d0) 0.001d0 +TEST-TOL6+)
(ASSERT-TO-TOLERANCE (GAMMA-QINV 1.0d0 1.0d0 1.0d0) 0.0d0 +DBL-EPSILON+))
......@@ -283,10 +283,8 @@
(assert-to-tolerance (tdist-Pinv 8.40941797784686861d-01 300.0d0) 1.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Pinv 9.32666983425369137d-01 300.0d0) 1.5d0 +test-tol6+)
(assert-to-tolerance (tdist-Pinv 9.76799239508425455d-01 300.0d0) 2.0d0 +test-tol6+)
(lisp-unit::assert-true
(let ((val (tdist-Pinv 1.0d0 300.0d0)))
(and (infinityp val) (plusp val))))
(assert-posinf (tdist-Pinv 1.0d0 300.0d0))
(assert-to-tolerance (tdist-Qinv 5.00000000000000000d-01 300.0d0) 0.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 4.99601390099057051d-1 300.0d0) 0.001d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 4.96013966979440912d-1 300.0d0) 0.01d0 +test-tol6+)
......@@ -298,9 +296,7 @@
(assert-to-tolerance (tdist-Qinv 8.279313677d-21 300.0d0) 10.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 1.93159812815803978d-57 300.0d0) 20.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 1.02557519997736154d-232 300.0d0) 100.0d0 +test-tol6+)
(lisp-unit::assert-true
(let ((val (tdist-Qinv 0.0d0 300.0d0)))
(and (infinityp val) (plusp val))))
(assert-posinf (tdist-Qinv 0.0d0 300.0d0))
(assert-to-tolerance (tdist-Pinv 4.99601390099057051d-01 300.0d0) -0.001d0 +test-tol6+)
(assert-to-tolerance (tdist-Pinv 4.96013966979440912d-01 300.0d0) -0.01d0 +test-tol6+)
......@@ -312,9 +308,7 @@
(assert-to-tolerance (tdist-Pinv 8.279313675556272534d-21 300.0d0) -10.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Pinv 1.93159812815803978d-57 300.0d0) -20.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Pinv 1.02557519997736154d-232 300.0d0) -100.0d0 +test-tol6+)
(lisp-unit::assert-true
(let ((val (tdist-Pinv 0.0d0 300.0d0)))
(and (infinityp val) (minusp val))))
(assert-neginf (tdist-Pinv 0.0d0 300.0d0))
(assert-to-tolerance (tdist-Qinv 5.00398609900942949d-1 300.0d0) -0.001d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 5.03986033020559088d-1 300.0d0) -0.01d0 +test-tol6+)
......@@ -323,6 +317,4 @@
(assert-to-tolerance (tdist-Qinv 8.40941797784686862d-1 300.0d0) -1.0d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 9.32666983425369137d-1 300.0d0) -1.5d0 +test-tol6+)
(assert-to-tolerance (tdist-Qinv 9.76799239508425455d-1 300.0d0) -2.0d0 +test-tol6+)
(lisp-unit::assert-true
(let ((val (tdist-Qinv 1.0d0 300.0d0)))
(and (infinityp val) (minusp val)))))
(assert-neginf (tdist-Qinv 1.0d0 300.0d0)))
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment