\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)\left(\frac{-1}{3} \cdot \left(im \cdot \left(im \cdot im\right)\right) + \left({im}^{5} \cdot \frac{-1}{60} - \left(im + im\right)\right)\right) \cdot \left(0.5 \cdot \cos re\right)double f(double re, double im) {
double r8523981 = 0.5;
double r8523982 = re;
double r8523983 = cos(r8523982);
double r8523984 = r8523981 * r8523983;
double r8523985 = 0.0;
double r8523986 = im;
double r8523987 = r8523985 - r8523986;
double r8523988 = exp(r8523987);
double r8523989 = exp(r8523986);
double r8523990 = r8523988 - r8523989;
double r8523991 = r8523984 * r8523990;
return r8523991;
}
double f(double re, double im) {
double r8523992 = -0.3333333333333333;
double r8523993 = im;
double r8523994 = r8523993 * r8523993;
double r8523995 = r8523993 * r8523994;
double r8523996 = r8523992 * r8523995;
double r8523997 = 5.0;
double r8523998 = pow(r8523993, r8523997);
double r8523999 = -0.016666666666666666;
double r8524000 = r8523998 * r8523999;
double r8524001 = r8523993 + r8523993;
double r8524002 = r8524000 - r8524001;
double r8524003 = r8523996 + r8524002;
double r8524004 = 0.5;
double r8524005 = re;
double r8524006 = cos(r8524005);
double r8524007 = r8524004 * r8524006;
double r8524008 = r8524003 * r8524007;
return r8524008;
}




Bits error versus re




Bits error versus im
Results
| Original | 58.0 |
|---|---|
| Target | 0.3 |
| Herbie | 0.7 |
Initial program 58.0
Taylor expanded around 0 0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2019151
(FPCore (re im)
:name "math.sin on complex, imaginary part"
:herbie-target
(if (< (fabs im) 1) (- (* (cos re) (+ (+ im (* (* (* 1/6 im) im) im)) (* (* (* (* (* 1/120 im) im) im) im) im)))) (* (* 0.5 (cos re)) (- (exp (- 0 im)) (exp im))))
(* (* 0.5 (cos re)) (- (exp (- 0 im)) (exp im))))