\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{x}{\frac{\left(\sqrt[3]{\sqrt[3]{y}} \cdot \sqrt[3]{\sqrt[3]{y}}\right) \cdot \sqrt[3]{\sqrt[3]{y}}}{\sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}}}} \cdot \frac{1}{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}}}}double f(double x, double y, double z, double t, double a, double b) {
double r20225001 = x;
double r20225002 = y;
double r20225003 = z;
double r20225004 = log(r20225003);
double r20225005 = r20225002 * r20225004;
double r20225006 = t;
double r20225007 = 1.0;
double r20225008 = r20225006 - r20225007;
double r20225009 = a;
double r20225010 = log(r20225009);
double r20225011 = r20225008 * r20225010;
double r20225012 = r20225005 + r20225011;
double r20225013 = b;
double r20225014 = r20225012 - r20225013;
double r20225015 = exp(r20225014);
double r20225016 = r20225001 * r20225015;
double r20225017 = r20225016 / r20225002;
return r20225017;
}
double f(double x, double y, double z, double t, double a, double b) {
double r20225018 = x;
double r20225019 = y;
double r20225020 = cbrt(r20225019);
double r20225021 = cbrt(r20225020);
double r20225022 = r20225021 * r20225021;
double r20225023 = r20225022 * r20225021;
double r20225024 = a;
double r20225025 = log(r20225024);
double r20225026 = t;
double r20225027 = 1.0;
double r20225028 = r20225026 - r20225027;
double r20225029 = r20225025 * r20225028;
double r20225030 = z;
double r20225031 = log(r20225030);
double r20225032 = r20225031 * r20225019;
double r20225033 = r20225029 + r20225032;
double r20225034 = b;
double r20225035 = r20225033 - r20225034;
double r20225036 = exp(r20225035);
double r20225037 = cbrt(r20225036);
double r20225038 = r20225023 / r20225037;
double r20225039 = r20225018 / r20225038;
double r20225040 = 1.0;
double r20225041 = r20225020 * r20225020;
double r20225042 = r20225037 * r20225037;
double r20225043 = r20225041 / r20225042;
double r20225044 = r20225040 / r20225043;
double r20225045 = r20225039 * r20225044;
return r20225045;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 2.0 |
|---|---|
| Target | 10.8 |
| Herbie | 1.1 |
Initial program 2.0
rmApplied associate-/l*1.9
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.9
Applied times-frac1.9
Applied *-un-lft-identity1.9
Applied times-frac1.1
rmApplied add-cube-cbrt1.1
Final simplification1.1
herbie shell --seed 2019174
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1.0))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))