



Bits error versus x
| Original | 19.4 |
|---|---|
| Target | 0.6 |
| Herbie | 0.6 |
Initial program 19.4
rmApplied frac-sub19.3
Simplified19.3
rmApplied flip--19.2
Applied associate-/l/19.2
Simplified0.7
rmApplied pow10.7
Applied pow10.7
Applied pow10.7
Applied pow-prod-down0.7
Applied pow-prod-down0.7
Simplified0.6
rmApplied add-cbrt-cube0.6
Final simplification0.6
herbie shell --seed 2018365 +o rules:numerics
(FPCore (x)
:name "2isqrt (example 3.6)"
:herbie-target
(/ 1 (+ (* (+ x 1) (sqrt x)) (* x (sqrt (+ x 1)))))
(- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))))
Time bar (total: 13.4s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 17.4b
Found 4 expressions with local error:
| 2.8b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 0.2b | (/ 1 (sqrt x)) |
| 0.1b | (/ 1 (sqrt (+ x 1))) |
| 0.0b | (sqrt (+ x 1)) |
4 calls. Slowest were:
| 13.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 2.0ms | (/ 1 (sqrt (+ x 1))) |
| 1.0ms | (/ 1 (sqrt x)) |
4 calls. Slowest were:
| 38.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 28.0ms | (/ 1 (sqrt (+ x 1))) |
| 17.0ms | (/ 1 (sqrt x)) |
| 16.0ms | (sqrt (+ x 1)) |
92 calls. Slowest were:
| 213.0ms | (fma (- (sqrt (- x 1))) (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1)))))) |
| 187.0ms | (fma (sqrt (/ 1 (sqrt x))) (sqrt (/ 1 (sqrt x))) (- (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1))))))) |
| 183.0ms | (fma (- (sqrt (+ (* x x) (- (* 1 1) (* x 1))))) (/ 1 (sqrt (+ (pow x 3) (pow 1 3)))) (* (sqrt (+ (* x x) (- (* 1 1) (* x 1)))) (/ 1 (sqrt (+ (pow x 3) (pow 1 3)))))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 16.8b
Found 4 expressions with local error:
| 2.9b | (- (sqrt (+ x 1)) (sqrt x)) |
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
| 0.1b | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
| 0.0b | (sqrt (+ x 1)) |
4 calls. Slowest were:
| 15.0ms | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
| 9.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 4.0ms | (* (sqrt x) (sqrt (+ x 1))) |
4 calls. Slowest were:
| 57.0ms | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
| 33.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 28.0ms | (* (sqrt x) (sqrt (+ x 1))) |
| 15.0ms | (sqrt (+ x 1)) |
106 calls. Slowest were:
| 509.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 116.0ms | (/ (* (* (- (sqrt (+ x 1)) (sqrt x)) (- (sqrt (+ x 1)) (sqrt x))) (- (sqrt (+ x 1)) (sqrt x))) (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1))))) |
| 83.0ms | (/ (* (* (- (sqrt (+ x 1)) (sqrt x)) (- (sqrt (+ x 1)) (sqrt x))) (- (sqrt (+ x 1)) (sqrt x))) (* (* (* (sqrt x) (sqrt (+ x 1))) (* (sqrt x) (sqrt (+ x 1)))) (* (sqrt x) (sqrt (+ x 1))))) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.3b
Found 4 expressions with local error:
| 0.3b | (/ 1 (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x)))) |
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
| 0.1b | (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 0.0b | (+ (sqrt (+ x 1)) (sqrt x)) |
4 calls. Slowest were:
| 26.0ms | (/ 1 (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x)))) |
| 25.0ms | (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 8.0ms | (* (sqrt x) (sqrt (+ x 1))) |
4 calls. Slowest were:
| 71.0ms | (/ 1 (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x)))) |
| 62.0ms | (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 23.0ms | (* (sqrt x) (sqrt (+ x 1))) |
| 21.0ms | (+ (sqrt (+ x 1)) (sqrt x)) |
71 calls. Slowest were:
| 419.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 193.0ms | (* (* (sqrt x) (sqrt (- (* x x) (* 1 1)))) (+ (pow (sqrt (+ x 1)) 3) (pow (sqrt x) 3))) |
| 178.0ms | (* (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) (* (* (+ (sqrt (+ x 1)) (sqrt x)) (+ (sqrt (+ x 1)) (sqrt x))) (+ (sqrt (+ x 1)) (sqrt x)))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.3b | (/ 1 (pow (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) 1)) |
| 0.1b | (fma (sqrt x) x (sqrt x)) |
| 0.0b | (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) |
| 0.0b | (sqrt (+ x 1)) |
4 calls. Slowest were:
| 1.0ms | (/ 1 (pow (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) 1)) |
| 1.0ms | (sqrt (+ x 1)) |
| 0.0ms | (fma (sqrt x) x (sqrt x)) |
4 calls. Slowest were:
| 61.0ms | (/ 1 (pow (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) 1)) |
| 45.0ms | (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) |
| 28.0ms | (fma (sqrt x) x (sqrt x)) |
| 11.0ms | (sqrt (+ x 1)) |
25 calls. Slowest were:
| 54.0ms | (- (+ (* 1/2 x) 1) (* 1/8 (pow x 2))) |
| 28.0ms | (pow (cbrt (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x)))) 1) |
| 27.0ms | (cbrt (pow (fma x (sqrt (+ x 1)) (fma (sqrt x) x (sqrt x))) 1)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0.0b
0% (0.6b remaining)
Error of 0.6b against oracle of 0.0b and baseline of 0.6b