


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Initial program 31.0
Simplified4.7
rmApplied add-exp-log4.7
Applied pow-exp4.7
Applied div-exp0.0
Final simplification0.0
herbie shell --seed 2019005 +o rules:numerics
(FPCore (x.re x.im y.re y.im)
:name "powComplex, real part"
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))
Time bar (total: 2.4m)Debug log
1 calls. Slowest were:
| 115.0ms | (* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 1.7b
Found 4 expressions with local error:
| 26.7b | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 0.2b | (* (atan2 x.im x.re) y.re) |
| 0.2b | (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) |
| 0.1b | (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)) |
4 calls. Slowest were:
| 6.0ms | (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) |
| 4.0ms | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 2.0ms | (* (atan2 x.im x.re) y.re) |
4 calls. Slowest were:
| 157.0ms | (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)) |
| 146.0ms | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 139.0ms | (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) |
| 13.0ms | (* (atan2 x.im x.re) y.re) |
108 calls. Slowest were:
| 235.0ms | (/ (exp (* y.im (atan2 x.im x.re))) (exp (* -1 (* y.re (log (/ -1 x.im)))))) |
| 151.0ms | (/ (* (* (pow (exp y.im) (atan2 x.im x.re)) (pow (exp y.im) (atan2 x.im x.re))) (pow (exp y.im) (atan2 x.im x.re))) (* (* (pow (hypot x.re x.im) y.re) (pow (hypot x.re x.im) y.re)) (pow (hypot x.re x.im) y.re))) |
| 135.0ms | (- 1 (+ (* y.im (* (atan2 x.im x.re) (* y.re (log x.im)))) (* 1/2 (* (pow y.im 2) (pow (log x.im) 2))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 26.7b | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 3.2b | (log (pow (hypot x.re x.im) y.re)) |
| 0.2b | (* (atan2 x.im x.re) y.re) |
| 0.2b | (* y.im (atan2 x.im x.re)) |
4 calls. Slowest were:
| 2.0ms | (* (atan2 x.im x.re) y.re) |
| 2.0ms | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 1.0ms | (log (pow (hypot x.re x.im) y.re)) |
4 calls. Slowest were:
| 144.0ms | (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 116.0ms | (log (pow (hypot x.re x.im) y.re)) |
| 23.0ms | (* (atan2 x.im x.re) y.re) |
| 14.0ms | (* y.im (atan2 x.im x.re)) |
27 calls. Slowest were:
| 166.0ms | (- 1 (+ (* y.im (* (atan2 x.im x.re) (* y.re (log x.im)))) (* 1/2 (* (pow y.im 2) (pow (log x.im) 2))))) |
| 61.0ms | (cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
| 35.0ms | (cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ 1 x.re))))) |
2 alts after pruning (1 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 29.4b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 29.4b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 26.7b | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 0.3b | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
4 calls. Slowest were:
| 18.0ms | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 11.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 7.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
4 calls. Slowest were:
| 175.0ms | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 85.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
| 31.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 25.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
27 calls. Slowest were:
| 181.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 166.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 124.0ms | (* -1 (* y.re (log (/ -1 x.re)))) |
2 alts after pruning (0 fresh and 2 done)
Merged error: 0.0b
0% (0.0b remaining)
Error of 0.0b against oracle of 0.0b and baseline of 0.0b