\frac{x0}{1 - x1} - x0\frac{\frac{\frac{x0}{\sqrt{1 - x1}}}{\sqrt{1 - x1}} \cdot \frac{x0}{1 - x1} - x0 \cdot x0}{\left(\sqrt[3]{x0 + \frac{x0}{1 - x1}} \cdot \sqrt[3]{x0 + \frac{x0}{1 - x1}}\right) \cdot \sqrt[3]{x0 + \frac{x0}{1 - x1}}}double f(double x0, double x1) {
double r128044 = x0;
double r128045 = 1.0;
double r128046 = x1;
double r128047 = r128045 - r128046;
double r128048 = r128044 / r128047;
double r128049 = r128048 - r128044;
return r128049;
}
double f(double x0, double x1) {
double r128050 = x0;
double r128051 = 1.0;
double r128052 = x1;
double r128053 = r128051 - r128052;
double r128054 = sqrt(r128053);
double r128055 = r128050 / r128054;
double r128056 = r128055 / r128054;
double r128057 = r128050 / r128053;
double r128058 = r128056 * r128057;
double r128059 = r128050 * r128050;
double r128060 = r128058 - r128059;
double r128061 = r128050 + r128057;
double r128062 = cbrt(r128061);
double r128063 = r128062 * r128062;
double r128064 = r128063 * r128062;
double r128065 = r128060 / r128064;
return r128065;
}




Bits error versus x0




Bits error versus x1
Results
| Original | 7.9 |
|---|---|
| Target | 0.3 |
| Herbie | 5.6 |
Initial program 7.9
rmApplied flip--7.3
Simplified7.3
rmApplied add-sqr-sqrt5.6
Applied associate-/r*5.6
rmApplied add-cube-cbrt5.6
Final simplification5.6
herbie shell --seed 2019199
(FPCore (x0 x1)
:name "(- (/ x0 (- 1 x1)) x0)"
:pre (or (and (== x0 1.855) (== x1 0.000209)) (and (== x0 2.985) (== x1 0.0186)))
:herbie-target
(/ (* x0 x1) (- 1.0 x1))
(- (/ x0 (- 1.0 x1)) x0))