x \cdot \left(\frac{y}{z} - \frac{t}{1 - z}\right)x \cdot \left(\frac{y}{z} - t \cdot \frac{1}{1 - z}\right)double f(double x, double y, double z, double t) {
double r468392 = x;
double r468393 = y;
double r468394 = z;
double r468395 = r468393 / r468394;
double r468396 = t;
double r468397 = 1.0;
double r468398 = r468397 - r468394;
double r468399 = r468396 / r468398;
double r468400 = r468395 - r468399;
double r468401 = r468392 * r468400;
return r468401;
}
double f(double x, double y, double z, double t) {
double r468402 = x;
double r468403 = y;
double r468404 = z;
double r468405 = r468403 / r468404;
double r468406 = t;
double r468407 = 1.0;
double r468408 = 1.0;
double r468409 = r468408 - r468404;
double r468410 = r468407 / r468409;
double r468411 = r468406 * r468410;
double r468412 = r468405 - r468411;
double r468413 = r468402 * r468412;
return r468413;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 4.8 |
|---|---|
| Target | 4.4 |
| Herbie | 4.8 |
Initial program 4.8
rmApplied div-inv4.8
Final simplification4.8
herbie shell --seed 2020047 +o rules:numerics
(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 z)))) -7.623226303312042e-196) (* x (- (/ y z) (* t (/ 1 (- 1 z))))) (if (< (* x (- (/ y z) (/ t (- 1 z)))) 1.4133944927702302e-211) (+ (/ (* y x) z) (- (/ (* t x) (- 1 z)))) (* x (- (/ y z) (* t (/ 1 (- 1 z)))))))
(* x (- (/ y z) (/ t (- 1 z)))))