\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 r180463 = x_re;
double r180464 = r180463 * r180463;
double r180465 = x_im;
double r180466 = r180465 * r180465;
double r180467 = r180464 - r180466;
double r180468 = r180467 * r180463;
double r180469 = r180463 * r180465;
double r180470 = r180465 * r180463;
double r180471 = r180469 + r180470;
double r180472 = r180471 * r180465;
double r180473 = r180468 - r180472;
return r180473;
}
double f(double x_re, double x_im) {
double r180474 = x_re;
double r180475 = 3.0;
double r180476 = pow(r180474, r180475);
double r180477 = x_im;
double r180478 = r180475 * r180477;
double r180479 = r180478 * r180474;
double r180480 = r180479 * r180477;
double r180481 = r180476 - r180480;
return r180481;
}




Bits error versus x.re




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