
(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 (/ z (* x x)))
(t_3
(/
1.0
(+
(/
(+
(/ a x)
(-
(/
(+
(/ b x)
(- (/ (* z (- t_2 (/ a x))) x) (/ 27464.7644705 (* x x))))
y)
t_2))
y)
(/ 1.0 x))))
(t_4 (- (+ x (/ z y)) (* a (/ x y)))))
(if (<= y -3.9e+145)
t_4
(if (<= y -1.8e+56)
t_3
(if (<= y 1.4e+41)
(fma
y
(/ (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t_1)
(/ t t_1))
(if (<= y 8.2e+103) 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 = z / (x * x);
double t_3 = 1.0 / ((((a / x) + ((((b / x) + (((z * (t_2 - (a / x))) / x) - (27464.7644705 / (x * x)))) / y) - t_2)) / y) + (1.0 / x));
double t_4 = (x + (z / y)) - (a * (x / y));
double tmp;
if (y <= -3.9e+145) {
tmp = t_4;
} else if (y <= -1.8e+56) {
tmp = t_3;
} else if (y <= 1.4e+41) {
tmp = fma(y, (fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616) / t_1), (t / t_1));
} else if (y <= 8.2e+103) {
tmp = 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(z / Float64(x * x)) t_3 = Float64(1.0 / Float64(Float64(Float64(Float64(a / x) + Float64(Float64(Float64(Float64(b / x) + Float64(Float64(Float64(z * Float64(t_2 - Float64(a / x))) / x) - Float64(27464.7644705 / Float64(x * x)))) / y) - t_2)) / y) + Float64(1.0 / x))) t_4 = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))) tmp = 0.0 if (y <= -3.9e+145) tmp = t_4; elseif (y <= -1.8e+56) tmp = t_3; elseif (y <= 1.4e+41) tmp = fma(y, Float64(fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616) / t_1), Float64(t / t_1)); elseif (y <= 8.2e+103) tmp = 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[(z / N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(1.0 / N[(N[(N[(N[(a / x), $MachinePrecision] + N[(N[(N[(N[(b / x), $MachinePrecision] + N[(N[(N[(z * N[(t$95$2 - N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] - N[(27464.7644705 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision] - t$95$2), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision] + N[(1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.9e+145], t$95$4, If[LessEqual[y, -1.8e+56], t$95$3, If[LessEqual[y, 1.4e+41], N[(y * N[(N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.2e+103], t$95$3, 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 := \frac{z}{x \cdot x}\\
t_3 := \frac{1}{\frac{\frac{a}{x} + \left(\frac{\frac{b}{x} + \left(\frac{z \cdot \left(t\_2 - \frac{a}{x}\right)}{x} - \frac{27464.7644705}{x \cdot x}\right)}{y} - t\_2\right)}{y} + \frac{1}{x}}\\
t_4 := \left(x + \frac{z}{y}\right) - a \cdot \frac{x}{y}\\
\mathbf{if}\;y \leq -3.9 \cdot 10^{+145}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;y \leq -1.8 \cdot 10^{+56}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq 1.4 \cdot 10^{+41}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right)}{t\_1}, \frac{t}{t\_1}\right)\\
\mathbf{elif}\;y \leq 8.2 \cdot 10^{+103}:\\
\;\;\;\;t\_3\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
if y < -3.8999999999999998e145 or 8.2000000000000003e103 < y Initial program 0.0%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6480.2
Applied rewrites80.2%
if -3.8999999999999998e145 < y < -1.79999999999999999e56 or 1.4e41 < y < 8.2000000000000003e103Initial program 9.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f649.5
Applied rewrites9.5%
Taylor expanded in y around -inf
Applied rewrites72.0%
if -1.79999999999999999e56 < y < 1.4e41Initial program 95.8%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites95.9%
Final simplification88.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (+ y a) b))
(t_2 (/ t (fma y (fma y t_1 c) i)))
(t_3
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_3 -1e-124)
t_2
(if (<= t_3 5e-16)
(/ (fma y (fma y 27464.7644705 230661.510616) t) (fma (* y y) t_1 i))
(if (<= t_3 INFINITY) t_2 (- (+ x (/ z y)) (* a (/ x y))))))))
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 = t / fma(y, fma(y, t_1, c), i);
double t_3 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_3 <= -1e-124) {
tmp = t_2;
} else if (t_3 <= 5e-16) {
tmp = fma(y, fma(y, 27464.7644705, 230661.510616), t) / fma((y * y), t_1, i);
} else if (t_3 <= ((double) INFINITY)) {
tmp = t_2;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, Float64(y + a), b) t_2 = Float64(t / fma(y, fma(y, t_1, c), i)) t_3 = 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)))))))) tmp = 0.0 if (t_3 <= -1e-124) tmp = t_2; elseif (t_3 <= 5e-16) tmp = Float64(fma(y, fma(y, 27464.7644705, 230661.510616), t) / fma(Float64(y * y), t_1, i)); elseif (t_3 <= Inf) tmp = t_2; else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); 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[(t / N[(y * N[(y * t$95$1 + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = 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]}, If[LessEqual[t$95$3, -1e-124], t$95$2, If[LessEqual[t$95$3, 5e-16], N[(N[(y * N[(y * 27464.7644705 + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(N[(y * y), $MachinePrecision] * t$95$1 + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, Infinity], t$95$2, N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, y + a, b\right)\\
t_2 := \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, t\_1, c\right), i\right)}\\
t_3 := \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)}\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{-124}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{-16}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, 27464.7644705, 230661.510616\right), t\right)}{\mathsf{fma}\left(y \cdot y, t\_1, i\right)}\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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)) < -9.99999999999999933e-125 or 5.0000000000000004e-16 < (/.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 96.7%
Taylor expanded in t 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
lower-+.f6479.2
Applied rewrites79.2%
if -9.99999999999999933e-125 < (/.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)) < 5.0000000000000004e-16Initial program 82.7%
Taylor expanded in c 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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6475.5
Applied rewrites75.5%
Taylor expanded in y around 0
Applied rewrites61.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification70.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (+ y a) b))
(t_2 (/ t (fma y (fma y t_1 c) i)))
(t_3
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_3 -1e-124)
t_2
(if (<= t_3 5e-16)
(/ (fma y 230661.510616 t) (fma (* y y) t_1 i))
(if (<= t_3 INFINITY) t_2 (- (+ x (/ z y)) (* a (/ x y))))))))
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 = t / fma(y, fma(y, t_1, c), i);
double t_3 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_3 <= -1e-124) {
tmp = t_2;
} else if (t_3 <= 5e-16) {
tmp = fma(y, 230661.510616, t) / fma((y * y), t_1, i);
} else if (t_3 <= ((double) INFINITY)) {
tmp = t_2;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, Float64(y + a), b) t_2 = Float64(t / fma(y, fma(y, t_1, c), i)) t_3 = 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)))))))) tmp = 0.0 if (t_3 <= -1e-124) tmp = t_2; elseif (t_3 <= 5e-16) tmp = Float64(fma(y, 230661.510616, t) / fma(Float64(y * y), t_1, i)); elseif (t_3 <= Inf) tmp = t_2; else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); 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[(t / N[(y * N[(y * t$95$1 + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = 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]}, If[LessEqual[t$95$3, -1e-124], t$95$2, If[LessEqual[t$95$3, 5e-16], N[(N[(y * 230661.510616 + t), $MachinePrecision] / N[(N[(y * y), $MachinePrecision] * t$95$1 + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, Infinity], t$95$2, N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, y + a, b\right)\\
t_2 := \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, t\_1, c\right), i\right)}\\
t_3 := \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)}\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{-124}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{-16}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, 230661.510616, t\right)}{\mathsf{fma}\left(y \cdot y, t\_1, i\right)}\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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)) < -9.99999999999999933e-125 or 5.0000000000000004e-16 < (/.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 96.7%
Taylor expanded in t 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
lower-+.f6479.2
Applied rewrites79.2%
if -9.99999999999999933e-125 < (/.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)) < 5.0000000000000004e-16Initial program 82.7%
Taylor expanded in c 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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6475.5
Applied rewrites75.5%
Taylor expanded in y around 0
Applied rewrites59.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification69.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (fma y (fma y (fma y (+ y a) b) c) i)))
(t_2
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_2 -1e-124)
t_1
(if (<= t_2 1e-100)
(/ (fma y (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t) i)
(if (<= t_2 INFINITY) t_1 (- (+ x (/ z y)) (* a (/ x y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_2 <= -1e-124) {
tmp = t_1;
} else if (t_2 <= 1e-100) {
tmp = fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / i;
} else if (t_2 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)) t_2 = 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)))))))) tmp = 0.0 if (t_2 <= -1e-124) tmp = t_1; elseif (t_2 <= 1e-100) tmp = Float64(fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / i); elseif (t_2 <= Inf) tmp = t_1; else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = 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]}, If[LessEqual[t$95$2, -1e-124], t$95$1, If[LessEqual[t$95$2, 1e-100], N[(N[(y * N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[t$95$2, Infinity], t$95$1, N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
t_2 := \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)}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10^{-100}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right), t\right)}{i}\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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)) < -9.99999999999999933e-125 or 1e-100 < (/.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 97.0%
Taylor expanded in t 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
lower-+.f6476.8
Applied rewrites76.8%
if -9.99999999999999933e-125 < (/.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)) < 1e-100Initial program 79.6%
Taylor expanded in i 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.f6457.4
Applied rewrites57.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification69.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (fma y (fma y (fma y (+ y a) b) c) i)))
(t_2
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_2 -1e-124)
t_1
(if (<= t_2 1e-100)
(fma y (/ 230661.510616 i) (/ t i))
(if (<= t_2 INFINITY) t_1 (- (+ x (/ z y)) (* a (/ x y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_2 <= -1e-124) {
tmp = t_1;
} else if (t_2 <= 1e-100) {
tmp = fma(y, (230661.510616 / i), (t / i));
} else if (t_2 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)) t_2 = 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)))))))) tmp = 0.0 if (t_2 <= -1e-124) tmp = t_1; elseif (t_2 <= 1e-100) tmp = fma(y, Float64(230661.510616 / i), Float64(t / i)); elseif (t_2 <= Inf) tmp = t_1; else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = 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]}, If[LessEqual[t$95$2, -1e-124], t$95$1, If[LessEqual[t$95$2, 1e-100], N[(y * N[(230661.510616 / i), $MachinePrecision] + N[(t / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, Infinity], t$95$1, N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
t_2 := \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)}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10^{-100}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{230661.510616}{i}, \frac{t}{i}\right)\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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)) < -9.99999999999999933e-125 or 1e-100 < (/.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 97.0%
Taylor expanded in t 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
lower-+.f6476.8
Applied rewrites76.8%
if -9.99999999999999933e-125 < (/.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)) < 1e-100Initial program 79.6%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites79.7%
Taylor expanded in y around 0
Applied rewrites63.4%
Taylor expanded in y around 0
Applied rewrites53.0%
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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification67.9%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (fma y (fma y (fma y (+ y a) b) c) i)))
(t_2
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_2 -1e-124)
t_1
(if (<= t_2 1e-100)
(fma y (/ 230661.510616 i) (/ t i))
(if (<= t_2 INFINITY) 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 = t / fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_2 <= -1e-124) {
tmp = t_1;
} else if (t_2 <= 1e-100) {
tmp = fma(y, (230661.510616 / i), (t / i));
} else if (t_2 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)) t_2 = 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)))))))) tmp = 0.0 if (t_2 <= -1e-124) tmp = t_1; elseif (t_2 <= 1e-100) tmp = fma(y, Float64(230661.510616 / i), Float64(t / i)); elseif (t_2 <= Inf) tmp = 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[(t / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = 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]}, If[LessEqual[t$95$2, -1e-124], t$95$1, If[LessEqual[t$95$2, 1e-100], N[(y * N[(230661.510616 / i), $MachinePrecision] + N[(t / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, Infinity], t$95$1, N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
t_2 := \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)}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10^{-100}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{230661.510616}{i}, \frac{t}{i}\right)\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;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)) < -9.99999999999999933e-125 or 1e-100 < (/.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 97.0%
Taylor expanded in t 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
lower-+.f6476.8
Applied rewrites76.8%
if -9.99999999999999933e-125 < (/.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)) < 1e-100Initial program 79.6%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites79.7%
Taylor expanded in y around 0
Applied rewrites63.4%
Taylor expanded in y around 0
Applied rewrites53.0%
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 t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites0.6%
Taylor expanded in y around inf
associate--l+N/A
div-subN/A
lower-+.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6460.1
Applied rewrites60.1%
Final simplification65.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_1 1e-100)
(fma y (/ 230661.510616 i) (/ t i))
(if (<= t_1 INFINITY)
(/ t (fma (fma y (+ y a) b) (* y 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 t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_1 <= 1e-100) {
tmp = fma(y, (230661.510616 / i), (t / i));
} else if (t_1 <= ((double) INFINITY)) {
tmp = t / fma(fma(y, (y + a), b), (y * y), i);
} else {
tmp = x + ((z - (x * a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = 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)))))))) tmp = 0.0 if (t_1 <= 1e-100) tmp = fma(y, Float64(230661.510616 / i), Float64(t / i)); elseif (t_1 <= Inf) tmp = Float64(t / fma(fma(y, Float64(y + a), b), Float64(y * 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_] := Block[{t$95$1 = 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]}, If[LessEqual[t$95$1, 1e-100], N[(y * N[(230661.510616 / i), $MachinePrecision] + N[(t / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(t / N[(N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] * N[(y * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \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)}\\
\mathbf{if}\;t\_1 \leq 10^{-100}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{230661.510616}{i}, \frac{t}{i}\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{t}{\mathsf{fma}\left(\mathsf{fma}\left(y, y + a, b\right), y \cdot 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)) < 1e-100Initial program 87.8%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites87.9%
Taylor expanded in y around 0
Applied rewrites58.9%
Taylor expanded in y around 0
Applied rewrites50.2%
if 1e-100 < (/.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 99.5%
Taylor expanded in c 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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6481.1
Applied rewrites81.1%
Taylor expanded in t around inf
Applied rewrites64.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 t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites0.6%
Taylor expanded in y around inf
associate--l+N/A
div-subN/A
lower-+.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6460.1
Applied rewrites60.1%
Final simplification56.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1
(+
t
(*
y
(+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))))
(if (<= (/ t_1 (+ i (* y (+ c (* y (+ b (* y (+ y a)))))))) INFINITY)
(/ t_1 (+ i (/ y (/ 1.0 (fma y (fma y (+ y a) b) c)))))
(- (+ x (/ z y)) (* a (/ x y))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))));
double tmp;
if ((t_1 / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = t_1 / (i + (y / (1.0 / fma(y, fma(y, (y + a), b), c))));
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) tmp = 0.0 if (Float64(t_1 / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(t_1 / Float64(i + Float64(y / Float64(1.0 / fma(y, fma(y, Float64(y + a), b), c))))); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = 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]}, If[LessEqual[N[(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]), $MachinePrecision], Infinity], N[(t$95$1 / N[(i + N[(y / N[(1.0 / N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)\\
\mathbf{if}\;\frac{t\_1}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t\_1}{i + \frac{y}{\frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right)}}}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-+.f64N/A
lower-/.f6491.0
Applied rewrites91.0%
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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification82.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1
(+
t
(*
y
(+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))))
(if (<= (/ t_1 (+ i (* y (+ c (* y (+ b (* y (+ y a)))))))) INFINITY)
(/ t_1 (+ i (fma y c (* y (* y (fma y (+ y a) b))))))
(- (+ x (/ z y)) (* a (/ x y))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))));
double tmp;
if ((t_1 / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = t_1 / (i + fma(y, c, (y * (y * fma(y, (y + a), b)))));
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) tmp = 0.0 if (Float64(t_1 / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(t_1 / Float64(i + fma(y, c, Float64(y * Float64(y * fma(y, Float64(y + a), b)))))); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = 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]}, If[LessEqual[N[(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]), $MachinePrecision], Infinity], N[(t$95$1 / N[(i + N[(y * c + N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)\\
\mathbf{if}\;\frac{t\_1}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t\_1}{i + \mathsf{fma}\left(y, c, y \cdot \left(y \cdot \mathsf{fma}\left(y, y + a, b\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f6490.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6490.9
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6490.9
Applied rewrites90.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification82.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))))
(if (<= t_1 INFINITY) t_1 (- (+ x (/ z y)) (* a (/ x y))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = (x + (z / y)) - (a * (x / y));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a))))))) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = (x + (z / y)) - (a * (x / y)) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = 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)))))))) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a))))))); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = (x + (z / y)) - (a * (x / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = 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]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \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)}\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification82.5%
(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 (fma x y z) 27464.7644705) 230661.510616) t)
(/ 1.0 (fma y (fma y (fma y (+ y a) b) c) i)))
(- (+ x (/ z y)) (* a (/ x 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, fma(x, y, z), 27464.7644705), 230661.510616), t) * (1.0 / fma(y, fma(y, fma(y, (y + a), b), c), i));
} else {
tmp = (x + (z / y)) - (a * (x / 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, fma(x, y, z), 27464.7644705), 230661.510616), t) * Float64(1.0 / fma(y, fma(y, fma(y, Float64(y + a), b), c), i))); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / 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 * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] * N[(1.0 / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $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(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right), t\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}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
Applied rewrites90.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification82.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 (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y y b) c) i))
(- (+ x (/ z y)) (* a (/ x 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, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, y, b), c), i);
} else {
tmp = (x + (z / y)) - (a * (x / 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, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, y, b), c), i)); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / 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 * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * y + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $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, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
Taylor expanded in a 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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6482.8
Applied rewrites82.8%
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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification77.2%
(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 y (fma y 27464.7644705 230661.510616) t) t_1)
(- (+ x (/ z y)) (* a (/ x 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(y, fma(y, 27464.7644705, 230661.510616), t) / t_1;
} else {
tmp = (x + (z / y)) - (a * (x / 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(y, fma(y, 27464.7644705, 230661.510616), t) / t_1); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / 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[(y * N[(y * 27464.7644705 + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / t$95$1), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $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(y, \mathsf{fma}\left(y, 27464.7644705, 230661.510616\right), t\right)}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6479.0
Applied rewrites79.0%
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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification74.8%
(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 y 230661.510616 t) t_1)
(- (+ x (/ z y)) (* a (/ x 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(y, 230661.510616, t) / t_1;
} else {
tmp = (x + (z / y)) - (a * (x / 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(y, 230661.510616, t) / t_1); else tmp = Float64(Float64(x + Float64(z / y)) - Float64(a * Float64(x / 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[(y * 230661.510616 + t), $MachinePrecision] / t$95$1), $MachinePrecision], N[(N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $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(y, 230661.510616, t\right)}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{z}{y}\right) - a \cdot \frac{x}{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 90.9%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6477.9
Applied rewrites77.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%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
+-commutativeN/A
distribute-lft-inN/A
lower-fma.f64N/A
lower-*.f640.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f640.0
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f640.0
Applied rewrites0.0%
Taylor expanded in y around inf
lower--.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.9
Applied rewrites66.9%
Final simplification74.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -3.5e+53)
t_1
(if (<= y 14.6)
(fma y (/ 230661.510616 i) (/ t i))
(if (<= y 2.4e+79) (/ (* y x) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.5e+53) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = fma(y, (230661.510616 / i), (t / i));
} else if (y <= 2.4e+79) {
tmp = (y * x) / a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -3.5e+53) tmp = t_1; elseif (y <= 14.6) tmp = fma(y, Float64(230661.510616 / i), Float64(t / i)); elseif (y <= 2.4e+79) tmp = Float64(Float64(y * x) / a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.5e+53], t$95$1, If[LessEqual[y, 14.6], N[(y * N[(230661.510616 / i), $MachinePrecision] + N[(t / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e+79], N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -3.5 \cdot 10^{+53}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 14.6:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{230661.510616}{i}, \frac{t}{i}\right)\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{+79}:\\
\;\;\;\;\frac{y \cdot x}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.50000000000000019e53 or 2.39999999999999986e79 < y Initial program 1.4%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites1.9%
Taylor expanded in y around inf
associate--l+N/A
div-subN/A
lower-+.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6462.7
Applied rewrites62.7%
if -3.50000000000000019e53 < y < 14.5999999999999996Initial program 96.9%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites97.0%
Taylor expanded in y around 0
Applied rewrites66.0%
Taylor expanded in y around 0
Applied rewrites56.5%
if 14.5999999999999996 < y < 2.39999999999999986e79Initial program 42.1%
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-*.f6414.3
Applied rewrites14.3%
Taylor expanded in y around inf
Applied rewrites33.2%
Final simplification56.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -3.5e+53)
t_1
(if (<= y 14.6)
(fma 230661.510616 (/ y i) (/ t i))
(if (<= y 2.4e+79) (/ (* y x) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.5e+53) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = fma(230661.510616, (y / i), (t / i));
} else if (y <= 2.4e+79) {
tmp = (y * x) / a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -3.5e+53) tmp = t_1; elseif (y <= 14.6) tmp = fma(230661.510616, Float64(y / i), Float64(t / i)); elseif (y <= 2.4e+79) tmp = Float64(Float64(y * x) / a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.5e+53], t$95$1, If[LessEqual[y, 14.6], N[(230661.510616 * N[(y / i), $MachinePrecision] + N[(t / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e+79], N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -3.5 \cdot 10^{+53}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 14.6:\\
\;\;\;\;\mathsf{fma}\left(230661.510616, \frac{y}{i}, \frac{t}{i}\right)\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{+79}:\\
\;\;\;\;\frac{y \cdot x}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.50000000000000019e53 or 2.39999999999999986e79 < y Initial program 1.4%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites1.9%
Taylor expanded in y around inf
associate--l+N/A
div-subN/A
lower-+.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6462.7
Applied rewrites62.7%
if -3.50000000000000019e53 < y < 14.5999999999999996Initial program 96.9%
Taylor expanded in c 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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6475.8
Applied rewrites75.8%
Taylor expanded in y around 0
Applied rewrites56.4%
if 14.5999999999999996 < y < 2.39999999999999986e79Initial program 42.1%
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-*.f6414.3
Applied rewrites14.3%
Taylor expanded in y around inf
Applied rewrites33.2%
Final simplification56.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -4e+53)
t_1
(if (<= y 14.6) (/ t i) (if (<= y 2.4e+79) (/ (* y x) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -4e+53) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = t / i;
} else if (y <= 2.4e+79) {
tmp = (y * x) / a;
} 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 + ((z - (x * a)) / y)
if (y <= (-4d+53)) then
tmp = t_1
else if (y <= 14.6d0) then
tmp = t / i
else if (y <= 2.4d+79) then
tmp = (y * x) / a
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 + ((z - (x * a)) / y);
double tmp;
if (y <= -4e+53) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = t / i;
} else if (y <= 2.4e+79) {
tmp = (y * x) / a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z - (x * a)) / y) tmp = 0 if y <= -4e+53: tmp = t_1 elif y <= 14.6: tmp = t / i elif y <= 2.4e+79: tmp = (y * x) / a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -4e+53) tmp = t_1; elseif (y <= 14.6) tmp = Float64(t / i); elseif (y <= 2.4e+79) tmp = Float64(Float64(y * x) / a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z - (x * a)) / y); tmp = 0.0; if (y <= -4e+53) tmp = t_1; elseif (y <= 14.6) tmp = t / i; elseif (y <= 2.4e+79) tmp = (y * x) / a; 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[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+53], t$95$1, If[LessEqual[y, 14.6], N[(t / i), $MachinePrecision], If[LessEqual[y, 2.4e+79], N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+53}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 14.6:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{+79}:\\
\;\;\;\;\frac{y \cdot x}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -4e53 or 2.39999999999999986e79 < y Initial program 1.4%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites1.9%
Taylor expanded in y around inf
associate--l+N/A
div-subN/A
lower-+.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6462.7
Applied rewrites62.7%
if -4e53 < y < 14.5999999999999996Initial program 96.9%
Taylor expanded in y around 0
lower-/.f6446.5
Applied rewrites46.5%
if 14.5999999999999996 < y < 2.39999999999999986e79Initial program 42.1%
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-*.f6414.3
Applied rewrites14.3%
Taylor expanded in y around inf
Applied rewrites33.2%
Final simplification50.8%
(FPCore (x y z t a b c i) :precision binary64 (let* ((t_1 (/ (* y x) a))) (if (<= y -7e+16) t_1 (if (<= y 14.6) (/ t i) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (y * x) / a;
double tmp;
if (y <= -7e+16) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = t / 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 = (y * x) / a
if (y <= (-7d+16)) then
tmp = t_1
else if (y <= 14.6d0) then
tmp = t / 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 = (y * x) / a;
double tmp;
if (y <= -7e+16) {
tmp = t_1;
} else if (y <= 14.6) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * x) / a tmp = 0 if y <= -7e+16: tmp = t_1 elif y <= 14.6: tmp = t / i else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * x) / a) tmp = 0.0 if (y <= -7e+16) tmp = t_1; elseif (y <= 14.6) tmp = Float64(t / i); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (y * x) / a; tmp = 0.0; if (y <= -7e+16) tmp = t_1; elseif (y <= 14.6) tmp = t / i; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(y * x), $MachinePrecision] / a), $MachinePrecision]}, If[LessEqual[y, -7e+16], t$95$1, If[LessEqual[y, 14.6], N[(t / i), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y \cdot x}{a}\\
\mathbf{if}\;y \leq -7 \cdot 10^{+16}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 14.6:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -7e16 or 14.5999999999999996 < y Initial program 12.1%
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-*.f645.6
Applied rewrites5.6%
Taylor expanded in y around inf
Applied rewrites16.5%
if -7e16 < y < 14.5999999999999996Initial program 99.6%
Taylor expanded in y around 0
lower-/.f6448.6
Applied rewrites48.6%
Final simplification33.7%
(FPCore (x y z t a b c i) :precision binary64 (if (<= i -1e-224) (/ t i) (if (<= i 4e-64) (/ t (* y c)) (/ t i))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (i <= -1e-224) {
tmp = t / i;
} else if (i <= 4e-64) {
tmp = t / (y * c);
} else {
tmp = t / i;
}
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 (i <= (-1d-224)) then
tmp = t / i
else if (i <= 4d-64) then
tmp = t / (y * c)
else
tmp = t / i
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 (i <= -1e-224) {
tmp = t / i;
} else if (i <= 4e-64) {
tmp = t / (y * c);
} else {
tmp = t / i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if i <= -1e-224: tmp = t / i elif i <= 4e-64: tmp = t / (y * c) else: tmp = t / i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (i <= -1e-224) tmp = Float64(t / i); elseif (i <= 4e-64) tmp = Float64(t / Float64(y * c)); else tmp = Float64(t / i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (i <= -1e-224) tmp = t / i; elseif (i <= 4e-64) tmp = t / (y * c); else tmp = t / i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[i, -1e-224], N[(t / i), $MachinePrecision], If[LessEqual[i, 4e-64], N[(t / N[(y * c), $MachinePrecision]), $MachinePrecision], N[(t / i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -1 \cdot 10^{-224}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;i \leq 4 \cdot 10^{-64}:\\
\;\;\;\;\frac{t}{y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{i}\\
\end{array}
\end{array}
if i < -1e-224 or 3.99999999999999986e-64 < i Initial program 55.0%
Taylor expanded in y around 0
lower-/.f6431.0
Applied rewrites31.0%
if -1e-224 < i < 3.99999999999999986e-64Initial program 73.0%
Taylor expanded in c 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-*.f6440.4
Applied rewrites40.4%
Taylor expanded in y around 0
Applied rewrites34.5%
Final simplification31.8%
(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 59.0%
Taylor expanded in y around 0
lower-/.f6427.6
Applied rewrites27.6%
(FPCore (x y z t a b c i) :precision binary64 (/ z a))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return z / a;
}
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 = z / a
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return z / a;
}
def code(x, y, z, t, a, b, c, i): return z / a
function code(x, y, z, t, a, b, c, i) return Float64(z / a) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = z / a; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(z / a), $MachinePrecision]
\begin{array}{l}
\\
\frac{z}{a}
\end{array}
Initial program 59.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-*.f646.8
Applied rewrites6.8%
Taylor expanded in z around inf
Applied rewrites5.9%
herbie shell --seed 2024219
(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)))