


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Initial program 31.8
Simplified6.8
rmApplied pow-to-exp6.8
Applied pow-exp6.8
Applied div-exp0.1
Final simplification0.1
herbie shell --seed 2019005 +o rules:numerics
(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: 3.9m)Debug log
1 calls. Slowest were:
| 103.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)))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 3.8b
Found 4 expressions with local error:
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 0.3b | (* (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:
| 12.0ms | (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) |
| 6.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 5.0ms | (* (atan2 x.im x.re) y.re) |
4 calls. Slowest were:
| 171.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 160.0ms | (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)) |
| 116.0ms | (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) |
| 17.0ms | (* (atan2 x.im x.re) y.re) |
108 calls. Slowest were:
| 224.0ms | (/ (exp (* y.im (atan2 x.im x.re))) (exp (* -1 (* y.re (log (/ -1 x.im)))))) |
| 142.0ms | (/ (exp (* y.im (atan2 x.im x.re))) (exp (* -1 (* y.re (log (/ 1 x.im)))))) |
| 134.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))) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 0.3b | (* (log (hypot x.re x.im)) y.re) |
| 0.3b | (* (atan2 x.im x.re) y.re) |
| 0.2b | (* y.im (atan2 x.im x.re)) |
4 calls. Slowest were:
| 3.0ms | (* (log (hypot x.re x.im)) y.re) |
| 2.0ms | (* (atan2 x.im x.re) y.re) |
| 2.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
4 calls. Slowest were:
| 119.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 98.0ms | (* (log (hypot x.re x.im)) y.re) |
| 18.0ms | (* (atan2 x.im x.re) y.re) |
| 12.0ms | (* y.im (atan2 x.im x.re)) |
20 calls. Slowest were:
| 121.0ms | (* -1 (* y.re (log (/ -1 x.re)))) |
| 104.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
| 86.0ms | (* -1 (* y.re (log (/ 1 x.re)))) |
5 alts after pruning (4 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 0.7b | (cbrt (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) |
4 calls. Slowest were:
| 4.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 4.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 4.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
4 calls. Slowest were:
| 202.0ms | (cbrt (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) |
| 131.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 123.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 101.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
18 calls. Slowest were:
| 352.0ms | (- (+ (* 1/3 (* y.im (* (exp (* 1/3 (- (log (/ 1 (pow (atan2 x.im x.re) 2))) (* 2 (log y.re))))) (log x.im)))) (exp (* 1/3 (+ (log (atan2 x.im x.re)) (log y.re))))) (* 7/81 (* (pow y.im 2) (* (pow (log x.im) 2) (exp (* 1/3 (+ (log (atan2 x.im x.re)) (log y.re)))))))) |
| 123.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
| 97.0ms | (pow (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) 1/3) |
6 alts after pruning (4 fresh and 2 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 37.9b | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 10.6b | (cbrt (* (* (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) |
4 calls. Slowest were:
| 26.0ms | (cbrt (* (* (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) |
| 4.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 4.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
4 calls. Slowest were:
| 148.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 132.0ms | (cbrt (* (* (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) |
| 123.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
| 101.0ms | (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) |
18 calls. Slowest were:
| 145.0ms | (cbrt (* (- (cos (- (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)) (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (cos (+ (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)) (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) (sin (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) |
| 122.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
| 94.0ms | (sin (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 0.0b
0% (0.0b remaining)
Error of 0.1b against oracle of 0.1b and baseline of 0.1b