double code(double x, double y, double z, double t) {
return ((double) (x - ((double) (((double) log(((double) (((double) (1.0 - y)) + ((double) (y * ((double) exp(z)))))))) / t))));
}
double code(double x, double y, double z, double t) {
double VAR;
if ((z <= -1.1970149978622335e-08)) {
VAR = ((double) (x + ((double) (((double) log(((double) (1.0 + ((double) (((double) (y * ((double) exp(z)))) - y)))))) * ((double) (-1.0 / t))))));
} else {
double VAR_1;
if ((z <= 4.0973833117696044e-84)) {
VAR_1 = ((double) (x - ((double) (((double) (1.0 * ((double) (y * ((double) (z / t)))))) + ((double) (((double) log(1.0)) / t))))));
} else {
VAR_1 = ((double) (x - ((double) (((double) log(((double) (1.0 + ((double) (y * ((double) (((double) (0.16666666666666666 * ((double) pow(z, 3.0)))) + ((double) (z + ((double) (z * ((double) (z * 0.5)))))))))))))) / t))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 25.0 |
|---|---|
| Target | 16.6 |
| Herbie | 8.3 |
if z < -1.19701499786223348e-8Initial program 12.3
Simplified12.3
rmApplied div-inv12.4
if -1.19701499786223348e-8 < z < 4.09738331176960443e-84Initial program 30.5
Simplified15.2
Taylor expanded around 0 6.4
Simplified5.3
if 4.09738331176960443e-84 < z Initial program 31.5
Simplified23.0
Taylor expanded around 0 14.1
Simplified14.1
Final simplification8.3
herbie shell --seed 2020191
(FPCore (x y z t)
:name "System.Random.MWC.Distributions:truncatedExp from mwc-random-0.13.3.2"
:precision binary64
:herbie-target
(if (< z -2.8874623088207947e+119) (- (- x (/ (/ (neg 0.5) (* y t)) (* z z))) (* (/ (neg 0.5) (* y t)) (/ (/ 2.0 z) (* z z)))) (- x (/ (log (+ 1.0 (* z y))) t)))
(- x (/ (log (+ (- 1.0 y) (* y (exp z)))) t)))