\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)\left({im}^{5} \cdot \frac{-1}{60} - \left(2 - \frac{-1}{3} \cdot \left(im \cdot im\right)\right) \cdot im\right) \cdot \left(0.5 \cdot \sin re\right)double f(double re, double im) {
double r11231956 = 0.5;
double r11231957 = re;
double r11231958 = sin(r11231957);
double r11231959 = r11231956 * r11231958;
double r11231960 = im;
double r11231961 = -r11231960;
double r11231962 = exp(r11231961);
double r11231963 = exp(r11231960);
double r11231964 = r11231962 - r11231963;
double r11231965 = r11231959 * r11231964;
return r11231965;
}
double f(double re, double im) {
double r11231966 = im;
double r11231967 = 5.0;
double r11231968 = pow(r11231966, r11231967);
double r11231969 = -0.016666666666666666;
double r11231970 = r11231968 * r11231969;
double r11231971 = 2.0;
double r11231972 = -0.3333333333333333;
double r11231973 = r11231966 * r11231966;
double r11231974 = r11231972 * r11231973;
double r11231975 = r11231971 - r11231974;
double r11231976 = r11231975 * r11231966;
double r11231977 = r11231970 - r11231976;
double r11231978 = 0.5;
double r11231979 = re;
double r11231980 = sin(r11231979);
double r11231981 = r11231978 * r11231980;
double r11231982 = r11231977 * r11231981;
return r11231982;
}




Bits error versus re




Bits error versus im
Results
| Original | 43.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.7 |
Initial program 43.7
Taylor expanded around 0 0.7
Simplified0.7
rmApplied pow10.7
Applied pow10.7
Applied pow10.7
Applied pow-prod-down0.7
Applied pow-prod-down0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2019172
(FPCore (re im)
:name "math.cos on complex, imaginary part"
:herbie-target
(if (< (fabs im) 1.0) (- (* (sin re) (+ (+ im (* (* (* 0.16666666666666666 im) im) im)) (* (* (* (* (* 0.008333333333333333 im) im) im) im) im)))) (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))
(* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))