\left(\left(\left(\sqrt{x + 1} - \sqrt{x}\right) + \left(\sqrt{y + 1} - \sqrt{y}\right)\right) + \left(\sqrt{z + 1} - \sqrt{z}\right)\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\left(\left(\frac{1}{\sqrt{x + 1} + \sqrt{x}} + \left(\sqrt{y + 1} - \sqrt{y}\right)\right) + \left(\sqrt{z + 1} - \sqrt{z}\right)\right) + \frac{1}{\sqrt{t + 1} + \sqrt{t}}double code(double x, double y, double z, double t) {
return ((double) (((double) (((double) (((double) (((double) sqrt(((double) (x + 1.0)))) - ((double) sqrt(x)))) + ((double) (((double) sqrt(((double) (y + 1.0)))) - ((double) sqrt(y)))))) + ((double) (((double) sqrt(((double) (z + 1.0)))) - ((double) sqrt(z)))))) + ((double) (((double) sqrt(((double) (t + 1.0)))) - ((double) sqrt(t))))));
}
double code(double x, double y, double z, double t) {
return ((double) (((double) (((double) (((double) (1.0 / ((double) (((double) sqrt(((double) (x + 1.0)))) + ((double) sqrt(x)))))) + ((double) (((double) sqrt(((double) (y + 1.0)))) - ((double) sqrt(y)))))) + ((double) (((double) sqrt(((double) (z + 1.0)))) - ((double) sqrt(z)))))) + ((double) (1.0 / ((double) (((double) sqrt(((double) (t + 1.0)))) + ((double) sqrt(t))))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 5.4 |
|---|---|
| Target | 1.5 |
| Herbie | 2.9 |
Initial program 5.4
rmApplied flip--5.4
Simplified5.1
Taylor expanded around 0 4.2
rmApplied flip--4.1
Simplified3.9
Taylor expanded around 0 2.9
Final simplification2.9
herbie shell --seed 2020169
(FPCore (x y z t)
:name "Main:z from "
:precision binary64
:herbie-target
(+ (+ (+ (/ 1.0 (+ (sqrt (+ x 1.0)) (sqrt x))) (/ 1.0 (+ (sqrt (+ y 1.0)) (sqrt y)))) (/ 1.0 (+ (sqrt (+ z 1.0)) (sqrt z)))) (- (sqrt (+ t 1.0)) (sqrt t)))
(+ (+ (+ (- (sqrt (+ x 1.0)) (sqrt x)) (- (sqrt (+ y 1.0)) (sqrt y))) (- (sqrt (+ z 1.0)) (sqrt z))) (- (sqrt (+ t 1.0)) (sqrt t))))