\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)\left(im \cdot -2 + \left({im}^{3} \cdot \frac{-1}{3} - {im}^{5} \cdot \frac{1}{60}\right)\right) \cdot \left(0.5 \cdot \sin re\right)double f(double re, double im) {
double r143735 = 0.5;
double r143736 = re;
double r143737 = sin(r143736);
double r143738 = r143735 * r143737;
double r143739 = im;
double r143740 = -r143739;
double r143741 = exp(r143740);
double r143742 = exp(r143739);
double r143743 = r143741 - r143742;
double r143744 = r143738 * r143743;
return r143744;
}
double f(double re, double im) {
double r143745 = im;
double r143746 = -2.0;
double r143747 = r143745 * r143746;
double r143748 = 3.0;
double r143749 = pow(r143745, r143748);
double r143750 = -0.3333333333333333;
double r143751 = r143749 * r143750;
double r143752 = 5.0;
double r143753 = pow(r143745, r143752);
double r143754 = 0.016666666666666666;
double r143755 = r143753 * r143754;
double r143756 = r143751 - r143755;
double r143757 = r143747 + r143756;
double r143758 = 0.5;
double r143759 = re;
double r143760 = sin(r143759);
double r143761 = r143758 * r143760;
double r143762 = r143757 * r143761;
return r143762;
}




Bits error versus re




Bits error versus im
Results
| Original | 43.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.8 |
Initial program 43.7
Taylor expanded around 0 0.8
Simplified0.8
rmApplied pow10.8
Applied pow10.8
Applied pow10.8
Applied pow-prod-down0.8
Applied pow-prod-down0.8
Simplified0.8
Final simplification0.8
herbie shell --seed 2019194
(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))))