


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
if y.re < -1.8506614466369815e+150Initial program 44.1
Simplified44.1
rmApplied add-sqr-sqrt44.1
Applied *-un-lft-identity44.1
Applied times-frac44.1
Simplified44.1
Simplified28.2
rmApplied associate-*l/28.1
Simplified28.1
Taylor expanded around -inf 13.3
Simplified13.3
if -1.8506614466369815e+150 < y.re < 6.560820932670249e+185Initial program 20.1
Simplified20.1
rmApplied add-sqr-sqrt20.1
Applied *-un-lft-identity20.1
Applied times-frac20.2
Simplified20.1
Simplified12.4
rmApplied associate-*l/12.2
Simplified12.2
rmApplied clear-num12.3
if 6.560820932670249e+185 < y.re Initial program 44.6
Simplified44.6
rmApplied add-sqr-sqrt44.6
Applied *-un-lft-identity44.6
Applied times-frac44.6
Simplified44.6
Simplified30.7
rmApplied associate-*l/30.7
Simplified30.7
rmApplied clear-num30.7
Taylor expanded around 0 12.2
Final simplification12.5
herbie shell --seed 2019022 +o rules:numerics
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, real part"
(/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))))
Time bar (total: 18.6s)Debug log
| 42.0ms | 270× | body | 80 | valid |
| 13.0ms | (/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 26.8b
Found 3 expressions with local error:
| 26.8b | (/ (fma x.re y.re (* x.im y.im)) (fma y.im y.im (* y.re y.re))) |
| 0.3b | (fma x.re y.re (* x.im y.im)) |
| 0.0b | (fma y.im y.im (* y.re y.re)) |
| 11× | add-cube-cbrt |
| 11× | *-un-lft-identity |
| 11× | add-sqr-sqrt |
| 9× | times-frac |
| 5× | add-exp-log |
| 5× | add-cbrt-cube |
| 3× | add-log-exp |
| 3× | associate-/l* |
| 3× | log1p-expm1-u |
| 3× | associate-/r* |
| 3× | pow1 |
| 3× | expm1-log1p-u |
| 2× | fma-udef |
| 1× | div-inv |
| 1× | div-exp |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 6.0ms | (/ (fma x.re y.re (* x.im y.im)) (fma y.im y.im (* y.re y.re))) |
| 0.0ms | (fma x.re y.re (* x.im y.im)) |
| 0.0ms | (fma y.im y.im (* y.re y.re)) |
| 88.0ms | (/ (fma x.re y.re (* x.im y.im)) (fma y.im y.im (* y.re y.re))) |
| 26.0ms | (fma y.im y.im (* y.re y.re)) |
| 26.0ms | (fma x.re y.re (* x.im y.im)) |
| 206.0ms | (/ (* (* (fma x.re y.re (* x.im y.im)) (fma x.re y.re (* x.im y.im))) (fma x.re y.re (* x.im y.im))) (* (* (fma y.im y.im (* y.re y.re)) (fma y.im y.im (* y.re y.re))) (fma y.im y.im (* y.re y.re)))) |
| 52.0ms | (/ (* (cbrt (fma x.re y.re (* x.im y.im))) (cbrt (fma x.re y.re (* x.im y.im)))) (* (cbrt (fma y.im y.im (* y.re y.re))) (cbrt (fma y.im y.im (* y.re y.re))))) |
| 20.0ms | (/ (* (cbrt (fma x.re y.re (* x.im y.im))) (cbrt (fma x.re y.re (* x.im y.im)))) (sqrt (fma y.im y.im (* y.re y.re)))) |
| 19.0ms | (/ (sqrt (fma x.re y.re (* x.im y.im))) (* (cbrt (fma y.im y.im (* y.re y.re))) (cbrt (fma y.im y.im (* y.re y.re))))) |
| 15.0ms | (/ (sqrt (fma x.re y.re (* x.im y.im))) (sqrt (fma y.im y.im (* y.re y.re)))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 16.4b
Found 4 expressions with local error:
| 16.2b | (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)) |
| 0.5b | (fma y.im x.im (* x.re y.re)) |
| 0.2b | (* (/ 1 (hypot y.im y.re)) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re))) |
| 0.0b | (hypot y.im y.re) |
| 20× | add-cube-cbrt |
| 20× | *-un-lft-identity |
| 20× | add-sqr-sqrt |
| 18× | times-frac |
| 16× | add-exp-log |
| 13× | associate-*r* |
| 11× | add-cbrt-cube |
| 6× | pow1 |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | prod-exp |
| 4× | associate-*l* |
| 4× | expm1-log1p-u |
| 3× | associate-/l* |
| 3× | div-inv |
| 3× | div-exp |
| 3× | associate-/r* |
| 2× | cbrt-unprod |
| 2× | rec-exp |
| 2× | cbrt-undiv |
| 1× | fma-udef |
| 1× | *-commutative |
| 1× | associate-*r/ |
| 1× | associate-*l/ |
| 1× | pow-prod-down |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | hypot-udef |
| 1× | frac-times |
| 10.0ms | (* (/ 1 (hypot y.im y.re)) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re))) |
| 3.0ms | (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)) |
| 0.0ms | (hypot y.im y.re) |
| 0.0ms | (fma y.im x.im (* x.re y.re)) |
| 134.0ms | (* (/ 1 (hypot y.im y.re)) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re))) |
| 63.0ms | (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)) |
| 36.0ms | (hypot y.im y.re) |
| 35.0ms | (fma y.im x.im (* x.re y.re)) |
| 579.0ms | (* (* (* (/ 1 (hypot y.im y.re)) (/ 1 (hypot y.im y.re))) (/ 1 (hypot y.im y.re))) (* (* (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re))) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)))) |
| 210.0ms | (* (/ 1 (hypot y.im y.re)) (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re))) |
| 162.0ms | (/ (* (* (fma y.im x.im (* x.re y.re)) (fma y.im x.im (* x.re y.re))) (fma y.im x.im (* x.re y.re))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re))) |
| 157.0ms | (* (* (* (/ 1 (hypot y.im y.re)) (/ 1 (hypot y.im y.re))) (/ 1 (hypot y.im y.re))) (/ (* (* (fma y.im x.im (* x.re y.re)) (fma y.im x.im (* x.re y.re))) (fma y.im x.im (* x.re y.re))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re)))) |
| 107.0ms | (+ (log (/ 1 (hypot y.im y.re))) (log (/ (fma y.im x.im (* x.re y.re)) (hypot y.im y.re)))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 7.6b
Found 4 expressions with local error:
| 16.2b | (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) |
| 0.3b | (fma y.re x.re (* x.im y.im)) |
| 0.1b | (/ (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) (hypot y.im y.re)) |
| 0.0b | (hypot y.im y.re) |
| 84× | times-frac |
| 54× | add-cube-cbrt |
| 54× | *-un-lft-identity |
| 54× | add-sqr-sqrt |
| 16× | associate-/l* |
| 11× | add-exp-log |
| 11× | add-cbrt-cube |
| 6× | div-inv |
| 6× | associate-/r* |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | div-exp |
| 4× | cbrt-undiv |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | associate-/l/ |
| 1× | fma-udef |
| 1× | hypot-udef |
| 6.0ms | (/ (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) (hypot y.im y.re)) |
| 4.0ms | (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) |
| 0.0ms | (fma y.re x.re (* x.im y.im)) |
| 0.0ms | (hypot y.im y.re) |
| 100.0ms | (/ (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) (hypot y.im y.re)) |
| 66.0ms | (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) |
| 34.0ms | (hypot y.im y.re) |
| 28.0ms | (fma y.re x.re (* x.im y.im)) |
| 362.0ms | (/ (sqrt (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re))) (sqrt (hypot y.im y.re))) |
| 292.0ms | (/ (/ (* (* (fma y.re x.re (* x.im y.im)) (fma y.re x.re (* x.im y.im))) (fma y.re x.re (* x.im y.im))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re))) |
| 209.0ms | (/ (* (* (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re)) (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re))) (/ (fma y.re x.re (* x.im y.im)) (hypot y.im y.re))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re))) |
| 177.0ms | (/ (/ (* (cbrt (fma y.re x.re (* x.im y.im))) (cbrt (fma y.re x.re (* x.im y.im)))) (* (cbrt (hypot y.im y.re)) (cbrt (hypot y.im y.re)))) 1) |
| 139.0ms | (/ (* (* (fma y.re x.re (* x.im y.im)) (fma y.re x.re (* x.im y.im))) (fma y.re x.re (* x.im y.im))) (* (* (hypot y.im y.re) (hypot y.im y.re)) (hypot y.im y.re))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 2.5b
Found 4 expressions with local error:
| 16.3b | (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im))) |
| 0.3b | (fma y.re x.re (* x.im y.im)) |
| 0.2b | (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) |
| 0.1b | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (hypot y.im y.re)) |
| 33× | times-frac |
| 29× | add-cube-cbrt |
| 29× | *-un-lft-identity |
| 29× | add-sqr-sqrt |
| 19× | associate-/r* |
| 16× | add-exp-log |
| 8× | associate-/l* |
| 8× | div-inv |
| 8× | add-cbrt-cube |
| 6× | div-exp |
| 5× | associate-/r/ |
| 5× | pow1 |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | rec-exp |
| 4× | expm1-log1p-u |
| 3× | frac-2neg |
| 3× | clear-num |
| 2× | cbrt-undiv |
| 1× | associate-/l/ |
| 1× | inv-pow |
| 1× | pow-flip |
| 1× | fma-udef |
| 9.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (hypot y.im y.re)) |
| 7.0ms | (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im))) |
| 3.0ms | (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) |
| 0.0ms | (fma y.re x.re (* x.im y.im)) |
| 154.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (hypot y.im y.re)) |
| 125.0ms | (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) |
| 81.0ms | (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im))) |
| 40.0ms | (fma y.re x.re (* x.im y.im)) |
| 406.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (cbrt (hypot y.im y.re))) |
| 375.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (sqrt (hypot y.im y.re))) |
| 340.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (cbrt (hypot y.im y.re))) |
| 312.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (sqrt (hypot y.im y.re))) |
| 271.0ms | (/ (/ 1 (/ (hypot y.im y.re) (fma y.re x.re (* x.im y.im)))) (hypot y.im y.re)) |
10 alts after pruning (7 fresh and 3 done)
Merged error: 2.5b
28.9% (10.0b remaining)
Error of 12.5b against oracle of 2.5b and baseline of 16.5b
| 1.1s | 8447× | body | 80 | valid |