



Bits error versus re




Bits error versus im
Results
| Original | 43.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.8 |
Initial program 43.6
Taylor expanded around 0 0.8
rmApplied distribute-neg-in0.8
Applied distribute-lft-in0.8
Simplified0.8
Final simplification0.8
herbie shell --seed 2019004
(FPCore (re im)
:name "math.cos on complex, imaginary part"
:herbie-target
(if (< (fabs im) 1) (- (* (sin re) (+ (+ im (* (* (* 1/6 im) im) im)) (* (* (* (* (* 1/120 im) im) im) im) im)))) (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))
(* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))
Time bar (total: 35.0s)Debug log
1 calls. Slowest were:
| 532.0ms | (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 44.0b
Found 2 expressions with local error:
| 6.9b | (- (exp (- im)) (exp im)) |
| 0.0b | (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))) |
2 calls. Slowest were:
| 27.0ms | (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))) |
| 13.0ms | (- (exp (- im)) (exp im)) |
2 calls. Slowest were:
| 120.0ms | (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))) |
| 62.0ms | (- (exp (- im)) (exp im)) |
24 calls. Slowest were:
| 157.0ms | (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))) |
| 154.0ms | (* 0.5 (* (sin re) (- (exp (- im)) (exp im)))) |
| 137.0ms | (* 0.5 (* (sin re) (- (exp (* -1 im)) (exp im)))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 0.2b | (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))) |
| 0.2b | (* 1/3 (pow im 3)) |
| 0.0b | (* 1/60 (pow im 5)) |
| 0.0b | (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))) |
4 calls. Slowest were:
| 19.0ms | (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))) |
| 9.0ms | (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))) |
| 1.0ms | (* 1/3 (pow im 3)) |
4 calls. Slowest were:
| 177.0ms | (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))) |
| 24.0ms | (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))) |
| 20.0ms | (* 1/60 (pow im 5)) |
| 15.0ms | (* 1/3 (pow im 3)) |
39 calls. Slowest were:
| 301.0ms | (* (* 0.5 (sin re)) (- (+ (* 1/60 (pow im 5)) (* 2 im)))) |
| 282.0ms | (* (* (* (* 0.5 (sin re)) (* 0.5 (sin re))) (* 0.5 (sin re))) (* (* (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 233.0ms | (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 0.2b | (* (+ (* im -2) (* -1/60 (pow im 5))) (* 0.5 (sin re))) |
| 0.2b | (* 1/3 (pow im 3)) |
| 0.1b | (* (* 0.5 (sin re)) (- (* 1/3 (pow im 3)))) |
| 0.0b | (* -1/60 (pow im 5)) |
4 calls. Slowest were:
| 13.0ms | (* (+ (* im -2) (* -1/60 (pow im 5))) (* 0.5 (sin re))) |
| 12.0ms | (* (* 0.5 (sin re)) (- (* 1/3 (pow im 3)))) |
| 2.0ms | (* -1/60 (pow im 5)) |
4 calls. Slowest were:
| 91.0ms | (* (+ (* im -2) (* -1/60 (pow im 5))) (* 0.5 (sin re))) |
| 83.0ms | (* (* 0.5 (sin re)) (- (* 1/3 (pow im 3)))) |
| 28.0ms | (* 1/3 (pow im 3)) |
| 26.0ms | (* -1/60 (pow im 5)) |
36 calls. Slowest were:
| 381.0ms | (+ (log (* 0.5 (sin re))) (log (- (* 1/3 (pow im 3))))) |
| 191.0ms | (* (* 0.5 (sin re)) (- (* 1/3 (pow im 3)))) |
| 153.0ms | (* (+ (pow (* im -2) 3) (pow (* -1/60 (pow im 5)) 3)) (* 0.5 (sin re))) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 14.1b | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 14.1b | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 14.1b | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 0.5b | (* (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))))) |
4 calls. Slowest were:
| 25.0ms | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 21.0ms | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 20.0ms | (* (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))))) |
4 calls. Slowest were:
| 366.0ms | (* (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im))))))) |
| 357.0ms | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 326.0ms | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
| 317.0ms | (cbrt (* (* 0.5 (sin re)) (- (+ (* 1/3 (pow im 3)) (+ (* 1/60 (pow im 5)) (* 2 im)))))) |
49 calls. Slowest were:
| 409.0ms | (+ (* (cbrt 0.5) (exp (* 1/3 (- (log (* 1/60 (sin re))) (* 5 (log (/ -1 im))))))) (* 20/3 (/ (* (cbrt 0.5) (exp (* 1/3 (- (log (* 1/60 (sin re))) (* 5 (log (/ -1 im))))))) (pow im 2)))) |
| 380.0ms | (+ (* 20/3 (/ (* (cbrt -0.5) (exp (* 1/3 (- (log (* 1/60 (sin re))) (* 5 (log (/ 1 im))))))) (pow im 2))) (* (cbrt -0.5) (exp (* 1/3 (- (log (* 1/60 (sin re))) (* 5 (log (/ 1 im)))))))) |
| 380.0ms | (- (+ (* (cbrt -0.5) (exp (* 1/3 (+ (log 2) (+ (log im) (log re)))))) (* 1/18 (* (pow im 2) (* (cbrt -0.5) (exp (* 1/3 (+ (log 2) (+ (log im) (log re))))))))) (* 1/18 (* (pow re 2) (* (cbrt -0.5) (exp (* 1/3 (+ (log 2) (+ (log im) (log re))))))))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 0.0b
0% (0.7b remaining)
Error of 0.8b against oracle of 0.1b and baseline of 0.8b