@@ -1,5 +1,7 @@ (in-package :randist) +(declaim (optimize (debug 3))) + (defun setup-cut-point-randist (p) (let* ((L 0) (J 0) @@ -28,7 +30,7 @@ (let ((M (length p))) #'(lambda () (let* ((U (random-uniform)) - (X (aref I (1- (floor (* M U)))))) + (X (aref I (floor (* M U))))) (tagbody start (if (<= U (aref Q (1- X))) @@ -37,7 +39,7 @@ (incf X) (go start))) end) - X))))) + (1- X)))))) cut-point.lisp @@ -34,5 +34,5 @@ make-discrete-random-var ;; Cut-point method for discrete distributions - ;; make-discrete-monotone-random-var + make-discrete-monotone-random-var )) packages.lisp 84eca18365245c027915fbad3313cc8169740b9a Leonardo Varuzza varuzza@gmail.com http://github.com/lvaruzza/cl-randist/commit/32486467102fb42487fbe7f2de9f2150a095e97b 32486467102fb42487fbe7f2de9f2150a095e97b 2008-02-24T08:33:49-08:00 2008-02-24T08:33:49-08:00 Fix cut-point 5ac1e960ce7d36f9fffd9d43d42b8de07daa4c30 Leonardo Varuzza varuzza@gmail.com