


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
if y.re < -1.7748707103617984e-55Initial program 37.3
Taylor expanded around -inf 0.2
Simplified0.2
rmApplied add-cbrt-cube0.2
if -1.7748707103617984e-55 < y.re < -4.860137094138367e-280 or 2.478594506845658e-159 < y.re Initial program 22.9
Taylor expanded around -inf 23.0
Simplified23.0
Taylor expanded around -inf 14.7
Simplified14.7
if -4.860137094138367e-280 < y.re < 2.478594506845658e-159Initial program 31.4
rmApplied add-cube-cbrt31.4
Final simplification10.8
herbie shell --seed 2019008
(FPCore (x.re x.im y.re y.im)
:name "powComplex, imaginary part"
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))
Time bar (total: 4.8m)Debug log
1 calls. Slowest were:
| 99.0ms | (* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 31.1b
Found 4 expressions with local error:
| 35.2b | (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 29.1b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 29.1b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 0.3b | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) |
4 calls. Slowest were:
| 18.0ms | (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 3.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) |
| 2.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
4 calls. Slowest were:
| 157.0ms | (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))) |
| 88.0ms | (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) |
| 22.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 18.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
27 calls. Slowest were:
| 143.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 108.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 79.0ms | (* -1 (* y.im (log (/ 1 x.re)))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 11.6b
Found 4 expressions with local error:
| 35.0b | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 29.1b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 0.3b | (* (log (- x.re)) y.im) |
| 0.3b | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 8.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 3.0ms | (* (log (- x.re)) y.im) |
| 2.0ms | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 115.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 72.0ms | (* (log (- x.re)) y.im) |
| 24.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
| 18.0ms | (* (atan2 x.im x.re) y.im) |
23 calls. Slowest were:
| 176.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 111.0ms | (sin (- (+ (* (atan2 x.im x.re) y.re) (* y.im (log -1))) (* y.im (log (/ 1 x.re))))) |
| 103.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 35.0b | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 35.0b | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 35.0b | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 29.1b | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
4 calls. Slowest were:
| 12.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 9.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 8.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
4 calls. Slowest were:
| 140.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 137.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 124.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 19.0ms | (sqrt (+ (* x.re x.re) (* x.im x.im))) |
21 calls. Slowest were:
| 134.0ms | (sqrt (- (* (* x.re x.re) (* x.re x.re)) (* (* x.im x.im) (* x.im x.im)))) |
| 102.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
| 88.0ms | (sin (- (+ (* (atan2 x.im x.re) y.re) (* y.im (log -1))) (* y.im (log (/ 1 x.re))))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 35.0b | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 0.3b | (* (log (- x.re)) y.im) |
| 0.3b | (* (log (- x.re)) y.re) |
| 0.3b | (* (atan2 x.im x.re) y.im) |
4 calls. Slowest were:
| 8.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 3.0ms | (* (log (- x.re)) y.im) |
| 3.0ms | (* (log (- x.re)) y.re) |
4 calls. Slowest were:
| 131.0ms | (sin (+ (* (log (- x.re)) y.im) (* (atan2 x.im x.re) y.re))) |
| 103.0ms | (* (log (- x.re)) y.re) |
| 76.0ms | (* (log (- x.re)) y.im) |
| 20.0ms | (* (atan2 x.im x.re) y.im) |
19 calls. Slowest were:
| 105.0ms | (* -1 (* y.im (log (/ -1 x.re)))) |
| 96.0ms | (* -1 (* y.re (log (/ -1 x.re)))) |
| 73.0ms | (sin (- (+ (* (atan2 x.im x.re) y.re) (* y.im (log -1))) (* y.im (log (/ 1 x.re))))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 6.0b
58% (5.8b remaining)
Error of 10.8b against oracle of 4.9b and baseline of 18.8b