\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)\cos re \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot -0.16666666666666666\right) - im\right)
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))
(FPCore (re im) :precision binary64 (* (cos re) (- (* (* im im) (* im -0.16666666666666666)) im)))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp(0.0 - im) - exp(im));
}
double code(double re, double im) {
return cos(re) * (((im * im) * (im * -0.16666666666666666)) - im);
}







Bits error versus re







Bits error versus im
Results
| Original | 58.0 |
|---|---|
| Target | 0.3 |
| Herbie | 0.9 |
| Alternative 1 | |
|---|---|
| Error | 1.3 |
| Cost | 6656 |
| Alternative 2 | |
|---|---|
| Error | 60.6 |
| Cost | 64 |
| Alternative 3 | |
|---|---|
| Error | 61.5 |
| Cost | 64 |

Initial program 58.0
Simplified58.0
Taylor expanded around 0 0.9
Simplified0.9
rmApplied unpow3_binary64_18490.9
Applied associate-*l*_binary64_17240.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2021044
(FPCore (re im)
:name "math.sin on complex, imaginary part"
:precision binary64
:herbie-target
(if (< (fabs im) 1.0) (- (* (cos re) (+ (+ im (* (* (* 0.16666666666666666 im) im) im)) (* (* (* (* (* 0.008333333333333333 im) im) im) im) im)))) (* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))
(* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))