\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.91893853320467001\right) + \frac{\left(\left(y + 7.93650079365100015 \cdot 10^{-4}\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.0833333333333329956}{x}\left(\left(x - 0.5\right) \cdot \log \left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right) + \left(\left(\left(x - 0.5\right) \cdot \log \left(\sqrt[3]{x}\right) - x\right) + 0.91893853320467001\right)\right) + \frac{\left(\left(y + 7.93650079365100015 \cdot 10^{-4}\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.0833333333333329956}{x}double f(double x, double y, double z) {
double r429025 = x;
double r429026 = 0.5;
double r429027 = r429025 - r429026;
double r429028 = log(r429025);
double r429029 = r429027 * r429028;
double r429030 = r429029 - r429025;
double r429031 = 0.91893853320467;
double r429032 = r429030 + r429031;
double r429033 = y;
double r429034 = 0.0007936500793651;
double r429035 = r429033 + r429034;
double r429036 = z;
double r429037 = r429035 * r429036;
double r429038 = 0.0027777777777778;
double r429039 = r429037 - r429038;
double r429040 = r429039 * r429036;
double r429041 = 0.083333333333333;
double r429042 = r429040 + r429041;
double r429043 = r429042 / r429025;
double r429044 = r429032 + r429043;
return r429044;
}
double f(double x, double y, double z) {
double r429045 = x;
double r429046 = 0.5;
double r429047 = r429045 - r429046;
double r429048 = cbrt(r429045);
double r429049 = r429048 * r429048;
double r429050 = log(r429049);
double r429051 = r429047 * r429050;
double r429052 = log(r429048);
double r429053 = r429047 * r429052;
double r429054 = r429053 - r429045;
double r429055 = 0.91893853320467;
double r429056 = r429054 + r429055;
double r429057 = r429051 + r429056;
double r429058 = y;
double r429059 = 0.0007936500793651;
double r429060 = r429058 + r429059;
double r429061 = z;
double r429062 = r429060 * r429061;
double r429063 = 0.0027777777777778;
double r429064 = r429062 - r429063;
double r429065 = r429064 * r429061;
double r429066 = 0.083333333333333;
double r429067 = r429065 + r429066;
double r429068 = r429067 / r429045;
double r429069 = r429057 + r429068;
return r429069;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 5.8 |
|---|---|
| Target | 1.0 |
| Herbie | 5.8 |
Initial program 5.8
rmApplied add-cube-cbrt5.8
Applied log-prod5.9
Applied distribute-lft-in5.9
Applied associate--l+5.8
Applied associate-+l+5.8
Final simplification5.8
herbie shell --seed 2020046 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:$slogFactorial from math-functions-0.1.5.2, B"
:precision binary64
: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)))