Commit e4ebbf49 authored by Jay Belanger's avatar Jay Belanger
Browse files

(math-random-digit): Rename to `math-random-three-digit-number'.

(math-random-digits): Don't depend on representation of integers.
parent 33f400c3
......@@ -615,7 +615,8 @@
;;; Avoid various pitfalls that may lurk in the built-in (random) function!
;;; Shuffling algorithm from Numerical Recipes, section 7.1.
(defvar math-random-last)
(defun math-random-digit ()
(defun math-random-three-digit-number ()
"Return a random three digit number."
(let (i)
(or (and (boundp 'var-RandSeed) (eq var-RandSeed math-last-RandSeed))
(math-init-random-base))
......@@ -635,17 +636,17 @@
;;; Produce an N-digit random integer.
(defun math-random-digits (n)
(cond ((<= n 6)
(math-scale-right (+ (* (math-random-digit) 1000) (math-random-digit))
(- 6 n)))
(t (let* ((slop (% (- 900003 n) 3))
(i (/ (+ n slop) 3))
(digs nil))
(while (> i 0)
(setq digs (cons (math-random-digit) digs)
i (1- i)))
(math-normalize (math-scale-right (cons 'bigpos digs)
slop))))))
"Produce a random N digit integer."
(let* ((slop (% (- 3 (% n 3)) 3))
(i (/ (+ n slop) 3))
(rnum 0))
(while (> i 0)
(setq rnum
(math-add
(math-random-three-digit-number)
(math-mul rnum 1000)))
(setq i (1- i)))
(math-normalize (math-scale-right rnum slop))))
;;; Produce a uniformly-distributed random float 0 <= N < 1.
(defun math-random-float ()
......
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