\left(0.5 \cdot \cos re\right) \cdot \left(e^{0.0 - im} - e^{im}\right)\left(0.5 \cdot \cos re\right) \cdot \left({im}^{3} \cdot \frac{-1}{3} - \left(\frac{1}{60} \cdot {im}^{5} + 2 \cdot im\right)\right)double f(double re, double im) {
double r129957 = 0.5;
double r129958 = re;
double r129959 = cos(r129958);
double r129960 = r129957 * r129959;
double r129961 = 0.0;
double r129962 = im;
double r129963 = r129961 - r129962;
double r129964 = exp(r129963);
double r129965 = exp(r129962);
double r129966 = r129964 - r129965;
double r129967 = r129960 * r129966;
return r129967;
}
double f(double re, double im) {
double r129968 = 0.5;
double r129969 = re;
double r129970 = cos(r129969);
double r129971 = r129968 * r129970;
double r129972 = im;
double r129973 = 3.0;
double r129974 = pow(r129972, r129973);
double r129975 = -0.3333333333333333;
double r129976 = r129974 * r129975;
double r129977 = 0.016666666666666666;
double r129978 = 5.0;
double r129979 = pow(r129972, r129978);
double r129980 = r129977 * r129979;
double r129981 = 2.0;
double r129982 = r129981 * r129972;
double r129983 = r129980 + r129982;
double r129984 = r129976 - r129983;
double r129985 = r129971 * r129984;
return r129985;
}




Bits error versus re




Bits error versus im
Results
| Original | 58.1 |
|---|---|
| Target | 0.2 |
| Herbie | 0.7 |
Initial program 58.1
Taylor expanded around 0 0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2019323
(FPCore (re im)
:name "math.sin on complex, imaginary part"
:precision binary64
:herbie-target
(if (< (fabs im) 1) (- (* (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))))