


Bits error versus u1



Bits error versus u2
Initial program 0.4
Initial simplification0.4
Taylor expanded around -inf 62.0
Simplified0.4
rmApplied add-cbrt-cube0.6
rmApplied add-sqr-sqrt0.6
Applied associate-*l*0.6
Final simplification0.6
herbie shell --seed 2018365 +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: 23.0s)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.3b | (* PI (* 2 u2)) |
| 0.2b | (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:
| 9.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 3.0ms | (pow (* -2 (log u1)) 0.5) |
| 2.0ms | (* PI (* 2 u2)) |
4 calls. Slowest were:
| 881.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 469.0ms | (fma (cos (* PI (* 2 u2))) (* (pow (* -2 (log u1)) 0.5) 1/6) 0.5) |
| 380.0ms | (pow (* -2 (log u1)) 0.5) |
| 88.0ms | (* PI (* 2 u2)) |
18 calls. Slowest were:
| 213.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)))) |
| 146.0ms | (pow (* -2 (- (log -1) (log (/ -1 u1)))) 0.5) |
| 81.0ms | (+ (* 1/6 (* (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5) (cos (* 2 (* u2 PI))))) 0.5) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 0.3b | (* PI (* 2 u2)) |
| 0.2b | (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:
| 6.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:
| 824.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 434.0ms | (fma (cos (* PI (* 2 u2))) (* (pow (* -2 (log u1)) 0.5) 1/6) 0.5) |
| 390.0ms | (pow (* -2 (log u1)) 0.5) |
| 70.0ms | (* PI (* 2 u2)) |
18 calls. Slowest were:
| 190.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)))) |
| 136.0ms | (pow (* -2 (- (log -1) (log (/ -1 u1)))) 0.5) |
| 89.0ms | (+ (* 1/6 (* (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5) (cos (* 2 (* u2 PI))))) 0.5) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.5b | (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:
| 29.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))) |
| 28.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:
| 1.0s | (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))) |
| 969.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
| 950.0ms | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (pow (* -2 (log u1)) 0.5) 1/6)) |
| 936.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
27 calls. Slowest were:
| 637.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))) |
| 82.0ms | (+ (log (* (pow (* -2 (log u1)) 0.5) 1/6)) (log (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 68.0ms | (* 1/6 (pow (* (pow (- (log -1) (log (/ -1 u1))) 1.0) (pow -2 1.0)) 0.5)) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.5b | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (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) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) |
| 0.4b | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
| 0.4b | (* (pow (* -2 (log u1)) 0.5) 1/6) |
4 calls. Slowest were:
| 74.0ms | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 66.0ms | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) |
| 10.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
4 calls. Slowest were:
| 1.1s | (cbrt (* (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) (* (pow (* -2 (log u1)) 0.5) 1/6))) |
| 995.0ms | (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6) |
| 983.0ms | (* (* (pow (* -2 (log u1)) 0.5) 1/6) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) |
| 973.0ms | (* (pow (* -2 (log u1)) 0.5) 1/6) |
30 calls. Slowest were:
| 738.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)) (* (* (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6))) (* (sqrt (pow (* -2 (log u1)) 0.5)) (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)))) |
| 174.0ms | (+ (log (* (pow (* -2 (log u1)) 0.5) 1/6)) (+ (log (sqrt (pow (* -2 (log u1)) 0.5))) (log (* (sqrt (pow (* -2 (log u1)) 0.5)) 1/6)))) |
| 151.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)) (* (* (* (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)))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 0.0b
0% (0.5b remaining)
Error of 0.6b against oracle of 0.0b and baseline of 0.6b