\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\begin{array}{l}
\mathbf{if}\;\frac{y \cdot \left(y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\right) + t}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i} \leq -8.574197889873893 \cdot 10^{-40}:\\
\;\;\;\;230661.510616 \cdot \frac{y}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(\frac{t}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(27464.7644705 \cdot \frac{y \cdot y}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(x \cdot \frac{{y}^{4}}{\left(y + a\right) \cdot {y}^{3} + \left(i + y \cdot \left(c + y \cdot b\right)\right)} + z \cdot \frac{{y}^{3}}{\left(y + a\right) \cdot {y}^{3} + \left(i + y \cdot \left(c + y \cdot b\right)\right)}\right)\right)\right)\\
\mathbf{elif}\;\frac{y \cdot \left(y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\right) + t}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i} \leq 2.3219598672624144 \cdot 10^{+220}:\\
\;\;\;\;\frac{y \cdot \left(y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\right) + t}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{elif}\;\frac{y \cdot \left(y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\right) + t}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i} \leq \infty:\\
\;\;\;\;230661.510616 \cdot \frac{y}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(\frac{t}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(27464.7644705 \cdot \frac{y \cdot y}{{y}^{4} + \left(a \cdot {y}^{3} + \left(b \cdot \left(y \cdot y\right) + \left(i + y \cdot c\right)\right)\right)} + \left(x \cdot \frac{{y}^{4}}{\left(y + a\right) \cdot {y}^{3} + \left(i + y \cdot \left(c + y \cdot b\right)\right)} + z \cdot \frac{{y}^{3}}{\left(y + a\right) \cdot {y}^{3} + \left(i + y \cdot \left(c + y \cdot b\right)\right)}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}(FPCore (x y z t a b c i) :precision binary64 (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
(* y (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
t)
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
-8.574197889873893e-40)
(+
(*
230661.510616
(/
y
(+ (pow y 4.0) (+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c)))))))
(+
(/
t
(+ (pow y 4.0) (+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c))))))
(+
(*
27464.7644705
(/
(* y y)
(+ (pow y 4.0) (+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c)))))))
(+
(*
x
(/ (pow y 4.0) (+ (* (+ y a) (pow y 3.0)) (+ i (* y (+ c (* y b)))))))
(*
z
(/
(pow y 3.0)
(+ (* (+ y a) (pow y 3.0)) (+ i (* y (+ c (* y b)))))))))))
(if (<=
(/
(+
(* y (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
t)
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
2.3219598672624144e+220)
(/
(+ (* y (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)) t)
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
(if (<=
(/
(+
(* y (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
t)
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
INFINITY)
(+
(*
230661.510616
(/
y
(+
(pow y 4.0)
(+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c)))))))
(+
(/
t
(+
(pow y 4.0)
(+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c))))))
(+
(*
27464.7644705
(/
(* y y)
(+
(pow y 4.0)
(+ (* a (pow y 3.0)) (+ (* b (* y y)) (+ i (* y c)))))))
(+
(*
x
(/
(pow y 4.0)
(+ (* (+ y a) (pow y 3.0)) (+ i (* y (+ c (* y b)))))))
(*
z
(/
(pow y 3.0)
(+ (* (+ y a) (pow y 3.0)) (+ i (* y (+ c (* y b)))))))))))
x))))double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((((y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)) + t) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)) <= -8.574197889873893e-40) {
tmp = (230661.510616 * (y / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c))))))) + ((t / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c)))))) + ((27464.7644705 * ((y * y) / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c))))))) + ((x * (pow(y, 4.0) / (((y + a) * pow(y, 3.0)) + (i + (y * (c + (y * b))))))) + (z * (pow(y, 3.0) / (((y + a) * pow(y, 3.0)) + (i + (y * (c + (y * b))))))))));
} else if ((((y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)) + t) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)) <= 2.3219598672624144e+220) {
tmp = ((y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)) + t) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if ((((y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)) + t) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)) <= ((double) INFINITY)) {
tmp = (230661.510616 * (y / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c))))))) + ((t / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c)))))) + ((27464.7644705 * ((y * y) / (pow(y, 4.0) + ((a * pow(y, 3.0)) + ((b * (y * y)) + (i + (y * c))))))) + ((x * (pow(y, 4.0) / (((y + a) * pow(y, 3.0)) + (i + (y * (c + (y * b))))))) + (z * (pow(y, 3.0) / (((y + a) * pow(y, 3.0)) + (i + (y * (c + (y * b))))))))));
} else {
tmp = x;
}
return tmp;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus i
Results
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) 471841060772561/17179869184) y) 7925469156333415/34359738368) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < -8.5741978898738933e-40 or 2.32195986726241441e220 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) 471841060772561/17179869184) y) 7925469156333415/34359738368) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 9.8
Taylor expanded around 0 11.1
Simplified11.1
rmApplied *-un-lft-identity_binary64_110111.1
Applied times-frac_binary64_11076.3
Simplified6.3
Simplified6.3
rmApplied *-un-lft-identity_binary64_11016.3
Applied times-frac_binary64_11071.6
Simplified1.6
Simplified1.6
if -8.5741978898738933e-40 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) 471841060772561/17179869184) y) 7925469156333415/34359738368) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < 2.32195986726241441e220Initial program 5.3
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) 471841060772561/17179869184) y) 7925469156333415/34359738368) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 64.0
Taylor expanded around inf 27.7
Final simplification13.6
herbie shell --seed 2020344
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
:precision binary64
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))