
(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 20 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
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* x y))))))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))))
(if (<= t_1 INFINITY) t_1 (+ x (/ z 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 + (x * y)))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = x + (z / 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 + (x * y)))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = x + (z / y);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = x + (z / 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(x * y)))))))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(x + Float64(z / 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 + (x * y)))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = x + (z / 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[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(x + N[(z / 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 + x \cdot y\right)\right)\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z}{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.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%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr0.0%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6470.6%
Simplified70.6%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6476.2%
Simplified76.2%
Final simplification84.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -2.6e+35)
t_1
(if (<= y 1.6e-132)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* x (* y y)))))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
(if (<= y 3.6e+62)
(*
(/ 1.0 (+ i (* y (+ c (* (+ y a) (* y y))))))
(+
t
(*
y
(+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* x y)))))))))
(if (<= y 2.2e+99) (* z (/ y b)) 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 / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 1.6e-132) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 3.6e+62) {
tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y))))))));
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} 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 / y)
if (y <= (-2.6d+35)) then
tmp = t_1
else if (y <= 1.6d-132) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)
else if (y <= 3.6d+62) then
tmp = (1.0d0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * (z + (x * y))))))))
else if (y <= 2.2d+99) then
tmp = z * (y / b)
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 / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 1.6e-132) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 3.6e+62) {
tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y))))))));
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -2.6e+35: tmp = t_1 elif y <= 1.6e-132: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) elif y <= 3.6e+62: tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) elif y <= 2.2e+99: tmp = z * (y / b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -2.6e+35) tmp = t_1; elseif (y <= 1.6e-132) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(x * Float64(y * y))))))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); elseif (y <= 3.6e+62) tmp = Float64(Float64(1.0 / Float64(i + Float64(y * Float64(c + Float64(Float64(y + a) * Float64(y * y)))))) * Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(x * y))))))))); elseif (y <= 2.2e+99) tmp = Float64(z * Float64(y / b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -2.6e+35) tmp = t_1; elseif (y <= 1.6e-132) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i); elseif (y <= 3.6e+62) tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))); elseif (y <= 2.2e+99) tmp = z * (y / b); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.6e+35], t$95$1, If[LessEqual[y, 1.6e-132], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.6e+62], N[(N[(1.0 / N[(i + N[(y * N[(c + N[(N[(y + a), $MachinePrecision] * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -2.6 \cdot 10^{+35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.6 \cdot 10^{-132}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + x \cdot \left(y \cdot y\right)\right)\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+62}:\\
\;\;\;\;\frac{1}{i + y \cdot \left(c + \left(y + a\right) \cdot \left(y \cdot y\right)\right)} \cdot \left(t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + x \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.60000000000000007e35 or 2.19999999999999978e99 < y Initial program 5.0%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.0%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6472.5%
Simplified72.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6478.0%
Simplified78.0%
if -2.60000000000000007e35 < y < 1.6000000000000001e-132Initial program 99.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6496.7%
Simplified96.7%
if 1.6000000000000001e-132 < y < 3.6e62Initial program 88.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr88.2%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6480.7%
Simplified80.7%
if 3.6e62 < y < 2.19999999999999978e99Initial program 0.7%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f640.7%
Simplified0.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr3.6%
Taylor expanded in b around inf
/-lowering-/.f6451.6%
Simplified51.6%
Final simplification84.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -2.6e+35)
t_1
(if (<= y 2e-134)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* x (* y y)))))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
(if (<= y 3.6e+62)
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* x y))))))))
(+ i (* y (+ c (* (+ y a) (* y y))))))
(if (<= y 2.2e+99) (* z (/ y b)) 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 / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 2e-134) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 3.6e+62) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) / (i + (y * (c + ((y + a) * (y * y)))));
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} 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 / y)
if (y <= (-2.6d+35)) then
tmp = t_1
else if (y <= 2d-134) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)
else if (y <= 3.6d+62) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * (z + (x * y)))))))) / (i + (y * (c + ((y + a) * (y * y)))))
else if (y <= 2.2d+99) then
tmp = z * (y / b)
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 / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 2e-134) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 3.6e+62) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) / (i + (y * (c + ((y + a) * (y * y)))));
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -2.6e+35: tmp = t_1 elif y <= 2e-134: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) elif y <= 3.6e+62: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) / (i + (y * (c + ((y + a) * (y * y))))) elif y <= 2.2e+99: tmp = z * (y / b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -2.6e+35) tmp = t_1; elseif (y <= 2e-134) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(x * Float64(y * y))))))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); elseif (y <= 3.6e+62) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(x * y)))))))) / Float64(i + Float64(y * Float64(c + Float64(Float64(y + a) * Float64(y * y)))))); elseif (y <= 2.2e+99) tmp = Float64(z * Float64(y / b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -2.6e+35) tmp = t_1; elseif (y <= 2e-134) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i); elseif (y <= 3.6e+62) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (x * y)))))))) / (i + (y * (c + ((y + a) * (y * y))))); elseif (y <= 2.2e+99) tmp = z * (y / b); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.6e+35], t$95$1, If[LessEqual[y, 2e-134], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.6e+62], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(N[(y + a), $MachinePrecision] * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -2.6 \cdot 10^{+35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2 \cdot 10^{-134}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + x \cdot \left(y \cdot y\right)\right)\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+62}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + x \cdot y\right)\right)\right)}{i + y \cdot \left(c + \left(y + a\right) \cdot \left(y \cdot y\right)\right)}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.60000000000000007e35 or 2.19999999999999978e99 < y Initial program 5.0%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.0%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6472.5%
Simplified72.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6478.0%
Simplified78.0%
if -2.60000000000000007e35 < y < 2.00000000000000008e-134Initial program 99.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6496.7%
Simplified96.7%
if 2.00000000000000008e-134 < y < 3.6e62Initial program 88.1%
Taylor expanded in b around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
Simplified80.5%
if 3.6e62 < y < 2.19999999999999978e99Initial program 0.7%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f640.7%
Simplified0.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr3.6%
Taylor expanded in b around inf
/-lowering-/.f6451.6%
Simplified51.6%
Final simplification84.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -1.3e+35)
t_1
(if (<= y 3.6e+62)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) i))
(if (<= y 2.2e+99) (* z (/ y b)) 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 / y);
double tmp;
if (y <= -1.3e+35) {
tmp = t_1;
} else if (y <= 3.6e+62) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} 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 / y)
if (y <= (-1.3d+35)) then
tmp = t_1
else if (y <= 3.6d+62) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i)
else if (y <= 2.2d+99) then
tmp = z * (y / b)
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 / y);
double tmp;
if (y <= -1.3e+35) {
tmp = t_1;
} else if (y <= 3.6e+62) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -1.3e+35: tmp = t_1 elif y <= 3.6e+62: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) elif y <= 2.2e+99: tmp = z * (y / b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -1.3e+35) tmp = t_1; elseif (y <= 3.6e+62) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); elseif (y <= 2.2e+99) tmp = Float64(z * Float64(y / b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -1.3e+35) tmp = t_1; elseif (y <= 3.6e+62) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i); elseif (y <= 2.2e+99) tmp = z * (y / b); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.3e+35], t$95$1, If[LessEqual[y, 3.6e+62], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -1.3 \cdot 10^{+35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+62}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.30000000000000003e35 or 2.19999999999999978e99 < y Initial program 5.0%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.0%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6472.5%
Simplified72.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6478.0%
Simplified78.0%
if -1.30000000000000003e35 < y < 3.6e62Initial program 96.0%
Taylor expanded in x around 0
*-lowering-*.f6488.6%
Simplified88.6%
if 3.6e62 < y < 2.19999999999999978e99Initial program 0.7%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f640.7%
Simplified0.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr3.6%
Taylor expanded in b around inf
/-lowering-/.f6451.6%
Simplified51.6%
Final simplification82.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -2.6e+35)
t_1
(if (<= y 3.4e+41)
(/
(+ t (* y (+ 230661.510616 (* x (* y (* y y))))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) 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 + (z / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 3.4e+41) {
tmp = (t + (y * (230661.510616 + (x * (y * (y * y)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y)
if (y <= (-2.6d+35)) then
tmp = t_1
else if (y <= 3.4d+41) then
tmp = (t + (y * (230661.510616d0 + (x * (y * (y * y)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y);
double tmp;
if (y <= -2.6e+35) {
tmp = t_1;
} else if (y <= 3.4e+41) {
tmp = (t + (y * (230661.510616 + (x * (y * (y * y)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -2.6e+35: tmp = t_1 elif y <= 3.4e+41: tmp = (t + (y * (230661.510616 + (x * (y * (y * y)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -2.6e+35) tmp = t_1; elseif (y <= 3.4e+41) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(x * Float64(y * Float64(y * y)))))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -2.6e+35) tmp = t_1; elseif (y <= 3.4e+41) tmp = (t + (y * (230661.510616 + (x * (y * (y * y)))))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.6e+35], t$95$1, If[LessEqual[y, 3.4e+41], N[(N[(t + N[(y * N[(230661.510616 + N[(x * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -2.6 \cdot 10^{+35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.4 \cdot 10^{+41}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + x \cdot \left(y \cdot \left(y \cdot y\right)\right)\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.60000000000000007e35 or 3.39999999999999998e41 < y Initial program 6.2%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr6.2%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.0%
Simplified71.0%
if -2.60000000000000007e35 < y < 3.39999999999999998e41Initial program 96.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.8%
Simplified88.8%
Final simplification80.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -3.9e+34)
t_1
(if (<= y 1.38e+41)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* x (* y y)))))))
(+ i (* y (+ c (* y b)))))
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 / y);
double tmp;
if (y <= -3.9e+34) {
tmp = t_1;
} else if (y <= 1.38e+41) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / (i + (y * (c + (y * b))));
} 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 / y)
if (y <= (-3.9d+34)) then
tmp = t_1
else if (y <= 1.38d+41) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (x * (y * y))))))) / (i + (y * (c + (y * b))))
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 / y);
double tmp;
if (y <= -3.9e+34) {
tmp = t_1;
} else if (y <= 1.38e+41) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / (i + (y * (c + (y * b))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -3.9e+34: tmp = t_1 elif y <= 1.38e+41: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / (i + (y * (c + (y * b)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -3.9e+34) tmp = t_1; elseif (y <= 1.38e+41) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(x * Float64(y * y))))))) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -3.9e+34) tmp = t_1; elseif (y <= 1.38e+41) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (x * (y * y))))))) / (i + (y * (c + (y * b)))); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.9e+34], t$95$1, If[LessEqual[y, 1.38e+41], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -3.9 \cdot 10^{+34}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.38 \cdot 10^{+41}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + x \cdot \left(y \cdot y\right)\right)\right)}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.90000000000000019e34 or 1.3800000000000001e41 < y Initial program 6.2%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr6.2%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.0%
Simplified71.0%
if -3.90000000000000019e34 < y < 1.3800000000000001e41Initial program 96.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6489.6%
Simplified89.6%
Taylor expanded in y around 0
*-lowering-*.f6484.4%
Simplified84.4%
Final simplification78.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -5.2e+32)
t_1
(if (<= y 1.06e+47)
(/
(+ t (* y (+ 230661.510616 (* y 27464.7644705))))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) 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 + (z / y);
double tmp;
if (y <= -5.2e+32) {
tmp = t_1;
} else if (y <= 1.06e+47) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y)
if (y <= (-5.2d+32)) then
tmp = t_1
else if (y <= 1.06d+47) then
tmp = (t + (y * (230661.510616d0 + (y * 27464.7644705d0)))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y);
double tmp;
if (y <= -5.2e+32) {
tmp = t_1;
} else if (y <= 1.06e+47) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -5.2e+32: tmp = t_1 elif y <= 1.06e+47: tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -5.2e+32) tmp = t_1; elseif (y <= 1.06e+47) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * 27464.7644705)))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -5.2e+32) tmp = t_1; elseif (y <= 1.06e+47) tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.2e+32], t$95$1, If[LessEqual[y, 1.06e+47], N[(N[(t + N[(y * N[(230661.510616 + N[(y * 27464.7644705), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -5.2 \cdot 10^{+32}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.06 \cdot 10^{+47}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot 27464.7644705\right)}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.2000000000000004e32 or 1.05999999999999996e47 < y Initial program 5.4%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.4%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.6%
Simplified66.6%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.5%
Simplified71.5%
if -5.2000000000000004e32 < y < 1.05999999999999996e47Initial program 96.7%
Taylor expanded in y around 0
*-commutativeN/A
*-lowering-*.f6482.7%
Simplified82.7%
Final simplification77.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -2.5e+34)
t_1
(if (<= y 9e+47)
(/
(+ t (* y 230661.510616))
(+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) 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 + (z / y);
double tmp;
if (y <= -2.5e+34) {
tmp = t_1;
} else if (y <= 9e+47) {
tmp = (t + (y * 230661.510616)) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y)
if (y <= (-2.5d+34)) then
tmp = t_1
else if (y <= 9d+47) then
tmp = (t + (y * 230661.510616d0)) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y);
double tmp;
if (y <= -2.5e+34) {
tmp = t_1;
} else if (y <= 9e+47) {
tmp = (t + (y * 230661.510616)) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -2.5e+34: tmp = t_1 elif y <= 9e+47: tmp = (t + (y * 230661.510616)) / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -2.5e+34) tmp = t_1; elseif (y <= 9e+47) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -2.5e+34) tmp = t_1; elseif (y <= 9e+47) tmp = (t + (y * 230661.510616)) / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.5e+34], t$95$1, If[LessEqual[y, 9e+47], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -2.5 \cdot 10^{+34}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 9 \cdot 10^{+47}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.4999999999999999e34 or 8.99999999999999958e47 < y Initial program 5.4%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.4%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.6%
Simplified66.6%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.5%
Simplified71.5%
if -2.4999999999999999e34 < y < 8.99999999999999958e47Initial program 96.7%
Taylor expanded in y around 0
*-commutativeN/A
*-lowering-*.f6481.2%
Simplified81.2%
Final simplification76.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -1.5e+35)
t_1
(if (<= y 2.8e+44)
(*
(/ 1.0 (+ i (* y (+ c (* (+ y a) (* y y))))))
(+ t (* y 230661.510616)))
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 / y);
double tmp;
if (y <= -1.5e+35) {
tmp = t_1;
} else if (y <= 2.8e+44) {
tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * 230661.510616));
} 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 / y)
if (y <= (-1.5d+35)) then
tmp = t_1
else if (y <= 2.8d+44) then
tmp = (1.0d0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * 230661.510616d0))
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 / y);
double tmp;
if (y <= -1.5e+35) {
tmp = t_1;
} else if (y <= 2.8e+44) {
tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * 230661.510616));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -1.5e+35: tmp = t_1 elif y <= 2.8e+44: tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * 230661.510616)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -1.5e+35) tmp = t_1; elseif (y <= 2.8e+44) tmp = Float64(Float64(1.0 / Float64(i + Float64(y * Float64(c + Float64(Float64(y + a) * Float64(y * y)))))) * Float64(t + Float64(y * 230661.510616))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -1.5e+35) tmp = t_1; elseif (y <= 2.8e+44) tmp = (1.0 / (i + (y * (c + ((y + a) * (y * y)))))) * (t + (y * 230661.510616)); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.5e+35], t$95$1, If[LessEqual[y, 2.8e+44], N[(N[(1.0 / N[(i + N[(y * N[(c + N[(N[(y + a), $MachinePrecision] * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -1.5 \cdot 10^{+35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.8 \cdot 10^{+44}:\\
\;\;\;\;\frac{1}{i + y \cdot \left(c + \left(y + a\right) \cdot \left(y \cdot y\right)\right)} \cdot \left(t + y \cdot 230661.510616\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.49999999999999995e35 or 2.8000000000000001e44 < y Initial program 5.4%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr5.4%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.6%
Simplified66.6%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.5%
Simplified71.5%
if -1.49999999999999995e35 < y < 2.8000000000000001e44Initial program 96.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr96.6%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6486.4%
Simplified86.4%
Taylor expanded in y around 0
*-lowering-*.f6473.3%
Simplified73.3%
Final simplification72.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -9.5e+33)
t_1
(if (<= y 3.5e+41)
(/ t (+ (* y (+ (* y (+ (* y (+ y a)) b)) c)) 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 + (z / y);
double tmp;
if (y <= -9.5e+33) {
tmp = t_1;
} else if (y <= 3.5e+41) {
tmp = t / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y)
if (y <= (-9.5d+33)) then
tmp = t_1
else if (y <= 3.5d+41) then
tmp = t / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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 + (z / y);
double tmp;
if (y <= -9.5e+33) {
tmp = t_1;
} else if (y <= 3.5e+41) {
tmp = t / ((y * ((y * ((y * (y + a)) + b)) + c)) + i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -9.5e+33: tmp = t_1 elif y <= 3.5e+41: tmp = t / ((y * ((y * ((y * (y + a)) + b)) + c)) + i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -9.5e+33) tmp = t_1; elseif (y <= 3.5e+41) tmp = Float64(t / Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(y + a)) + b)) + c)) + i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -9.5e+33) tmp = t_1; elseif (y <= 3.5e+41) tmp = t / ((y * ((y * ((y * (y + a)) + b)) + c)) + 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.5e+33], t$95$1, If[LessEqual[y, 3.5e+41], N[(t / N[(N[(y * N[(N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -9.5 \cdot 10^{+33}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{+41}:\\
\;\;\;\;\frac{t}{y \cdot \left(y \cdot \left(y \cdot \left(y + a\right) + b\right) + c\right) + i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -9.5000000000000003e33 or 3.4999999999999999e41 < y Initial program 6.2%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr6.2%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6471.0%
Simplified71.0%
if -9.5000000000000003e33 < y < 3.4999999999999999e41Initial program 96.7%
Taylor expanded in y around 0
Simplified68.9%
Final simplification69.9%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -6e-14)
t_1
(if (<= y 1.4e+40)
(* t (/ -1.0 (- (- 0.0 i) (* y (+ c (* y b))))))
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 / y);
double tmp;
if (y <= -6e-14) {
tmp = t_1;
} else if (y <= 1.4e+40) {
tmp = t * (-1.0 / ((0.0 - i) - (y * (c + (y * b)))));
} 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 / y)
if (y <= (-6d-14)) then
tmp = t_1
else if (y <= 1.4d+40) then
tmp = t * ((-1.0d0) / ((0.0d0 - i) - (y * (c + (y * b)))))
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 / y);
double tmp;
if (y <= -6e-14) {
tmp = t_1;
} else if (y <= 1.4e+40) {
tmp = t * (-1.0 / ((0.0 - i) - (y * (c + (y * b)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -6e-14: tmp = t_1 elif y <= 1.4e+40: tmp = t * (-1.0 / ((0.0 - i) - (y * (c + (y * b))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -6e-14) tmp = t_1; elseif (y <= 1.4e+40) tmp = Float64(t * Float64(-1.0 / Float64(Float64(0.0 - i) - Float64(y * Float64(c + Float64(y * b)))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -6e-14) tmp = t_1; elseif (y <= 1.4e+40) tmp = t * (-1.0 / ((0.0 - i) - (y * (c + (y * b))))); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6e-14], t$95$1, If[LessEqual[y, 1.4e+40], N[(t * N[(-1.0 / N[(N[(0.0 - i), $MachinePrecision] - N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -6 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.4 \cdot 10^{+40}:\\
\;\;\;\;t \cdot \frac{-1}{\left(0 - i\right) - y \cdot \left(c + y \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.9999999999999997e-14 or 1.4000000000000001e40 < y Initial program 11.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr11.9%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.2%
Simplified62.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6466.9%
Simplified66.9%
if -5.9999999999999997e-14 < y < 1.4000000000000001e40Initial program 96.5%
Taylor expanded in t around -inf
Simplified93.4%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6471.4%
Simplified71.4%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6470.0%
Simplified70.0%
Final simplification68.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ t (* y 230661.510616))) (t_2 (+ x (/ z y))))
(if (<= y -8e-14)
t_2
(if (<= y 2.6e-122) (/ t_1 i) (if (<= y 5.3e+38) (/ (/ t_1 y) c) t_2)))))
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);
double t_2 = x + (z / y);
double tmp;
if (y <= -8e-14) {
tmp = t_2;
} else if (y <= 2.6e-122) {
tmp = t_1 / i;
} else if (y <= 5.3e+38) {
tmp = (t_1 / y) / c;
} else {
tmp = t_2;
}
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) :: t_2
real(8) :: tmp
t_1 = t + (y * 230661.510616d0)
t_2 = x + (z / y)
if (y <= (-8d-14)) then
tmp = t_2
else if (y <= 2.6d-122) then
tmp = t_1 / i
else if (y <= 5.3d+38) then
tmp = (t_1 / y) / c
else
tmp = t_2
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 = t + (y * 230661.510616);
double t_2 = x + (z / y);
double tmp;
if (y <= -8e-14) {
tmp = t_2;
} else if (y <= 2.6e-122) {
tmp = t_1 / i;
} else if (y <= 5.3e+38) {
tmp = (t_1 / y) / c;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t + (y * 230661.510616) t_2 = x + (z / y) tmp = 0 if y <= -8e-14: tmp = t_2 elif y <= 2.6e-122: tmp = t_1 / i elif y <= 5.3e+38: tmp = (t_1 / y) / c else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t + Float64(y * 230661.510616)) t_2 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -8e-14) tmp = t_2; elseif (y <= 2.6e-122) tmp = Float64(t_1 / i); elseif (y <= 5.3e+38) tmp = Float64(Float64(t_1 / y) / c); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t + (y * 230661.510616); t_2 = x + (z / y); tmp = 0.0; if (y <= -8e-14) tmp = t_2; elseif (y <= 2.6e-122) tmp = t_1 / i; elseif (y <= 5.3e+38) tmp = (t_1 / y) / c; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8e-14], t$95$2, If[LessEqual[y, 2.6e-122], N[(t$95$1 / i), $MachinePrecision], If[LessEqual[y, 5.3e+38], N[(N[(t$95$1 / y), $MachinePrecision] / c), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + y \cdot 230661.510616\\
t_2 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -8 \cdot 10^{-14}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 2.6 \cdot 10^{-122}:\\
\;\;\;\;\frac{t\_1}{i}\\
\mathbf{elif}\;y \leq 5.3 \cdot 10^{+38}:\\
\;\;\;\;\frac{\frac{t\_1}{y}}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -7.99999999999999999e-14 or 5.30000000000000024e38 < y Initial program 11.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr11.9%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.2%
Simplified62.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6466.9%
Simplified66.9%
if -7.99999999999999999e-14 < y < 2.59999999999999975e-122Initial program 99.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6497.6%
Simplified97.6%
Taylor expanded in i around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.5%
Simplified67.5%
Taylor expanded in y around 0
Simplified67.5%
if 2.59999999999999975e-122 < y < 5.30000000000000024e38Initial program 88.8%
Taylor expanded in t around -inf
Simplified78.7%
Taylor expanded in c around -inf
/-lowering-/.f64N/A
Simplified35.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6429.2%
Simplified29.2%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.6e-14)
x
(if (<= y 4.4e+30)
(/ t i)
(if (<= y 2.4e+99) (* z (/ y b)) (if (<= y 2.8e+164) (/ z y) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.6e-14) {
tmp = x;
} else if (y <= 4.4e+30) {
tmp = t / i;
} else if (y <= 2.4e+99) {
tmp = z * (y / b);
} else if (y <= 2.8e+164) {
tmp = z / y;
} else {
tmp = x;
}
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 <= (-2.6d-14)) then
tmp = x
else if (y <= 4.4d+30) then
tmp = t / i
else if (y <= 2.4d+99) then
tmp = z * (y / b)
else if (y <= 2.8d+164) then
tmp = z / y
else
tmp = x
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 <= -2.6e-14) {
tmp = x;
} else if (y <= 4.4e+30) {
tmp = t / i;
} else if (y <= 2.4e+99) {
tmp = z * (y / b);
} else if (y <= 2.8e+164) {
tmp = z / y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -2.6e-14: tmp = x elif y <= 4.4e+30: tmp = t / i elif y <= 2.4e+99: tmp = z * (y / b) elif y <= 2.8e+164: tmp = z / y else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.6e-14) tmp = x; elseif (y <= 4.4e+30) tmp = Float64(t / i); elseif (y <= 2.4e+99) tmp = Float64(z * Float64(y / b)); elseif (y <= 2.8e+164) tmp = Float64(z / y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -2.6e-14) tmp = x; elseif (y <= 4.4e+30) tmp = t / i; elseif (y <= 2.4e+99) tmp = z * (y / b); elseif (y <= 2.8e+164) tmp = z / y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.6e-14], x, If[LessEqual[y, 4.4e+30], N[(t / i), $MachinePrecision], If[LessEqual[y, 2.4e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.8e+164], N[(z / y), $MachinePrecision], x]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.6 \cdot 10^{-14}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 4.4 \cdot 10^{+30}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{elif}\;y \leq 2.8 \cdot 10^{+164}:\\
\;\;\;\;\frac{z}{y}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -2.59999999999999997e-14 or 2.8000000000000002e164 < y Initial program 12.6%
Taylor expanded in y around inf
Simplified54.0%
if -2.59999999999999997e-14 < y < 4.4e30Initial program 97.2%
Taylor expanded in y around 0
/-lowering-/.f6447.7%
Simplified47.7%
if 4.4e30 < y < 2.4000000000000001e99Initial program 20.7%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f644.5%
Simplified4.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr6.6%
Taylor expanded in b around inf
/-lowering-/.f6441.4%
Simplified41.4%
if 2.4000000000000001e99 < y < 2.8000000000000002e164Initial program 0.4%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f640.4%
Simplified0.4%
Taylor expanded in y around inf
/-lowering-/.f6461.8%
Simplified61.8%
Final simplification50.6%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -8.5e-14)
x
(if (<= y 1.8e-80)
(/ t i)
(if (<= y 3.4e+99) (/ z a) (if (<= y 7.5e+164) (/ z y) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -8.5e-14) {
tmp = x;
} else if (y <= 1.8e-80) {
tmp = t / i;
} else if (y <= 3.4e+99) {
tmp = z / a;
} else if (y <= 7.5e+164) {
tmp = z / y;
} else {
tmp = x;
}
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 <= (-8.5d-14)) then
tmp = x
else if (y <= 1.8d-80) then
tmp = t / i
else if (y <= 3.4d+99) then
tmp = z / a
else if (y <= 7.5d+164) then
tmp = z / y
else
tmp = x
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 <= -8.5e-14) {
tmp = x;
} else if (y <= 1.8e-80) {
tmp = t / i;
} else if (y <= 3.4e+99) {
tmp = z / a;
} else if (y <= 7.5e+164) {
tmp = z / y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -8.5e-14: tmp = x elif y <= 1.8e-80: tmp = t / i elif y <= 3.4e+99: tmp = z / a elif y <= 7.5e+164: tmp = z / y else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -8.5e-14) tmp = x; elseif (y <= 1.8e-80) tmp = Float64(t / i); elseif (y <= 3.4e+99) tmp = Float64(z / a); elseif (y <= 7.5e+164) tmp = Float64(z / y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -8.5e-14) tmp = x; elseif (y <= 1.8e-80) tmp = t / i; elseif (y <= 3.4e+99) tmp = z / a; elseif (y <= 7.5e+164) tmp = z / y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -8.5e-14], x, If[LessEqual[y, 1.8e-80], N[(t / i), $MachinePrecision], If[LessEqual[y, 3.4e+99], N[(z / a), $MachinePrecision], If[LessEqual[y, 7.5e+164], N[(z / y), $MachinePrecision], x]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8.5 \cdot 10^{-14}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-80}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;y \leq 3.4 \cdot 10^{+99}:\\
\;\;\;\;\frac{z}{a}\\
\mathbf{elif}\;y \leq 7.5 \cdot 10^{+164}:\\
\;\;\;\;\frac{z}{y}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -8.50000000000000038e-14 or 7.49999999999999976e164 < y Initial program 12.6%
Taylor expanded in y around inf
Simplified54.0%
if -8.50000000000000038e-14 < y < 1.8e-80Initial program 99.7%
Taylor expanded in y around 0
/-lowering-/.f6458.3%
Simplified58.3%
if 1.8e-80 < y < 3.39999999999999984e99Initial program 63.5%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.0%
Simplified17.0%
Taylor expanded in a around inf
/-lowering-/.f6418.1%
Simplified18.1%
if 3.39999999999999984e99 < y < 7.49999999999999976e164Initial program 0.4%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f640.4%
Simplified0.4%
Taylor expanded in y around inf
/-lowering-/.f6461.8%
Simplified61.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -8.5e-14)
t_1
(if (<= y 3.5e+41) (* t (/ -1.0 (- (* y (- 0.0 c)) 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 + (z / y);
double tmp;
if (y <= -8.5e-14) {
tmp = t_1;
} else if (y <= 3.5e+41) {
tmp = t * (-1.0 / ((y * (0.0 - c)) - 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 + (z / y)
if (y <= (-8.5d-14)) then
tmp = t_1
else if (y <= 3.5d+41) then
tmp = t * ((-1.0d0) / ((y * (0.0d0 - c)) - 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 + (z / y);
double tmp;
if (y <= -8.5e-14) {
tmp = t_1;
} else if (y <= 3.5e+41) {
tmp = t * (-1.0 / ((y * (0.0 - c)) - i));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -8.5e-14: tmp = t_1 elif y <= 3.5e+41: tmp = t * (-1.0 / ((y * (0.0 - c)) - i)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -8.5e-14) tmp = t_1; elseif (y <= 3.5e+41) tmp = Float64(t * Float64(-1.0 / Float64(Float64(y * Float64(0.0 - c)) - i))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -8.5e-14) tmp = t_1; elseif (y <= 3.5e+41) tmp = t * (-1.0 / ((y * (0.0 - c)) - 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.5e-14], t$95$1, If[LessEqual[y, 3.5e+41], N[(t * N[(-1.0 / N[(N[(y * N[(0.0 - c), $MachinePrecision]), $MachinePrecision] - i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -8.5 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{+41}:\\
\;\;\;\;t \cdot \frac{-1}{y \cdot \left(0 - c\right) - i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.50000000000000038e-14 or 3.4999999999999999e41 < y Initial program 11.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr11.9%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.2%
Simplified62.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6466.9%
Simplified66.9%
if -8.50000000000000038e-14 < y < 3.4999999999999999e41Initial program 96.5%
Taylor expanded in t around -inf
Simplified93.4%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6471.4%
Simplified71.4%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6464.1%
Simplified64.1%
Final simplification65.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -6e-14)
t_1
(if (<= y 8.5e-8)
(/ (+ t (* y 230661.510616)) i)
(if (<= y 2.2e+99) (* z (/ y b)) 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 / y);
double tmp;
if (y <= -6e-14) {
tmp = t_1;
} else if (y <= 8.5e-8) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} 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 / y)
if (y <= (-6d-14)) then
tmp = t_1
else if (y <= 8.5d-8) then
tmp = (t + (y * 230661.510616d0)) / i
else if (y <= 2.2d+99) then
tmp = z * (y / b)
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 / y);
double tmp;
if (y <= -6e-14) {
tmp = t_1;
} else if (y <= 8.5e-8) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -6e-14: tmp = t_1 elif y <= 8.5e-8: tmp = (t + (y * 230661.510616)) / i elif y <= 2.2e+99: tmp = z * (y / b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -6e-14) tmp = t_1; elseif (y <= 8.5e-8) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); elseif (y <= 2.2e+99) tmp = Float64(z * Float64(y / b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -6e-14) tmp = t_1; elseif (y <= 8.5e-8) tmp = (t + (y * 230661.510616)) / i; elseif (y <= 2.2e+99) tmp = z * (y / b); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6e-14], t$95$1, If[LessEqual[y, 8.5e-8], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -6 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{-8}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.9999999999999997e-14 or 2.19999999999999978e99 < y Initial program 11.5%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr11.5%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6467.8%
Simplified67.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6473.0%
Simplified73.0%
if -5.9999999999999997e-14 < y < 8.49999999999999935e-8Initial program 99.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.9%
Simplified93.9%
Taylor expanded in i around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.9%
Simplified58.9%
Taylor expanded in y around 0
Simplified58.1%
if 8.49999999999999935e-8 < y < 2.19999999999999978e99Initial program 43.0%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6411.7%
Simplified11.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr16.4%
Taylor expanded in b around inf
/-lowering-/.f6425.7%
Simplified25.7%
Final simplification61.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ z y))))
(if (<= y -4.8e-14)
t_1
(if (<= y 4.7e+30) (/ t i) (if (<= y 2.2e+99) (* z (/ y b)) 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 / y);
double tmp;
if (y <= -4.8e-14) {
tmp = t_1;
} else if (y <= 4.7e+30) {
tmp = t / i;
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} 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 / y)
if (y <= (-4.8d-14)) then
tmp = t_1
else if (y <= 4.7d+30) then
tmp = t / i
else if (y <= 2.2d+99) then
tmp = z * (y / b)
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 / y);
double tmp;
if (y <= -4.8e-14) {
tmp = t_1;
} else if (y <= 4.7e+30) {
tmp = t / i;
} else if (y <= 2.2e+99) {
tmp = z * (y / b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + (z / y) tmp = 0 if y <= -4.8e-14: tmp = t_1 elif y <= 4.7e+30: tmp = t / i elif y <= 2.2e+99: tmp = z * (y / b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(z / y)) tmp = 0.0 if (y <= -4.8e-14) tmp = t_1; elseif (y <= 4.7e+30) tmp = Float64(t / i); elseif (y <= 2.2e+99) tmp = Float64(z * Float64(y / b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + (z / y); tmp = 0.0; if (y <= -4.8e-14) tmp = t_1; elseif (y <= 4.7e+30) tmp = t / i; elseif (y <= 2.2e+99) tmp = z * (y / b); 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[(z / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4.8e-14], t$95$1, If[LessEqual[y, 4.7e+30], N[(t / i), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z * N[(y / b), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z}{y}\\
\mathbf{if}\;y \leq -4.8 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.7 \cdot 10^{+30}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;z \cdot \frac{y}{b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -4.8e-14 or 2.19999999999999978e99 < y Initial program 11.5%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr11.5%
Taylor expanded in y around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6467.8%
Simplified67.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f6473.0%
Simplified73.0%
if -4.8e-14 < y < 4.6999999999999999e30Initial program 97.2%
Taylor expanded in y around 0
/-lowering-/.f6447.7%
Simplified47.7%
if 4.6999999999999999e30 < y < 2.19999999999999978e99Initial program 20.7%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f644.5%
Simplified4.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr6.6%
Taylor expanded in b around inf
/-lowering-/.f6441.4%
Simplified41.4%
Final simplification58.9%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y -9.5e-18) x (if (<= y 1.8e-80) (/ t i) (if (<= y 2.2e+99) (/ z a) x))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -9.5e-18) {
tmp = x;
} else if (y <= 1.8e-80) {
tmp = t / i;
} else if (y <= 2.2e+99) {
tmp = z / a;
} else {
tmp = x;
}
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 <= (-9.5d-18)) then
tmp = x
else if (y <= 1.8d-80) then
tmp = t / i
else if (y <= 2.2d+99) then
tmp = z / a
else
tmp = x
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 <= -9.5e-18) {
tmp = x;
} else if (y <= 1.8e-80) {
tmp = t / i;
} else if (y <= 2.2e+99) {
tmp = z / a;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -9.5e-18: tmp = x elif y <= 1.8e-80: tmp = t / i elif y <= 2.2e+99: tmp = z / a else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -9.5e-18) tmp = x; elseif (y <= 1.8e-80) tmp = Float64(t / i); elseif (y <= 2.2e+99) tmp = Float64(z / a); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -9.5e-18) tmp = x; elseif (y <= 1.8e-80) tmp = t / i; elseif (y <= 2.2e+99) tmp = z / a; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -9.5e-18], x, If[LessEqual[y, 1.8e-80], N[(t / i), $MachinePrecision], If[LessEqual[y, 2.2e+99], N[(z / a), $MachinePrecision], x]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.5 \cdot 10^{-18}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-80}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+99}:\\
\;\;\;\;\frac{z}{a}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -9.5000000000000003e-18 or 2.19999999999999978e99 < y Initial program 11.5%
Taylor expanded in y around inf
Simplified50.1%
if -9.5000000000000003e-18 < y < 1.8e-80Initial program 99.7%
Taylor expanded in y around 0
/-lowering-/.f6458.3%
Simplified58.3%
if 1.8e-80 < y < 2.19999999999999978e99Initial program 63.5%
Taylor expanded in z around inf
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.0%
Simplified17.0%
Taylor expanded in a around inf
/-lowering-/.f6418.1%
Simplified18.1%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y -5.8e-14) x (if (<= y 4.6e-5) (/ t i) x)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -5.8e-14) {
tmp = x;
} else if (y <= 4.6e-5) {
tmp = t / i;
} else {
tmp = x;
}
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 <= (-5.8d-14)) then
tmp = x
else if (y <= 4.6d-5) then
tmp = t / i
else
tmp = x
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 <= -5.8e-14) {
tmp = x;
} else if (y <= 4.6e-5) {
tmp = t / i;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -5.8e-14: tmp = x elif y <= 4.6e-5: tmp = t / i else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -5.8e-14) tmp = x; elseif (y <= 4.6e-5) tmp = Float64(t / i); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -5.8e-14) tmp = x; elseif (y <= 4.6e-5) tmp = t / i; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -5.8e-14], x, If[LessEqual[y, 4.6e-5], N[(t / i), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.8 \cdot 10^{-14}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 4.6 \cdot 10^{-5}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -5.8000000000000005e-14 or 4.6e-5 < y Initial program 15.4%
Taylor expanded in y around inf
Simplified42.7%
if -5.8000000000000005e-14 < y < 4.6e-5Initial program 99.6%
Taylor expanded in y around 0
/-lowering-/.f6450.7%
Simplified50.7%
(FPCore (x y z t a b c i) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return x;
}
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
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return x;
}
def code(x, y, z, t, a, b, c, i): return x
function code(x, y, z, t, a, b, c, i) return x end
function tmp = code(x, y, z, t, a, b, c, i) tmp = x; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 53.6%
Taylor expanded in y around inf
Simplified24.9%
herbie shell --seed 2024161
(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)))