



Bits error versus x
Results
| Original | 19.8 |
|---|---|
| Target | 0.6 |
| Herbie | 0.4 |
Initial program 19.8
rmApplied frac-sub19.8
Simplified19.8
rmApplied flip--19.6
Applied associate-/l/19.6
Simplified0.7
rmApplied associate-/r*0.4
rmApplied div-inv0.4
Final simplification0.4
herbie shell --seed 2019002 +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: 18.5s)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.3b
Found 3 expressions with local error:
| 3.3b | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 0.2b | (/ 1 (sqrt x)) |
| 0.1b | (/ 1 (sqrt (+ x 1))) |
3 calls. Slowest were:
| 14.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 3.0ms | (/ 1 (sqrt (+ x 1))) |
| 2.0ms | (/ 1 (sqrt x)) |
3 calls. Slowest were:
| 42.0ms | (- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))) |
| 26.0ms | (/ 1 (sqrt x)) |
| 12.0ms | (/ 1 (sqrt (+ x 1))) |
83 calls. Slowest were:
| 212.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)))))) |
| 203.0ms | (fma (- (sqrt (- x 1))) (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1)))))) |
| 182.0ms | (fma (- (sqrt (- x 1))) (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (sqrt (- x 1)) (/ 1 (sqrt (- (* x x) (* 1 1)))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 20.7b
Found 3 expressions with local error:
| 3.3b | (- (sqrt (+ x 1)) (sqrt x)) |
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
| 0.1b | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
3 calls. Slowest were:
| 36.0ms | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
| 13.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 9.0ms | (* (sqrt x) (sqrt (+ x 1))) |
3 calls. Slowest were:
| 93.0ms | (/ (- (sqrt (+ x 1)) (sqrt x)) (* (sqrt x) (sqrt (+ x 1)))) |
| 37.0ms | (- (sqrt (+ x 1)) (sqrt x)) |
| 35.0ms | (* (sqrt x) (sqrt (+ x 1))) |
97 calls. Slowest were:
| 604.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 105.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))))) |
| 103.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))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.7b
Found 4 expressions with local error:
| 0.6b | (/ 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:
| 19.0ms | (/ 1 (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x)))) |
| 18.0ms | (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 5.0ms | (+ (sqrt (+ x 1)) (sqrt x)) |
4 calls. Slowest were:
| 77.0ms | (/ 1 (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x)))) |
| 64.0ms | (* (* (sqrt x) (sqrt (+ x 1))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 30.0ms | (+ (sqrt (+ x 1)) (sqrt x)) |
| 29.0ms | (* (sqrt x) (sqrt (+ x 1))) |
71 calls. Slowest were:
| 508.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 224.0ms | (* (* (sqrt x) (sqrt (- (* x x) (* 1 1)))) (+ (pow (sqrt (+ x 1)) 3) (pow (sqrt x) 3))) |
| 206.0ms | (* (* (sqrt x) (sqrt (- (* x x) (* 1 1)))) (- (* (sqrt (+ x 1)) (sqrt (+ x 1))) (* (sqrt x) (sqrt x)))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.3b | (* (sqrt x) (sqrt (+ x 1))) |
| 0.1b | (/ (/ 1 (* (sqrt x) (sqrt (+ x 1)))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 0.1b | (/ 1 (* (sqrt x) (sqrt (+ x 1)))) |
| 0.0b | (+ (sqrt (+ x 1)) (sqrt x)) |
4 calls. Slowest were:
| 13.0ms | (/ (/ 1 (* (sqrt x) (sqrt (+ x 1)))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 8.0ms | (+ (sqrt (+ x 1)) (sqrt x)) |
| 4.0ms | (* (sqrt x) (sqrt (+ x 1))) |
4 calls. Slowest were:
| 90.0ms | (/ (/ 1 (* (sqrt x) (sqrt (+ x 1)))) (+ (sqrt (+ x 1)) (sqrt x))) |
| 37.0ms | (* (sqrt x) (sqrt (+ x 1))) |
| 24.0ms | (/ 1 (* (sqrt x) (sqrt (+ x 1)))) |
| 20.0ms | (+ (sqrt (+ x 1)) (sqrt x)) |
102 calls. Slowest were:
| 512.0ms | (* (* (* (sqrt x) (sqrt x)) (sqrt x)) (* (* (sqrt (+ x 1)) (sqrt (+ x 1))) (sqrt (+ x 1)))) |
| 157.0ms | (- (+ (* 3/8 (/ 1 (pow x 3))) (/ 1 x)) (* 1/2 (/ 1 (pow x 2)))) |
| 107.0ms | (- (log (/ 1 (* (sqrt x) (sqrt (+ x 1))))) (log (+ (sqrt (+ x 1)) (sqrt x)))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0b
0% (0.4b remaining)
Error of 0.4b against oracle of 0.0b and baseline of 0.4b