


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
if y.im < 9.816734328325423e+110Initial program 22.7
rmApplied add-sqr-sqrt22.7
Applied associate-/r*22.6
rmApplied *-un-lft-identity22.6
Applied associate-/l*22.7
if 9.816734328325423e+110 < y.im Initial program 40.4
rmApplied add-sqr-sqrt40.4
Applied associate-/r*40.4
Taylor expanded around 0 39.6
Simplified39.6
Final simplification25.5
herbie shell --seed 2019016
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, imaginary part"
(/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))))
Time bar (total: 14.9s)Debug log
| 274× | (pre true 80) |
| 274× | (body real 80) |
| 45.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 26.5b
Found 3 expressions with local error:
| 26.5b | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 1.2b | (- (* x.im y.re) (* x.re y.im)) |
| 0.0b | (+ (* y.re y.re) (* y.im y.im)) |
| 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 |
| 7.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 4.0ms | (+ (* y.re y.re) (* y.im y.im)) |
| 3.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 92.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 52.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 24.0ms | (+ (* y.re y.re) (* y.im y.im)) |
| 261.0ms | (- (pow (* x.im y.re) 3) (pow (* x.re y.im) 3)) |
| 204.0ms | (- (* (* x.im y.re) (* x.im y.re)) (* (* x.re y.im) (* x.re y.im))) |
| 180.0ms | (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im)))) |
| 100.0ms | (/ (* (* (- (* x.im y.re) (* x.re y.im)) (- (* x.im y.re) (* x.re y.im))) (- (* x.im y.re) (* x.re y.im))) (* (* (+ (* y.re y.re) (* y.im y.im)) (+ (* y.re y.re) (* y.im y.im))) (+ (* y.re y.re) (* y.im y.im)))) |
| 71.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 26.0b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 16.5b | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 1.2b | (- (* x.im y.re) (* x.re y.im)) |
| 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 |
| 14.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 4.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 96.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 46.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 32.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 31.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 271.0ms | (- (pow (* x.im y.re) 3) (pow (* x.re y.im) 3)) |
| 252.0ms | (- (* (* x.im y.re) (* x.im y.re)) (* (* x.re y.im) (* x.re y.im))) |
| 219.0ms | (/ (* (* (- (* x.im y.re) (* x.re y.im)) (- (* x.im y.re) (* x.re y.im))) (- (* x.im y.re) (* x.re y.im))) (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im))))) |
| 143.0ms | (sqrt (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im))))) |
| 131.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 22.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 16.5b | (/ (- (* y.re x.im) (* y.im x.re)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 1.2b | (- (* y.re x.im) (* y.im x.re)) |
| 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 | (/ (- (* y.re x.im) (* y.im x.re)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 8.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 4.0ms | (- (* y.re x.im) (* y.im x.re)) |
| 4.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 73.0ms | (/ (- (* y.re x.im) (* y.im x.re)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 45.0ms | (- (* y.re x.im) (* y.im x.re)) |
| 29.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 20.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 285.0ms | (- (pow (* y.re x.im) 3) (pow (* y.im x.re) 3)) |
| 252.0ms | (- (* (* y.re x.im) (* y.re x.im)) (* (* y.im x.re) (* y.im x.re))) |
| 223.0ms | (/ (* (* (- (* y.re x.im) (* y.im x.re)) (- (* y.re x.im) (* y.im x.re))) (- (* y.re x.im) (* y.im x.re))) (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im))))) |
| 180.0ms | (sqrt (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im))))) |
| 141.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 22.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.2b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 16.6b | (/ (sqrt (+ (* y.re y.re) (* y.im y.im))) (- (* x.im y.re) (* x.re y.im))) |
| 1.2b | (- (* x.im y.re) (* x.re y.im)) |
| 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 |
| 12.0ms | (/ (sqrt (+ (* y.re y.re) (* y.im y.im))) (- (* x.im y.re) (* x.re y.im))) |
| 3.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 3.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 3.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 61.0ms | (/ (sqrt (+ (* y.re y.re) (* y.im y.im))) (- (* x.im y.re) (* x.re y.im))) |
| 39.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 23.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 19.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 149.0ms | (/ (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (* (* (- (* x.im y.re) (* x.re y.im)) (- (* x.im y.re) (* x.re y.im))) (- (* x.im y.re) (* x.re y.im)))) |
| 123.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
| 116.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
| 103.0ms | (+ (* (* x.im y.re) (* x.im y.re)) (+ (* (* x.re y.im) (* x.re y.im)) (* (* x.im y.re) (* x.re y.im)))) |
| 96.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
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) |