\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 r123332 = 0.5;
double r123333 = re;
double r123334 = cos(r123333);
double r123335 = r123332 * r123334;
double r123336 = 0.0;
double r123337 = im;
double r123338 = r123336 - r123337;
double r123339 = exp(r123338);
double r123340 = exp(r123337);
double r123341 = r123339 - r123340;
double r123342 = r123335 * r123341;
return r123342;
}
double f(double re, double im) {
double r123343 = 0.5;
double r123344 = re;
double r123345 = cos(r123344);
double r123346 = r123343 * r123345;
double r123347 = im;
double r123348 = 3.0;
double r123349 = pow(r123347, r123348);
double r123350 = -0.3333333333333333;
double r123351 = r123349 * r123350;
double r123352 = 0.016666666666666666;
double r123353 = 5.0;
double r123354 = pow(r123347, r123353);
double r123355 = r123352 * r123354;
double r123356 = 2.0;
double r123357 = r123356 * r123347;
double r123358 = r123355 + r123357;
double r123359 = r123351 - r123358;
double r123360 = r123346 * r123359;
return r123360;
}




Bits error versus re




Bits error versus im
Results
| Original | 58.0 |
|---|---|
| Target | 0.2 |
| Herbie | 0.8 |
Initial program 58.0
Taylor expanded around 0 0.8
Simplified0.8
Final simplification0.8
herbie shell --seed 2019325
(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))))