



Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
Results
| Original | 25.7 |
|---|---|
| Target | 0.5 |
| Herbie | 25.5 |
if d < 9.816734328325423e+110Initial program 22.7
rmApplied add-sqr-sqrt22.7
Applied associate-/r*22.6
rmApplied clear-num22.7
if 9.816734328325423e+110 < d Initial program 40.4
rmApplied add-sqr-sqrt40.4
Applied associate-/r*40.4
Taylor expanded around -inf 40.4
Taylor expanded around 0 39.6
Simplified39.6
Final simplification25.5
herbie shell --seed 2019016
(FPCore (a b c d)
:name "Complex division, imag part"
:herbie-target
(if (< (fabs d) (fabs c)) (/ (- b (* a (/ d c))) (+ c (* d (/ d c)))) (/ (+ (- a) (* b (/ c d))) (+ d (* c (/ c d)))))
(/ (- (* b c) (* a d)) (+ (* c c) (* d d))))
Time bar (total: 15.6s)Debug log
| 274× | (pre true 80) |
| 274× | (body real 80) |
| 83.0ms | (/ (- (* b c) (* a d)) (+ (* c c) (* d d))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 26.5b
Found 3 expressions with local error:
| 26.5b | (/ (- (* b c) (* a d)) (+ (* c c) (* d d))) |
| 1.2b | (- (* b c) (* a d)) |
| 0.0b | (+ (* c c) (* d d)) |
| 11× | add-cube-cbrt |
| 11× | *-un-lft-identity |
| 11× | add-sqr-sqrt |
| 9× | times-frac |
| 7× | add-log-exp |
| 5× | add-exp-log |
| 5× | add-cbrt-cube |
| 3× | associate-/l* |
| 3× | associate-/r* |
| 3× | pow1 |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-/r/ |
| 2× | flip3-- |
| 2× | flip3-+ |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | sum-log |
| 1× | clear-num |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 13.0ms | (/ (- (* b c) (* a d)) (+ (* c c) (* d d))) |
| 3.0ms | (+ (* c c) (* d d)) |
| 2.0ms | (- (* b c) (* a d)) |
| 72.0ms | (/ (- (* b c) (* a d)) (+ (* c c) (* d d))) |
| 35.0ms | (- (* b c) (* a d)) |
| 15.0ms | (+ (* c c) (* d d)) |
| 340.0ms | (- (* (* b c) (* b c)) (* (* a d) (* a d))) |
| 313.0ms | (- (pow (* b c) 3) (pow (* a d) 3)) |
| 164.0ms | (+ (* (* c c) (* c c)) (- (* (* d d) (* d d)) (* (* c c) (* d d)))) |
| 103.0ms | (/ (* (* (- (* b c) (* a d)) (- (* b c) (* a d))) (- (* b c) (* a d))) (* (* (+ (* c c) (* d d)) (+ (* c c) (* d d))) (+ (* c c) (* d d)))) |
| 71.0ms | (/ (- (* b c) (* a d)) (cbrt (+ (* c c) (* d d)))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 26.0b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 16.5b | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 1.2b | (- (* b c) (* a d)) |
| 23× | add-sqr-sqrt |
| 21× | add-cube-cbrt |
| 21× | *-un-lft-identity |
| 18× | sqrt-prod |
| 18× | times-frac |
| 6× | add-log-exp |
| 6× | add-exp-log |
| 6× | sqrt-div |
| 6× | add-cbrt-cube |
| 6× | associate-/r* |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | flip-+ |
| 3× | flip3-+ |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | associate-/r/ |
| 2× | pow1/2 |
| 2× | flip3-- |
| 2× | rem-sqrt-square |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 8.0ms | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 2.0ms | (sqrt (+ (* c c) (* d d))) |
| 2.0ms | (- (* b c) (* a d)) |
| 2.0ms | (sqrt (+ (* c c) (* d d))) |
| 123.0ms | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 50.0ms | (- (* b c) (* a d)) |
| 24.0ms | (sqrt (+ (* c c) (* d d))) |
| 19.0ms | (sqrt (+ (* c c) (* d d))) |
| 326.0ms | (- (pow (* b c) 3) (pow (* a d) 3)) |
| 265.0ms | (- (* (* b c) (* b c)) (* (* a d) (* a d))) |
| 244.0ms | (/ (* (* (- (* b c) (* a d)) (- (* b c) (* a d))) (- (* b c) (* a d))) (* (* (sqrt (+ (* c c) (* d d))) (sqrt (+ (* c c) (* d d)))) (sqrt (+ (* c c) (* d d))))) |
| 182.0ms | (sqrt (+ (* (* c c) (* c c)) (- (* (* d d) (* d d)) (* (* c c) (* d d))))) |
| 156.0ms | (sqrt (- (* (* c c) (* c c)) (* (* d d) (* d d)))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 22.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 16.5b | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 1.2b | (- (* b c) (* a d)) |
| 23× | add-sqr-sqrt |
| 21× | add-cube-cbrt |
| 21× | *-un-lft-identity |
| 18× | sqrt-prod |
| 18× | times-frac |
| 6× | add-log-exp |
| 6× | add-exp-log |
| 6× | sqrt-div |
| 6× | add-cbrt-cube |
| 6× | associate-/r* |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | flip-+ |
| 3× | flip3-+ |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | associate-/r/ |
| 2× | pow1/2 |
| 2× | flip3-- |
| 2× | rem-sqrt-square |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 15.0ms | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 4.0ms | (sqrt (+ (* c c) (* d d))) |
| 4.0ms | (sqrt (+ (* c c) (* d d))) |
| 4.0ms | (- (* b c) (* a d)) |
| 101.0ms | (/ (- (* b c) (* a d)) (sqrt (+ (* c c) (* d d)))) |
| 31.0ms | (- (* b c) (* a d)) |
| 30.0ms | (sqrt (+ (* c c) (* d d))) |
| 30.0ms | (sqrt (+ (* c c) (* d d))) |
| 296.0ms | (- (pow (* b c) 3) (pow (* a d) 3)) |
| 285.0ms | (- (* (* b c) (* b c)) (* (* a d) (* a d))) |
| 234.0ms | (/ (* (* (- (* b c) (* a d)) (- (* b c) (* a d))) (- (* b c) (* a d))) (* (* (sqrt (+ (* c c) (* d d))) (sqrt (+ (* c c) (* d d)))) (sqrt (+ (* c c) (* d d))))) |
| 136.0ms | (sqrt (- (* (* c c) (* c c)) (* (* d d) (* d d)))) |
| 127.0ms | (sqrt (+ (* (* c c) (* c c)) (- (* (* d d) (* d d)) (* (* c c) (* d d))))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 22.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 28.2b | (sqrt (+ (* c c) (* d d))) |
| 16.6b | (/ (sqrt (+ (* c c) (* d d))) (- (* b c) (* a d))) |
| 1.2b | (- (* b c) (* a d)) |
| 23× | add-sqr-sqrt |
| 21× | add-cube-cbrt |
| 21× | *-un-lft-identity |
| 18× | sqrt-prod |
| 18× | times-frac |
| 6× | add-log-exp |
| 6× | associate-/l* |
| 6× | add-exp-log |
| 6× | sqrt-div |
| 6× | add-cbrt-cube |
| 4× | pow1 |
| 3× | flip-+ |
| 3× | flip3-+ |
| 3× | associate-/r* |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | associate-/r/ |
| 2× | pow1/2 |
| 2× | flip3-- |
| 2× | rem-sqrt-square |
| 1× | div-inv |
| 1× | div-exp |
| 1× | diff-log |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 15.0ms | (/ (sqrt (+ (* c c) (* d d))) (- (* b c) (* a d))) |
| 4.0ms | (sqrt (+ (* c c) (* d d))) |
| 4.0ms | (- (* b c) (* a d)) |
| 4.0ms | (sqrt (+ (* c c) (* d d))) |
| 84.0ms | (/ (sqrt (+ (* c c) (* d d))) (- (* b c) (* a d))) |
| 33.0ms | (- (* b c) (* a d)) |
| 31.0ms | (sqrt (+ (* c c) (* d d))) |
| 30.0ms | (sqrt (+ (* c c) (* d d))) |
| 201.0ms | (/ (* (* (sqrt (+ (* c c) (* d d))) (sqrt (+ (* c c) (* d d)))) (sqrt (+ (* c c) (* d d)))) (* (* (- (* b c) (* a d)) (- (* b c) (* a d))) (- (* b c) (* a d)))) |
| 174.0ms | (sqrt (- (* (* c c) (* c c)) (* (* d d) (* d d)))) |
| 147.0ms | (+ (* (* b c) (* b c)) (+ (* (* a d) (* a d)) (* (* b c) (* a d)))) |
| 141.0ms | (sqrt (- (* (* c c) (* c c)) (* (* d d) (* d d)))) |
| 114.0ms | (sqrt (- (* (* c c) (* c c)) (* (* d d) (* d d)))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 22.9b
5.5% (2.4b remaining)
Error of 25.5b against oracle of 23.1b and baseline of 25.7b
| 8500× | (pre true 80) |
| 8500× | (body real 80) |