



Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.7 |
| Herbie | 19.9 |
Initial program 19.7
rmApplied add-sqr-sqrt19.9
rmApplied frac-sub19.9
Simplified19.9
rmApplied add-log-exp19.9
rmApplied add-sqr-sqrt19.9
Applied log-prod19.9
Final simplification19.9
herbie shell --seed 2019004 +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: 27.2s)Debug log
1 calls. Slowest were:
| 6.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 21.2b
Found 4 expressions with local error:
| 3.2b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 0.3b | (/ 1 (sqrt x)) |
| 0.2b | (/ 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))) |
| 2.0ms | (/ 1 (sqrt x)) |
4 calls. Slowest were:
| 41.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 23.0ms | (sqrt (+ x 1)) |
| 22.0ms | (/ 1 (sqrt (+ x 1))) |
| 19.0ms | (/ 1 (sqrt x)) |
92 calls. Slowest were:
| 200.0ms | (fma (- (sqrt (- x 1))) (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1)))))) |
| 180.0ms | (fma (sqrt (/ 1 (sqrt x))) (sqrt (/ 1 (sqrt x))) (- (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1))))))) |
| 178.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)))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 20.6b
Found 4 expressions with local error:
| 3.2b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 3.2b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 0.3b | (/ 1 (sqrt x)) |
| 0.3b | (/ 1 (sqrt x)) |
4 calls. Slowest were:
| 9.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 7.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 1.0ms | (/ 1 (sqrt x)) |
4 calls. Slowest were:
| 37.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 36.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 11.0ms | (/ 1 (sqrt x)) |
| 11.0ms | (/ 1 (sqrt x)) |
140 calls. Slowest were:
| 179.0ms | (fma (- (sqrt (- x 1))) (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1)))))) |
| 174.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)))))) |
| 169.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)))))) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 20.6b
Found 4 expressions with local error:
| 3.2b | (- (sqrt (+ x 1)) (sqrt x)) |
| 3.2b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
| 0.3b | (/ 1 (sqrt x)) |
4 calls. Slowest were:
| 14.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 13.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 9.0ms | (* (sqrt x) (sqrt (+ x 1))) |
4 calls. Slowest were:
| 37.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 33.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 22.0ms | (* (sqrt x) (sqrt (+ x 1))) |
| 12.0ms | (/ 1 (sqrt x)) |
138 calls. Slowest were:
| 718.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 162.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)))))) |
| 159.0ms | (fma (* (cbrt (/ 1 (sqrt x))) (cbrt (/ 1 (sqrt x)))) (cbrt (/ 1 (sqrt x))) (- (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1))))))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 20.6b
Found 4 expressions with local error:
| 3.2b | (- (sqrt (+ x 1)) (sqrt x)) |
| 3.2b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 2.6b | (log (exp (- (sqrt (+ x 1)) (sqrt x)))) |
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
4 calls. Slowest were:
| 12.0ms | (log (exp (- (sqrt (+ x 1)) (sqrt x)))) |
| 7.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 6.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
4 calls. Slowest were:
| 29.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 27.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 24.0ms | (log (exp (- (sqrt (+ x 1)) (sqrt x)))) |
| 21.0ms | (* (sqrt x) (sqrt (+ x 1))) |
164 calls. Slowest were:
| 375.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 106.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)))))) |
| 106.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)))))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 20.6b
0% (0.7b remaining)
Error of 19.9b against oracle of 19.2b and baseline of 19.9b