
(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)))
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);
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static 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);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\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}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(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)))
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);
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static 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);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\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}
\end{array}
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (fma y (fma y (+ y a) b) c) i))
(t_2 (+ x (/ (- z (* x a)) y)))
(t_3 (/ t t_1))
(t_4
(fma y (/ z (* y y)) (fma x (/ (* y y) (+ b (* y (+ y a)))) t_3))))
(if (<= y -1.32e+154)
t_2
(if (<= y -1.28e+51)
t_4
(if (<= y 2.7e+59)
(fma
(/ (fma y (fma y z 27464.7644705) 230661.510616) t_1)
y
(fma x (/ (* y (* y (* y y))) t_1) t_3))
(if (<= y 1.25e+154) t_4 t_2))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = x + ((z - (x * a)) / y);
double t_3 = t / t_1;
double t_4 = fma(y, (z / (y * y)), fma(x, ((y * y) / (b + (y * (y + a)))), t_3));
double tmp;
if (y <= -1.32e+154) {
tmp = t_2;
} else if (y <= -1.28e+51) {
tmp = t_4;
} else if (y <= 2.7e+59) {
tmp = fma((fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), y, fma(x, ((y * (y * (y * y))) / t_1), t_3));
} else if (y <= 1.25e+154) {
tmp = t_4;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, fma(y, fma(y, Float64(y + a), b), c), i) t_2 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) t_3 = Float64(t / t_1) t_4 = fma(y, Float64(z / Float64(y * y)), fma(x, Float64(Float64(y * y) / Float64(b + Float64(y * Float64(y + a)))), t_3)) tmp = 0.0 if (y <= -1.32e+154) tmp = t_2; elseif (y <= -1.28e+51) tmp = t_4; elseif (y <= 2.7e+59) tmp = fma(Float64(fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), y, fma(x, Float64(Float64(y * Float64(y * Float64(y * y))) / t_1), t_3)); elseif (y <= 1.25e+154) tmp = t_4; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t / t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(y * N[(z / N[(y * y), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(y * y), $MachinePrecision] / N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.32e+154], t$95$2, If[LessEqual[y, -1.28e+51], t$95$4, If[LessEqual[y, 2.7e+59], N[(N[(N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] * y + N[(x * N[(N[(y * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+154], t$95$4, t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)\\
t_2 := x + \frac{z - x \cdot a}{y}\\
t_3 := \frac{t}{t\_1}\\
t_4 := \mathsf{fma}\left(y, \frac{z}{y \cdot y}, \mathsf{fma}\left(x, \frac{y \cdot y}{b + y \cdot \left(y + a\right)}, t\_3\right)\right)\\
\mathbf{if}\;y \leq -1.32 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.28 \cdot 10^{+51}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;y \leq 2.7 \cdot 10^{+59}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right)}{t\_1}, y, \mathsf{fma}\left(x, \frac{y \cdot \left(y \cdot \left(y \cdot y\right)\right)}{t\_1}, t\_3\right)\right)\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+154}:\\
\;\;\;\;t\_4\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -1.31999999999999998e154 or 1.25000000000000001e154 < y Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6476.2
Applied rewrites76.2%
if -1.31999999999999998e154 < y < -1.27999999999999993e51 or 2.7000000000000001e59 < y < 1.25000000000000001e154Initial program 5.0%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites5.0%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6417.7
Applied rewrites17.7%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6457.0
Applied rewrites57.0%
Taylor expanded in y around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6477.5
Applied rewrites77.5%
if -1.27999999999999993e51 < y < 2.7000000000000001e59Initial program 93.1%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites93.8%
Applied rewrites93.8%
Final simplification86.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ i (* y (+ c (* y (+ b (* y (+ y a)))))))))
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
t_1)
INFINITY)
(/
(+
t
(/ y (/ 1.0 (fma y (fma y (fma x y z) 27464.7644705) 230661.510616))))
t_1)
(+ x (/ (- z (* x a)) y)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * (b + (y * (y + a))))));
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / t_1) <= ((double) INFINITY)) {
tmp = (t + (y / (1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / t_1;
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a))))))) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / t_1) <= Inf) tmp = Float64(Float64(t + Float64(y / Float64(1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / t_1); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], Infinity], N[(N[(t + N[(y / N[(1.0 / N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)\\
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{t\_1} \leq \infty:\\
\;\;\;\;\frac{t + \frac{y}{\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right)}}}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
Applied rewrites91.5%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification80.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ i (* y (+ c (* y (+ b (* y (+ y a)))))))))
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
t_1)
INFINITY)
(/
(fma (fma y (fma y x z) 27464.7644705) (* y y) (fma y 230661.510616 t))
t_1)
(+ x (/ (- z (* x a)) y)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * (b + (y * (y + a))))));
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / t_1) <= ((double) INFINITY)) {
tmp = fma(fma(y, fma(y, x, z), 27464.7644705), (y * y), fma(y, 230661.510616, t)) / t_1;
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a))))))) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / t_1) <= Inf) tmp = Float64(fma(fma(y, fma(y, x, z), 27464.7644705), Float64(y * y), fma(y, 230661.510616, t)) / t_1); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], Infinity], N[(N[(N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] * N[(y * y), $MachinePrecision] + N[(y * 230661.510616 + t), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)\\
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{t\_1} \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), y \cdot y, \mathsf{fma}\left(y, 230661.510616, t\right)\right)}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
flip-+N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites80.9%
Applied rewrites91.5%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification80.3%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(*
(fma (fma y (fma y x z) 27464.7644705) (* y y) (fma y 230661.510616 t))
(/ 1.0 (fma y (fma y (fma y (+ y a) b) c) i)))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = fma(fma(y, fma(y, x, z), 27464.7644705), (y * y), fma(y, 230661.510616, t)) * (1.0 / fma(y, fma(y, fma(y, (y + a), b), c), i));
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(fma(fma(y, fma(y, x, z), 27464.7644705), Float64(y * y), fma(y, 230661.510616, t)) * Float64(1.0 / fma(y, fma(y, fma(y, Float64(y + a), b), c), i))); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] * N[(y * y), $MachinePrecision] + N[(y * 230661.510616 + t), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), y \cdot y, \mathsf{fma}\left(y, 230661.510616, t\right)\right) \cdot \frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
flip-+N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites80.9%
Applied rewrites91.2%
lift-fma.f64N/A
lift-fma.f64N/A
lift-fma.f64N/A
lift-fma.f64N/A
distribute-rgt-inN/A
*-commutativeN/A
associate-+l+N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
lift-fma.f64N/A
lower-fma.f6491.3
Applied rewrites91.3%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification80.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (fma y (fma y (+ y a) b) c) i))
(t_2 (+ b (* y (+ y a))))
(t_3 (fma x (/ (* y y) t_2) (/ t t_1)))
(t_4 (+ x (/ (- z (* x a)) y))))
(if (<= y -1.32e+154)
t_4
(if (<= y -4.2e+27)
(fma y (/ z (* y y)) t_3)
(if (<= y 140000000.0)
(/
(+
t
(/
y
(/ 1.0 (fma y (fma y (fma x y z) 27464.7644705) 230661.510616))))
(+ i (* y (+ c (* y t_2)))))
(if (<= y 5e+138)
(fma y (/ (fma y (fma y z 27464.7644705) 230661.510616) t_1) t_3)
t_4))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = b + (y * (y + a));
double t_3 = fma(x, ((y * y) / t_2), (t / t_1));
double t_4 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -1.32e+154) {
tmp = t_4;
} else if (y <= -4.2e+27) {
tmp = fma(y, (z / (y * y)), t_3);
} else if (y <= 140000000.0) {
tmp = (t + (y / (1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / (i + (y * (c + (y * t_2))));
} else if (y <= 5e+138) {
tmp = fma(y, (fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), t_3);
} else {
tmp = t_4;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, fma(y, fma(y, Float64(y + a), b), c), i) t_2 = Float64(b + Float64(y * Float64(y + a))) t_3 = fma(x, Float64(Float64(y * y) / t_2), Float64(t / t_1)) t_4 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -1.32e+154) tmp = t_4; elseif (y <= -4.2e+27) tmp = fma(y, Float64(z / Float64(y * y)), t_3); elseif (y <= 140000000.0) tmp = Float64(Float64(t + Float64(y / Float64(1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / Float64(i + Float64(y * Float64(c + Float64(y * t_2))))); elseif (y <= 5e+138) tmp = fma(y, Float64(fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), t_3); else tmp = t_4; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$2 = N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x * N[(N[(y * y), $MachinePrecision] / t$95$2), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.32e+154], t$95$4, If[LessEqual[y, -4.2e+27], N[(y * N[(z / N[(y * y), $MachinePrecision]), $MachinePrecision] + t$95$3), $MachinePrecision], If[LessEqual[y, 140000000.0], N[(N[(t + N[(y / N[(1.0 / N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e+138], N[(y * N[(N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + t$95$3), $MachinePrecision], t$95$4]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)\\
t_2 := b + y \cdot \left(y + a\right)\\
t_3 := \mathsf{fma}\left(x, \frac{y \cdot y}{t\_2}, \frac{t}{t\_1}\right)\\
t_4 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -1.32 \cdot 10^{+154}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;y \leq -4.2 \cdot 10^{+27}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{z}{y \cdot y}, t\_3\right)\\
\mathbf{elif}\;y \leq 140000000:\\
\;\;\;\;\frac{t + \frac{y}{\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right)}}}{i + y \cdot \left(c + y \cdot t\_2\right)}\\
\mathbf{elif}\;y \leq 5 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right)}{t\_1}, t\_3\right)\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
if y < -1.31999999999999998e154 or 5.00000000000000016e138 < y Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6476.2
Applied rewrites76.2%
if -1.31999999999999998e154 < y < -4.19999999999999989e27Initial program 0.9%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites4.0%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f644.4
Applied rewrites4.4%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6442.2
Applied rewrites42.2%
Taylor expanded in y around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6463.1
Applied rewrites63.1%
if -4.19999999999999989e27 < y < 1.4e8Initial program 98.2%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
Applied rewrites98.2%
if 1.4e8 < y < 5.00000000000000016e138Initial program 36.7%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites40.6%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6454.5
Applied rewrites54.5%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6478.0
Applied rewrites78.0%
Final simplification86.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (+ y a) b))
(t_2 (fma y (fma y t_1 c) i))
(t_3 (+ b (* y (+ y a))))
(t_4 (+ x (/ (- z (* x a)) y)))
(t_5 (/ t t_2)))
(if (<= y -1.32e+154)
t_4
(if (<= y -4.2e+27)
(fma y (/ z (* y y)) (fma x (/ (* y y) t_3) t_5))
(if (<= y 140000000.0)
(/
(+
t
(/
y
(/ 1.0 (fma y (fma y (fma x y z) 27464.7644705) 230661.510616))))
(+ i (* y (+ c (* y t_3)))))
(if (<= y 5e+138)
(fma
(/ (fma y (fma y z 27464.7644705) 230661.510616) t_2)
y
(fma y (* x (/ y t_1)) t_5))
t_4))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, (y + a), b);
double t_2 = fma(y, fma(y, t_1, c), i);
double t_3 = b + (y * (y + a));
double t_4 = x + ((z - (x * a)) / y);
double t_5 = t / t_2;
double tmp;
if (y <= -1.32e+154) {
tmp = t_4;
} else if (y <= -4.2e+27) {
tmp = fma(y, (z / (y * y)), fma(x, ((y * y) / t_3), t_5));
} else if (y <= 140000000.0) {
tmp = (t + (y / (1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / (i + (y * (c + (y * t_3))));
} else if (y <= 5e+138) {
tmp = fma((fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_2), y, fma(y, (x * (y / t_1)), t_5));
} else {
tmp = t_4;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, Float64(y + a), b) t_2 = fma(y, fma(y, t_1, c), i) t_3 = Float64(b + Float64(y * Float64(y + a))) t_4 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) t_5 = Float64(t / t_2) tmp = 0.0 if (y <= -1.32e+154) tmp = t_4; elseif (y <= -4.2e+27) tmp = fma(y, Float64(z / Float64(y * y)), fma(x, Float64(Float64(y * y) / t_3), t_5)); elseif (y <= 140000000.0) tmp = Float64(Float64(t + Float64(y / Float64(1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / Float64(i + Float64(y * Float64(c + Float64(y * t_3))))); elseif (y <= 5e+138) tmp = fma(Float64(fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_2), y, fma(y, Float64(x * Float64(y / t_1)), t_5)); else tmp = t_4; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$2 = N[(y * N[(y * t$95$1 + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$3 = N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(t / t$95$2), $MachinePrecision]}, If[LessEqual[y, -1.32e+154], t$95$4, If[LessEqual[y, -4.2e+27], N[(y * N[(z / N[(y * y), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(y * y), $MachinePrecision] / t$95$3), $MachinePrecision] + t$95$5), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 140000000.0], N[(N[(t + N[(y / N[(1.0 / N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e+138], N[(N[(N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$2), $MachinePrecision] * y + N[(y * N[(x * N[(y / t$95$1), $MachinePrecision]), $MachinePrecision] + t$95$5), $MachinePrecision]), $MachinePrecision], t$95$4]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, y + a, b\right)\\
t_2 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, t\_1, c\right), i\right)\\
t_3 := b + y \cdot \left(y + a\right)\\
t_4 := x + \frac{z - x \cdot a}{y}\\
t_5 := \frac{t}{t\_2}\\
\mathbf{if}\;y \leq -1.32 \cdot 10^{+154}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;y \leq -4.2 \cdot 10^{+27}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{z}{y \cdot y}, \mathsf{fma}\left(x, \frac{y \cdot y}{t\_3}, t\_5\right)\right)\\
\mathbf{elif}\;y \leq 140000000:\\
\;\;\;\;\frac{t + \frac{y}{\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right)}}}{i + y \cdot \left(c + y \cdot t\_3\right)}\\
\mathbf{elif}\;y \leq 5 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right)}{t\_2}, y, \mathsf{fma}\left(y, x \cdot \frac{y}{t\_1}, t\_5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
if y < -1.31999999999999998e154 or 5.00000000000000016e138 < y Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6476.2
Applied rewrites76.2%
if -1.31999999999999998e154 < y < -4.19999999999999989e27Initial program 0.9%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites4.0%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f644.4
Applied rewrites4.4%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6442.2
Applied rewrites42.2%
Taylor expanded in y around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6463.1
Applied rewrites63.1%
if -4.19999999999999989e27 < y < 1.4e8Initial program 98.2%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
Applied rewrites98.2%
if 1.4e8 < y < 5.00000000000000016e138Initial program 36.7%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites40.6%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6454.5
Applied rewrites54.5%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6478.0
Applied rewrites78.0%
Applied rewrites77.8%
Final simplification86.1%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(*
(/ 1.0 (fma y (fma y (fma y (+ y a) b) c) i))
(fma y (fma y (fma y (fma x y z) 27464.7644705) 230661.510616) t))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = (1.0 / fma(y, fma(y, fma(y, (y + a), b), c), i)) * fma(y, fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616), t);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(Float64(1.0 / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)) * fma(y, fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616), t)); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(1.0 / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)} \cdot \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right), t\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Applied rewrites91.2%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification80.2%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6485.3
Applied rewrites85.3%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification76.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(*
(/ 1.0 (fma y (fma y (fma y (+ y a) b) c) i))
(fma y (fma y 27464.7644705 230661.510616) t))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = (1.0 / fma(y, fma(y, fma(y, (y + a), b), c), i)) * fma(y, fma(y, 27464.7644705, 230661.510616), t);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(Float64(1.0 / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)) * fma(y, fma(y, 27464.7644705, 230661.510616), t)); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(1.0 / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * 27464.7644705 + 230661.510616), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)} \cdot \mathsf{fma}\left(y, \mathsf{fma}\left(y, 27464.7644705, 230661.510616\right), t\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
flip-+N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites80.9%
Applied rewrites91.2%
Taylor expanded in y around 0
Applied rewrites76.9%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification71.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (fma y 230661.510616 t) (fma (fma y (fma y (+ y a) b) c) y i))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = fma(y, 230661.510616, t) / fma(fma(y, fma(y, (y + a), b), c), y, i);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(fma(y, 230661.510616, t) / fma(fma(y, fma(y, Float64(y + a), b), c), y, i)); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y * 230661.510616 + t), $MachinePrecision] / N[(N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, 230661.510616, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-+.f64N/A
*-commutativeN/A
lift-*.f64N/A
+-commutativeN/A
lift-+.f64N/A
lower-fma.f6476.3
Applied rewrites76.3%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification71.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ b (* y (+ y a))))
(t_2
(fma
y
(/ z (* y y))
(fma x (/ (* y y) t_1) (/ t (fma y (fma y (fma y (+ y a) b) c) i)))))
(t_3 (+ x (/ (- z (* x a)) y))))
(if (<= y -1.32e+154)
t_3
(if (<= y -4.2e+27)
t_2
(if (<= y 3e+62)
(/
(+
t
(/
y
(/ 1.0 (fma y (fma y (fma x y z) 27464.7644705) 230661.510616))))
(+ i (* y (+ c (* y t_1)))))
(if (<= y 1.25e+154) t_2 t_3))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = b + (y * (y + a));
double t_2 = fma(y, (z / (y * y)), fma(x, ((y * y) / t_1), (t / fma(y, fma(y, fma(y, (y + a), b), c), i))));
double t_3 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -1.32e+154) {
tmp = t_3;
} else if (y <= -4.2e+27) {
tmp = t_2;
} else if (y <= 3e+62) {
tmp = (t + (y / (1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / (i + (y * (c + (y * t_1))));
} else if (y <= 1.25e+154) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(b + Float64(y * Float64(y + a))) t_2 = fma(y, Float64(z / Float64(y * y)), fma(x, Float64(Float64(y * y) / t_1), Float64(t / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)))) t_3 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -1.32e+154) tmp = t_3; elseif (y <= -4.2e+27) tmp = t_2; elseif (y <= 3e+62) tmp = Float64(Float64(t + Float64(y / Float64(1.0 / fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616)))) / Float64(i + Float64(y * Float64(c + Float64(y * t_1))))); elseif (y <= 1.25e+154) tmp = t_2; else tmp = t_3; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(y * N[(z / N[(y * y), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(y * y), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.32e+154], t$95$3, If[LessEqual[y, -4.2e+27], t$95$2, If[LessEqual[y, 3e+62], N[(N[(t + N[(y / N[(1.0 / N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+154], t$95$2, t$95$3]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b + y \cdot \left(y + a\right)\\
t_2 := \mathsf{fma}\left(y, \frac{z}{y \cdot y}, \mathsf{fma}\left(x, \frac{y \cdot y}{t\_1}, \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\right)\right)\\
t_3 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -1.32 \cdot 10^{+154}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq -4.2 \cdot 10^{+27}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 3 \cdot 10^{+62}:\\
\;\;\;\;\frac{t + \frac{y}{\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right)}}}{i + y \cdot \left(c + y \cdot t\_1\right)}\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if y < -1.31999999999999998e154 or 1.25000000000000001e154 < y Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6476.2
Applied rewrites76.2%
if -1.31999999999999998e154 < y < -4.19999999999999989e27 or 3e62 < y < 1.25000000000000001e154Initial program 2.8%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites5.0%
Taylor expanded in i around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6414.9
Applied rewrites14.9%
Taylor expanded in c around 0
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f6454.6
Applied rewrites54.6%
Taylor expanded in y around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6472.1
Applied rewrites72.1%
if -4.19999999999999989e27 < y < 3e62Initial program 95.0%
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
Applied rewrites95.0%
Final simplification86.1%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (fma y 230661.510616 t) (fma (fma y b c) y i))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = fma(y, 230661.510616, t) / fma(fma(y, b, c), y, i);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(fma(y, 230661.510616, t) / fma(fma(y, b, c), y, i)); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y * 230661.510616 + t), $MachinePrecision] / N[(N[(y * b + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, 230661.510616, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(y, b, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
Taylor expanded in y around 0
lower-*.f6472.7
Applied rewrites72.7%
lift-*.f64N/A
lift-+.f64N/A
lower-fma.f6472.7
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6472.7
Applied rewrites72.7%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification69.3%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (fma y 230661.510616 t) (+ i (* y c)))
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = fma(y, 230661.510616, t) / (i + (y * c));
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(fma(y, 230661.510616, t) / Float64(i + Float64(y * c))); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y * 230661.510616 + t), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, 230661.510616, t\right)}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
Taylor expanded in y around 0
lower-*.f6467.4
Applied rewrites67.4%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification66.1%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (+ t (* y 230661.510616)) i)
(+ x (/ (- z (* x a)) y))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Double.POSITIVE_INFINITY) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if ((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= math.inf: tmp = (t + (y * 230661.510616)) / i else: tmp = x + ((z - (x * a)) / y) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); else tmp = Float64(x + Float64(Float64(z - Float64(x * a)) / y)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Inf) tmp = (t + (y * 230661.510616)) / i; else tmp = x + ((z - (x * a)) / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z - x \cdot a}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
Taylor expanded in i around inf
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6444.9
Applied rewrites44.9%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6464.3
Applied rewrites64.3%
Final simplification52.9%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (+ t (* y 230661.510616)) i)
(* y (/ x a))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = y * (x / a);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Double.POSITIVE_INFINITY) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = y * (x / a);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if ((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= math.inf: tmp = (t + (y * 230661.510616)) / i else: tmp = y * (x / a) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); else tmp = Float64(y * Float64(x / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Inf) tmp = (t + (y * 230661.510616)) / i; else tmp = y * (x / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(y * N[(x / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{x}{a}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
Taylor expanded in i around inf
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6444.9
Applied rewrites44.9%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites0.1%
Taylor expanded in a around inf
lower-/.f6410.6
Applied rewrites10.6%
Taylor expanded in y around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-*.f64N/A
associate-+r+N/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
Applied rewrites15.2%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f6419.3
Applied rewrites19.3%
Final simplification34.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ (+ t (* y 230661.510616)) i)
(/ (* y x) a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = (y * x) / a;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Double.POSITIVE_INFINITY) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = (y * x) / a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if ((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= math.inf: tmp = (t + (y * 230661.510616)) / i else: tmp = (y * x) / a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); else tmp = Float64(Float64(y * x) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Inf) tmp = (t + (y * 230661.510616)) / i; else tmp = (y * x) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot x}{a}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6476.3
Applied rewrites76.3%
Taylor expanded in i around inf
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6444.9
Applied rewrites44.9%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in a around inf
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f640.3
Applied rewrites0.3%
Taylor expanded in y around inf
lower-/.f64N/A
lower-*.f6410.1
Applied rewrites10.1%
Final simplification30.6%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ t i)
(/ (* y x) a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = t / i;
} else {
tmp = (y * x) / a;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Double.POSITIVE_INFINITY) {
tmp = t / i;
} else {
tmp = (y * x) / a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if ((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= math.inf: tmp = t / i else: tmp = (y * x) / a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(t / i); else tmp = Float64(Float64(y * x) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Inf) tmp = t / i; else tmp = (y * x) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(t / i), $MachinePrecision], N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot x}{a}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 91.4%
Taylor expanded in y around 0
lower-/.f6440.4
Applied rewrites40.4%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in a around inf
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f640.3
Applied rewrites0.3%
Taylor expanded in y around inf
lower-/.f64N/A
lower-*.f6410.1
Applied rewrites10.1%
Final simplification28.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (- x (/ (* x a) y))))
(if (<= y -4.4e+109)
t_1
(if (<= y -4e-64)
(* y (/ x a))
(if (<= y 0.065) (/ (+ t (* y 230661.510616)) i) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x - ((x * a) / y);
double tmp;
if (y <= -4.4e+109) {
tmp = t_1;
} else if (y <= -4e-64) {
tmp = y * (x / a);
} else if (y <= 0.065) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((x * a) / y)
if (y <= (-4.4d+109)) then
tmp = t_1
else if (y <= (-4d-64)) then
tmp = y * (x / a)
else if (y <= 0.065d0) then
tmp = (t + (y * 230661.510616d0)) / i
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x - ((x * a) / y);
double tmp;
if (y <= -4.4e+109) {
tmp = t_1;
} else if (y <= -4e-64) {
tmp = y * (x / a);
} else if (y <= 0.065) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x - ((x * a) / y) tmp = 0 if y <= -4.4e+109: tmp = t_1 elif y <= -4e-64: tmp = y * (x / a) elif y <= 0.065: tmp = (t + (y * 230661.510616)) / i else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x - Float64(Float64(x * a) / y)) tmp = 0.0 if (y <= -4.4e+109) tmp = t_1; elseif (y <= -4e-64) tmp = Float64(y * Float64(x / a)); elseif (y <= 0.065) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x - ((x * a) / y); tmp = 0.0; if (y <= -4.4e+109) tmp = t_1; elseif (y <= -4e-64) tmp = y * (x / a); elseif (y <= 0.065) tmp = (t + (y * 230661.510616)) / i; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4.4e+109], t$95$1, If[LessEqual[y, -4e-64], N[(y * N[(x / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.065], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{x \cdot a}{y}\\
\mathbf{if}\;y \leq -4.4 \cdot 10^{+109}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -4 \cdot 10^{-64}:\\
\;\;\;\;y \cdot \frac{x}{a}\\
\mathbf{elif}\;y \leq 0.065:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -4.3999999999999998e109 or 0.065000000000000002 < y Initial program 10.5%
Taylor expanded in x around inf
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f644.3
Applied rewrites4.3%
Taylor expanded in y around inf
mul-1-negN/A
lower-+.f64N/A
lower-neg.f64N/A
lower-/.f64N/A
lower-*.f6453.1
Applied rewrites53.1%
if -4.3999999999999998e109 < y < -3.99999999999999986e-64Initial program 45.7%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites45.8%
Taylor expanded in a around inf
lower-/.f6423.9
Applied rewrites23.9%
Taylor expanded in y around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-*.f64N/A
associate-+r+N/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
Applied rewrites14.8%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f6420.0
Applied rewrites20.0%
if -3.99999999999999986e-64 < y < 0.065000000000000002Initial program 99.7%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6495.6
Applied rewrites95.6%
Taylor expanded in i around inf
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6461.4
Applied rewrites61.4%
Final simplification51.1%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 1.7e-29) (/ t i) (* z (/ 1.0 a))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.7e-29) {
tmp = t / i;
} else {
tmp = z * (1.0 / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 1.7d-29) then
tmp = t / i
else
tmp = z * (1.0d0 / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.7e-29) {
tmp = t / i;
} else {
tmp = z * (1.0 / a);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 1.7e-29: tmp = t / i else: tmp = z * (1.0 / a) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 1.7e-29) tmp = Float64(t / i); else tmp = Float64(z * Float64(1.0 / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 1.7e-29) tmp = t / i; else tmp = z * (1.0 / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 1.7e-29], N[(t / i), $MachinePrecision], N[(z * N[(1.0 / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.7 \cdot 10^{-29}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{1}{a}\\
\end{array}
\end{array}
if y < 1.69999999999999986e-29Initial program 66.8%
Taylor expanded in y around 0
lower-/.f6433.1
Applied rewrites33.1%
if 1.69999999999999986e-29 < y Initial program 21.0%
Taylor expanded in a around inf
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f647.4
Applied rewrites7.4%
Taylor expanded in z around inf
lower-/.f6412.3
Applied rewrites12.3%
clear-numN/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6412.3
Applied rewrites12.3%
Final simplification27.2%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 1.7e-29) (/ t i) (/ z a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.7e-29) {
tmp = t / i;
} else {
tmp = z / a;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 1.7d-29) then
tmp = t / i
else
tmp = z / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.7e-29) {
tmp = t / i;
} else {
tmp = z / a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 1.7e-29: tmp = t / i else: tmp = z / a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 1.7e-29) tmp = Float64(t / i); else tmp = Float64(z / a); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 1.7e-29) tmp = t / i; else tmp = z / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 1.7e-29], N[(t / i), $MachinePrecision], N[(z / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.7 \cdot 10^{-29}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{z}{a}\\
\end{array}
\end{array}
if y < 1.69999999999999986e-29Initial program 66.8%
Taylor expanded in y around 0
lower-/.f6433.1
Applied rewrites33.1%
if 1.69999999999999986e-29 < y Initial program 21.0%
Taylor expanded in a around inf
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f647.4
Applied rewrites7.4%
Taylor expanded in z around inf
lower-/.f6412.3
Applied rewrites12.3%
(FPCore (x y z t a b c i) :precision binary64 (/ t i))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return t / i;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = t / i
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return t / i;
}
def code(x, y, z, t, a, b, c, i): return t / i
function code(x, y, z, t, a, b, c, i) return Float64(t / i) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = t / i; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(t / i), $MachinePrecision]
\begin{array}{l}
\\
\frac{t}{i}
\end{array}
Initial program 53.9%
Taylor expanded in y around 0
lower-/.f6425.2
Applied rewrites25.2%
herbie shell --seed 2024212
(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)))