
(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 17 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 (if (or (<= y -1.8e-31) (not (<= y 3.9e-170))) (+ x (/ (- (/ t y) y) (* z 3.0))) (+ x (* 0.3333333333333333 (/ (/ t z) y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.8e-31) || !(y <= 3.9e-170)) {
tmp = x + (((t / y) - y) / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
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 <= (-1.8d-31)) .or. (.not. (y <= 3.9d-170))) then
tmp = x + (((t / y) - y) / (z * 3.0d0))
else
tmp = x + (0.3333333333333333d0 * ((t / z) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.8e-31) || !(y <= 3.9e-170)) {
tmp = x + (((t / y) - y) / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -1.8e-31) or not (y <= 3.9e-170): tmp = x + (((t / y) - y) / (z * 3.0)) else: tmp = x + (0.3333333333333333 * ((t / z) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -1.8e-31) || !(y <= 3.9e-170)) tmp = Float64(x + Float64(Float64(Float64(t / y) - y) / Float64(z * 3.0))); else tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(t / z) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -1.8e-31) || ~((y <= 3.9e-170))) tmp = x + (((t / y) - y) / (z * 3.0)); else tmp = x + (0.3333333333333333 * ((t / z) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -1.8e-31], N[Not[LessEqual[y, 3.9e-170]], $MachinePrecision]], N[(x + N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.8 \cdot 10^{-31} \lor \neg \left(y \leq 3.9 \cdot 10^{-170}\right):\\
\;\;\;\;x + \frac{\frac{t}{y} - y}{z \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\end{array}
\end{array}
if y < -1.80000000000000002e-31 or 3.90000000000000021e-170 < y Initial program 98.7%
sub-neg98.7%
associate-+l+98.7%
+-commutative98.7%
remove-double-neg98.7%
distribute-frac-neg98.7%
distribute-neg-in98.7%
remove-double-neg98.7%
sub-neg98.7%
neg-mul-198.7%
times-frac99.2%
distribute-frac-neg99.2%
neg-mul-199.2%
*-commutative99.2%
associate-/l*99.1%
*-commutative99.1%
Simplified99.7%
*-commutative99.7%
clear-num99.7%
div-inv99.7%
metadata-eval99.7%
un-div-inv99.8%
Applied egg-rr99.8%
if -1.80000000000000002e-31 < y < 3.90000000000000021e-170Initial program 91.9%
sub-neg91.9%
associate-+l+91.9%
+-commutative91.9%
remove-double-neg91.9%
distribute-frac-neg91.9%
distribute-neg-in91.9%
remove-double-neg91.9%
sub-neg91.9%
neg-mul-191.9%
times-frac88.9%
distribute-frac-neg88.9%
neg-mul-188.9%
*-commutative88.9%
associate-/l*88.9%
*-commutative88.9%
Simplified88.9%
Taylor expanded in t around inf 91.8%
div-inv90.0%
*-commutative90.0%
Applied egg-rr90.0%
associate-*r/91.8%
*-rgt-identity91.8%
associate-/r*98.7%
Simplified98.7%
Final simplification99.4%
(FPCore (x y z t) :precision binary64 (if (<= (+ (/ t (* y (* z 3.0))) (- x (/ y (* z 3.0)))) 1e+282) (+ (+ x (/ t (* z (* y 3.0)))) (/ y (* z -3.0))) (+ x (* (/ 1.0 (* z 3.0)) (- (/ t y) y)))))
double code(double x, double y, double z, double t) {
double tmp;
if (((t / (y * (z * 3.0))) + (x - (y / (z * 3.0)))) <= 1e+282) {
tmp = (x + (t / (z * (y * 3.0)))) + (y / (z * -3.0));
} else {
tmp = x + ((1.0 / (z * 3.0)) * ((t / y) - y));
}
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 (((t / (y * (z * 3.0d0))) + (x - (y / (z * 3.0d0)))) <= 1d+282) then
tmp = (x + (t / (z * (y * 3.0d0)))) + (y / (z * (-3.0d0)))
else
tmp = x + ((1.0d0 / (z * 3.0d0)) * ((t / y) - y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (((t / (y * (z * 3.0))) + (x - (y / (z * 3.0)))) <= 1e+282) {
tmp = (x + (t / (z * (y * 3.0)))) + (y / (z * -3.0));
} else {
tmp = x + ((1.0 / (z * 3.0)) * ((t / y) - y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if ((t / (y * (z * 3.0))) + (x - (y / (z * 3.0)))) <= 1e+282: tmp = (x + (t / (z * (y * 3.0)))) + (y / (z * -3.0)) else: tmp = x + ((1.0 / (z * 3.0)) * ((t / y) - y)) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(Float64(t / Float64(y * Float64(z * 3.0))) + Float64(x - Float64(y / Float64(z * 3.0)))) <= 1e+282) tmp = Float64(Float64(x + Float64(t / Float64(z * Float64(y * 3.0)))) + Float64(y / Float64(z * -3.0))); else tmp = Float64(x + Float64(Float64(1.0 / Float64(z * 3.0)) * Float64(Float64(t / y) - y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (((t / (y * (z * 3.0))) + (x - (y / (z * 3.0)))) <= 1e+282) tmp = (x + (t / (z * (y * 3.0)))) + (y / (z * -3.0)); else tmp = x + ((1.0 / (z * 3.0)) * ((t / y) - y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(N[(t / N[(y * N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e+282], N[(N[(x + N[(t / N[(z * N[(y * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(1.0 / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] * N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t}{y \cdot \left(z \cdot 3\right)} + \left(x - \frac{y}{z \cdot 3}\right) \leq 10^{+282}:\\
\;\;\;\;\left(x + \frac{t}{z \cdot \left(y \cdot 3\right)}\right) + \frac{y}{z \cdot -3}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{1}{z \cdot 3} \cdot \left(\frac{t}{y} - y\right)\\
\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))) < 1.00000000000000003e282Initial program 98.5%
+-commutative98.5%
associate-+r-98.5%
sub-neg98.5%
associate-*l*98.5%
*-commutative98.5%
distribute-frac-neg298.5%
distribute-rgt-neg-in98.5%
metadata-eval98.5%
Simplified98.5%
if 1.00000000000000003e282 < (+.f64 (-.f64 x (/.f64 y (*.f64 z #s(literal 3 binary64)))) (/.f64 t (*.f64 (*.f64 z #s(literal 3 binary64)) y))) Initial program 86.0%
sub-neg86.0%
associate-+l+86.0%
+-commutative86.0%
remove-double-neg86.0%
distribute-frac-neg86.0%
distribute-neg-in86.0%
remove-double-neg86.0%
sub-neg86.0%
neg-mul-186.0%
times-frac97.7%
distribute-frac-neg97.7%
neg-mul-197.7%
*-commutative97.7%
associate-/l*97.7%
*-commutative97.7%
Simplified100.0%
clear-num100.0%
div-inv100.0%
metadata-eval100.0%
inv-pow100.0%
Applied egg-rr100.0%
unpow-1100.0%
Simplified100.0%
Final simplification98.8%
(FPCore (x y z t)
:precision binary64
(if (<= y -2.9e+33)
(* y (/ -0.3333333333333333 z))
(if (<= y 1.7e-44)
(* 0.3333333333333333 (/ (/ t z) y))
(if (<= y 12000.0)
x
(if (<= y 1.08e+42)
(* t (/ (/ 0.3333333333333333 y) z))
(/ (* y -0.3333333333333333) z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -2.9e+33) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 1.7e-44) {
tmp = 0.3333333333333333 * ((t / z) / y);
} else if (y <= 12000.0) {
tmp = x;
} else if (y <= 1.08e+42) {
tmp = t * ((0.3333333333333333 / y) / z);
} else {
tmp = (y * -0.3333333333333333) / 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 <= (-2.9d+33)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 1.7d-44) then
tmp = 0.3333333333333333d0 * ((t / z) / y)
else if (y <= 12000.0d0) then
tmp = x
else if (y <= 1.08d+42) then
tmp = t * ((0.3333333333333333d0 / y) / z)
else
tmp = (y * (-0.3333333333333333d0)) / z
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -2.9e+33) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 1.7e-44) {
tmp = 0.3333333333333333 * ((t / z) / y);
} else if (y <= 12000.0) {
tmp = x;
} else if (y <= 1.08e+42) {
tmp = t * ((0.3333333333333333 / y) / z);
} else {
tmp = (y * -0.3333333333333333) / z;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -2.9e+33: tmp = y * (-0.3333333333333333 / z) elif y <= 1.7e-44: tmp = 0.3333333333333333 * ((t / z) / y) elif y <= 12000.0: tmp = x elif y <= 1.08e+42: tmp = t * ((0.3333333333333333 / y) / z) else: tmp = (y * -0.3333333333333333) / z return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -2.9e+33) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 1.7e-44) tmp = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)); elseif (y <= 12000.0) tmp = x; elseif (y <= 1.08e+42) tmp = Float64(t * Float64(Float64(0.3333333333333333 / y) / z)); else tmp = Float64(Float64(y * -0.3333333333333333) / z); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -2.9e+33) tmp = y * (-0.3333333333333333 / z); elseif (y <= 1.7e-44) tmp = 0.3333333333333333 * ((t / z) / y); elseif (y <= 12000.0) tmp = x; elseif (y <= 1.08e+42) tmp = t * ((0.3333333333333333 / y) / z); else tmp = (y * -0.3333333333333333) / z; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -2.9e+33], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.7e-44], N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 12000.0], x, If[LessEqual[y, 1.08e+42], N[(t * N[(N[(0.3333333333333333 / y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.9 \cdot 10^{+33}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{-44}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\mathbf{elif}\;y \leq 12000:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.08 \cdot 10^{+42}:\\
\;\;\;\;t \cdot \frac{\frac{0.3333333333333333}{y}}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -2.90000000000000025e33Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-lft-identity99.7%
Simplified99.7%
Taylor expanded in y around inf 73.7%
*-commutative73.7%
associate-*l/73.7%
associate-/l*73.8%
Simplified73.8%
if -2.90000000000000025e33 < y < 1.70000000000000008e-44Initial program 93.6%
+-commutative93.6%
associate-+r-93.6%
sub-neg93.6%
associate-*l*93.6%
*-commutative93.6%
distribute-frac-neg293.6%
distribute-rgt-neg-in93.6%
metadata-eval93.6%
Simplified93.6%
*-un-lft-identity93.6%
times-frac93.6%
Applied egg-rr93.6%
associate-*r/93.6%
associate-*l/93.6%
*-lft-identity93.6%
Simplified93.6%
Taylor expanded in t around inf 60.5%
associate-/l/63.4%
Simplified63.4%
if 1.70000000000000008e-44 < y < 12000Initial program 100.0%
+-commutative100.0%
associate-+r-100.0%
sub-neg100.0%
associate-*l*100.0%
*-commutative100.0%
distribute-frac-neg2100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in z around inf 65.5%
if 12000 < y < 1.08e42Initial program 88.9%
sub-neg88.9%
associate-+l+88.9%
+-commutative88.9%
remove-double-neg88.9%
distribute-frac-neg88.9%
distribute-neg-in88.9%
remove-double-neg88.9%
sub-neg88.9%
neg-mul-188.9%
times-frac88.9%
distribute-frac-neg88.9%
neg-mul-188.9%
*-commutative88.9%
associate-/l*88.9%
*-commutative88.9%
Simplified99.8%
Taylor expanded in t around inf 78.2%
Taylor expanded in t around inf 78.0%
associate-*r/78.2%
metadata-eval78.2%
Simplified78.2%
Taylor expanded in y around 0 57.2%
associate-/r*57.2%
Simplified57.2%
if 1.08e42 < y Initial program 99.9%
+-commutative99.9%
associate-+r-99.9%
sub-neg99.9%
associate-*l*99.9%
*-commutative99.9%
distribute-frac-neg299.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
*-un-lft-identity99.9%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.8%
associate-*l/99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 78.2%
*-commutative78.2%
associate-*l/78.2%
Applied egg-rr78.2%
(FPCore (x y z t)
:precision binary64
(if (<= y -4.1e+34)
(* y (/ -0.3333333333333333 z))
(if (<= y 4.8e-44)
(* 0.3333333333333333 (/ (/ t z) y))
(if (<= y 7400.0)
x
(if (<= y 1.25e+42)
(* 0.3333333333333333 (/ t (* y z)))
(/ (* y -0.3333333333333333) z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -4.1e+34) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 4.8e-44) {
tmp = 0.3333333333333333 * ((t / z) / y);
} else if (y <= 7400.0) {
tmp = x;
} else if (y <= 1.25e+42) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = (y * -0.3333333333333333) / 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.1d+34)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 4.8d-44) then
tmp = 0.3333333333333333d0 * ((t / z) / y)
else if (y <= 7400.0d0) then
tmp = x
else if (y <= 1.25d+42) then
tmp = 0.3333333333333333d0 * (t / (y * z))
else
tmp = (y * (-0.3333333333333333d0)) / z
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -4.1e+34) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 4.8e-44) {
tmp = 0.3333333333333333 * ((t / z) / y);
} else if (y <= 7400.0) {
tmp = x;
} else if (y <= 1.25e+42) {
tmp = 0.3333333333333333 * (t / (y * z));
} else {
tmp = (y * -0.3333333333333333) / z;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -4.1e+34: tmp = y * (-0.3333333333333333 / z) elif y <= 4.8e-44: tmp = 0.3333333333333333 * ((t / z) / y) elif y <= 7400.0: tmp = x elif y <= 1.25e+42: tmp = 0.3333333333333333 * (t / (y * z)) else: tmp = (y * -0.3333333333333333) / z return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -4.1e+34) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 4.8e-44) tmp = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)); elseif (y <= 7400.0) tmp = x; elseif (y <= 1.25e+42) tmp = Float64(0.3333333333333333 * Float64(t / Float64(y * z))); else tmp = Float64(Float64(y * -0.3333333333333333) / z); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -4.1e+34) tmp = y * (-0.3333333333333333 / z); elseif (y <= 4.8e-44) tmp = 0.3333333333333333 * ((t / z) / y); elseif (y <= 7400.0) tmp = x; elseif (y <= 1.25e+42) tmp = 0.3333333333333333 * (t / (y * z)); else tmp = (y * -0.3333333333333333) / z; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -4.1e+34], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e-44], N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7400.0], x, If[LessEqual[y, 1.25e+42], N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.1 \cdot 10^{+34}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-44}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\mathbf{elif}\;y \leq 7400:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+42}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -4.0999999999999998e34Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-lft-identity99.7%
Simplified99.7%
Taylor expanded in y around inf 73.7%
*-commutative73.7%
associate-*l/73.7%
associate-/l*73.8%
Simplified73.8%
if -4.0999999999999998e34 < y < 4.80000000000000017e-44Initial program 93.6%
+-commutative93.6%
associate-+r-93.6%
sub-neg93.6%
associate-*l*93.6%
*-commutative93.6%
distribute-frac-neg293.6%
distribute-rgt-neg-in93.6%
metadata-eval93.6%
Simplified93.6%
*-un-lft-identity93.6%
times-frac93.6%
Applied egg-rr93.6%
associate-*r/93.6%
associate-*l/93.6%
*-lft-identity93.6%
Simplified93.6%
Taylor expanded in t around inf 60.5%
associate-/l/63.4%
Simplified63.4%
if 4.80000000000000017e-44 < y < 7400Initial program 100.0%
+-commutative100.0%
associate-+r-100.0%
sub-neg100.0%
associate-*l*100.0%
*-commutative100.0%
distribute-frac-neg2100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in z around inf 65.5%
if 7400 < y < 1.25000000000000002e42Initial program 88.9%
+-commutative88.9%
associate-+r-88.9%
sub-neg88.9%
associate-*l*88.9%
*-commutative88.9%
distribute-frac-neg288.9%
distribute-rgt-neg-in88.9%
metadata-eval88.9%
Simplified88.9%
*-un-lft-identity88.9%
times-frac88.7%
Applied egg-rr88.7%
associate-*r/88.7%
associate-*l/88.7%
*-lft-identity88.7%
Simplified88.7%
Taylor expanded in t around inf 57.2%
if 1.25000000000000002e42 < y Initial program 99.9%
+-commutative99.9%
associate-+r-99.9%
sub-neg99.9%
associate-*l*99.9%
*-commutative99.9%
distribute-frac-neg299.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
*-un-lft-identity99.9%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.8%
associate-*l/99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 78.2%
*-commutative78.2%
associate-*l/78.2%
Applied egg-rr78.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ t (* y z)))))
(if (<= y -5.4e+34)
(* y (/ -0.3333333333333333 z))
(if (<= y 2.7e-45)
t_1
(if (<= y 19500.0)
x
(if (<= y 1e+42) t_1 (/ (* y -0.3333333333333333) z)))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * (t / (y * z));
double tmp;
if (y <= -5.4e+34) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 2.7e-45) {
tmp = t_1;
} else if (y <= 19500.0) {
tmp = x;
} else if (y <= 1e+42) {
tmp = t_1;
} else {
tmp = (y * -0.3333333333333333) / 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 = 0.3333333333333333d0 * (t / (y * z))
if (y <= (-5.4d+34)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 2.7d-45) then
tmp = t_1
else if (y <= 19500.0d0) then
tmp = x
else if (y <= 1d+42) then
tmp = t_1
else
tmp = (y * (-0.3333333333333333d0)) / z
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * (t / (y * z));
double tmp;
if (y <= -5.4e+34) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 2.7e-45) {
tmp = t_1;
} else if (y <= 19500.0) {
tmp = x;
} else if (y <= 1e+42) {
tmp = t_1;
} else {
tmp = (y * -0.3333333333333333) / z;
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * (t / (y * z)) tmp = 0 if y <= -5.4e+34: tmp = y * (-0.3333333333333333 / z) elif y <= 2.7e-45: tmp = t_1 elif y <= 19500.0: tmp = x elif y <= 1e+42: tmp = t_1 else: tmp = (y * -0.3333333333333333) / z return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(t / Float64(y * z))) tmp = 0.0 if (y <= -5.4e+34) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 2.7e-45) tmp = t_1; elseif (y <= 19500.0) tmp = x; elseif (y <= 1e+42) tmp = t_1; else tmp = Float64(Float64(y * -0.3333333333333333) / z); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * (t / (y * z)); tmp = 0.0; if (y <= -5.4e+34) tmp = y * (-0.3333333333333333 / z); elseif (y <= 2.7e-45) tmp = t_1; elseif (y <= 19500.0) tmp = x; elseif (y <= 1e+42) tmp = t_1; else tmp = (y * -0.3333333333333333) / z; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.4e+34], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.7e-45], t$95$1, If[LessEqual[y, 19500.0], x, If[LessEqual[y, 1e+42], t$95$1, N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\mathbf{if}\;y \leq -5.4 \cdot 10^{+34}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 2.7 \cdot 10^{-45}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 19500:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 10^{+42}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -5.4000000000000001e34Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-lft-identity99.7%
Simplified99.7%
Taylor expanded in y around inf 73.7%
*-commutative73.7%
associate-*l/73.7%
associate-/l*73.8%
Simplified73.8%
if -5.4000000000000001e34 < y < 2.69999999999999985e-45 or 19500 < y < 1.00000000000000004e42Initial program 93.3%
+-commutative93.3%
associate-+r-93.3%
sub-neg93.3%
associate-*l*93.3%
*-commutative93.3%
distribute-frac-neg293.3%
distribute-rgt-neg-in93.3%
metadata-eval93.3%
Simplified93.3%
*-un-lft-identity93.3%
times-frac93.3%
Applied egg-rr93.3%
associate-*r/93.3%
associate-*l/93.3%
*-lft-identity93.3%
Simplified93.3%
Taylor expanded in t around inf 60.3%
if 2.69999999999999985e-45 < y < 19500Initial program 100.0%
+-commutative100.0%
associate-+r-100.0%
sub-neg100.0%
associate-*l*100.0%
*-commutative100.0%
distribute-frac-neg2100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in z around inf 65.5%
if 1.00000000000000004e42 < y Initial program 99.9%
+-commutative99.9%
associate-+r-99.9%
sub-neg99.9%
associate-*l*99.9%
*-commutative99.9%
distribute-frac-neg299.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
*-un-lft-identity99.9%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.8%
associate-*l/99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 78.2%
*-commutative78.2%
associate-*l/78.2%
Applied egg-rr78.2%
(FPCore (x y z t) :precision binary64 (if (or (<= y -1.65e-30) (not (<= y 1.45e-169))) (+ x (* 0.3333333333333333 (/ (- (/ t y) y) z))) (+ x (* 0.3333333333333333 (/ (/ t z) y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.65e-30) || !(y <= 1.45e-169)) {
tmp = x + (0.3333333333333333 * (((t / y) - y) / z));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
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 <= (-1.65d-30)) .or. (.not. (y <= 1.45d-169))) then
tmp = x + (0.3333333333333333d0 * (((t / y) - y) / z))
else
tmp = x + (0.3333333333333333d0 * ((t / z) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.65e-30) || !(y <= 1.45e-169)) {
tmp = x + (0.3333333333333333 * (((t / y) - y) / z));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -1.65e-30) or not (y <= 1.45e-169): tmp = x + (0.3333333333333333 * (((t / y) - y) / z)) else: tmp = x + (0.3333333333333333 * ((t / z) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -1.65e-30) || !(y <= 1.45e-169)) tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(Float64(t / y) - y) / z))); else tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(t / z) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -1.65e-30) || ~((y <= 1.45e-169))) tmp = x + (0.3333333333333333 * (((t / y) - y) / z)); else tmp = x + (0.3333333333333333 * ((t / z) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -1.65e-30], N[Not[LessEqual[y, 1.45e-169]], $MachinePrecision]], N[(x + N[(0.3333333333333333 * N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.65 \cdot 10^{-30} \lor \neg \left(y \leq 1.45 \cdot 10^{-169}\right):\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{\frac{t}{y} - y}{z}\\
\mathbf{else}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\end{array}
\end{array}
if y < -1.6500000000000001e-30 or 1.4500000000000001e-169 < y Initial program 98.7%
sub-neg98.7%
associate-+l+98.7%
+-commutative98.7%
remove-double-neg98.7%
distribute-frac-neg98.7%
distribute-neg-in98.7%
remove-double-neg98.7%
sub-neg98.7%
neg-mul-198.7%
times-frac99.2%
distribute-frac-neg99.2%
neg-mul-199.2%
*-commutative99.2%
associate-/l*99.1%
*-commutative99.1%
Simplified99.7%
Taylor expanded in z around 0 99.7%
if -1.6500000000000001e-30 < y < 1.4500000000000001e-169Initial program 91.9%
sub-neg91.9%
associate-+l+91.9%
+-commutative91.9%
remove-double-neg91.9%
distribute-frac-neg91.9%
distribute-neg-in91.9%
remove-double-neg91.9%
sub-neg91.9%
neg-mul-191.9%
times-frac88.9%
distribute-frac-neg88.9%
neg-mul-188.9%
*-commutative88.9%
associate-/l*88.9%
*-commutative88.9%
Simplified88.9%
Taylor expanded in t around inf 91.8%
div-inv90.0%
*-commutative90.0%
Applied egg-rr90.0%
associate-*r/91.8%
*-rgt-identity91.8%
associate-/r*98.7%
Simplified98.7%
Final simplification99.3%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -3.25e-31)
(+ x (* t_1 (/ 0.3333333333333333 z)))
(if (<= y 6.9e-171)
(+ x (* 0.3333333333333333 (/ (/ t z) 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 <= -3.25e-31) {
tmp = x + (t_1 * (0.3333333333333333 / z));
} else if (y <= 6.9e-171) {
tmp = x + (0.3333333333333333 * ((t / z) / 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 <= (-3.25d-31)) then
tmp = x + (t_1 * (0.3333333333333333d0 / z))
else if (y <= 6.9d-171) then
tmp = x + (0.3333333333333333d0 * ((t / z) / 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 <= -3.25e-31) {
tmp = x + (t_1 * (0.3333333333333333 / z));
} else if (y <= 6.9e-171) {
tmp = x + (0.3333333333333333 * ((t / z) / 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 <= -3.25e-31: tmp = x + (t_1 * (0.3333333333333333 / z)) elif y <= 6.9e-171: tmp = x + (0.3333333333333333 * ((t / z) / 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 <= -3.25e-31) tmp = Float64(x + Float64(t_1 * Float64(0.3333333333333333 / z))); elseif (y <= 6.9e-171) tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(t / z) / 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 <= -3.25e-31) tmp = x + (t_1 * (0.3333333333333333 / z)); elseif (y <= 6.9e-171) tmp = x + (0.3333333333333333 * ((t / z) / 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, -3.25e-31], N[(x + N[(t$95$1 * N[(0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.9e-171], N[(x + N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $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 -3.25 \cdot 10^{-31}:\\
\;\;\;\;x + t\_1 \cdot \frac{0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 6.9 \cdot 10^{-171}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\mathbf{else}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{t\_1}{z}\\
\end{array}
\end{array}
if y < -3.24999999999999983e-31Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
if -3.24999999999999983e-31 < y < 6.89999999999999979e-171Initial program 91.9%
sub-neg91.9%
associate-+l+91.9%
+-commutative91.9%
remove-double-neg91.9%
distribute-frac-neg91.9%
distribute-neg-in91.9%
remove-double-neg91.9%
sub-neg91.9%
neg-mul-191.9%
times-frac88.9%
distribute-frac-neg88.9%
neg-mul-188.9%
*-commutative88.9%
associate-/l*88.9%
*-commutative88.9%
Simplified88.9%
Taylor expanded in t around inf 91.8%
div-inv90.0%
*-commutative90.0%
Applied egg-rr90.0%
associate-*r/91.8%
*-rgt-identity91.8%
associate-/r*98.7%
Simplified98.7%
if 6.89999999999999979e-171 < y Initial program 97.9%
sub-neg97.9%
associate-+l+97.9%
+-commutative97.9%
remove-double-neg97.9%
distribute-frac-neg97.9%
distribute-neg-in97.9%
remove-double-neg97.9%
sub-neg97.9%
neg-mul-197.9%
times-frac98.9%
distribute-frac-neg98.9%
neg-mul-198.9%
*-commutative98.9%
associate-/l*98.7%
*-commutative98.7%
Simplified99.6%
Taylor expanded in z around 0 99.7%
Final simplification99.4%
(FPCore (x y z t) :precision binary64 (if (or (<= y -3.4e+34) (not (<= y 6.2e+43))) (- x (/ y (* z 3.0))) (+ x (* 0.3333333333333333 (/ (/ t z) y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.4e+34) || !(y <= 6.2e+43)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
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 <= (-3.4d+34)) .or. (.not. (y <= 6.2d+43))) then
tmp = x - (y / (z * 3.0d0))
else
tmp = x + (0.3333333333333333d0 * ((t / z) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.4e+34) || !(y <= 6.2e+43)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * ((t / z) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -3.4e+34) or not (y <= 6.2e+43): tmp = x - (y / (z * 3.0)) else: tmp = x + (0.3333333333333333 * ((t / z) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -3.4e+34) || !(y <= 6.2e+43)) tmp = Float64(x - Float64(y / Float64(z * 3.0))); else tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(t / z) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -3.4e+34) || ~((y <= 6.2e+43))) tmp = x - (y / (z * 3.0)); else tmp = x + (0.3333333333333333 * ((t / z) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -3.4e+34], N[Not[LessEqual[y, 6.2e+43]], $MachinePrecision]], N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.4 \cdot 10^{+34} \lor \neg \left(y \leq 6.2 \cdot 10^{+43}\right):\\
\;\;\;\;x - \frac{y}{z \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\end{array}
\end{array}
if y < -3.3999999999999999e34 or 6.2000000000000003e43 < y Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
*-commutative99.7%
clear-num99.6%
div-inv99.7%
metadata-eval99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in t around 0 94.7%
neg-mul-194.7%
Simplified94.7%
if -3.3999999999999999e34 < y < 6.2000000000000003e43Initial program 93.8%
sub-neg93.8%
associate-+l+93.8%
+-commutative93.8%
remove-double-neg93.8%
distribute-frac-neg93.8%
distribute-neg-in93.8%
remove-double-neg93.8%
sub-neg93.8%
neg-mul-193.8%
times-frac92.6%
distribute-frac-neg92.6%
neg-mul-192.6%
*-commutative92.6%
associate-/l*92.6%
*-commutative92.6%
Simplified93.2%
Taylor expanded in t around inf 89.7%
div-inv88.6%
*-commutative88.6%
Applied egg-rr88.6%
associate-*r/89.7%
*-rgt-identity89.7%
associate-/r*94.5%
Simplified94.5%
Final simplification94.6%
(FPCore (x y z t) :precision binary64 (if (or (<= y -2.8e+33) (not (<= y 1e+42))) (- x (/ y (* z 3.0))) (+ x (* 0.3333333333333333 (/ t (* y z))))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.8e+33) || !(y <= 1e+42)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * (t / (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 <= (-2.8d+33)) .or. (.not. (y <= 1d+42))) then
tmp = x - (y / (z * 3.0d0))
else
tmp = x + (0.3333333333333333d0 * (t / (y * z)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.8e+33) || !(y <= 1e+42)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = x + (0.3333333333333333 * (t / (y * z)));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -2.8e+33) or not (y <= 1e+42): tmp = x - (y / (z * 3.0)) else: tmp = x + (0.3333333333333333 * (t / (y * z))) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -2.8e+33) || !(y <= 1e+42)) tmp = Float64(x - Float64(y / Float64(z * 3.0))); else tmp = Float64(x + Float64(0.3333333333333333 * Float64(t / Float64(y * z)))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -2.8e+33) || ~((y <= 1e+42))) tmp = x - (y / (z * 3.0)); else tmp = x + (0.3333333333333333 * (t / (y * z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -2.8e+33], N[Not[LessEqual[y, 1e+42]], $MachinePrecision]], N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(0.3333333333333333 * N[(t / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.8 \cdot 10^{+33} \lor \neg \left(y \leq 10^{+42}\right):\\
\;\;\;\;x - \frac{y}{z \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{t}{y \cdot z}\\
\end{array}
\end{array}
if y < -2.8000000000000001e33 or 1.00000000000000004e42 < y Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
*-commutative99.7%
clear-num99.6%
div-inv99.7%
metadata-eval99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in t around 0 94.7%
neg-mul-194.7%
Simplified94.7%
if -2.8000000000000001e33 < y < 1.00000000000000004e42Initial program 93.8%
sub-neg93.8%
associate-+l+93.8%
+-commutative93.8%
remove-double-neg93.8%
distribute-frac-neg93.8%
distribute-neg-in93.8%
remove-double-neg93.8%
sub-neg93.8%
neg-mul-193.8%
times-frac92.6%
distribute-frac-neg92.6%
neg-mul-192.6%
*-commutative92.6%
associate-/l*92.6%
*-commutative92.6%
Simplified93.2%
Taylor expanded in t around inf 89.7%
Final simplification91.9%
(FPCore (x y z t) :precision binary64 (if (or (<= y -2.2e-31) (not (<= y 4.5e-44))) (- x (/ y (* z 3.0))) (/ (* 0.3333333333333333 (/ t z)) y)))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.2e-31) || !(y <= 4.5e-44)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = (0.3333333333333333 * (t / z)) / y;
}
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 <= (-2.2d-31)) .or. (.not. (y <= 4.5d-44))) then
tmp = x - (y / (z * 3.0d0))
else
tmp = (0.3333333333333333d0 * (t / z)) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.2e-31) || !(y <= 4.5e-44)) {
tmp = x - (y / (z * 3.0));
} else {
tmp = (0.3333333333333333 * (t / z)) / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -2.2e-31) or not (y <= 4.5e-44): tmp = x - (y / (z * 3.0)) else: tmp = (0.3333333333333333 * (t / z)) / y return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -2.2e-31) || !(y <= 4.5e-44)) tmp = Float64(x - Float64(y / Float64(z * 3.0))); else tmp = Float64(Float64(0.3333333333333333 * Float64(t / z)) / y); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -2.2e-31) || ~((y <= 4.5e-44))) tmp = x - (y / (z * 3.0)); else tmp = (0.3333333333333333 * (t / z)) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -2.2e-31], N[Not[LessEqual[y, 4.5e-44]], $MachinePrecision]], N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.3333333333333333 * N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.2 \cdot 10^{-31} \lor \neg \left(y \leq 4.5 \cdot 10^{-44}\right):\\
\;\;\;\;x - \frac{y}{z \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.3333333333333333 \cdot \frac{t}{z}}{y}\\
\end{array}
\end{array}
if y < -2.2000000000000001e-31 or 4.4999999999999999e-44 < y Initial program 99.2%
sub-neg99.2%
associate-+l+99.2%
+-commutative99.2%
remove-double-neg99.2%
distribute-frac-neg99.2%
distribute-neg-in99.2%
remove-double-neg99.2%
sub-neg99.2%
neg-mul-199.2%
times-frac99.1%
distribute-frac-neg99.1%
neg-mul-199.1%
*-commutative99.1%
associate-/l*99.0%
*-commutative99.0%
Simplified99.7%
*-commutative99.7%
clear-num99.7%
div-inv99.7%
metadata-eval99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in t around 0 88.5%
neg-mul-188.5%
Simplified88.5%
if -2.2000000000000001e-31 < y < 4.4999999999999999e-44Initial program 93.1%
+-commutative93.1%
associate-+r-93.1%
sub-neg93.1%
associate-*l*93.1%
*-commutative93.1%
distribute-frac-neg293.1%
distribute-rgt-neg-in93.1%
metadata-eval93.1%
Simplified93.1%
Taylor expanded in y around 0 95.7%
Taylor expanded in t around inf 66.1%
Final simplification78.3%
(FPCore (x y z t) :precision binary64 (if (or (<= y -2.3e-29) (not (<= y 3.6e-45))) (+ x (* y (/ -0.3333333333333333 z))) (* 0.3333333333333333 (/ (/ t z) y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.3e-29) || !(y <= 3.6e-45)) {
tmp = x + (y * (-0.3333333333333333 / z));
} else {
tmp = 0.3333333333333333 * ((t / z) / y);
}
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 <= (-2.3d-29)) .or. (.not. (y <= 3.6d-45))) then
tmp = x + (y * ((-0.3333333333333333d0) / z))
else
tmp = 0.3333333333333333d0 * ((t / z) / y)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.3e-29) || !(y <= 3.6e-45)) {
tmp = x + (y * (-0.3333333333333333 / z));
} else {
tmp = 0.3333333333333333 * ((t / z) / y);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -2.3e-29) or not (y <= 3.6e-45): tmp = x + (y * (-0.3333333333333333 / z)) else: tmp = 0.3333333333333333 * ((t / z) / y) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -2.3e-29) || !(y <= 3.6e-45)) tmp = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))); else tmp = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -2.3e-29) || ~((y <= 3.6e-45))) tmp = x + (y * (-0.3333333333333333 / z)); else tmp = 0.3333333333333333 * ((t / z) / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -2.3e-29], N[Not[LessEqual[y, 3.6e-45]], $MachinePrecision]], N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.3 \cdot 10^{-29} \lor \neg \left(y \leq 3.6 \cdot 10^{-45}\right):\\
\;\;\;\;x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\end{array}
\end{array}
if y < -2.29999999999999991e-29 or 3.60000000000000001e-45 < y Initial program 99.2%
sub-neg99.2%
associate-+l+99.2%
+-commutative99.2%
remove-double-neg99.2%
distribute-frac-neg99.2%
distribute-neg-in99.2%
remove-double-neg99.2%
sub-neg99.2%
neg-mul-199.2%
times-frac99.1%
distribute-frac-neg99.1%
neg-mul-199.1%
*-commutative99.1%
associate-/l*99.0%
*-commutative99.0%
Simplified99.7%
Taylor expanded in t around 0 88.3%
associate-*r/88.3%
*-commutative88.3%
associate-*r/88.3%
Simplified88.3%
if -2.29999999999999991e-29 < y < 3.60000000000000001e-45Initial program 93.1%
+-commutative93.1%
associate-+r-93.1%
sub-neg93.1%
associate-*l*93.1%
*-commutative93.1%
distribute-frac-neg293.1%
distribute-rgt-neg-in93.1%
metadata-eval93.1%
Simplified93.1%
*-un-lft-identity93.1%
times-frac93.1%
Applied egg-rr93.1%
associate-*r/93.1%
associate-*l/93.1%
*-lft-identity93.1%
Simplified93.1%
Taylor expanded in t around inf 62.9%
associate-/l/66.1%
Simplified66.1%
Final simplification78.2%
(FPCore (x y z t)
:precision binary64
(if (<= y -1.55e-27)
(+ x (* y (/ -0.3333333333333333 z)))
(if (<= y 5.5e-43)
(/ (* 0.3333333333333333 (/ t z)) y)
(- x (* 0.3333333333333333 (/ y z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.55e-27) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 5.5e-43) {
tmp = (0.3333333333333333 * (t / z)) / y;
} 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 <= (-1.55d-27)) then
tmp = x + (y * ((-0.3333333333333333d0) / z))
else if (y <= 5.5d-43) then
tmp = (0.3333333333333333d0 * (t / z)) / y
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 <= -1.55e-27) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 5.5e-43) {
tmp = (0.3333333333333333 * (t / z)) / y;
} else {
tmp = x - (0.3333333333333333 * (y / z));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -1.55e-27: tmp = x + (y * (-0.3333333333333333 / z)) elif y <= 5.5e-43: tmp = (0.3333333333333333 * (t / z)) / y else: tmp = x - (0.3333333333333333 * (y / z)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -1.55e-27) tmp = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))); elseif (y <= 5.5e-43) tmp = Float64(Float64(0.3333333333333333 * Float64(t / z)) / y); 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 <= -1.55e-27) tmp = x + (y * (-0.3333333333333333 / z)); elseif (y <= 5.5e-43) tmp = (0.3333333333333333 * (t / z)) / y; else tmp = x - (0.3333333333333333 * (y / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -1.55e-27], N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5.5e-43], N[(N[(0.3333333333333333 * N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.55 \cdot 10^{-27}:\\
\;\;\;\;x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 5.5 \cdot 10^{-43}:\\
\;\;\;\;\frac{0.3333333333333333 \cdot \frac{t}{z}}{y}\\
\mathbf{else}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\end{array}
\end{array}
if y < -1.5499999999999999e-27Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 91.1%
associate-*r/91.1%
*-commutative91.1%
associate-*r/91.2%
Simplified91.2%
if -1.5499999999999999e-27 < y < 5.50000000000000013e-43Initial program 93.1%
+-commutative93.1%
associate-+r-93.1%
sub-neg93.1%
associate-*l*93.1%
*-commutative93.1%
distribute-frac-neg293.1%
distribute-rgt-neg-in93.1%
metadata-eval93.1%
Simplified93.1%
Taylor expanded in y around 0 95.7%
Taylor expanded in t around inf 66.1%
if 5.50000000000000013e-43 < y Initial program 98.5%
+-commutative98.5%
associate-+r-98.5%
sub-neg98.5%
associate-*l*98.5%
*-commutative98.5%
distribute-frac-neg298.5%
distribute-rgt-neg-in98.5%
metadata-eval98.5%
Simplified98.5%
Taylor expanded in t around 0 85.7%
metadata-eval85.7%
cancel-sign-sub-inv85.7%
Simplified85.7%
(FPCore (x y z t)
:precision binary64
(if (<= y -2e-30)
(+ x (* y (/ -0.3333333333333333 z)))
(if (<= y 6.6e-46)
(* 0.3333333333333333 (/ (/ t z) y))
(- x (* 0.3333333333333333 (/ y z))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -2e-30) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 6.6e-46) {
tmp = 0.3333333333333333 * ((t / z) / y);
} 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 <= (-2d-30)) then
tmp = x + (y * ((-0.3333333333333333d0) / z))
else if (y <= 6.6d-46) then
tmp = 0.3333333333333333d0 * ((t / z) / y)
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 <= -2e-30) {
tmp = x + (y * (-0.3333333333333333 / z));
} else if (y <= 6.6e-46) {
tmp = 0.3333333333333333 * ((t / z) / y);
} else {
tmp = x - (0.3333333333333333 * (y / z));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -2e-30: tmp = x + (y * (-0.3333333333333333 / z)) elif y <= 6.6e-46: tmp = 0.3333333333333333 * ((t / z) / y) else: tmp = x - (0.3333333333333333 * (y / z)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -2e-30) tmp = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))); elseif (y <= 6.6e-46) tmp = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)); 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 <= -2e-30) tmp = x + (y * (-0.3333333333333333 / z)); elseif (y <= 6.6e-46) tmp = 0.3333333333333333 * ((t / z) / y); else tmp = x - (0.3333333333333333 * (y / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -2e-30], N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.6e-46], N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2 \cdot 10^{-30}:\\
\;\;\;\;x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 6.6 \cdot 10^{-46}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\mathbf{else}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\end{array}
\end{array}
if y < -2e-30Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 91.1%
associate-*r/91.1%
*-commutative91.1%
associate-*r/91.2%
Simplified91.2%
if -2e-30 < y < 6.60000000000000027e-46Initial program 93.1%
+-commutative93.1%
associate-+r-93.1%
sub-neg93.1%
associate-*l*93.1%
*-commutative93.1%
distribute-frac-neg293.1%
distribute-rgt-neg-in93.1%
metadata-eval93.1%
Simplified93.1%
*-un-lft-identity93.1%
times-frac93.1%
Applied egg-rr93.1%
associate-*r/93.1%
associate-*l/93.1%
*-lft-identity93.1%
Simplified93.1%
Taylor expanded in t around inf 62.9%
associate-/l/66.1%
Simplified66.1%
if 6.60000000000000027e-46 < y Initial program 98.5%
+-commutative98.5%
associate-+r-98.5%
sub-neg98.5%
associate-*l*98.5%
*-commutative98.5%
distribute-frac-neg298.5%
distribute-rgt-neg-in98.5%
metadata-eval98.5%
Simplified98.5%
Taylor expanded in t around 0 85.7%
metadata-eval85.7%
cancel-sign-sub-inv85.7%
Simplified85.7%
(FPCore (x y z t) :precision binary64 (if (or (<= y -1.4e+47) (not (<= y 1.65e+22))) (* -0.3333333333333333 (/ y z)) x))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.4e+47) || !(y <= 1.65e+22)) {
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 ((y <= (-1.4d+47)) .or. (.not. (y <= 1.65d+22))) 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 ((y <= -1.4e+47) || !(y <= 1.65e+22)) {
tmp = -0.3333333333333333 * (y / z);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -1.4e+47) or not (y <= 1.65e+22): tmp = -0.3333333333333333 * (y / z) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -1.4e+47) || !(y <= 1.65e+22)) 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 ((y <= -1.4e+47) || ~((y <= 1.65e+22))) tmp = -0.3333333333333333 * (y / z); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -1.4e+47], N[Not[LessEqual[y, 1.65e+22]], $MachinePrecision]], N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.4 \cdot 10^{+47} \lor \neg \left(y \leq 1.65 \cdot 10^{+22}\right):\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -1.39999999999999994e47 or 1.6499999999999999e22 < y Initial program 99.0%
+-commutative99.0%
associate-+r-99.0%
sub-neg99.0%
associate-*l*99.0%
*-commutative99.0%
distribute-frac-neg299.0%
distribute-rgt-neg-in99.0%
metadata-eval99.0%
Simplified99.0%
*-un-lft-identity99.0%
times-frac98.8%
Applied egg-rr98.8%
associate-*r/98.9%
associate-*l/98.9%
*-lft-identity98.9%
Simplified98.9%
Taylor expanded in y around inf 74.1%
if -1.39999999999999994e47 < y < 1.6499999999999999e22Initial program 94.3%
+-commutative94.3%
associate-+r-94.3%
sub-neg94.3%
associate-*l*94.3%
*-commutative94.3%
distribute-frac-neg294.3%
distribute-rgt-neg-in94.3%
metadata-eval94.3%
Simplified94.3%
Taylor expanded in z around inf 37.4%
Final simplification53.9%
(FPCore (x y z t) :precision binary64 (if (<= y -2.15e+51) (* y (/ -0.3333333333333333 z)) (if (<= y 48000000000000.0) x (/ (* y -0.3333333333333333) z))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -2.15e+51) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 48000000000000.0) {
tmp = x;
} else {
tmp = (y * -0.3333333333333333) / 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 <= (-2.15d+51)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 48000000000000.0d0) then
tmp = x
else
tmp = (y * (-0.3333333333333333d0)) / z
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -2.15e+51) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 48000000000000.0) {
tmp = x;
} else {
tmp = (y * -0.3333333333333333) / z;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -2.15e+51: tmp = y * (-0.3333333333333333 / z) elif y <= 48000000000000.0: tmp = x else: tmp = (y * -0.3333333333333333) / z return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -2.15e+51) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 48000000000000.0) tmp = x; else tmp = Float64(Float64(y * -0.3333333333333333) / z); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -2.15e+51) tmp = y * (-0.3333333333333333 / z); elseif (y <= 48000000000000.0) tmp = x; else tmp = (y * -0.3333333333333333) / z; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -2.15e+51], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 48000000000000.0], x, N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.15 \cdot 10^{+51}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 48000000000000:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -2.1499999999999999e51Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-lft-identity99.7%
Simplified99.7%
Taylor expanded in y around inf 74.5%
*-commutative74.5%
associate-*l/74.5%
associate-/l*74.6%
Simplified74.6%
if -2.1499999999999999e51 < y < 4.8e13Initial program 94.3%
+-commutative94.3%
associate-+r-94.3%
sub-neg94.3%
associate-*l*94.3%
*-commutative94.3%
distribute-frac-neg294.3%
distribute-rgt-neg-in94.3%
metadata-eval94.3%
Simplified94.3%
Taylor expanded in z around inf 37.4%
if 4.8e13 < y Initial program 98.2%
+-commutative98.2%
associate-+r-98.2%
sub-neg98.2%
associate-*l*98.2%
*-commutative98.2%
distribute-frac-neg298.2%
distribute-rgt-neg-in98.2%
metadata-eval98.2%
Simplified98.2%
*-un-lft-identity98.2%
times-frac98.0%
Applied egg-rr98.0%
associate-*r/98.1%
associate-*l/98.1%
*-lft-identity98.1%
Simplified98.1%
Taylor expanded in y around inf 73.6%
*-commutative73.6%
associate-*l/73.7%
Applied egg-rr73.7%
(FPCore (x y z t) :precision binary64 (if (<= y -4.8e+51) (* y (/ -0.3333333333333333 z)) (if (<= y 145000000000.0) x (* -0.3333333333333333 (/ y z)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -4.8e+51) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 145000000000.0) {
tmp = x;
} else {
tmp = -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.8d+51)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 145000000000.0d0) then
tmp = x
else
tmp = (-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.8e+51) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 145000000000.0) {
tmp = x;
} else {
tmp = -0.3333333333333333 * (y / z);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -4.8e+51: tmp = y * (-0.3333333333333333 / z) elif y <= 145000000000.0: tmp = x else: tmp = -0.3333333333333333 * (y / z) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -4.8e+51) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 145000000000.0) tmp = x; else tmp = Float64(-0.3333333333333333 * Float64(y / z)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -4.8e+51) tmp = y * (-0.3333333333333333 / z); elseif (y <= 145000000000.0) tmp = x; else tmp = -0.3333333333333333 * (y / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -4.8e+51], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 145000000000.0], x, N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.8 \cdot 10^{+51}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 145000000000:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{z}\\
\end{array}
\end{array}
if y < -4.7999999999999997e51Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-lft-identity99.7%
Simplified99.7%
Taylor expanded in y around inf 74.5%
*-commutative74.5%
associate-*l/74.5%
associate-/l*74.6%
Simplified74.6%
if -4.7999999999999997e51 < y < 1.45e11Initial program 94.3%
+-commutative94.3%
associate-+r-94.3%
sub-neg94.3%
associate-*l*94.3%
*-commutative94.3%
distribute-frac-neg294.3%
distribute-rgt-neg-in94.3%
metadata-eval94.3%
Simplified94.3%
Taylor expanded in z around inf 37.4%
if 1.45e11 < y Initial program 98.2%
+-commutative98.2%
associate-+r-98.2%
sub-neg98.2%
associate-*l*98.2%
*-commutative98.2%
distribute-frac-neg298.2%
distribute-rgt-neg-in98.2%
metadata-eval98.2%
Simplified98.2%
*-un-lft-identity98.2%
times-frac98.0%
Applied egg-rr98.0%
associate-*r/98.1%
associate-*l/98.1%
*-lft-identity98.1%
Simplified98.1%
Taylor expanded in y around inf 73.6%
(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 96.4%
+-commutative96.4%
associate-+r-96.4%
sub-neg96.4%
associate-*l*96.4%
*-commutative96.4%
distribute-frac-neg296.4%
distribute-rgt-neg-in96.4%
metadata-eval96.4%
Simplified96.4%
Taylor expanded in z around inf 29.6%
(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 2024106
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:alt
(+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y))
(+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))