\frac{x0}{1 - x1} - x0\begin{array}{l}
\mathbf{if}\;x1 \le 0.0023027651367187492:\\
\;\;\;\;\frac{x0 \cdot \frac{\frac{{\left(\frac{{x0}^{3}}{{\left(1 - x1\right)}^{6}}\right)}^{3} - {\left({x0}^{3}\right)}^{3}}{\frac{{x0}^{6}}{{\left(\sqrt{1} + \sqrt{x1}\right)}^{12} \cdot {\left(\sqrt{1} - \sqrt{x1}\right)}^{12}} + \left({x0}^{6} + {\left(\frac{x0}{1 - x1}\right)}^{6}\right)}}{\frac{x0}{\frac{{\left(1 - x1\right)}^{4}}{x0}} + x0 \cdot \left(x0 + \frac{x0}{\left(1 - x1\right) \cdot \left(1 - x1\right)}\right)}}{x0 + \frac{x0}{1 - x1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x0 \cdot \frac{\frac{{\left(\frac{{x0}^{3}}{{\left(\sqrt{1} + \sqrt{x1}\right)}^{6} \cdot {\left(\sqrt{1} - \sqrt{x1}\right)}^{6}}\right)}^{3} - {\left({x0}^{3}\right)}^{3}}{\left({x0}^{6} + {\left(\frac{x0}{1 - x1}\right)}^{6}\right) + \frac{{x0}^{6}}{{\left(1 - x1\right)}^{12}}}}{\frac{x0}{\frac{{\left(1 - x1\right)}^{4}}{x0}} + x0 \cdot \left(x0 + \frac{x0}{\left(1 - x1\right) \cdot \left(1 - x1\right)}\right)}}{x0 + \frac{x0}{1 - x1}}\\
\end{array}double code(double x0, double x1) {
return ((double) ((x0 / ((double) (1.0 - x1))) - x0));
}
double code(double x0, double x1) {
double VAR;
if ((x1 <= 0.002302765136718749)) {
VAR = (((double) (x0 * ((((double) (((double) pow((((double) pow(x0, 3.0)) / ((double) pow(((double) (1.0 - x1)), 6.0))), 3.0)) - ((double) pow(((double) pow(x0, 3.0)), 3.0)))) / ((double) ((((double) pow(x0, 6.0)) / ((double) (((double) pow(((double) (((double) sqrt(1.0)) + ((double) sqrt(x1)))), 12.0)) * ((double) pow(((double) (((double) sqrt(1.0)) - ((double) sqrt(x1)))), 12.0))))) + ((double) (((double) pow(x0, 6.0)) + ((double) pow((x0 / ((double) (1.0 - x1))), 6.0))))))) / ((double) ((x0 / (((double) pow(((double) (1.0 - x1)), 4.0)) / x0)) + ((double) (x0 * ((double) (x0 + (x0 / ((double) (((double) (1.0 - x1)) * ((double) (1.0 - x1)))))))))))))) / ((double) (x0 + (x0 / ((double) (1.0 - x1))))));
} else {
VAR = (((double) (x0 * ((((double) (((double) pow((((double) pow(x0, 3.0)) / ((double) (((double) pow(((double) (((double) sqrt(1.0)) + ((double) sqrt(x1)))), 6.0)) * ((double) pow(((double) (((double) sqrt(1.0)) - ((double) sqrt(x1)))), 6.0))))), 3.0)) - ((double) pow(((double) pow(x0, 3.0)), 3.0)))) / ((double) (((double) (((double) pow(x0, 6.0)) + ((double) pow((x0 / ((double) (1.0 - x1))), 6.0)))) + (((double) pow(x0, 6.0)) / ((double) pow(((double) (1.0 - x1)), 12.0)))))) / ((double) ((x0 / (((double) pow(((double) (1.0 - x1)), 4.0)) / x0)) + ((double) (x0 * ((double) (x0 + (x0 / ((double) (((double) (1.0 - x1)) * ((double) (1.0 - x1)))))))))))))) / ((double) (x0 + (x0 / ((double) (1.0 - x1))))));
}
return VAR;
}




Bits error versus x0




Bits error versus x1
Results
| Original | 8.3 |
|---|---|
| Target | 0.5 |
| Herbie | 3.5 |
if x1 < 0.0023027651367187492Initial program 11.3
rmApplied flip--11.4
Simplified9.1
Simplified9.1
rmApplied flip3--7.8
Simplified7.8
Simplified7.8
rmApplied flip3--6.3
Simplified6.3
rmApplied add-sqr-sqrt6.3
Applied add-sqr-sqrt6.3
Applied difference-of-squares6.2
Applied unpow-prod-down6.2
if 0.0023027651367187492 < x1 Initial program 5.5
rmApplied flip--4.0
Simplified4.7
Simplified4.7
rmApplied flip3--4.4
Simplified4.4
Simplified4.5
rmApplied flip3--4.4
Simplified4.3
rmApplied add-sqr-sqrt4.3
Applied add-sqr-sqrt4.3
Applied difference-of-squares4.3
Applied unpow-prod-down1.0
Final simplification3.5
herbie shell --seed 2020182
(FPCore (x0 x1)
:name "(- (/ x0 (- 1 x1)) x0)"
:precision binary64
: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))