\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(\log \left(\sqrt{y}\right) \cdot x + \left(\log \left(\sqrt{y}\right) \cdot x + \left(\left(\log 1 - y \cdot 1\right) \cdot z - \left(\left(z \cdot 0.5\right) \cdot y\right) \cdot y\right)\right)\right) - tdouble f(double x, double y, double z, double t) {
double r8916988 = x;
double r8916989 = y;
double r8916990 = log(r8916989);
double r8916991 = r8916988 * r8916990;
double r8916992 = z;
double r8916993 = 1.0;
double r8916994 = r8916993 - r8916989;
double r8916995 = log(r8916994);
double r8916996 = r8916992 * r8916995;
double r8916997 = r8916991 + r8916996;
double r8916998 = t;
double r8916999 = r8916997 - r8916998;
return r8916999;
}
double f(double x, double y, double z, double t) {
double r8917000 = y;
double r8917001 = sqrt(r8917000);
double r8917002 = log(r8917001);
double r8917003 = x;
double r8917004 = r8917002 * r8917003;
double r8917005 = 1.0;
double r8917006 = log(r8917005);
double r8917007 = r8917000 * r8917005;
double r8917008 = r8917006 - r8917007;
double r8917009 = z;
double r8917010 = r8917008 * r8917009;
double r8917011 = 0.5;
double r8917012 = r8917009 * r8917011;
double r8917013 = r8917012 * r8917000;
double r8917014 = r8917013 * r8917000;
double r8917015 = r8917010 - r8917014;
double r8917016 = r8917004 + r8917015;
double r8917017 = r8917004 + r8917016;
double r8917018 = t;
double r8917019 = r8917017 - r8917018;
return r8917019;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 9.6
Taylor expanded around 0 0.4
Simplified0.4
rmApplied add-sqr-sqrt0.4
Applied log-prod0.4
Applied distribute-rgt-in0.4
Applied associate-+l+0.4
Taylor expanded around 0 0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019179
(FPCore (x y z t)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, B"
:herbie-target
(- (* (- z) (+ (+ (* 0.5 (* y y)) y) (* (/ 0.3333333333333333 (* 1.0 (* 1.0 1.0))) (* y (* y y))))) (- t (* x (log y))))
(- (+ (* x (log y)) (* z (log (- 1.0 y)))) t))