Lines Matching refs:fl
17 (let-values ([(p q s?) (cond [(p . fl< . 0.5) (values p (fl- 1.0 p) #f)]
18 [else (values (fl- 1.0 p) p #t)])])
20 (define r (fl/ p q))
21 (define g (fl* r (fl+ n 1.0)))
27 (cond [(u . fl< . f) k]
28 [(k . fl> . 110.0) (reject)]
29 [else (let ([k (fl+ k 1.0)])
30 (loop k (fl* f (fl- (fl/ g k) r)) (fl- u f)))]))))
31 (if s? (fl- n k) k)))))
36 (let-values ([(p q s?) (cond [(p . fl< . 0.5) (values p (fl- 1.0 p) #f)]
37 [else (values (fl- 1.0 p) p #t)])])
39 (define m (flfloor (fl* (fl+ n 1.0) p)))
41 (define b (fl+ 1.15 (fl* 2.53 σ)))
42 (define a (+ -0.0873 (fl* 0.0248 b) (fl* 0.01 p)))
43 (define c (fl+ 0.5 (fl* n p)))
44 (define α (fl* σ (fl+ 2.83 (fl/ 5.1 b))))
45 (define vr (fl- 0.92 (fl/ 4.2 b)))
51 (define u (fl- (random) 0.5))
52 (define us (fl- 0.5 (flabs u)))
53 (define k (flfloor (fl+ c (fl* u (fl+ b (fl/ (fl* 2.0 a) us))))))
54 (cond [(or (k . fl< . 0.0) (k . fl> . n)) (loop)]
55 [(and (us . fl>= . 0.07) (v . fl<= . vr)) k]
57 (let ([v (fl* v (fl/ α (fl+ b (fl/ a (fl* us us)))))])
59 (fllog-factorial (fl- n m))
61 (- (fllog-factorial (fl- n k)))
62 (fl* (fl- k m) (fllog (fl/ p q)))))
63 (cond [((fllog v) . fl<= . h) k]
65 (if s? (fl- n k) k)))))
69 (define q (fl- 1.0 p))
70 (define μ (fl- (fl* (fl+ n 1.0) p) 0.5))
72 (define γ (fl/ (fl- q p) σ))
77 (define k (flround (fl+ μ (fl* σ (fl+ z (fl/ (fl* γ (fl- (fl* z z) 1.0)) 6.0))))))
78 (if (and (k . fl>= . 0.0) (k . fl<= . n)) k (loop))))))
84 (define q (fl- 1.0 p))
85 (fl/ (fl* 0.4784 (fl+ (fl* p p) (fl* q q))) (flsqrt (* n p q))))
90 [(or (not (integer? n)) (n . fl< . 0.0) (p . fl< . 0.0) (p . fl> . 1.0))
92 [(or (fl= n 0.0) (fl= p 0.0))
94 [(fl= p 1.0)
96 [(and (n . fl> . 1e8)
97 ((flbinomial-normal-appx-error-bound n p) . fl< . (flexp -10.0)))
99 [((fl* n (flmin p (fl- 1.0 p))) . fl>= . 10.0)