\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\begin{array}{l}
\mathbf{if}\;x \le -4.113171276926967275782769070286056694835 \cdot 10^{70} \lor \neg \left(x \le 1.783473391079485331997464252300643745944 \cdot 10^{47}\right):\\
\;\;\;\;\left(\frac{y}{{x}^{2}} + 4.16438922227999963610045597306452691555 \cdot x\right) - 110.1139242984810948655649553984403610229\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x, 4.16438922227999963610045597306452691555, 78.69949241540000173245061887428164482117\right), x, 137.5194164160000127594685181975364685059\right), x, y\right), x, z\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x + 43.3400022514000013984514225739985704422, x, 263.5050747210000281484099105000495910645\right), x, 313.3992158940000081202015280723571777344\right), x, 47.06687660600000100430406746454536914825\right)} \cdot \left(x - 2\right)\\
\end{array}double f(double x, double y, double z) {
double r411372 = x;
double r411373 = 2.0;
double r411374 = r411372 - r411373;
double r411375 = 4.16438922228;
double r411376 = r411372 * r411375;
double r411377 = 78.6994924154;
double r411378 = r411376 + r411377;
double r411379 = r411378 * r411372;
double r411380 = 137.519416416;
double r411381 = r411379 + r411380;
double r411382 = r411381 * r411372;
double r411383 = y;
double r411384 = r411382 + r411383;
double r411385 = r411384 * r411372;
double r411386 = z;
double r411387 = r411385 + r411386;
double r411388 = r411374 * r411387;
double r411389 = 43.3400022514;
double r411390 = r411372 + r411389;
double r411391 = r411390 * r411372;
double r411392 = 263.505074721;
double r411393 = r411391 + r411392;
double r411394 = r411393 * r411372;
double r411395 = 313.399215894;
double r411396 = r411394 + r411395;
double r411397 = r411396 * r411372;
double r411398 = 47.066876606;
double r411399 = r411397 + r411398;
double r411400 = r411388 / r411399;
return r411400;
}
double f(double x, double y, double z) {
double r411401 = x;
double r411402 = -4.113171276926967e+70;
bool r411403 = r411401 <= r411402;
double r411404 = 1.7834733910794853e+47;
bool r411405 = r411401 <= r411404;
double r411406 = !r411405;
bool r411407 = r411403 || r411406;
double r411408 = y;
double r411409 = 2.0;
double r411410 = pow(r411401, r411409);
double r411411 = r411408 / r411410;
double r411412 = 4.16438922228;
double r411413 = r411412 * r411401;
double r411414 = r411411 + r411413;
double r411415 = 110.1139242984811;
double r411416 = r411414 - r411415;
double r411417 = 78.6994924154;
double r411418 = fma(r411401, r411412, r411417);
double r411419 = 137.519416416;
double r411420 = fma(r411418, r411401, r411419);
double r411421 = fma(r411420, r411401, r411408);
double r411422 = z;
double r411423 = fma(r411421, r411401, r411422);
double r411424 = 43.3400022514;
double r411425 = r411401 + r411424;
double r411426 = 263.505074721;
double r411427 = fma(r411425, r411401, r411426);
double r411428 = 313.399215894;
double r411429 = fma(r411427, r411401, r411428);
double r411430 = 47.066876606;
double r411431 = fma(r411429, r411401, r411430);
double r411432 = r411423 / r411431;
double r411433 = 2.0;
double r411434 = r411401 - r411433;
double r411435 = r411432 * r411434;
double r411436 = r411407 ? r411416 : r411435;
return r411436;
}




Bits error versus x




Bits error versus y




Bits error versus z
| Original | 27.1 |
|---|---|
| Target | 0.6 |
| Herbie | 0.7 |
if x < -4.113171276926967e+70 or 1.7834733910794853e+47 < x Initial program 62.6
Simplified59.5
rmApplied *-un-lft-identity59.5
Applied add-sqr-sqrt59.5
Applied times-frac59.5
Applied associate-/r*59.5
Simplified59.5
Taylor expanded around inf 0.5
if -4.113171276926967e+70 < x < 1.7834733910794853e+47Initial program 2.1
Simplified0.9
rmApplied clear-num0.9
rmApplied div-inv0.9
Applied add-sqr-sqrt0.9
Applied times-frac0.9
Simplified0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2019353 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C"
:precision binary64
:herbie-target
(if (< x -3.326128725870005e+62) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811) (if (< x 9.429991714554673e+55) (* (/ (- x 2) 1) (/ (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z) (+ (* (+ (+ (* 263.505074721 x) (+ (* 43.3400022514 (* x x)) (* x (* x x)))) 313.399215894) x) 47.066876606))) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811)))
(/ (* (- x 2) (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z)) (+ (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x) 47.066876606)))