double code(double x0, double x1) {
return ((double) (((double) (x0 / ((double) (1.0 - x1)))) - x0));
}
double code(double x0, double x1) {
double VAR;
if ((x1 <= 0.00021208908081054686)) {
VAR = ((double) (((double) (x0 * ((double) (((double) log(((double) cbrt(((double) exp(((double) (((double) (x0 / ((double) (1.0 + ((double) (x1 * ((double) (x1 - 2.0)))))))) - x0)))))))) + ((double) (((double) log(((double) cbrt(((double) exp(((double) (((double) (x0 / ((double) (1.0 + ((double) (x1 * ((double) (x1 - 2.0)))))))) - x0)))))))) * 2.0)))))) / ((double) (x0 + ((double) (x0 / ((double) (1.0 - x1))))))));
} else {
VAR = ((double) (((double) (x0 * ((double) log(((double) exp(((double) (((double) (x0 / ((double) (1.0 + ((double) (x1 * ((double) (x1 - 2.0)))))))) - x0)))))))) / ((double) (x0 + ((double) (1.0 / ((double) (((double) (1.0 - x1)) / x0))))))));
}
return VAR;
}




Bits error versus x0




Bits error versus x1
Results
| Original | 8.4 |
|---|---|
| Target | 0.5 |
| Herbie | 3.3 |
if x1 < 2.12089080810546861e-4Initial program 11.3
rmApplied flip--11.4
Simplified9.1
Simplified9.1
Taylor expanded around 0 11.3
Simplified9.1
rmApplied add-log-exp9.1
Applied add-log-exp9.1
Applied diff-log8.9
Simplified8.9
rmApplied add-cube-cbrt10.6
Applied log-prod9.9
Simplified6.6
if 2.12089080810546861e-4 < x1 Initial program 5.5
rmApplied flip--4.0
Simplified4.7
Simplified4.7
Taylor expanded around 0 2.8
Simplified2.8
rmApplied add-log-exp2.8
Applied add-log-exp2.8
Applied diff-log3.9
Simplified1.0
rmApplied clear-num0
Final simplification3.3
herbie shell --seed 2020191
(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))