core.clj

(ns pi.core
  (:gen-class))

(set! *unchecked-math* :warn-on-boxed)

(defn gen-point [] [(- (* 2.0 (rand)) 1.0) (- (* 2.0 (rand)) 1.0)])

(defn test-point [[^double x ^double y]] (< (+ (* x x) (* y y)) 1.0))

(defn gen-and-test-point []
	(let [x (- (* 2.0 (rand)) 1.0)
		y (- (* 2.0 (rand)) 1.0)]
		(if (< (+ (* x x) (* y y)) 1.0) 1 0)))

(defn sample1
	[n]
		(transduce (map (fn [i] (if (test-point (gen-point)) 1 0))) + (range n)))

(defn sample2
	[n]
	(transduce (map (fn [i] (gen-and-test-point))) + (range n)))

(defn estimate-pi 
	[n, sample] (* 4.0 (/ (sample n) n)))

(defn -main
  "I don't do a whole lot ... yet."
  [& args]
  (time (estimate-pi 10000000 sample1))
  (time (estimate-pi 10000000 sample2)))