\left(0.5 \cdot \cos re\right) \cdot \left(e^{0.0 - im} - e^{im}\right)0.5 \cdot \left(\left({im}^{3} \cdot \frac{-1}{3} - \left(\left(im + im\right) + {im}^{5} \cdot \frac{1}{60}\right)\right) \cdot \cos re\right)double f(double re, double im) {
double r272415 = 0.5;
double r272416 = re;
double r272417 = cos(r272416);
double r272418 = r272415 * r272417;
double r272419 = 0.0;
double r272420 = im;
double r272421 = r272419 - r272420;
double r272422 = exp(r272421);
double r272423 = exp(r272420);
double r272424 = r272422 - r272423;
double r272425 = r272418 * r272424;
return r272425;
}
double f(double re, double im) {
double r272426 = 0.5;
double r272427 = im;
double r272428 = 3.0;
double r272429 = pow(r272427, r272428);
double r272430 = -0.3333333333333333;
double r272431 = r272429 * r272430;
double r272432 = r272427 + r272427;
double r272433 = 5.0;
double r272434 = pow(r272427, r272433);
double r272435 = 0.016666666666666666;
double r272436 = r272434 * r272435;
double r272437 = r272432 + r272436;
double r272438 = r272431 - r272437;
double r272439 = re;
double r272440 = cos(r272439);
double r272441 = r272438 * r272440;
double r272442 = r272426 * r272441;
return r272442;
}




Bits error versus re




Bits error versus im
Results
| Original | 58.0 |
|---|---|
| Target | 0.3 |
| Herbie | 0.8 |
Initial program 58.0
Simplified58.0
Taylor expanded around 0 0.8
Simplified0.8
Final simplification0.8
herbie shell --seed 2019179
(FPCore (re im)
:name "math.sin on complex, imaginary part"
:herbie-target
(if (< (fabs im) 1.0) (- (* (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))))