\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)\left(\left(im \cdot \left(im \cdot im\right)\right) \cdot \frac{-1}{3} - \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 r4810999 = 0.5;
double r4811000 = re;
double r4811001 = cos(r4811000);
double r4811002 = r4810999 * r4811001;
double r4811003 = 0.0;
double r4811004 = im;
double r4811005 = r4811003 - r4811004;
double r4811006 = exp(r4811005);
double r4811007 = exp(r4811004);
double r4811008 = r4811006 - r4811007;
double r4811009 = r4811002 * r4811008;
return r4811009;
}
double f(double re, double im) {
double r4811010 = im;
double r4811011 = r4811010 * r4811010;
double r4811012 = r4811010 * r4811011;
double r4811013 = -0.3333333333333333;
double r4811014 = r4811012 * r4811013;
double r4811015 = 5.0;
double r4811016 = pow(r4811010, r4811015);
double r4811017 = 0.016666666666666666;
double r4811018 = r4811016 * r4811017;
double r4811019 = r4811010 + r4811010;
double r4811020 = r4811018 + r4811019;
double r4811021 = r4811014 - r4811020;
double r4811022 = 0.5;
double r4811023 = re;
double r4811024 = cos(r4811023);
double r4811025 = r4811022 * r4811024;
double r4811026 = r4811021 * r4811025;
return r4811026;
}




Bits error versus re




Bits error versus im
Results
| Original | 57.8 |
|---|---|
| Target | 0.3 |
| Herbie | 0.9 |
Initial program 57.8
Taylor expanded around 0 0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2019154
(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))))