


Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 28.7
Initial simplification28.7
rmApplied flip--28.8
Applied associate-/l/28.8
Simplified0.5
rmApplied associate-/l*0.3
Final simplification0.3
herbie shell --seed 2018360
(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: 19.5s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 29.1b
Found 4 expressions with local error:
| 28.9b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.2b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 0.2b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 24.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 10.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 66.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 51.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 46.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 25.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 272.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 182.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 181.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.4b
Found 4 expressions with local error:
| 0.4b | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 0.2b | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 0.2b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.2b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 29.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 23.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 15.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
4 calls. Slowest were:
| 160.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 77.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 51.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 25.0ms | (- (* b b) (* (* c a) 4)) |
33 calls. Slowest were:
| 908.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 815.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (- (* b b) (* (* c a) 4))) b) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 264.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (- (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.2b | (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a) |
| 0.2b | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 0.2b | (/ (* c -4) (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a)) |
| 0.2b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
4 calls. Slowest were:
| 16.0ms | (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a) |
| 14.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 8.0ms | (/ (* c -4) (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a)) |
4 calls. Slowest were:
| 125.0ms | (/ (* c -4) (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a)) |
| 92.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 90.0ms | (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a) |
| 45.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
39 calls. Slowest were:
| 417.0ms | (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a) |
| 407.0ms | (/ -4 (sqrt (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a))) |
| 400.0ms | (sqrt (/ (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) a)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.4b | (/ (* (* c -4) a) (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b))) |
| 0.2b | (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)) |
| 0.2b | (- (* b b) (* (* c a) 4)) |
| 0.2b | (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) |
4 calls. Slowest were:
| 11.0ms | (/ (* (* c -4) a) (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b))) |
| 9.0ms | (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) |
| 7.0ms | (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)) |
4 calls. Slowest were:
| 240.0ms | (/ (* (* c -4) a) (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b))) |
| 84.0ms | (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) |
| 84.0ms | (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)) |
| 40.0ms | (- (* b b) (* (* c a) 4)) |
44 calls. Slowest were:
| 427.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)) (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b))) (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)))) |
| 270.0ms | (sqrt (+ (* (* b b) (* b b)) (+ (* (* (* c a) 4) (* (* c a) 4)) (* (* b b) (* (* c a) 4))))) |
| 227.0ms | (- (log (* (* c -4) a)) (log (+ (* (* 2 a) (sqrt (- (* b b) (* (* c a) 4)))) (* (* 2 a) b)))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 0b
0% (0.3b remaining)
Error of 0.3b against oracle of 0.0b and baseline of 0.3b