
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ x (/ (/ (- (/ t y) y) z) 3.0))))
(if (<= y -8.5e-52)
t_1
(if (<= y 4.8e-17) (+ x (/ (/ t (* z 3.0)) y)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x + ((((t / y) - y) / z) / 3.0);
double tmp;
if (y <= -8.5e-52) {
tmp = t_1;
} else if (y <= 4.8e-17) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((((t / y) - y) / z) / 3.0d0)
if (y <= (-8.5d-52)) then
tmp = t_1
else if (y <= 4.8d-17) then
tmp = x + ((t / (z * 3.0d0)) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x + ((((t / y) - y) / z) / 3.0);
double tmp;
if (y <= -8.5e-52) {
tmp = t_1;
} else if (y <= 4.8e-17) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x + ((((t / y) - y) / z) / 3.0) tmp = 0 if y <= -8.5e-52: tmp = t_1 elif y <= 4.8e-17: tmp = x + ((t / (z * 3.0)) / y) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x + Float64(Float64(Float64(Float64(t / y) - y) / z) / 3.0)) tmp = 0.0 if (y <= -8.5e-52) tmp = t_1; elseif (y <= 4.8e-17) tmp = Float64(x + Float64(Float64(t / Float64(z * 3.0)) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x + ((((t / y) - y) / z) / 3.0); tmp = 0.0; if (y <= -8.5e-52) tmp = t_1; elseif (y <= 4.8e-17) tmp = x + ((t / (z * 3.0)) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x + N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.5e-52], t$95$1, If[LessEqual[y, 4.8e-17], N[(x + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{\frac{\frac{t}{y} - y}{z}}{3}\\
\mathbf{if}\;y \leq -8.5 \cdot 10^{-52}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-17}:\\
\;\;\;\;x + \frac{\frac{t}{z \cdot 3}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.50000000000000006e-52 or 4.79999999999999973e-17 < y Initial program 97.8%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.7%
Applied egg-rr99.7%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f6499.9%
Applied egg-rr99.9%
if -8.50000000000000006e-52 < y < 4.79999999999999973e-17Initial program 92.1%
Taylor expanded in x around inf
Simplified90.7%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6497.2%
Applied egg-rr97.2%
Final simplification98.7%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (+ (- x (/ y (* z 3.0))) (/ t (* y (* z 3.0)))))) (if (<= t_1 1e+307) t_1 (+ x (* (/ (- (/ t y) y) -3.0) (/ -1.0 z))))))
double code(double x, double y, double z, double t) {
double t_1 = (x - (y / (z * 3.0))) + (t / (y * (z * 3.0)));
double tmp;
if (t_1 <= 1e+307) {
tmp = t_1;
} else {
tmp = x + ((((t / y) - y) / -3.0) * (-1.0 / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x - (y / (z * 3.0d0))) + (t / (y * (z * 3.0d0)))
if (t_1 <= 1d+307) then
tmp = t_1
else
tmp = x + ((((t / y) - y) / (-3.0d0)) * ((-1.0d0) / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x - (y / (z * 3.0))) + (t / (y * (z * 3.0)));
double tmp;
if (t_1 <= 1e+307) {
tmp = t_1;
} else {
tmp = x + ((((t / y) - y) / -3.0) * (-1.0 / z));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x - (y / (z * 3.0))) + (t / (y * (z * 3.0))) tmp = 0 if t_1 <= 1e+307: tmp = t_1 else: tmp = x + ((((t / y) - y) / -3.0) * (-1.0 / z)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(y * Float64(z * 3.0)))) tmp = 0.0 if (t_1 <= 1e+307) tmp = t_1; else tmp = Float64(x + Float64(Float64(Float64(Float64(t / y) - y) / -3.0) * Float64(-1.0 / z))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x - (y / (z * 3.0))) + (t / (y * (z * 3.0))); tmp = 0.0; if (t_1 <= 1e+307) tmp = t_1; else tmp = x + ((((t / y) - y) / -3.0) * (-1.0 / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(y * N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 1e+307], t$95$1, N[(x + N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / -3.0), $MachinePrecision] * N[(-1.0 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{y \cdot \left(z \cdot 3\right)}\\
\mathbf{if}\;t\_1 \leq 10^{+307}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\frac{t}{y} - y}{-3} \cdot \frac{-1}{z}\\
\end{array}
\end{array}
if (+.f64 (-.f64 x (/.f64 y (*.f64 z #s(literal 3 binary64)))) (/.f64 t (*.f64 (*.f64 z #s(literal 3 binary64)) y))) < 9.99999999999999986e306Initial program 98.1%
if 9.99999999999999986e306 < (+.f64 (-.f64 x (/.f64 y (*.f64 z #s(literal 3 binary64)))) (/.f64 t (*.f64 (*.f64 z #s(literal 3 binary64)) y))) Initial program 82.3%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.9%
Simplified99.9%
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
metadata-evalN/A
distribute-rgt-neg-inN/A
associate-*r/N/A
*-commutativeN/A
distribute-lft-neg-inN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
metadata-evalN/A
/-lowering-/.f64100.0%
Applied egg-rr100.0%
Final simplification98.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -2.6e-24)
(+ x (/ t_1 (* z 3.0)))
(if (<= y 1.8e-30)
(+ x (/ (/ t (* z 3.0)) y))
(- x (* -0.3333333333333333 (/ t_1 z)))))))
double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -2.6e-24) {
tmp = x + (t_1 / (z * 3.0));
} else if (y <= 1.8e-30) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = x - (-0.3333333333333333 * (t_1 / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (t / y) - y
if (y <= (-2.6d-24)) then
tmp = x + (t_1 / (z * 3.0d0))
else if (y <= 1.8d-30) then
tmp = x + ((t / (z * 3.0d0)) / y)
else
tmp = x - ((-0.3333333333333333d0) * (t_1 / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -2.6e-24) {
tmp = x + (t_1 / (z * 3.0));
} else if (y <= 1.8e-30) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = x - (-0.3333333333333333 * (t_1 / z));
}
return tmp;
}
def code(x, y, z, t): t_1 = (t / y) - y tmp = 0 if y <= -2.6e-24: tmp = x + (t_1 / (z * 3.0)) elif y <= 1.8e-30: tmp = x + ((t / (z * 3.0)) / y) else: tmp = x - (-0.3333333333333333 * (t_1 / z)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(t / y) - y) tmp = 0.0 if (y <= -2.6e-24) tmp = Float64(x + Float64(t_1 / Float64(z * 3.0))); elseif (y <= 1.8e-30) tmp = Float64(x + Float64(Float64(t / Float64(z * 3.0)) / y)); else tmp = Float64(x - Float64(-0.3333333333333333 * Float64(t_1 / z))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (t / y) - y; tmp = 0.0; if (y <= -2.6e-24) tmp = x + (t_1 / (z * 3.0)); elseif (y <= 1.8e-30) tmp = x + ((t / (z * 3.0)) / y); else tmp = x - (-0.3333333333333333 * (t_1 / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[y, -2.6e-24], N[(x + N[(t$95$1 / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.8e-30], N[(x + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(x - N[(-0.3333333333333333 * N[(t$95$1 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{y} - y\\
\mathbf{if}\;y \leq -2.6 \cdot 10^{-24}:\\
\;\;\;\;x + \frac{t\_1}{z \cdot 3}\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-30}:\\
\;\;\;\;x + \frac{\frac{t}{z \cdot 3}}{y}\\
\mathbf{else}:\\
\;\;\;\;x - -0.3333333333333333 \cdot \frac{t\_1}{z}\\
\end{array}
\end{array}
if y < -2.6e-24Initial program 97.0%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.9%
Applied egg-rr99.9%
if -2.6e-24 < y < 1.8000000000000002e-30Initial program 92.3%
Taylor expanded in x around inf
Simplified90.9%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6497.2%
Applied egg-rr97.2%
if 1.8000000000000002e-30 < y Initial program 98.4%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.7%
Simplified99.7%
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
metadata-evalN/A
distribute-rgt-neg-inN/A
associate-*r/N/A
distribute-rgt-neg-inN/A
times-fracN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f6499.7%
Applied egg-rr99.7%
Final simplification98.7%
(FPCore (x y z t)
:precision binary64
(if (<= y -8e-52)
(+ x (* (- y (/ t y)) (/ -0.3333333333333333 z)))
(if (<= y 7e-35)
(+ x (/ (/ t (* z 3.0)) y))
(- x (* -0.3333333333333333 (/ (- (/ t y) y) z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -8e-52) {
tmp = x + ((y - (t / y)) * (-0.3333333333333333 / z));
} else if (y <= 7e-35) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = x - (-0.3333333333333333 * (((t / y) - y) / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-8d-52)) then
tmp = x + ((y - (t / y)) * ((-0.3333333333333333d0) / z))
else if (y <= 7d-35) then
tmp = x + ((t / (z * 3.0d0)) / y)
else
tmp = x - ((-0.3333333333333333d0) * (((t / y) - y) / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -8e-52) {
tmp = x + ((y - (t / y)) * (-0.3333333333333333 / z));
} else if (y <= 7e-35) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = x - (-0.3333333333333333 * (((t / y) - y) / z));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -8e-52: tmp = x + ((y - (t / y)) * (-0.3333333333333333 / z)) elif y <= 7e-35: tmp = x + ((t / (z * 3.0)) / y) else: tmp = x - (-0.3333333333333333 * (((t / y) - y) / z)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -8e-52) tmp = Float64(x + Float64(Float64(y - Float64(t / y)) * Float64(-0.3333333333333333 / z))); elseif (y <= 7e-35) tmp = Float64(x + Float64(Float64(t / Float64(z * 3.0)) / y)); else tmp = Float64(x - Float64(-0.3333333333333333 * Float64(Float64(Float64(t / y) - y) / z))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -8e-52) tmp = x + ((y - (t / y)) * (-0.3333333333333333 / z)); elseif (y <= 7e-35) tmp = x + ((t / (z * 3.0)) / y); else tmp = x - (-0.3333333333333333 * (((t / y) - y) / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -8e-52], N[(x + N[(N[(y - N[(t / y), $MachinePrecision]), $MachinePrecision] * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7e-35], N[(x + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(x - N[(-0.3333333333333333 * N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8 \cdot 10^{-52}:\\
\;\;\;\;x + \left(y - \frac{t}{y}\right) \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 7 \cdot 10^{-35}:\\
\;\;\;\;x + \frac{\frac{t}{z \cdot 3}}{y}\\
\mathbf{else}:\\
\;\;\;\;x - -0.3333333333333333 \cdot \frac{\frac{t}{y} - y}{z}\\
\end{array}
\end{array}
if y < -8.0000000000000001e-52Initial program 97.2%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.8%
Simplified99.8%
if -8.0000000000000001e-52 < y < 6.99999999999999992e-35Initial program 92.0%
Taylor expanded in x around inf
Simplified90.5%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6497.1%
Applied egg-rr97.1%
if 6.99999999999999992e-35 < y Initial program 98.4%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.7%
Simplified99.7%
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
metadata-evalN/A
distribute-rgt-neg-inN/A
associate-*r/N/A
distribute-rgt-neg-inN/A
times-fracN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f6499.7%
Applied egg-rr99.7%
Final simplification98.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ x (* (- y (/ t y)) (/ -0.3333333333333333 z)))))
(if (<= y -8e-52)
t_1
(if (<= y 2.65e-39) (+ x (/ (/ t (* z 3.0)) y)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x + ((y - (t / y)) * (-0.3333333333333333 / z));
double tmp;
if (y <= -8e-52) {
tmp = t_1;
} else if (y <= 2.65e-39) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((y - (t / y)) * ((-0.3333333333333333d0) / z))
if (y <= (-8d-52)) then
tmp = t_1
else if (y <= 2.65d-39) then
tmp = x + ((t / (z * 3.0d0)) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x + ((y - (t / y)) * (-0.3333333333333333 / z));
double tmp;
if (y <= -8e-52) {
tmp = t_1;
} else if (y <= 2.65e-39) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x + ((y - (t / y)) * (-0.3333333333333333 / z)) tmp = 0 if y <= -8e-52: tmp = t_1 elif y <= 2.65e-39: tmp = x + ((t / (z * 3.0)) / y) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x + Float64(Float64(y - Float64(t / y)) * Float64(-0.3333333333333333 / z))) tmp = 0.0 if (y <= -8e-52) tmp = t_1; elseif (y <= 2.65e-39) tmp = Float64(x + Float64(Float64(t / Float64(z * 3.0)) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x + ((y - (t / y)) * (-0.3333333333333333 / z)); tmp = 0.0; if (y <= -8e-52) tmp = t_1; elseif (y <= 2.65e-39) tmp = x + ((t / (z * 3.0)) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x + N[(N[(y - N[(t / y), $MachinePrecision]), $MachinePrecision] * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8e-52], t$95$1, If[LessEqual[y, 2.65e-39], N[(x + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(y - \frac{t}{y}\right) \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{if}\;y \leq -8 \cdot 10^{-52}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.65 \cdot 10^{-39}:\\
\;\;\;\;x + \frac{\frac{t}{z \cdot 3}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.0000000000000001e-52 or 2.65000000000000002e-39 < y Initial program 97.8%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.8%
Simplified99.8%
if -8.0000000000000001e-52 < y < 2.65000000000000002e-39Initial program 92.0%
Taylor expanded in x around inf
Simplified90.5%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6497.1%
Applied egg-rr97.1%
Final simplification98.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ (/ y z) 3.0))))
(if (<= y -6.6e-24)
t_1
(if (<= y 1.1e+39) (+ x (/ (/ t (* z 3.0)) y)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 1.1e+39) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-6.6d-24)) then
tmp = t_1
else if (y <= 1.1d+39) then
tmp = x + ((t / (z * 3.0d0)) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 1.1e+39) {
tmp = x + ((t / (z * 3.0)) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -6.6e-24: tmp = t_1 elif y <= 1.1e+39: tmp = x + ((t / (z * 3.0)) / y) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -6.6e-24) tmp = t_1; elseif (y <= 1.1e+39) tmp = Float64(x + Float64(Float64(t / Float64(z * 3.0)) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -6.6e-24) tmp = t_1; elseif (y <= 1.1e+39) tmp = x + ((t / (z * 3.0)) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.6e-24], t$95$1, If[LessEqual[y, 1.1e+39], N[(x + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -6.6 \cdot 10^{-24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{+39}:\\
\;\;\;\;x + \frac{\frac{t}{z \cdot 3}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.59999999999999968e-24 or 1.1000000000000001e39 < y Initial program 97.6%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
Taylor expanded in y around inf
Simplified93.0%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6493.1%
Applied egg-rr93.1%
if -6.59999999999999968e-24 < y < 1.1000000000000001e39Initial program 93.2%
Taylor expanded in x around inf
Simplified88.8%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6493.7%
Applied egg-rr93.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ (/ y z) 3.0))))
(if (<= y -6.6e-24)
t_1
(if (<= y 2.6e+39) (+ x (/ (/ t z) (* y 3.0))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 2.6e+39) {
tmp = x + ((t / z) / (y * 3.0));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-6.6d-24)) then
tmp = t_1
else if (y <= 2.6d+39) then
tmp = x + ((t / z) / (y * 3.0d0))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 2.6e+39) {
tmp = x + ((t / z) / (y * 3.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -6.6e-24: tmp = t_1 elif y <= 2.6e+39: tmp = x + ((t / z) / (y * 3.0)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -6.6e-24) tmp = t_1; elseif (y <= 2.6e+39) tmp = Float64(x + Float64(Float64(t / z) / Float64(y * 3.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -6.6e-24) tmp = t_1; elseif (y <= 2.6e+39) tmp = x + ((t / z) / (y * 3.0)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.6e-24], t$95$1, If[LessEqual[y, 2.6e+39], N[(x + N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -6.6 \cdot 10^{-24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.6 \cdot 10^{+39}:\\
\;\;\;\;x + \frac{\frac{t}{z}}{y \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.59999999999999968e-24 or 2.6e39 < y Initial program 97.6%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
Taylor expanded in y around inf
Simplified93.0%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6493.1%
Applied egg-rr93.1%
if -6.59999999999999968e-24 < y < 2.6e39Initial program 93.2%
Taylor expanded in x around inf
Simplified88.8%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6493.7%
Applied egg-rr93.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ (/ y z) 3.0))))
(if (<= y -4.2e-24)
t_1
(if (<= y 1.1e+39) (+ x (/ t (* y (* z 3.0)))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -4.2e-24) {
tmp = t_1;
} else if (y <= 1.1e+39) {
tmp = x + (t / (y * (z * 3.0)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-4.2d-24)) then
tmp = t_1
else if (y <= 1.1d+39) then
tmp = x + (t / (y * (z * 3.0d0)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -4.2e-24) {
tmp = t_1;
} else if (y <= 1.1e+39) {
tmp = x + (t / (y * (z * 3.0)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -4.2e-24: tmp = t_1 elif y <= 1.1e+39: tmp = x + (t / (y * (z * 3.0))) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -4.2e-24) tmp = t_1; elseif (y <= 1.1e+39) tmp = Float64(x + Float64(t / Float64(y * Float64(z * 3.0)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -4.2e-24) tmp = t_1; elseif (y <= 1.1e+39) tmp = x + (t / (y * (z * 3.0))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4.2e-24], t$95$1, If[LessEqual[y, 1.1e+39], N[(x + N[(t / N[(y * N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -4.2 \cdot 10^{-24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{+39}:\\
\;\;\;\;x + \frac{t}{y \cdot \left(z \cdot 3\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -4.1999999999999999e-24 or 1.1000000000000001e39 < y Initial program 97.6%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
Taylor expanded in y around inf
Simplified93.0%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6493.1%
Applied egg-rr93.1%
if -4.1999999999999999e-24 < y < 1.1000000000000001e39Initial program 93.2%
Taylor expanded in x around inf
Simplified88.8%
Final simplification91.0%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (- x (/ (/ y z) 3.0)))) (if (<= y -9e-51) t_1 (if (<= y 4.2e-87) (/ (/ (/ t 3.0) z) y) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -9e-51) {
tmp = t_1;
} else if (y <= 4.2e-87) {
tmp = ((t / 3.0) / z) / y;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-9d-51)) then
tmp = t_1
else if (y <= 4.2d-87) then
tmp = ((t / 3.0d0) / z) / y
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -9e-51) {
tmp = t_1;
} else if (y <= 4.2e-87) {
tmp = ((t / 3.0) / z) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -9e-51: tmp = t_1 elif y <= 4.2e-87: tmp = ((t / 3.0) / z) / y else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -9e-51) tmp = t_1; elseif (y <= 4.2e-87) tmp = Float64(Float64(Float64(t / 3.0) / z) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -9e-51) tmp = t_1; elseif (y <= 4.2e-87) tmp = ((t / 3.0) / z) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9e-51], t$95$1, If[LessEqual[y, 4.2e-87], N[(N[(N[(t / 3.0), $MachinePrecision] / z), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -9 \cdot 10^{-51}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.2 \cdot 10^{-87}:\\
\;\;\;\;\frac{\frac{\frac{t}{3}}{z}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.99999999999999948e-51 or 4.20000000000000014e-87 < y Initial program 97.4%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6498.0%
Applied egg-rr98.0%
Taylor expanded in y around inf
Simplified85.5%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6485.6%
Applied egg-rr85.6%
if -8.99999999999999948e-51 < y < 4.20000000000000014e-87Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
*-commutativeN/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
div-invN/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6473.9%
Applied egg-rr73.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ (/ y z) 3.0))))
(if (<= y -5.6e-45)
t_1
(if (<= y 1.22e-89) (/ 0.3333333333333333 (* y (/ z t))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -5.6e-45) {
tmp = t_1;
} else if (y <= 1.22e-89) {
tmp = 0.3333333333333333 / (y * (z / t));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-5.6d-45)) then
tmp = t_1
else if (y <= 1.22d-89) then
tmp = 0.3333333333333333d0 / (y * (z / t))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -5.6e-45) {
tmp = t_1;
} else if (y <= 1.22e-89) {
tmp = 0.3333333333333333 / (y * (z / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -5.6e-45: tmp = t_1 elif y <= 1.22e-89: tmp = 0.3333333333333333 / (y * (z / t)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -5.6e-45) tmp = t_1; elseif (y <= 1.22e-89) tmp = Float64(0.3333333333333333 / Float64(y * Float64(z / t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -5.6e-45) tmp = t_1; elseif (y <= 1.22e-89) tmp = 0.3333333333333333 / (y * (z / t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.6e-45], t$95$1, If[LessEqual[y, 1.22e-89], N[(0.3333333333333333 / N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -5.6 \cdot 10^{-45}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.22 \cdot 10^{-89}:\\
\;\;\;\;\frac{0.3333333333333333}{y \cdot \frac{z}{t}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.6000000000000003e-45 or 1.22e-89 < y Initial program 97.4%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6498.0%
Applied egg-rr98.0%
Taylor expanded in y around inf
Simplified85.5%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6485.6%
Applied egg-rr85.6%
if -5.6000000000000003e-45 < y < 1.22e-89Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
associate-/l*N/A
clear-numN/A
frac-timesN/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6466.0%
Applied egg-rr66.0%
clear-numN/A
un-div-invN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.9%
Applied egg-rr73.9%
Final simplification81.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ (/ y z) 3.0))))
(if (<= y -1.3e-50)
t_1
(if (<= y 1.5e-87) (* 0.3333333333333333 (/ t (* y z))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -1.3e-50) {
tmp = t_1;
} else if (y <= 1.5e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y / z) / 3.0d0)
if (y <= (-1.3d-50)) then
tmp = t_1
else if (y <= 1.5d-87) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - ((y / z) / 3.0);
double tmp;
if (y <= -1.3e-50) {
tmp = t_1;
} else if (y <= 1.5e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - ((y / z) / 3.0) tmp = 0 if y <= -1.3e-50: tmp = t_1 elif y <= 1.5e-87: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(Float64(y / z) / 3.0)) tmp = 0.0 if (y <= -1.3e-50) tmp = t_1; elseif (y <= 1.5e-87) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - ((y / z) / 3.0); tmp = 0.0; if (y <= -1.3e-50) tmp = t_1; elseif (y <= 1.5e-87) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[(y / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.3e-50], t$95$1, If[LessEqual[y, 1.5e-87], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\frac{y}{z}}{3}\\
\mathbf{if}\;y \leq -1.3 \cdot 10^{-50}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.5 \cdot 10^{-87}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.3000000000000001e-50 or 1.50000000000000008e-87 < y Initial program 97.4%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6498.0%
Applied egg-rr98.0%
Taylor expanded in y around inf
Simplified85.5%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6485.6%
Applied egg-rr85.6%
if -1.3000000000000001e-50 < y < 1.50000000000000008e-87Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
associate-/l*N/A
div-invN/A
associate-*l*N/A
inv-powN/A
clear-numN/A
inv-powN/A
pow-prod-downN/A
inv-powN/A
associate-/l/N/A
clear-numN/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
Final simplification78.6%
(FPCore (x y z t)
:precision binary64
(if (<= y -7.5e-41)
(- x (/ y (* z 3.0)))
(if (<= y 1.15e-87)
(* 0.3333333333333333 (/ t (* y z)))
(+ x (* -0.3333333333333333 (/ y z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -7.5e-41) {
tmp = x - (y / (z * 3.0));
} else if (y <= 1.15e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = x + (-0.3333333333333333 * (y / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-7.5d-41)) then
tmp = x - (y / (z * 3.0d0))
else if (y <= 1.15d-87) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = x + ((-0.3333333333333333d0) * (y / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -7.5e-41) {
tmp = x - (y / (z * 3.0));
} else if (y <= 1.15e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = x + (-0.3333333333333333 * (y / z));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -7.5e-41: tmp = x - (y / (z * 3.0)) elif y <= 1.15e-87: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = x + (-0.3333333333333333 * (y / z)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -7.5e-41) tmp = Float64(x - Float64(y / Float64(z * 3.0))); elseif (y <= 1.15e-87) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = Float64(x + Float64(-0.3333333333333333 * Float64(y / z))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -7.5e-41) tmp = x - (y / (z * 3.0)); elseif (y <= 1.15e-87) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = x + (-0.3333333333333333 * (y / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -7.5e-41], N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.15e-87], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.5 \cdot 10^{-41}:\\
\;\;\;\;x - \frac{y}{z \cdot 3}\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{-87}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;x + -0.3333333333333333 \cdot \frac{y}{z}\\
\end{array}
\end{array}
if y < -7.50000000000000049e-41Initial program 97.2%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6499.9%
Applied egg-rr99.9%
Taylor expanded in y around inf
Simplified89.8%
if -7.50000000000000049e-41 < y < 1.1500000000000001e-87Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
associate-/l*N/A
div-invN/A
associate-*l*N/A
inv-powN/A
clear-numN/A
inv-powN/A
pow-prod-downN/A
inv-powN/A
associate-/l/N/A
clear-numN/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
if 1.1500000000000001e-87 < y Initial program 97.6%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6496.7%
Simplified96.7%
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
metadata-evalN/A
distribute-rgt-neg-inN/A
associate-*r/N/A
distribute-rgt-neg-inN/A
times-fracN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f6496.7%
Applied egg-rr96.7%
Taylor expanded in y around inf
Simplified82.3%
Final simplification78.5%
(FPCore (x y z t)
:precision binary64
(if (<= y -4.4e-45)
(+ x (* y (/ -0.3333333333333333 z)))
(if (<= y 1.7e-87)
(* 0.3333333333333333 (/ t (* y z)))
(+ x (* -0.3333333333333333 (/ y z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -4.4e-45) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 1.7e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = x + (-0.3333333333333333 * (y / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-4.4d-45)) then
tmp = x + (y * ((-0.3333333333333333d0) / z))
else if (y <= 1.7d-87) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = x + ((-0.3333333333333333d0) * (y / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -4.4e-45) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 1.7e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = x + (-0.3333333333333333 * (y / z));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -4.4e-45: tmp = x + (y * (-0.3333333333333333 / z)) elif y <= 1.7e-87: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = x + (-0.3333333333333333 * (y / z)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -4.4e-45) tmp = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))); elseif (y <= 1.7e-87) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = Float64(x + Float64(-0.3333333333333333 * Float64(y / z))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -4.4e-45) tmp = x + (y * (-0.3333333333333333 / z)); elseif (y <= 1.7e-87) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = x + (-0.3333333333333333 * (y / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -4.4e-45], N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.7e-87], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.4 \cdot 10^{-45}:\\
\;\;\;\;x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{-87}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;x + -0.3333333333333333 \cdot \frac{y}{z}\\
\end{array}
\end{array}
if y < -4.39999999999999987e-45Initial program 97.2%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.8%
Simplified99.8%
Taylor expanded in y around inf
Simplified89.7%
if -4.39999999999999987e-45 < y < 1.6999999999999999e-87Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
associate-/l*N/A
div-invN/A
associate-*l*N/A
inv-powN/A
clear-numN/A
inv-powN/A
pow-prod-downN/A
inv-powN/A
associate-/l/N/A
clear-numN/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
if 1.6999999999999999e-87 < y Initial program 97.6%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6496.7%
Simplified96.7%
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
metadata-evalN/A
distribute-rgt-neg-inN/A
associate-*r/N/A
distribute-rgt-neg-inN/A
times-fracN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f6496.7%
Applied egg-rr96.7%
Taylor expanded in y around inf
Simplified82.3%
Final simplification78.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ x (* y (/ -0.3333333333333333 z)))))
(if (<= y -3.4e-48)
t_1
(if (<= y 2.85e-87) (* 0.3333333333333333 (/ t (* y z))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x + (y * (-0.3333333333333333 / z));
double tmp;
if (y <= -3.4e-48) {
tmp = t_1;
} else if (y <= 2.85e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x + (y * ((-0.3333333333333333d0) / z))
if (y <= (-3.4d-48)) then
tmp = t_1
else if (y <= 2.85d-87) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x + (y * (-0.3333333333333333 / z));
double tmp;
if (y <= -3.4e-48) {
tmp = t_1;
} else if (y <= 2.85e-87) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x + (y * (-0.3333333333333333 / z)) tmp = 0 if y <= -3.4e-48: tmp = t_1 elif y <= 2.85e-87: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))) tmp = 0.0 if (y <= -3.4e-48) tmp = t_1; elseif (y <= 2.85e-87) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x + (y * (-0.3333333333333333 / z)); tmp = 0.0; if (y <= -3.4e-48) tmp = t_1; elseif (y <= 2.85e-87) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.4e-48], t$95$1, If[LessEqual[y, 2.85e-87], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{if}\;y \leq -3.4 \cdot 10^{-48}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.85 \cdot 10^{-87}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.40000000000000028e-48 or 2.85e-87 < y Initial program 97.4%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6498.0%
Simplified98.0%
Taylor expanded in y around inf
Simplified85.5%
if -3.40000000000000028e-48 < y < 2.85e-87Initial program 91.9%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6488.7%
Applied egg-rr88.7%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6473.8%
Simplified73.8%
associate-/l*N/A
div-invN/A
associate-*l*N/A
inv-powN/A
clear-numN/A
inv-powN/A
pow-prod-downN/A
inv-powN/A
associate-/l/N/A
clear-numN/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
Final simplification78.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (/ y z) -3.0)))
(if (<= y -6.6e-24)
t_1
(if (<= y 3.5e+39) (* 0.3333333333333333 (/ t (* y z))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (y / z) / -3.0;
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 3.5e+39) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (y / z) / (-3.0d0)
if (y <= (-6.6d-24)) then
tmp = t_1
else if (y <= 3.5d+39) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (y / z) / -3.0;
double tmp;
if (y <= -6.6e-24) {
tmp = t_1;
} else if (y <= 3.5e+39) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (y / z) / -3.0 tmp = 0 if y <= -6.6e-24: tmp = t_1 elif y <= 3.5e+39: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(y / z) / -3.0) tmp = 0.0 if (y <= -6.6e-24) tmp = t_1; elseif (y <= 3.5e+39) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (y / z) / -3.0; tmp = 0.0; if (y <= -6.6e-24) tmp = t_1; elseif (y <= 3.5e+39) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y / z), $MachinePrecision] / -3.0), $MachinePrecision]}, If[LessEqual[y, -6.6e-24], t$95$1, If[LessEqual[y, 3.5e+39], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\frac{y}{z}}{-3}\\
\mathbf{if}\;y \leq -6.6 \cdot 10^{-24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{+39}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.59999999999999968e-24 or 3.5000000000000002e39 < y Initial program 97.6%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6499.8%
Simplified99.8%
Taylor expanded in y around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6468.9%
Simplified68.9%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6468.9%
Applied egg-rr68.9%
clear-numN/A
un-div-invN/A
div-invN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
metadata-eval69.0%
Applied egg-rr69.0%
if -6.59999999999999968e-24 < y < 3.5000000000000002e39Initial program 93.2%
associate-+l-N/A
--lowering--.f64N/A
*-commutativeN/A
associate-/r*N/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6489.3%
Applied egg-rr89.3%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6465.4%
Simplified65.4%
associate-/l*N/A
div-invN/A
associate-*l*N/A
inv-powN/A
clear-numN/A
inv-powN/A
pow-prod-downN/A
inv-powN/A
associate-/l/N/A
clear-numN/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6460.4%
Applied egg-rr60.4%
Final simplification64.8%
(FPCore (x y z t) :precision binary64 (if (<= z -2e+40) x (if (<= z 7e+169) (/ (/ y z) -3.0) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2e+40) {
tmp = x;
} else if (z <= 7e+169) {
tmp = (y / z) / -3.0;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-2d+40)) then
tmp = x
else if (z <= 7d+169) then
tmp = (y / z) / (-3.0d0)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2e+40) {
tmp = x;
} else if (z <= 7e+169) {
tmp = (y / z) / -3.0;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2e+40: tmp = x elif z <= 7e+169: tmp = (y / z) / -3.0 else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2e+40) tmp = x; elseif (z <= 7e+169) tmp = Float64(Float64(y / z) / -3.0); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -2e+40) tmp = x; elseif (z <= 7e+169) tmp = (y / z) / -3.0; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -2e+40], x, If[LessEqual[z, 7e+169], N[(N[(y / z), $MachinePrecision] / -3.0), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2 \cdot 10^{+40}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 7 \cdot 10^{+169}:\\
\;\;\;\;\frac{\frac{y}{z}}{-3}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -2.00000000000000006e40 or 7.00000000000000038e169 < z Initial program 98.9%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6489.1%
Simplified89.1%
Taylor expanded in x around inf
Simplified59.3%
if -2.00000000000000006e40 < z < 7.00000000000000038e169Initial program 93.3%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6498.1%
Simplified98.1%
Taylor expanded in y around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.4%
Simplified53.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6453.4%
Applied egg-rr53.4%
clear-numN/A
un-div-invN/A
div-invN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
metadata-eval53.5%
Applied egg-rr53.5%
(FPCore (x y z t) :precision binary64 (if (<= z -1e+39) x (if (<= z 5.5e+170) (/ (* y -0.3333333333333333) z) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1e+39) {
tmp = x;
} else if (z <= 5.5e+170) {
tmp = (y * -0.3333333333333333) / z;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-1d+39)) then
tmp = x
else if (z <= 5.5d+170) then
tmp = (y * (-0.3333333333333333d0)) / z
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1e+39) {
tmp = x;
} else if (z <= 5.5e+170) {
tmp = (y * -0.3333333333333333) / z;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -1e+39: tmp = x elif z <= 5.5e+170: tmp = (y * -0.3333333333333333) / z else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -1e+39) tmp = x; elseif (z <= 5.5e+170) tmp = Float64(Float64(y * -0.3333333333333333) / z); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -1e+39) tmp = x; elseif (z <= 5.5e+170) tmp = (y * -0.3333333333333333) / z; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -1e+39], x, If[LessEqual[z, 5.5e+170], N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1 \cdot 10^{+39}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 5.5 \cdot 10^{+170}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -9.9999999999999994e38 or 5.4999999999999999e170 < z Initial program 98.9%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6489.1%
Simplified89.1%
Taylor expanded in x around inf
Simplified59.3%
if -9.9999999999999994e38 < z < 5.4999999999999999e170Initial program 93.3%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6498.1%
Simplified98.1%
Taylor expanded in y around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.4%
Simplified53.4%
Final simplification55.6%
(FPCore (x y z t) :precision binary64 (if (<= z -1.04e+40) x (if (<= z 6.6e+169) (* -0.3333333333333333 (/ y z)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1.04e+40) {
tmp = x;
} else if (z <= 6.6e+169) {
tmp = -0.3333333333333333 * (y / z);
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-1.04d+40)) then
tmp = x
else if (z <= 6.6d+169) then
tmp = (-0.3333333333333333d0) * (y / z)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1.04e+40) {
tmp = x;
} else if (z <= 6.6e+169) {
tmp = -0.3333333333333333 * (y / z);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -1.04e+40: tmp = x elif z <= 6.6e+169: tmp = -0.3333333333333333 * (y / z) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -1.04e+40) tmp = x; elseif (z <= 6.6e+169) tmp = Float64(-0.3333333333333333 * Float64(y / z)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -1.04e+40) tmp = x; elseif (z <= 6.6e+169) tmp = -0.3333333333333333 * (y / z); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -1.04e+40], x, If[LessEqual[z, 6.6e+169], N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.04 \cdot 10^{+40}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 6.6 \cdot 10^{+169}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -1.04e40 or 6.5999999999999994e169 < z Initial program 98.9%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6489.1%
Simplified89.1%
Taylor expanded in x around inf
Simplified59.3%
if -1.04e40 < z < 6.5999999999999994e169Initial program 93.3%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6498.1%
Simplified98.1%
Taylor expanded in y around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.4%
Simplified53.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6453.4%
Applied egg-rr53.4%
Final simplification55.6%
(FPCore (x y z t) :precision binary64 (if (<= z -7.8e+39) x (if (<= z 6.6e+169) (* y (/ -0.3333333333333333 z)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -7.8e+39) {
tmp = x;
} else if (z <= 6.6e+169) {
tmp = y * (-0.3333333333333333 / z);
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-7.8d+39)) then
tmp = x
else if (z <= 6.6d+169) then
tmp = y * ((-0.3333333333333333d0) / z)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -7.8e+39) {
tmp = x;
} else if (z <= 6.6e+169) {
tmp = y * (-0.3333333333333333 / z);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -7.8e+39: tmp = x elif z <= 6.6e+169: tmp = y * (-0.3333333333333333 / z) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -7.8e+39) tmp = x; elseif (z <= 6.6e+169) tmp = Float64(y * Float64(-0.3333333333333333 / z)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -7.8e+39) tmp = x; elseif (z <= 6.6e+169) tmp = y * (-0.3333333333333333 / z); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -7.8e+39], x, If[LessEqual[z, 6.6e+169], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -7.8 \cdot 10^{+39}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 6.6 \cdot 10^{+169}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -7.8000000000000002e39 or 6.5999999999999994e169 < z Initial program 98.9%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6489.1%
Simplified89.1%
Taylor expanded in x around inf
Simplified59.3%
if -7.8000000000000002e39 < z < 6.5999999999999994e169Initial program 93.3%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6498.1%
Simplified98.1%
Taylor expanded in y around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.4%
Simplified53.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6453.4%
Applied egg-rr53.4%
(FPCore (x y z t) :precision binary64 x)
double code(double x, double y, double z, double t) {
return x;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x
end function
public static double code(double x, double y, double z, double t) {
return x;
}
def code(x, y, z, t): return x
function code(x, y, z, t) return x end
function tmp = code(x, y, z, t) tmp = x; end
code[x_, y_, z_, t_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 95.4%
sub-negN/A
associate-+l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
unsub-negN/A
neg-mul-1N/A
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
distribute-neg-fracN/A
neg-mul-1N/A
times-fracN/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f6494.7%
Simplified94.7%
Taylor expanded in x around inf
Simplified28.3%
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y)))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + ((t / (z * 3.0d0)) / y)
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y)
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(Float64(t / Float64(z * 3.0)) / y)) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{\frac{t}{z \cdot 3}}{y}
\end{array}
herbie shell --seed 2024138
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:alt
(! :herbie-platform default (+ (- x (/ y (* z 3))) (/ (/ t (* z 3)) y)))
(+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))