\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.91893853320467\right) + \frac{\left(\left(y + 0.0007936500793651\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.083333333333333}{x}\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)}{x} + \left(\left(0.91893853320467 + \left(\left(\sqrt{x} - \sqrt{0.5}\right) \cdot \left(\left(\sqrt{x} + \sqrt{0.5}\right) \cdot \log \left(\sqrt{x}\right)\right) - x\right)\right) + \left(x - 0.5\right) \cdot \log \left(\sqrt{x}\right)\right)double f(double x, double y, double z) {
double r21787839 = x;
double r21787840 = 0.5;
double r21787841 = r21787839 - r21787840;
double r21787842 = log(r21787839);
double r21787843 = r21787841 * r21787842;
double r21787844 = r21787843 - r21787839;
double r21787845 = 0.91893853320467;
double r21787846 = r21787844 + r21787845;
double r21787847 = y;
double r21787848 = 0.0007936500793651;
double r21787849 = r21787847 + r21787848;
double r21787850 = z;
double r21787851 = r21787849 * r21787850;
double r21787852 = 0.0027777777777778;
double r21787853 = r21787851 - r21787852;
double r21787854 = r21787853 * r21787850;
double r21787855 = 0.083333333333333;
double r21787856 = r21787854 + r21787855;
double r21787857 = r21787856 / r21787839;
double r21787858 = r21787846 + r21787857;
return r21787858;
}
double f(double x, double y, double z) {
double r21787859 = 0.083333333333333;
double r21787860 = z;
double r21787861 = y;
double r21787862 = 0.0007936500793651;
double r21787863 = r21787861 + r21787862;
double r21787864 = r21787860 * r21787863;
double r21787865 = 0.0027777777777778;
double r21787866 = r21787864 - r21787865;
double r21787867 = r21787860 * r21787866;
double r21787868 = r21787859 + r21787867;
double r21787869 = x;
double r21787870 = r21787868 / r21787869;
double r21787871 = 0.91893853320467;
double r21787872 = sqrt(r21787869);
double r21787873 = 0.5;
double r21787874 = sqrt(r21787873);
double r21787875 = r21787872 - r21787874;
double r21787876 = r21787872 + r21787874;
double r21787877 = log(r21787872);
double r21787878 = r21787876 * r21787877;
double r21787879 = r21787875 * r21787878;
double r21787880 = r21787879 - r21787869;
double r21787881 = r21787871 + r21787880;
double r21787882 = r21787869 - r21787873;
double r21787883 = r21787882 * r21787877;
double r21787884 = r21787881 + r21787883;
double r21787885 = r21787870 + r21787884;
return r21787885;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 5.9 |
|---|---|
| Target | 1.1 |
| Herbie | 6.0 |
Initial program 5.9
rmApplied add-sqr-sqrt5.9
Applied log-prod5.9
Applied distribute-rgt-in5.9
Applied associate--l+5.9
Applied associate-+l+5.9
rmApplied add-sqr-sqrt5.9
Applied add-sqr-sqrt6.0
Applied difference-of-squares6.0
Applied associate-*r*6.0
Final simplification6.0
herbie shell --seed 2019162
(FPCore (x y z)
:name "Numeric.SpecFunctions:$slogFactorial from math-functions-0.1.5.2, B"
:herbie-target
(+ (+ (+ (* (- x 0.5) (log x)) (- 0.91893853320467 x)) (/ 0.083333333333333 x)) (* (/ z x) (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))
(+ (+ (- (* (- x 0.5) (log x)) x) 0.91893853320467) (/ (+ (* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z) 0.083333333333333) x)))