x \cdot \left(\frac{y}{z} - \frac{t}{1 - z}\right)\left(x \cdot \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{z} \cdot \sqrt[3]{z}}\right) \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{z}} - x \cdot \frac{t}{1 - z}(FPCore (x y z t) :precision binary64 (* x (- (/ y z) (/ t (- 1.0 z)))))
(FPCore (x y z t) :precision binary64 (- (* (* x (/ (* (cbrt y) (cbrt y)) (* (cbrt z) (cbrt z)))) (/ (cbrt y) (cbrt z))) (* x (/ t (- 1.0 z)))))
double code(double x, double y, double z, double t) {
return x * ((y / z) - (t / (1.0 - z)));
}
double code(double x, double y, double z, double t) {
return ((x * ((cbrt(y) * cbrt(y)) / (cbrt(z) * cbrt(z)))) * (cbrt(y) / cbrt(z))) - (x * (t / (1.0 - z)));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 4.5 |
|---|---|
| Target | 4.4 |
| Herbie | 1.9 |
Initial program 4.5
rmApplied sub-neg_binary644.5
Applied distribute-rgt-in_binary644.5
Simplified4.5
Simplified4.5
rmApplied add-cube-cbrt_binary645.0
Applied add-cube-cbrt_binary645.1
Applied times-frac_binary645.1
Applied associate-*r*_binary641.9
Final simplification1.9
herbie shell --seed 2020273
(FPCore (x y z t)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, C"
:precision binary64
:herbie-target
(if (< (* x (- (/ y z) (/ t (- 1.0 z)))) -7.623226303312042e-196) (* x (- (/ y z) (* t (/ 1.0 (- 1.0 z))))) (if (< (* x (- (/ y z) (/ t (- 1.0 z)))) 1.4133944927702302e-211) (+ (/ (* y x) z) (- (/ (* t x) (- 1.0 z)))) (* x (- (/ y z) (* t (/ 1.0 (- 1.0 z)))))))
(* x (- (/ y z) (/ t (- 1.0 z)))))