


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
if y.re < -2.0145832239889755e-11Initial program 37.1
Taylor expanded around 0 0.1
if -2.0145832239889755e-11 < y.re Initial program 27.6
Taylor expanded around 0 24.9
Taylor expanded around -inf 0.2
Simplified0.2
Final simplification0.1
herbie shell --seed 2018362
(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: 15.4s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 30.7b
Found 4 expressions with local error:
| 29.3b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 29.3b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 24.9b | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 0.3b | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 9.0ms | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 3.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 2.0ms | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 164.0ms | (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 28.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 20.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 14.0ms | (* (atan2 x.im x.re) y.im) |
27 calls. Slowest were:
| 151.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 142.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 67.0ms | (cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 17.7b
Found 4 expressions with local error:
| 29.3b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 0.3b | (* (atan2 x.im x.re) y.im) |
| 0.2b | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
| 0.0b | (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im)) |
4 calls. Slowest were:
| 9.0ms | (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im)) |
| 5.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
| 2.0ms | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 105.0ms | (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im)) |
| 80.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) |
| 26.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 15.0ms | (* (atan2 x.im x.re) y.im) |
23 calls. Slowest were:
| 183.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 116.0ms | (- (+ (* y.re (log (/ -1 x.re))) (* y.im (atan2 x.im x.re)))) |
| 89.0ms | (* -1 (* y.re (log (/ 1 x.re)))) |
3 alts after pruning (2 fresh and 1 done)
Merged error: 0.1b
Found 3 expressions with local error:
| 0.3b | (* (atan2 x.im x.re) y.im) |
| 0.3b | (* (log (- x.re)) y.re) |
| 0.0b | (- (* (log (- x.re)) y.re) (* (atan2 x.im x.re) y.im)) |
3 calls. Slowest were:
| 10.0ms | (- (* (log (- x.re)) y.re) (* (atan2 x.im x.re) y.im)) |
| 6.0ms | (* (log (- x.re)) y.re) |
| 5.0ms | (* (atan2 x.im x.re) y.im) |
3 calls. Slowest were:
| 151.0ms | (- (* (log (- x.re)) y.re) (* (atan2 x.im x.re) y.im)) |
| 78.0ms | (* (log (- x.re)) y.re) |
| 14.0ms | (* (atan2 x.im x.re) y.im) |
14 calls. Slowest were:
| 89.0ms | (- (+ (* y.re (log (/ -1 x.re))) (* y.im (atan2 x.im x.re)))) |
| 79.0ms | (- (* y.re (log -1)) (+ (* y.im (atan2 x.im x.re)) (* y.re (log (/ 1 x.re))))) |
| 69.0ms | (- (+ (* y.re (log -1)) (* y.re (log x.re))) (* y.im (atan2 x.im x.re))) |
3 alts after pruning (1 fresh and 2 done)
Merged error: 0.1b
Found 3 expressions with local error:
| 0.3b | (* (atan2 x.im x.re) y.im) |
| 0.1b | (* y.re (log x.im)) |
| 0.0b | (- (* y.re (log x.im)) (* (atan2 x.im x.re) y.im)) |
3 calls. Slowest were:
| 6.0ms | (- (* y.re (log x.im)) (* (atan2 x.im x.re) y.im)) |
| 5.0ms | (* (atan2 x.im x.re) y.im) |
| 2.0ms | (* y.re (log x.im)) |
3 calls. Slowest were:
| 99.0ms | (- (* y.re (log x.im)) (* (atan2 x.im x.re) y.im)) |
| 41.0ms | (* y.re (log x.im)) |
| 22.0ms | (* (atan2 x.im x.re) y.im) |
15 calls. Slowest were:
| 189.0ms | (- (* y.re (log -1)) (+ (* y.re (log (/ -1 x.im))) (* y.im (atan2 x.im x.re)))) |
| 106.0ms | (* y.re (- (log -1) (log (/ -1 x.im)))) |
| 62.0ms | (* -1 (* y.re (log (/ 1 x.im)))) |
3 alts after pruning (0 fresh and 3 done)
Merged error: 0.1b
98.2% (0.1b remaining)
Error of 0.1b against oracle of 0.0b and baseline of 7.1b