\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
{x.re}^{3} - \left(3 \cdot x.im\right) \cdot \left(x.re \cdot x.im\right)double f(double x_re, double x_im) {
double r808 = x_re;
double r809 = r808 * r808;
double r810 = x_im;
double r811 = r810 * r810;
double r812 = r809 - r811;
double r813 = r812 * r808;
double r814 = r808 * r810;
double r815 = r810 * r808;
double r816 = r814 + r815;
double r817 = r816 * r810;
double r818 = r813 - r817;
return r818;
}
double f(double x_re, double x_im) {
double r819 = x_re;
double r820 = 3.0;
double r821 = pow(r819, r820);
double r822 = x_im;
double r823 = r820 * r822;
double r824 = r819 * r822;
double r825 = r823 * r824;
double r826 = r821 - r825;
return r826;
}




Bits error versus x.re




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