Changeset 34 for trunk/monte-carlo.lisp

Timestamp:

03/11/08 01:34:13 (8 months ago)

Author:

lhealy

Message:

The classes/types in the different contexts are now gathered together
in one place, in *type-names* for the types and in *data-class-name*
for data classes, populated by #'add-data-class. Both defdata and
defmfun-all use the table and so mapping between various names is
consistent. The data class names are now different, *-double-float
and *-single-float replaces *-double and *-single. The regression
tests give the same results as before.

Files:

• trunk/monte-carlo.lisp (modified) (4 diffs)

trunk/monte-carlo.lisp

@@ -1 +1 @@
 ;; Monte Carlo Integration
 ;; Liam Healy Sat Feb  3 2007 - 17:42
-;; Time-stamp: <2008-02-17 18:29:59EST monte-carlo.lisp>
+;; Time-stamp: <2008-03-09 19:29:43EDT monte-carlo.lisp>
 ;; $Id$
 
@@ -292 +292 @@
 (defun random-walk-plain-example (&optional (nsamples 500000))
   (letm ((ws (monte-carlo-plain 3))
-         (lower (vector-double #(0.0d0 0.0d0 0.0d0)))
-         (upper (vector-double (vector pi pi pi)))
+         (lower (vector-double-float #(0.0d0 0.0d0 0.0d0)))
+         (upper (vector-double-float (vector pi pi pi)))
          (rng (random-number-generator *mt19937* 0)))
     (monte-carlo-integrate-plain monte-carlo-g lower upper nsamples rng ws)))
@@ -299 +299 @@
 (defun random-walk-miser-example (&optional (nsamples 500000))
   (letm ((ws (monte-carlo-miser 3))
-         (lower (vector-double #(0.0d0 0.0d0 0.0d0)))
-         (upper (vector-double (vector pi pi pi)))
+         (lower (vector-double-float #(0.0d0 0.0d0 0.0d0)))
+         (upper (vector-double-float (vector pi pi pi)))
          (rng (random-number-generator *mt19937* 0)))
     (monte-carlo-integrate-miser monte-carlo-g lower upper nsamples rng ws)))
@@ -306 +306 @@
 (defun random-walk-vegas-example (&optional (nsamples 500000))
   (letm ((ws (monte-carlo-vegas 3))
-         (lower (vector-double #(0.0d0 0.0d0 0.0d0)))
-         (upper (vector-double (vector pi pi pi)))
+         (lower (vector-double-float #(0.0d0 0.0d0 0.0d0)))
+         (upper (vector-double-float (vector pi pi pi)))
          (rng (random-number-generator *mt19937* 0)))
     (monte-carlo-integrate-vegas monte-carlo-g lower upper nsamples rng ws)))