\frac{x0}{1 - x1} - x0
\frac{e^{\log \left(\left(x0 \cdot x0\right) \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(-x1\right) \cdot -2\right)\right)}}{x0 + \frac{x0}{1 - x1}}
(FPCore (x0 x1) :precision binary64 (- (/ x0 (- 1.0 x1)) x0))
(FPCore (x0 x1) :precision binary64 (/ (exp (log (* (* x0 x0) (expm1 (* (log1p (- x1)) -2.0))))) (+ x0 (/ x0 (- 1.0 x1)))))
double code(double x0, double x1) {
return (x0 / (1.0 - x1)) - x0;
}
double code(double x0, double x1) {
return exp(log((x0 * x0) * expm1(log1p(-x1) * -2.0))) / (x0 + (x0 / (1.0 - x1)));
}




Bits error versus x0




Bits error versus x1
Results
| Original | 8.4 |
|---|---|
| Target | 0.5 |
| Herbie | 0 |
Initial program 8.4
rmApplied pow1_binary648.4
Applied sqr-pow_binary647.5
Simplified7.5
Simplified7.5
rmApplied flip--_binary647.7
Simplified7.8
Simplified7.7
Taylor expanded around 0 7.3
Simplified0.5
rmApplied add-exp-log_binary640
Final simplification0
herbie shell --seed 2021211
(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))