Commit ccce7a89 authored by Raymond Toy's avatar Raymond Toy
Browse files

Forgot to rename incomplete-gamma functions

!4 forgot to rename the incomplete-gamma functions in qd-bessel.lisp.
parent 3f2625e0
......@@ -37,6 +37,7 @@
:maintainer "See <http://www.common-lisp.net/project/oct>"
:licence "MIT"
:version "2013.11.26" ; Just use the date
:in-order-to ((test-op (test-op "oct/tests")))
:components
((:file "qd-package")
(:file "qd-rep" :depends-on ("qd-package"))
......@@ -82,13 +83,6 @@
; (defmethod perform ((op test-op) (c (eql (asdf:find-system :oct))))
; (oos 'test-op 'rt:do-tests))
(asdf:defsystem rt
:components
((:module "rt"
:components
((:file "rt-package")
(:file "rt")))))
(asdf:defsystem oct/tests
:depends-on (oct rt)
:version "2013.11.26" ; Just use the date
......@@ -97,7 +91,7 @@
(:file "qd-test")
(:file "rt-tests")))
(defmethod perform ((op test-op) (c (eql (asdf:find-system :oct))))
(defmethod perform ((op test-op) (c (eql (asdf:find-system :oct/tests))))
(let* ((rt-package (find-package "RT"))
(pass (funcall (intern "DO-TESTS" rt-package)))
(unexpected-failures (symbol-value (find-symbol "*UNEXPECTED-FAILURES*" rt-package)))
......
......@@ -140,7 +140,7 @@
(format t "sqrt2 = ~S~%" sqrt2)
(format t "exp/p/sqrt2 = ~S~%" exp/p/sqrt2))
(do* ((k 0 (1+ k))
(bk (/ (incomplete-gamma 1/2 p)
(bk (/ (incomplete-gamma-lower 1/2 p)
2 sqrt2 (sqrt p))
(- (/ (* bk (- k 1/2)) 2 p)
(/ exp/p/sqrt2 (ash 1 (+ k 1)))))
......@@ -265,7 +265,7 @@
(defun alpha (n z)
(let ((n (float n (realpart z))))
(/ (incomplete-gamma (1+ n) (/ z 2))
(/ (incomplete-gamma-lower (1+ n) (/ z 2))
(expt z (1+ n)))))
(defun alpha-iter (n z alpha-old)
......@@ -280,8 +280,8 @@
(defun beta (n z)
(let ((n (float n (realpart z))))
(/ (- (incomplete-gamma (1+ n) (/ z 2))
(incomplete-gamma (1+ n) (/ z -2)))
(/ (- (incomplete-gamma-lower (1+ n) (/ z 2))
(incomplete-gamma-lower (1+ n) (/ z -2)))
(expt z (1+ n)))))
(defun beta-iter (n z old-beta)
......@@ -444,7 +444,7 @@
(defun big-i (n theta z v)
(let* ((a (- 1 v n n)))
(* (expt z (- a))
(incomplete-gamma-tail a (* theta z)))))
(incomplete-gamma-upper a (* theta z)))))
(defun sum-big-ia (big-n v z)
(let ((big-n-1/2 (+ big-n 1/2))
......
(asdf:defsystem rt
:components
((:module "rt"
:components
((:file "rt-package")
(:file "rt")))))
(require :asdf)
(pushnew "./" asdf:*central-registry*)
(require :oct)
(asdf:oos 'asdf:load-op :oct :force t )
(dribble "test.log")
(asdf:oos 'asdf:test-op :oct/tests :force t)
(asdf:test-system :oct :force t)
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