


Bits error versus u1



Bits error versus u2
Initial program 0.4
Initial simplification0.4
rmApplied add-sqr-sqrt0.6
Applied associate-*l*0.6
Final simplification0.6
herbie shell --seed 2018362 +o rules:numerics
(FPCore (u1 u2)
:name "normal distribution"
:pre (and (<= 0 u1 1) (<= 0 u2 1))
(+ (* (* (/ 1 6) (pow (* -2 (log u1)) 0.5)) (cos (* (* 2 PI) u2))) 0.5))
Time bar (total: 20.7s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 0.4b
Found 4 expressions with local error:
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 0.2b | (* PI (* 2 u2)) |
| 0.1b | (pow (* -2 (log u1)) 0.5) |
| 0.0b | (fma (cos (* PI (* 2 u2))) (* (pow (* -2 (log u1)) 0.5) 1/6) 0.5) |
4 calls. Slowest were:
| 5.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 1.0ms | (pow (* -2 (log u1)) 0.5) |
| 1.0ms | (* PI (* 2 u2)) |
4 calls. Slowest were:
| 775.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 378.0ms | (fma (cos (* PI (* 2 u2))) (* (pow (* -2 (log u1)) 0.5) 1/6) 0.5) |
| 349.0ms | (pow (* -2 (log u1)) 0.5) |
| 49.0ms | (* PI (* 2 u2)) |
18 calls. Slowest were:
| 160.0ms | (- (+ 0.5 (* 1/6 (pow (* (pow -2 1.0) (pow (log u1) 1.0)) 0.5))) (* 1/3 (* (* (pow u2 2) (pow PI 2)) (pow (* (pow -2 1.0) (pow (log u1) 1.0)) 0.5)))) |
| 87.0ms | (+ (* 1/6 (* (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5) (cos (* 2 (* u2 PI))))) 0.5) |
| 83.0ms | (pow (* -2 (- (log -1) (log (/ -1 u1)))) 0.5) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.6b | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 0.5b | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) |
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
4 calls. Slowest were:
| 42.0ms | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 29.0ms | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) |
| 8.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
4 calls. Slowest were:
| 962.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 940.0ms | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 883.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 824.0ms | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) |
27 calls. Slowest were:
| 742.0ms | (* (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 64.0ms | (+ (log (* (pow (* -2 (log u1)) 0.5) 1/6)) (log (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 61.0ms | (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.4b | (expm1 (log1p (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 0.2b | (* PI (* 2 u2)) |
| 0.2b | (log1p (* (pow (* -2 (log u1)) 0.5) 1/6)) |
4 calls. Slowest were:
| 8.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 2.0ms | (* PI (* 2 u2)) |
| 1.0ms | (log1p (* (pow (* -2 (log u1)) 0.5) 1/6)) |
4 calls. Slowest were:
| 1.2s | (expm1 (log1p (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 1.1s | (log1p (* (pow (* -2 (log u1)) 0.5) 1/6)) |
| 935.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 48.0ms | (* PI (* 2 u2)) |
17 calls. Slowest were:
| 72.0ms | (log (+ (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) 1)) |
| 58.0ms | (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) |
| 47.0ms | (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.4b | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
| 0.4b | (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)) |
| 0.2b | (* PI (* 2 u2)) |
| 0.1b | (pow (* -2 (log u1)) 0.5) |
4 calls. Slowest were:
| 14.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)) |
| 7.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
| 1.0ms | (pow (* -2 (log u1)) 0.5) |
4 calls. Slowest were:
| 919.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
| 886.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)) |
| 346.0ms | (pow (* -2 (log u1)) 0.5) |
| 83.0ms | (* PI (* 2 u2)) |
30 calls. Slowest were:
| 502.0ms | (* (* (* (sqrt (pow (* -2 (log u1)) 0.5)) (sqrt (pow (* -2 (log u1)) 0.5))) (sqrt (pow (* -2 (log u1)) 0.5))) (* (* (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) |
| 127.0ms | (pow (* -2 (- (log -1) (log (/ -1 u1)))) 0.5) |
| 70.0ms | (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 0b
0% (0.5b remaining)
Error of 0.6b against oracle of 0.0b and baseline of 0.6b