\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(\left(3 \cdot x.im\right) \cdot x.re\right) \cdot x.imdouble f(double x_re, double x_im) {
double r226135 = x_re;
double r226136 = r226135 * r226135;
double r226137 = x_im;
double r226138 = r226137 * r226137;
double r226139 = r226136 - r226138;
double r226140 = r226139 * r226135;
double r226141 = r226135 * r226137;
double r226142 = r226137 * r226135;
double r226143 = r226141 + r226142;
double r226144 = r226143 * r226137;
double r226145 = r226140 - r226144;
return r226145;
}
double f(double x_re, double x_im) {
double r226146 = x_re;
double r226147 = 3.0;
double r226148 = pow(r226146, r226147);
double r226149 = x_im;
double r226150 = r226147 * r226149;
double r226151 = r226150 * r226146;
double r226152 = r226151 * r226149;
double r226153 = r226148 - r226152;
return r226153;
}




Bits error versus x.re




Bits error versus x.im
Results
| Original | 7.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
Initial program 7.4
Simplified0.2
rmApplied associate-*r*0.2
rmApplied associate-*r*0.2
Final simplification0.2
herbie shell --seed 2020021
(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)))