\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(\log \left(\sqrt{x}\right) \cdot \left(x - 0.5\right) + \left(\log \left(\sqrt{x}\right) \cdot \left(x - 0.5\right) - x\right)\right) + 0.91893853320467001\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 r426401 = x;
double r426402 = 0.5;
double r426403 = r426401 - r426402;
double r426404 = log(r426401);
double r426405 = r426403 * r426404;
double r426406 = r426405 - r426401;
double r426407 = 0.91893853320467;
double r426408 = r426406 + r426407;
double r426409 = y;
double r426410 = 0.0007936500793651;
double r426411 = r426409 + r426410;
double r426412 = z;
double r426413 = r426411 * r426412;
double r426414 = 0.0027777777777778;
double r426415 = r426413 - r426414;
double r426416 = r426415 * r426412;
double r426417 = 0.083333333333333;
double r426418 = r426416 + r426417;
double r426419 = r426418 / r426401;
double r426420 = r426408 + r426419;
return r426420;
}
double f(double x, double y, double z) {
double r426421 = x;
double r426422 = sqrt(r426421);
double r426423 = log(r426422);
double r426424 = 0.5;
double r426425 = r426421 - r426424;
double r426426 = r426423 * r426425;
double r426427 = r426426 - r426421;
double r426428 = r426426 + r426427;
double r426429 = 0.91893853320467;
double r426430 = r426428 + r426429;
double r426431 = y;
double r426432 = 0.0007936500793651;
double r426433 = r426431 + r426432;
double r426434 = z;
double r426435 = r426433 * r426434;
double r426436 = 0.0027777777777778;
double r426437 = r426435 - r426436;
double r426438 = r426437 * r426434;
double r426439 = 0.083333333333333;
double r426440 = r426438 + r426439;
double r426441 = r426440 / r426421;
double r426442 = r426430 + r426441;
return r426442;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.0 |
|---|---|
| Target | 1.1 |
| Herbie | 6.0 |
Initial program 6.0
rmApplied add-sqr-sqrt6.0
Applied log-prod6.0
Applied distribute-rgt-in6.0
Applied associate--l+6.0
Final simplification6.0
herbie shell --seed 2020047
(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)))