\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} - 3 \cdot \left(x.im \cdot \left(x.re \cdot x.im\right)\right)double f(double x_re, double x_im) {
double r223171 = x_re;
double r223172 = r223171 * r223171;
double r223173 = x_im;
double r223174 = r223173 * r223173;
double r223175 = r223172 - r223174;
double r223176 = r223175 * r223171;
double r223177 = r223171 * r223173;
double r223178 = r223173 * r223171;
double r223179 = r223177 + r223178;
double r223180 = r223179 * r223173;
double r223181 = r223176 - r223180;
return r223181;
}
double f(double x_re, double x_im) {
double r223182 = x_re;
double r223183 = 3.0;
double r223184 = pow(r223182, r223183);
double r223185 = x_im;
double r223186 = r223182 * r223185;
double r223187 = r223185 * r223186;
double r223188 = r223183 * r223187;
double r223189 = r223184 - r223188;
return r223189;
}




Bits error versus x.re




Bits error versus x.im
Results
| Original | 7.0 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 7.0
Simplified0.2
Final simplification0.2
herbie shell --seed 2020024
(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)))