\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)\left(-2 \cdot \left(re \cdot im\right)\right) \cdot 0.5
double f(double re, double im) {
double r9110052 = 0.5;
double r9110053 = re;
double r9110054 = sin(r9110053);
double r9110055 = r9110052 * r9110054;
double r9110056 = im;
double r9110057 = -r9110056;
double r9110058 = exp(r9110057);
double r9110059 = exp(r9110056);
double r9110060 = r9110058 - r9110059;
double r9110061 = r9110055 * r9110060;
return r9110061;
}
double f(double re, double im) {
double r9110062 = -2.0;
double r9110063 = re;
double r9110064 = im;
double r9110065 = r9110063 * r9110064;
double r9110066 = r9110062 * r9110065;
double r9110067 = 0.5;
double r9110068 = r9110066 * r9110067;
return r9110068;
}




Bits error versus re




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