\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.im + \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.re
3 \cdot \left(\left(x.im \cdot x.re\right) \cdot x.re\right) - {x.im}^{3}double f(double x_re, double x_im) {
double r181818 = x_re;
double r181819 = r181818 * r181818;
double r181820 = x_im;
double r181821 = r181820 * r181820;
double r181822 = r181819 - r181821;
double r181823 = r181822 * r181820;
double r181824 = r181818 * r181820;
double r181825 = r181820 * r181818;
double r181826 = r181824 + r181825;
double r181827 = r181826 * r181818;
double r181828 = r181823 + r181827;
return r181828;
}
double f(double x_re, double x_im) {
double r181829 = 3.0;
double r181830 = x_im;
double r181831 = x_re;
double r181832 = r181830 * r181831;
double r181833 = r181832 * r181831;
double r181834 = r181829 * r181833;
double r181835 = pow(r181830, r181829);
double r181836 = r181834 - r181835;
return r181836;
}




Bits error versus x.re




Bits error versus x.im
Results
| Original | 7.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 7.7
Simplified7.7
rmApplied associate-*r*0.2
Final simplification0.2
herbie shell --seed 2020036
(FPCore (x.re x.im)
:name "math.cube on complex, imaginary part"
:precision binary64
:herbie-target
(+ (* (* x.re x.im) (* 2 x.re)) (* (* x.im (- x.re x.im)) (+ x.re x.im)))
(+ (* (- (* x.re x.re) (* x.im x.im)) x.im) (* (+ (* x.re x.im) (* x.im x.re)) x.re)))