0
@@ -1,5 +1,6 @@
0
 (in-package :randist)
0
 
0
+(declaim (optimize (speed 3) (safety 1)))
0
 
0
 ;;    The negative binomial distribution has the form,
0
 ;;    prob(k) =  Gamma(n + k)/(Gamma(n) Gamma(k + 1))  p^n (1-p)^k 
0
@@ -16,5 +17,8 @@
0
 
0
 
0
 (defun random-negative-binomial (p n)
0
+  (declare (type double-float p)
0
+	   (integer n))
0
   (let ((X (random-gamma (coerce n 'double-float) 1d0)))
0
-    (random-poisson (* X (/ (- 1d0 p) p)))))
0
\ No newline at end of file
0
+    (random-poisson (* X (/ (- 1d0 p) p)))))
0
+

0
@@ -27,10 +27,6 @@
0
 p(x) dx = (a/b) / (x/b)^(a+1) dx     for x >= b"
0
 
0
   (declare (type double-float a b))
0
-<<<<<<< HEAD:pareto.lisp
0
-  (let* ((x (random-uniform))
0
-=======
0
   (let* ((x (random-pos))
0
->>>>>>> experimental:pareto.lisp
0
    (z (expt x (/ -1d0 a))))
0
-    (* b z)))
0
\ No newline at end of file
0
+    (* b z)))