\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)\left(0.5 \cdot \cos re\right) \cdot \left(\left(\left(\left(im \cdot im\right) \cdot \frac{-1}{3}\right) \cdot im - im \cdot 2\right) - \frac{1}{60} \cdot {im}^{5}\right)double f(double re, double im) {
double r6730574 = 0.5;
double r6730575 = re;
double r6730576 = cos(r6730575);
double r6730577 = r6730574 * r6730576;
double r6730578 = 0.0;
double r6730579 = im;
double r6730580 = r6730578 - r6730579;
double r6730581 = exp(r6730580);
double r6730582 = exp(r6730579);
double r6730583 = r6730581 - r6730582;
double r6730584 = r6730577 * r6730583;
return r6730584;
}
double f(double re, double im) {
double r6730585 = 0.5;
double r6730586 = re;
double r6730587 = cos(r6730586);
double r6730588 = r6730585 * r6730587;
double r6730589 = im;
double r6730590 = r6730589 * r6730589;
double r6730591 = -0.3333333333333333;
double r6730592 = r6730590 * r6730591;
double r6730593 = r6730592 * r6730589;
double r6730594 = 2.0;
double r6730595 = r6730589 * r6730594;
double r6730596 = r6730593 - r6730595;
double r6730597 = 0.016666666666666666;
double r6730598 = 5.0;
double r6730599 = pow(r6730589, r6730598);
double r6730600 = r6730597 * r6730599;
double r6730601 = r6730596 - r6730600;
double r6730602 = r6730588 * r6730601;
return r6730602;
}




Bits error versus re




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