\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 r10227823 = 0.5;
double r10227824 = re;
double r10227825 = sin(r10227824);
double r10227826 = r10227823 * r10227825;
double r10227827 = im;
double r10227828 = -r10227827;
double r10227829 = exp(r10227828);
double r10227830 = exp(r10227827);
double r10227831 = r10227829 - r10227830;
double r10227832 = r10227826 * r10227831;
return r10227832;
}
double f(double re, double im) {
double r10227833 = im;
double r10227834 = 5.0;
double r10227835 = pow(r10227833, r10227834);
double r10227836 = -0.016666666666666666;
double r10227837 = r10227835 * r10227836;
double r10227838 = 2.0;
double r10227839 = -0.3333333333333333;
double r10227840 = r10227833 * r10227833;
double r10227841 = r10227839 * r10227840;
double r10227842 = r10227838 - r10227841;
double r10227843 = r10227842 * r10227833;
double r10227844 = r10227837 - r10227843;
double r10227845 = 0.5;
double r10227846 = re;
double r10227847 = sin(r10227846);
double r10227848 = r10227845 * r10227847;
double r10227849 = r10227844 * r10227848;
return r10227849;
}




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))))