


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 28.7
rmApplied flip-+28.8
Applied associate-/l/28.8
Simplified0.4
rmApplied times-frac0.5
Simplified0.5
Simplified0.5
rmApplied pow10.5
Applied pow10.5
Applied pow-prod-down0.5
Simplified0.3
rmApplied div-inv0.4
Final simplification0.4
herbie shell --seed 2019007 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, narrow range"
:pre (and (< 1.0536712127723509e-08 a 94906265.62425156) (< 1.0536712127723509e-08 b 94906265.62425156) (< 1.0536712127723509e-08 c 94906265.62425156))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 43.5s)Debug log
1 calls. Slowest were:
| 87.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 28.6b
Found 4 expressions with local error:
| 28.5b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.2b | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 0.2b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 22.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 6.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 3.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 76.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 46.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 45.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 44.0ms | (- (* b b) (* (* 4 a) c)) |
40 calls. Slowest were:
| 385.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 326.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 172.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 0.4b
Found 4 expressions with local error:
| 0.4b | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 0.3b | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 0.2b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.2b | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 20.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 13.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 10.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
4 calls. Slowest were:
| 135.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 87.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 60.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 32.0ms | (- (* b b) (* (* 4 a) c)) |
143 calls. Slowest were:
| 1.2s | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 1.1s | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 409.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.3b | (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b))))) |
| 0.3b | (* (/ 2 a) (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b)))))) |
| 0.2b | (- (- b) (sqrt (fma -4 (* c a) (* b b)))) |
| 0.2b | (fma -4 (* c a) (* b b)) |
4 calls. Slowest were:
| 16.0ms | (* (/ 2 a) (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b)))))) |
| 12.0ms | (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b))))) |
| 8.0ms | (- (- b) (sqrt (fma -4 (* c a) (* b b)))) |
4 calls. Slowest were:
| 174.0ms | (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b))))) |
| 97.0ms | (* (/ 2 a) (/ (* c a) (- (- b) (sqrt (fma -4 (* c a) (* b b)))))) |
| 56.0ms | (- (- b) (sqrt (fma -4 (* c a) (* b b)))) |
| 47.0ms | (fma -4 (* c a) (* b b)) |
113 calls. Slowest were:
| 509.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma -4 (* c a) (* b b))) (sqrt (fma -4 (* c a) (* b b)))) (* (- b) (sqrt (fma -4 (* c a) (* b b)))))) |
| 410.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma -4 (* c a) (* b b))) (sqrt (fma -4 (* c a) (* b b)))) (* (- b) (sqrt (fma -4 (* c a) (* b b)))))) |
| 233.0ms | (+ (log (/ 2 a)) (- (log (* c a)) (log (- (- b) (sqrt (fma -4 (* c a) (* b b))))))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.2b | (- (- b) (sqrt (fma (* c a) -4 (* b b)))) |
| 0.2b | (/ (* 2 c) (- (- b) (sqrt (fma (* c a) -4 (* b b))))) |
| 0.2b | (fma (* c a) -4 (* b b)) |
| 0.1b | (sqrt (fma (* c a) -4 (* b b))) |
4 calls. Slowest were:
| 7.0ms | (/ (* 2 c) (- (- b) (sqrt (fma (* c a) -4 (* b b))))) |
| 5.0ms | (- (- b) (sqrt (fma (* c a) -4 (* b b)))) |
| 1.0ms | (sqrt (fma (* c a) -4 (* b b))) |
4 calls. Slowest were:
| 145.0ms | (/ (* 2 c) (- (- b) (sqrt (fma (* c a) -4 (* b b))))) |
| 57.0ms | (- (- b) (sqrt (fma (* c a) -4 (* b b)))) |
| 42.0ms | (sqrt (fma (* c a) -4 (* b b))) |
| 38.0ms | (fma (* c a) -4 (* b b)) |
94 calls. Slowest were:
| 390.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* c a) -4 (* b b))) (sqrt (fma (* c a) -4 (* b b)))) (* (- b) (sqrt (fma (* c a) -4 (* b b)))))) |
| 80.0ms | (fma (- (sqrt (cbrt (fma (* c a) -4 (* b b))))) (sqrt (* (cbrt (fma (* c a) -4 (* b b))) (cbrt (fma (* c a) -4 (* b b))))) (* (sqrt (cbrt (fma (* c a) -4 (* b b)))) (sqrt (* (cbrt (fma (* c a) -4 (* b b))) (cbrt (fma (* c a) -4 (* b b))))))) |
| 75.0ms | (fma 1 (- b) (- (* (sqrt (fma (* c a) -4 (* b b))) 1))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 0b
0% (0.4b remaining)
Error of 0.4b against oracle of 0.0b and baseline of 0.4b