
(FPCore (x y z t a) :precision binary64 (- (+ x y) (/ (* (- z t) y) (- a t))))
double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = (x + y) - (((z - t) * y) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
def code(x, y, z, t, a): return (x + y) - (((z - t) * y) / (a - t))
function code(x, y, z, t, a) return Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = (x + y) - (((z - t) * y) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (- (+ x y) (/ (* (- z t) y) (- a t))))
double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = (x + y) - (((z - t) * y) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
def code(x, y, z, t, a): return (x + y) - (((z - t) * y) / (a - t))
function code(x, y, z, t, a) return Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = (x + y) - (((z - t) * y) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}
\end{array}
(FPCore (x y z t a)
:precision binary64
(if (<= t -3.4e+102)
(- x (* y (/ (- a z) t)))
(if (<= t 2.8e+47)
(+ (+ x y) (/ (- z t) (/ (- t a) y)))
(+ (- x (* a (/ y t))) (* y (/ z t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -3.4e+102) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 2.8e+47) {
tmp = (x + y) + ((z - t) / ((t - a) / y));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-3.4d+102)) then
tmp = x - (y * ((a - z) / t))
else if (t <= 2.8d+47) then
tmp = (x + y) + ((z - t) / ((t - a) / y))
else
tmp = (x - (a * (y / t))) + (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -3.4e+102) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 2.8e+47) {
tmp = (x + y) + ((z - t) / ((t - a) / y));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -3.4e+102: tmp = x - (y * ((a - z) / t)) elif t <= 2.8e+47: tmp = (x + y) + ((z - t) / ((t - a) / y)) else: tmp = (x - (a * (y / t))) + (y * (z / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -3.4e+102) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= 2.8e+47) tmp = Float64(Float64(x + y) + Float64(Float64(z - t) / Float64(Float64(t - a) / y))); else tmp = Float64(Float64(x - Float64(a * Float64(y / t))) + Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -3.4e+102) tmp = x - (y * ((a - z) / t)); elseif (t <= 2.8e+47) tmp = (x + y) + ((z - t) / ((t - a) / y)); else tmp = (x - (a * (y / t))) + (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -3.4e+102], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.8e+47], N[(N[(x + y), $MachinePrecision] + N[(N[(z - t), $MachinePrecision] / N[(N[(t - a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x - N[(a * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -3.4 \cdot 10^{+102}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq 2.8 \cdot 10^{+47}:\\
\;\;\;\;\left(x + y\right) + \frac{z - t}{\frac{t - a}{y}}\\
\mathbf{else}:\\
\;\;\;\;\left(x - a \cdot \frac{y}{t}\right) + y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if t < -3.4e102Initial program 44.6%
Taylor expanded in y around 0 44.6%
associate-*l/56.4%
Simplified56.4%
*-commutative56.4%
clear-num56.8%
un-div-inv56.4%
Applied egg-rr56.4%
Taylor expanded in t around -inf 67.9%
mul-1-neg67.9%
unsub-neg67.9%
*-commutative67.9%
cancel-sign-sub-inv67.9%
neg-mul-167.9%
distribute-rgt-in68.0%
associate-/l*90.7%
neg-mul-190.7%
sub-neg90.7%
Simplified90.7%
if -3.4e102 < t < 2.79999999999999988e47Initial program 90.9%
Taylor expanded in y around 0 90.9%
associate-*l/94.1%
Simplified94.1%
*-commutative94.1%
clear-num94.1%
un-div-inv94.7%
Applied egg-rr94.7%
if 2.79999999999999988e47 < t Initial program 53.6%
Taylor expanded in y around 0 53.6%
associate-*l/62.0%
Simplified62.0%
Taylor expanded in t around inf 76.7%
sub-neg76.7%
mul-1-neg76.7%
unsub-neg76.7%
associate-/l*81.7%
mul-1-neg81.7%
remove-double-neg81.7%
associate-/l*86.8%
Simplified86.8%
Final simplification92.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* y (/ (- a z) t)))))
(if (<= a -2.6e+49)
(+ x y)
(if (<= a -2.2e-63)
t_1
(if (<= a -9.5e-67)
x
(if (<= a -1.86e-227)
(+ x (* y (/ z t)))
(if (<= a 6.5e-90)
(- x (/ (* y (- a z)) t))
(if (or (<= a 2.35e+26)
(and (not (<= a 2.2e+132)) (<= a 2.55e+154)))
t_1
(+ x y)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y * ((a - z) / t));
double tmp;
if (a <= -2.6e+49) {
tmp = x + y;
} else if (a <= -2.2e-63) {
tmp = t_1;
} else if (a <= -9.5e-67) {
tmp = x;
} else if (a <= -1.86e-227) {
tmp = x + (y * (z / t));
} else if (a <= 6.5e-90) {
tmp = x - ((y * (a - z)) / t);
} else if ((a <= 2.35e+26) || (!(a <= 2.2e+132) && (a <= 2.55e+154))) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y * ((a - z) / t))
if (a <= (-2.6d+49)) then
tmp = x + y
else if (a <= (-2.2d-63)) then
tmp = t_1
else if (a <= (-9.5d-67)) then
tmp = x
else if (a <= (-1.86d-227)) then
tmp = x + (y * (z / t))
else if (a <= 6.5d-90) then
tmp = x - ((y * (a - z)) / t)
else if ((a <= 2.35d+26) .or. (.not. (a <= 2.2d+132)) .and. (a <= 2.55d+154)) then
tmp = t_1
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y * ((a - z) / t));
double tmp;
if (a <= -2.6e+49) {
tmp = x + y;
} else if (a <= -2.2e-63) {
tmp = t_1;
} else if (a <= -9.5e-67) {
tmp = x;
} else if (a <= -1.86e-227) {
tmp = x + (y * (z / t));
} else if (a <= 6.5e-90) {
tmp = x - ((y * (a - z)) / t);
} else if ((a <= 2.35e+26) || (!(a <= 2.2e+132) && (a <= 2.55e+154))) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y * ((a - z) / t)) tmp = 0 if a <= -2.6e+49: tmp = x + y elif a <= -2.2e-63: tmp = t_1 elif a <= -9.5e-67: tmp = x elif a <= -1.86e-227: tmp = x + (y * (z / t)) elif a <= 6.5e-90: tmp = x - ((y * (a - z)) / t) elif (a <= 2.35e+26) or (not (a <= 2.2e+132) and (a <= 2.55e+154)): tmp = t_1 else: tmp = x + y return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y * Float64(Float64(a - z) / t))) tmp = 0.0 if (a <= -2.6e+49) tmp = Float64(x + y); elseif (a <= -2.2e-63) tmp = t_1; elseif (a <= -9.5e-67) tmp = x; elseif (a <= -1.86e-227) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= 6.5e-90) tmp = Float64(x - Float64(Float64(y * Float64(a - z)) / t)); elseif ((a <= 2.35e+26) || (!(a <= 2.2e+132) && (a <= 2.55e+154))) tmp = t_1; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y * ((a - z) / t)); tmp = 0.0; if (a <= -2.6e+49) tmp = x + y; elseif (a <= -2.2e-63) tmp = t_1; elseif (a <= -9.5e-67) tmp = x; elseif (a <= -1.86e-227) tmp = x + (y * (z / t)); elseif (a <= 6.5e-90) tmp = x - ((y * (a - z)) / t); elseif ((a <= 2.35e+26) || (~((a <= 2.2e+132)) && (a <= 2.55e+154))) tmp = t_1; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.6e+49], N[(x + y), $MachinePrecision], If[LessEqual[a, -2.2e-63], t$95$1, If[LessEqual[a, -9.5e-67], x, If[LessEqual[a, -1.86e-227], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 6.5e-90], N[(x - N[(N[(y * N[(a - z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[a, 2.35e+26], And[N[Not[LessEqual[a, 2.2e+132]], $MachinePrecision], LessEqual[a, 2.55e+154]]], t$95$1, N[(x + y), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - y \cdot \frac{a - z}{t}\\
\mathbf{if}\;a \leq -2.6 \cdot 10^{+49}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq -2.2 \cdot 10^{-63}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq -9.5 \cdot 10^{-67}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -1.86 \cdot 10^{-227}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq 6.5 \cdot 10^{-90}:\\
\;\;\;\;x - \frac{y \cdot \left(a - z\right)}{t}\\
\mathbf{elif}\;a \leq 2.35 \cdot 10^{+26} \lor \neg \left(a \leq 2.2 \cdot 10^{+132}\right) \land a \leq 2.55 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -2.59999999999999989e49 or 2.3499999999999999e26 < a < 2.19999999999999989e132 or 2.55e154 < a Initial program 76.7%
Taylor expanded in a around inf 77.1%
+-commutative77.1%
Simplified77.1%
if -2.59999999999999989e49 < a < -2.2e-63 or 6.4999999999999996e-90 < a < 2.3499999999999999e26 or 2.19999999999999989e132 < a < 2.55e154Initial program 71.4%
Taylor expanded in y around 0 71.4%
associate-*l/72.1%
Simplified72.1%
*-commutative72.1%
clear-num72.1%
un-div-inv72.0%
Applied egg-rr72.0%
Taylor expanded in t around -inf 66.4%
mul-1-neg66.4%
unsub-neg66.4%
*-commutative66.4%
cancel-sign-sub-inv66.4%
neg-mul-166.4%
distribute-rgt-in66.4%
associate-/l*73.7%
neg-mul-173.7%
sub-neg73.7%
Simplified73.7%
if -2.2e-63 < a < -9.4999999999999994e-67Initial program 100.0%
Taylor expanded in x around inf 100.0%
if -9.4999999999999994e-67 < a < -1.86000000000000011e-227Initial program 74.0%
Taylor expanded in y around 0 74.0%
associate-*l/82.3%
Simplified82.3%
*-commutative82.3%
clear-num82.2%
un-div-inv82.4%
Applied egg-rr82.4%
Taylor expanded in t around -inf 70.7%
mul-1-neg70.7%
unsub-neg70.7%
*-commutative70.7%
cancel-sign-sub-inv70.7%
neg-mul-170.7%
distribute-rgt-in70.7%
associate-/l*76.5%
neg-mul-176.5%
sub-neg76.5%
Simplified76.5%
Taylor expanded in a around 0 79.1%
neg-mul-179.1%
distribute-neg-frac279.1%
Simplified79.1%
sub-neg79.1%
+-commutative79.1%
distribute-lft-neg-in79.1%
add-sqr-sqrt63.8%
sqrt-unprod63.3%
sqr-neg63.3%
sqrt-unprod12.4%
add-sqr-sqrt60.6%
add-sqr-sqrt29.9%
sqrt-unprod66.7%
sqr-neg66.7%
sqrt-unprod42.6%
add-sqr-sqrt79.1%
Applied egg-rr79.1%
if -1.86000000000000011e-227 < a < 6.4999999999999996e-90Initial program 76.4%
Taylor expanded in t around inf 83.7%
associate--l+83.7%
distribute-lft-out--83.7%
div-sub83.7%
mul-1-neg83.7%
unsub-neg83.7%
*-commutative83.7%
distribute-lft-out--83.7%
Simplified83.7%
Final simplification78.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* z (/ y (- t a)))))
(if (<= z -1.4e+186)
t_1
(if (<= z -6.6e-187)
(+ x y)
(if (<= z -2.8e-293)
x
(if (<= z 1.6e+35)
(+ x y)
(if (<= z 3.5e+135)
t_1
(if (<= z 9.7e+154)
x
(if (or (<= z 6.5e+206) (not (<= z 1e+213)))
t_1
(+ x y))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = z * (y / (t - a));
double tmp;
if (z <= -1.4e+186) {
tmp = t_1;
} else if (z <= -6.6e-187) {
tmp = x + y;
} else if (z <= -2.8e-293) {
tmp = x;
} else if (z <= 1.6e+35) {
tmp = x + y;
} else if (z <= 3.5e+135) {
tmp = t_1;
} else if (z <= 9.7e+154) {
tmp = x;
} else if ((z <= 6.5e+206) || !(z <= 1e+213)) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = z * (y / (t - a))
if (z <= (-1.4d+186)) then
tmp = t_1
else if (z <= (-6.6d-187)) then
tmp = x + y
else if (z <= (-2.8d-293)) then
tmp = x
else if (z <= 1.6d+35) then
tmp = x + y
else if (z <= 3.5d+135) then
tmp = t_1
else if (z <= 9.7d+154) then
tmp = x
else if ((z <= 6.5d+206) .or. (.not. (z <= 1d+213))) then
tmp = t_1
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = z * (y / (t - a));
double tmp;
if (z <= -1.4e+186) {
tmp = t_1;
} else if (z <= -6.6e-187) {
tmp = x + y;
} else if (z <= -2.8e-293) {
tmp = x;
} else if (z <= 1.6e+35) {
tmp = x + y;
} else if (z <= 3.5e+135) {
tmp = t_1;
} else if (z <= 9.7e+154) {
tmp = x;
} else if ((z <= 6.5e+206) || !(z <= 1e+213)) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = z * (y / (t - a)) tmp = 0 if z <= -1.4e+186: tmp = t_1 elif z <= -6.6e-187: tmp = x + y elif z <= -2.8e-293: tmp = x elif z <= 1.6e+35: tmp = x + y elif z <= 3.5e+135: tmp = t_1 elif z <= 9.7e+154: tmp = x elif (z <= 6.5e+206) or not (z <= 1e+213): tmp = t_1 else: tmp = x + y return tmp
function code(x, y, z, t, a) t_1 = Float64(z * Float64(y / Float64(t - a))) tmp = 0.0 if (z <= -1.4e+186) tmp = t_1; elseif (z <= -6.6e-187) tmp = Float64(x + y); elseif (z <= -2.8e-293) tmp = x; elseif (z <= 1.6e+35) tmp = Float64(x + y); elseif (z <= 3.5e+135) tmp = t_1; elseif (z <= 9.7e+154) tmp = x; elseif ((z <= 6.5e+206) || !(z <= 1e+213)) tmp = t_1; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = z * (y / (t - a)); tmp = 0.0; if (z <= -1.4e+186) tmp = t_1; elseif (z <= -6.6e-187) tmp = x + y; elseif (z <= -2.8e-293) tmp = x; elseif (z <= 1.6e+35) tmp = x + y; elseif (z <= 3.5e+135) tmp = t_1; elseif (z <= 9.7e+154) tmp = x; elseif ((z <= 6.5e+206) || ~((z <= 1e+213))) tmp = t_1; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(z * N[(y / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.4e+186], t$95$1, If[LessEqual[z, -6.6e-187], N[(x + y), $MachinePrecision], If[LessEqual[z, -2.8e-293], x, If[LessEqual[z, 1.6e+35], N[(x + y), $MachinePrecision], If[LessEqual[z, 3.5e+135], t$95$1, If[LessEqual[z, 9.7e+154], x, If[Or[LessEqual[z, 6.5e+206], N[Not[LessEqual[z, 1e+213]], $MachinePrecision]], t$95$1, N[(x + y), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \frac{y}{t - a}\\
\mathbf{if}\;z \leq -1.4 \cdot 10^{+186}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -6.6 \cdot 10^{-187}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;z \leq -2.8 \cdot 10^{-293}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 1.6 \cdot 10^{+35}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{+135}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 9.7 \cdot 10^{+154}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 6.5 \cdot 10^{+206} \lor \neg \left(z \leq 10^{+213}\right):\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if z < -1.40000000000000009e186 or 1.59999999999999991e35 < z < 3.5000000000000003e135 or 9.69999999999999987e154 < z < 6.4999999999999995e206 or 9.99999999999999984e212 < z Initial program 78.0%
sub-neg78.0%
+-commutative78.0%
distribute-frac-neg78.0%
distribute-rgt-neg-out78.0%
associate-/l*91.8%
fma-define91.9%
distribute-frac-neg91.9%
distribute-neg-frac291.9%
sub-neg91.9%
distribute-neg-in91.9%
remove-double-neg91.9%
+-commutative91.9%
sub-neg91.9%
Simplified91.9%
Taylor expanded in z around inf 58.1%
associate-/l*64.5%
Simplified64.5%
associate-*r/58.1%
clear-num58.1%
Applied egg-rr58.1%
clear-num58.1%
*-commutative58.1%
associate-/l*74.7%
Applied egg-rr74.7%
if -1.40000000000000009e186 < z < -6.6e-187 or -2.80000000000000025e-293 < z < 1.59999999999999991e35 or 6.4999999999999995e206 < z < 9.99999999999999984e212Initial program 74.9%
Taylor expanded in a around inf 66.1%
+-commutative66.1%
Simplified66.1%
if -6.6e-187 < z < -2.80000000000000025e-293 or 3.5000000000000003e135 < z < 9.69999999999999987e154Initial program 70.6%
Taylor expanded in x around inf 93.8%
Final simplification70.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* y (/ z (- t a)))))
(if (<= z -1.12e+163)
t_1
(if (<= z -1.46e+119)
(- (+ x y) y)
(if (<= z -2.2e-187)
(+ x y)
(if (<= z -4.6e-295)
x
(if (or (<= z 1.95e+37) (and (not (<= z 8.5e+205)) (<= z 9.5e+212)))
(+ x y)
t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * (z / (t - a));
double tmp;
if (z <= -1.12e+163) {
tmp = t_1;
} else if (z <= -1.46e+119) {
tmp = (x + y) - y;
} else if (z <= -2.2e-187) {
tmp = x + y;
} else if (z <= -4.6e-295) {
tmp = x;
} else if ((z <= 1.95e+37) || (!(z <= 8.5e+205) && (z <= 9.5e+212))) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = y * (z / (t - a))
if (z <= (-1.12d+163)) then
tmp = t_1
else if (z <= (-1.46d+119)) then
tmp = (x + y) - y
else if (z <= (-2.2d-187)) then
tmp = x + y
else if (z <= (-4.6d-295)) then
tmp = x
else if ((z <= 1.95d+37) .or. (.not. (z <= 8.5d+205)) .and. (z <= 9.5d+212)) then
tmp = x + y
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = y * (z / (t - a));
double tmp;
if (z <= -1.12e+163) {
tmp = t_1;
} else if (z <= -1.46e+119) {
tmp = (x + y) - y;
} else if (z <= -2.2e-187) {
tmp = x + y;
} else if (z <= -4.6e-295) {
tmp = x;
} else if ((z <= 1.95e+37) || (!(z <= 8.5e+205) && (z <= 9.5e+212))) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * (z / (t - a)) tmp = 0 if z <= -1.12e+163: tmp = t_1 elif z <= -1.46e+119: tmp = (x + y) - y elif z <= -2.2e-187: tmp = x + y elif z <= -4.6e-295: tmp = x elif (z <= 1.95e+37) or (not (z <= 8.5e+205) and (z <= 9.5e+212)): tmp = x + y else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(z / Float64(t - a))) tmp = 0.0 if (z <= -1.12e+163) tmp = t_1; elseif (z <= -1.46e+119) tmp = Float64(Float64(x + y) - y); elseif (z <= -2.2e-187) tmp = Float64(x + y); elseif (z <= -4.6e-295) tmp = x; elseif ((z <= 1.95e+37) || (!(z <= 8.5e+205) && (z <= 9.5e+212))) tmp = Float64(x + y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * (z / (t - a)); tmp = 0.0; if (z <= -1.12e+163) tmp = t_1; elseif (z <= -1.46e+119) tmp = (x + y) - y; elseif (z <= -2.2e-187) tmp = x + y; elseif (z <= -4.6e-295) tmp = x; elseif ((z <= 1.95e+37) || (~((z <= 8.5e+205)) && (z <= 9.5e+212))) tmp = x + y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.12e+163], t$95$1, If[LessEqual[z, -1.46e+119], N[(N[(x + y), $MachinePrecision] - y), $MachinePrecision], If[LessEqual[z, -2.2e-187], N[(x + y), $MachinePrecision], If[LessEqual[z, -4.6e-295], x, If[Or[LessEqual[z, 1.95e+37], And[N[Not[LessEqual[z, 8.5e+205]], $MachinePrecision], LessEqual[z, 9.5e+212]]], N[(x + y), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z}{t - a}\\
\mathbf{if}\;z \leq -1.12 \cdot 10^{+163}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -1.46 \cdot 10^{+119}:\\
\;\;\;\;\left(x + y\right) - y\\
\mathbf{elif}\;z \leq -2.2 \cdot 10^{-187}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;z \leq -4.6 \cdot 10^{-295}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 1.95 \cdot 10^{+37} \lor \neg \left(z \leq 8.5 \cdot 10^{+205}\right) \land z \leq 9.5 \cdot 10^{+212}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.11999999999999996e163 or 1.9499999999999999e37 < z < 8.49999999999999997e205 or 9.4999999999999993e212 < z Initial program 77.6%
sub-neg77.6%
+-commutative77.6%
distribute-frac-neg77.6%
distribute-rgt-neg-out77.6%
associate-/l*91.3%
fma-define91.4%
distribute-frac-neg91.4%
distribute-neg-frac291.4%
sub-neg91.4%
distribute-neg-in91.4%
remove-double-neg91.4%
+-commutative91.4%
sub-neg91.4%
Simplified91.4%
Taylor expanded in z around inf 57.1%
associate-/l*62.8%
Simplified62.8%
if -1.11999999999999996e163 < z < -1.4600000000000001e119Initial program 90.7%
Taylor expanded in t around inf 70.5%
if -1.4600000000000001e119 < z < -2.20000000000000008e-187 or -4.6e-295 < z < 1.9499999999999999e37 or 8.49999999999999997e205 < z < 9.4999999999999993e212Initial program 74.1%
Taylor expanded in a around inf 67.1%
+-commutative67.1%
Simplified67.1%
if -2.20000000000000008e-187 < z < -4.6e-295Initial program 67.7%
Taylor expanded in x around inf 93.1%
Final simplification67.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (+ x y) (/ (* y z) a))))
(if (<= t -4.6e-26)
(- x (* y (/ (- a z) t)))
(if (<= t -1.05e-208)
t_1
(if (<= t -1.8e-274)
(* y (+ (/ (- z t) (- t a)) 1.0))
(if (<= t 2.3e-73) t_1 (+ x (* y (- (/ z t) (/ a t))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - ((y * z) / a);
double tmp;
if (t <= -4.6e-26) {
tmp = x - (y * ((a - z) / t));
} else if (t <= -1.05e-208) {
tmp = t_1;
} else if (t <= -1.8e-274) {
tmp = y * (((z - t) / (t - a)) + 1.0);
} else if (t <= 2.3e-73) {
tmp = t_1;
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = (x + y) - ((y * z) / a)
if (t <= (-4.6d-26)) then
tmp = x - (y * ((a - z) / t))
else if (t <= (-1.05d-208)) then
tmp = t_1
else if (t <= (-1.8d-274)) then
tmp = y * (((z - t) / (t - a)) + 1.0d0)
else if (t <= 2.3d-73) then
tmp = t_1
else
tmp = x + (y * ((z / t) - (a / t)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - ((y * z) / a);
double tmp;
if (t <= -4.6e-26) {
tmp = x - (y * ((a - z) / t));
} else if (t <= -1.05e-208) {
tmp = t_1;
} else if (t <= -1.8e-274) {
tmp = y * (((z - t) / (t - a)) + 1.0);
} else if (t <= 2.3e-73) {
tmp = t_1;
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x + y) - ((y * z) / a) tmp = 0 if t <= -4.6e-26: tmp = x - (y * ((a - z) / t)) elif t <= -1.05e-208: tmp = t_1 elif t <= -1.8e-274: tmp = y * (((z - t) / (t - a)) + 1.0) elif t <= 2.3e-73: tmp = t_1 else: tmp = x + (y * ((z / t) - (a / t))) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(x + y) - Float64(Float64(y * z) / a)) tmp = 0.0 if (t <= -4.6e-26) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= -1.05e-208) tmp = t_1; elseif (t <= -1.8e-274) tmp = Float64(y * Float64(Float64(Float64(z - t) / Float64(t - a)) + 1.0)); elseif (t <= 2.3e-73) tmp = t_1; else tmp = Float64(x + Float64(y * Float64(Float64(z / t) - Float64(a / t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (x + y) - ((y * z) / a); tmp = 0.0; if (t <= -4.6e-26) tmp = x - (y * ((a - z) / t)); elseif (t <= -1.05e-208) tmp = t_1; elseif (t <= -1.8e-274) tmp = y * (((z - t) / (t - a)) + 1.0); elseif (t <= 2.3e-73) tmp = t_1; else tmp = x + (y * ((z / t) - (a / t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(x + y), $MachinePrecision] - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.6e-26], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -1.05e-208], t$95$1, If[LessEqual[t, -1.8e-274], N[(y * N[(N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.3e-73], t$95$1, N[(x + N[(y * N[(N[(z / t), $MachinePrecision] - N[(a / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x + y\right) - \frac{y \cdot z}{a}\\
\mathbf{if}\;t \leq -4.6 \cdot 10^{-26}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq -1.05 \cdot 10^{-208}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -1.8 \cdot 10^{-274}:\\
\;\;\;\;y \cdot \left(\frac{z - t}{t - a} + 1\right)\\
\mathbf{elif}\;t \leq 2.3 \cdot 10^{-73}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \left(\frac{z}{t} - \frac{a}{t}\right)\\
\end{array}
\end{array}
if t < -4.60000000000000018e-26Initial program 60.3%
Taylor expanded in y around 0 60.3%
associate-*l/67.8%
Simplified67.8%
*-commutative67.8%
clear-num68.1%
un-div-inv67.9%
Applied egg-rr67.9%
Taylor expanded in t around -inf 73.6%
mul-1-neg73.6%
unsub-neg73.6%
*-commutative73.6%
cancel-sign-sub-inv73.6%
neg-mul-173.6%
distribute-rgt-in73.6%
associate-/l*88.1%
neg-mul-188.1%
sub-neg88.1%
Simplified88.1%
if -4.60000000000000018e-26 < t < -1.05000000000000006e-208 or -1.79999999999999991e-274 < t < 2.29999999999999988e-73Initial program 96.1%
Taylor expanded in t around 0 83.6%
if -1.05000000000000006e-208 < t < -1.79999999999999991e-274Initial program 94.0%
Taylor expanded in x around 0 70.0%
sub-neg70.0%
*-rgt-identity70.0%
associate-*r/70.0%
distribute-rgt-neg-in70.0%
mul-1-neg70.0%
distribute-lft-in70.0%
mul-1-neg70.0%
unsub-neg70.0%
Simplified70.0%
if 2.29999999999999988e-73 < t Initial program 58.4%
Taylor expanded in y around 0 58.4%
associate-*l/67.6%
Simplified67.6%
*-commutative67.6%
clear-num67.5%
un-div-inv67.5%
Applied egg-rr67.5%
Taylor expanded in t around -inf 71.3%
mul-1-neg71.3%
unsub-neg71.3%
*-commutative71.3%
cancel-sign-sub-inv71.3%
neg-mul-171.3%
distribute-rgt-in71.4%
associate-/l*78.4%
neg-mul-178.4%
sub-neg78.4%
Simplified78.4%
div-sub78.4%
Applied egg-rr78.4%
Final simplification82.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (+ x y) (/ (* y z) a))) (t_2 (- x (* y (/ (- a z) t)))))
(if (<= t -1.7e-21)
t_2
(if (<= t -1.05e-208)
t_1
(if (<= t -6e-274)
(* y (+ (/ (- z t) (- t a)) 1.0))
(if (<= t 4.2e-76) t_1 t_2))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - ((y * z) / a);
double t_2 = x - (y * ((a - z) / t));
double tmp;
if (t <= -1.7e-21) {
tmp = t_2;
} else if (t <= -1.05e-208) {
tmp = t_1;
} else if (t <= -6e-274) {
tmp = y * (((z - t) / (t - a)) + 1.0);
} else if (t <= 4.2e-76) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (x + y) - ((y * z) / a)
t_2 = x - (y * ((a - z) / t))
if (t <= (-1.7d-21)) then
tmp = t_2
else if (t <= (-1.05d-208)) then
tmp = t_1
else if (t <= (-6d-274)) then
tmp = y * (((z - t) / (t - a)) + 1.0d0)
else if (t <= 4.2d-76) then
tmp = t_1
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - ((y * z) / a);
double t_2 = x - (y * ((a - z) / t));
double tmp;
if (t <= -1.7e-21) {
tmp = t_2;
} else if (t <= -1.05e-208) {
tmp = t_1;
} else if (t <= -6e-274) {
tmp = y * (((z - t) / (t - a)) + 1.0);
} else if (t <= 4.2e-76) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x + y) - ((y * z) / a) t_2 = x - (y * ((a - z) / t)) tmp = 0 if t <= -1.7e-21: tmp = t_2 elif t <= -1.05e-208: tmp = t_1 elif t <= -6e-274: tmp = y * (((z - t) / (t - a)) + 1.0) elif t <= 4.2e-76: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(x + y) - Float64(Float64(y * z) / a)) t_2 = Float64(x - Float64(y * Float64(Float64(a - z) / t))) tmp = 0.0 if (t <= -1.7e-21) tmp = t_2; elseif (t <= -1.05e-208) tmp = t_1; elseif (t <= -6e-274) tmp = Float64(y * Float64(Float64(Float64(z - t) / Float64(t - a)) + 1.0)); elseif (t <= 4.2e-76) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (x + y) - ((y * z) / a); t_2 = x - (y * ((a - z) / t)); tmp = 0.0; if (t <= -1.7e-21) tmp = t_2; elseif (t <= -1.05e-208) tmp = t_1; elseif (t <= -6e-274) tmp = y * (((z - t) / (t - a)) + 1.0); elseif (t <= 4.2e-76) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(x + y), $MachinePrecision] - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.7e-21], t$95$2, If[LessEqual[t, -1.05e-208], t$95$1, If[LessEqual[t, -6e-274], N[(y * N[(N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.2e-76], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x + y\right) - \frac{y \cdot z}{a}\\
t_2 := x - y \cdot \frac{a - z}{t}\\
\mathbf{if}\;t \leq -1.7 \cdot 10^{-21}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq -1.05 \cdot 10^{-208}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -6 \cdot 10^{-274}:\\
\;\;\;\;y \cdot \left(\frac{z - t}{t - a} + 1\right)\\
\mathbf{elif}\;t \leq 4.2 \cdot 10^{-76}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if t < -1.7e-21 or 4.19999999999999985e-76 < t Initial program 59.3%
Taylor expanded in y around 0 59.3%
associate-*l/67.7%
Simplified67.7%
*-commutative67.7%
clear-num67.8%
un-div-inv67.7%
Applied egg-rr67.7%
Taylor expanded in t around -inf 72.4%
mul-1-neg72.4%
unsub-neg72.4%
*-commutative72.4%
cancel-sign-sub-inv72.4%
neg-mul-172.4%
distribute-rgt-in72.4%
associate-/l*82.9%
neg-mul-182.9%
sub-neg82.9%
Simplified82.9%
if -1.7e-21 < t < -1.05000000000000006e-208 or -5.99999999999999954e-274 < t < 4.19999999999999985e-76Initial program 96.1%
Taylor expanded in t around 0 83.6%
if -1.05000000000000006e-208 < t < -5.99999999999999954e-274Initial program 94.0%
Taylor expanded in x around 0 70.0%
sub-neg70.0%
*-rgt-identity70.0%
associate-*r/70.0%
distribute-rgt-neg-in70.0%
mul-1-neg70.0%
distribute-lft-in70.0%
mul-1-neg70.0%
unsub-neg70.0%
Simplified70.0%
Final simplification82.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* y (/ (- a z) t)))))
(if (<= t -2.8e+59)
t_1
(if (<= t 4.1e-79)
(- (+ x y) (/ (* y z) a))
(if (<= t 1.75e+15) (+ x (* y (/ z t))) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y * ((a - z) / t));
double tmp;
if (t <= -2.8e+59) {
tmp = t_1;
} else if (t <= 4.1e-79) {
tmp = (x + y) - ((y * z) / a);
} else if (t <= 1.75e+15) {
tmp = x + (y * (z / t));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y * ((a - z) / t))
if (t <= (-2.8d+59)) then
tmp = t_1
else if (t <= 4.1d-79) then
tmp = (x + y) - ((y * z) / a)
else if (t <= 1.75d+15) then
tmp = x + (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 a) {
double t_1 = x - (y * ((a - z) / t));
double tmp;
if (t <= -2.8e+59) {
tmp = t_1;
} else if (t <= 4.1e-79) {
tmp = (x + y) - ((y * z) / a);
} else if (t <= 1.75e+15) {
tmp = x + (y * (z / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y * ((a - z) / t)) tmp = 0 if t <= -2.8e+59: tmp = t_1 elif t <= 4.1e-79: tmp = (x + y) - ((y * z) / a) elif t <= 1.75e+15: tmp = x + (y * (z / t)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y * Float64(Float64(a - z) / t))) tmp = 0.0 if (t <= -2.8e+59) tmp = t_1; elseif (t <= 4.1e-79) tmp = Float64(Float64(x + y) - Float64(Float64(y * z) / a)); elseif (t <= 1.75e+15) tmp = Float64(x + Float64(y * Float64(z / t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y * ((a - z) / t)); tmp = 0.0; if (t <= -2.8e+59) tmp = t_1; elseif (t <= 4.1e-79) tmp = (x + y) - ((y * z) / a); elseif (t <= 1.75e+15) tmp = x + (y * (z / t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.8e+59], t$95$1, If[LessEqual[t, 4.1e-79], N[(N[(x + y), $MachinePrecision] - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.75e+15], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - y \cdot \frac{a - z}{t}\\
\mathbf{if}\;t \leq -2.8 \cdot 10^{+59}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 4.1 \cdot 10^{-79}:\\
\;\;\;\;\left(x + y\right) - \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 1.75 \cdot 10^{+15}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -2.7999999999999998e59 or 1.75e15 < t Initial program 53.0%
Taylor expanded in y around 0 53.0%
associate-*l/62.5%
Simplified62.5%
*-commutative62.5%
clear-num62.7%
un-div-inv62.6%
Applied egg-rr62.6%
Taylor expanded in t around -inf 74.2%
mul-1-neg74.2%
unsub-neg74.2%
*-commutative74.2%
cancel-sign-sub-inv74.2%
neg-mul-174.2%
distribute-rgt-in74.4%
associate-/l*88.1%
neg-mul-188.1%
sub-neg88.1%
Simplified88.1%
if -2.7999999999999998e59 < t < 4.09999999999999994e-79Initial program 94.8%
Taylor expanded in t around 0 77.2%
if 4.09999999999999994e-79 < t < 1.75e15Initial program 73.2%
Taylor expanded in y around 0 73.2%
associate-*l/81.9%
Simplified81.9%
*-commutative81.9%
clear-num81.9%
un-div-inv82.0%
Applied egg-rr82.0%
Taylor expanded in t around -inf 60.4%
mul-1-neg60.4%
unsub-neg60.4%
*-commutative60.4%
cancel-sign-sub-inv60.4%
neg-mul-160.4%
distribute-rgt-in60.4%
associate-/l*60.3%
neg-mul-160.3%
sub-neg60.3%
Simplified60.3%
Taylor expanded in a around 0 75.7%
neg-mul-175.7%
distribute-neg-frac275.7%
Simplified75.7%
sub-neg75.7%
+-commutative75.7%
distribute-lft-neg-in75.7%
add-sqr-sqrt54.2%
sqrt-unprod60.0%
sqr-neg60.0%
sqrt-unprod16.5%
add-sqr-sqrt44.5%
add-sqr-sqrt0.0%
sqrt-unprod75.7%
sqr-neg75.7%
sqrt-unprod75.5%
add-sqr-sqrt75.7%
Applied egg-rr75.7%
Final simplification81.7%
(FPCore (x y z t a)
:precision binary64
(if (<= t -1.7e+101)
(- x (* y (/ (- a z) t)))
(if (<= t 2.4e+47)
(- (+ x y) (* (- z t) (/ y (- a t))))
(+ (- x (* a (/ y t))) (* y (/ z t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.7e+101) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 2.4e+47) {
tmp = (x + y) - ((z - t) * (y / (a - t)));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-1.7d+101)) then
tmp = x - (y * ((a - z) / t))
else if (t <= 2.4d+47) then
tmp = (x + y) - ((z - t) * (y / (a - t)))
else
tmp = (x - (a * (y / t))) + (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.7e+101) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 2.4e+47) {
tmp = (x + y) - ((z - t) * (y / (a - t)));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.7e+101: tmp = x - (y * ((a - z) / t)) elif t <= 2.4e+47: tmp = (x + y) - ((z - t) * (y / (a - t))) else: tmp = (x - (a * (y / t))) + (y * (z / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.7e+101) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= 2.4e+47) tmp = Float64(Float64(x + y) - Float64(Float64(z - t) * Float64(y / Float64(a - t)))); else tmp = Float64(Float64(x - Float64(a * Float64(y / t))) + Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.7e+101) tmp = x - (y * ((a - z) / t)); elseif (t <= 2.4e+47) tmp = (x + y) - ((z - t) * (y / (a - t))); else tmp = (x - (a * (y / t))) + (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.7e+101], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.4e+47], N[(N[(x + y), $MachinePrecision] - N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x - N[(a * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.7 \cdot 10^{+101}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{+47}:\\
\;\;\;\;\left(x + y\right) - \left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{else}:\\
\;\;\;\;\left(x - a \cdot \frac{y}{t}\right) + y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if t < -1.70000000000000009e101Initial program 44.6%
Taylor expanded in y around 0 44.6%
associate-*l/56.4%
Simplified56.4%
*-commutative56.4%
clear-num56.8%
un-div-inv56.4%
Applied egg-rr56.4%
Taylor expanded in t around -inf 67.9%
mul-1-neg67.9%
unsub-neg67.9%
*-commutative67.9%
cancel-sign-sub-inv67.9%
neg-mul-167.9%
distribute-rgt-in68.0%
associate-/l*90.7%
neg-mul-190.7%
sub-neg90.7%
Simplified90.7%
if -1.70000000000000009e101 < t < 2.40000000000000019e47Initial program 90.9%
Taylor expanded in y around 0 90.9%
associate-*l/94.1%
Simplified94.1%
if 2.40000000000000019e47 < t Initial program 53.6%
Taylor expanded in y around 0 53.6%
associate-*l/62.0%
Simplified62.0%
Taylor expanded in t around inf 76.7%
sub-neg76.7%
mul-1-neg76.7%
unsub-neg76.7%
associate-/l*81.7%
mul-1-neg81.7%
remove-double-neg81.7%
associate-/l*86.8%
Simplified86.8%
Final simplification92.0%
(FPCore (x y z t a)
:precision binary64
(if (<= t -2.2e+102)
(- x (* y (/ (- a z) t)))
(if (<= t 3.9e-33)
(+ (+ x y) (* z (/ y (- t a))))
(+ (- x (* a (/ y t))) (* y (/ z t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.2e+102) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (z * (y / (t - a)));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-2.2d+102)) then
tmp = x - (y * ((a - z) / t))
else if (t <= 3.9d-33) then
tmp = (x + y) + (z * (y / (t - a)))
else
tmp = (x - (a * (y / t))) + (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.2e+102) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (z * (y / (t - a)));
} else {
tmp = (x - (a * (y / t))) + (y * (z / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.2e+102: tmp = x - (y * ((a - z) / t)) elif t <= 3.9e-33: tmp = (x + y) + (z * (y / (t - a))) else: tmp = (x - (a * (y / t))) + (y * (z / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.2e+102) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= 3.9e-33) tmp = Float64(Float64(x + y) + Float64(z * Float64(y / Float64(t - a)))); else tmp = Float64(Float64(x - Float64(a * Float64(y / t))) + Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -2.2e+102) tmp = x - (y * ((a - z) / t)); elseif (t <= 3.9e-33) tmp = (x + y) + (z * (y / (t - a))); else tmp = (x - (a * (y / t))) + (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.2e+102], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.9e-33], N[(N[(x + y), $MachinePrecision] + N[(z * N[(y / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x - N[(a * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.2 \cdot 10^{+102}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq 3.9 \cdot 10^{-33}:\\
\;\;\;\;\left(x + y\right) + z \cdot \frac{y}{t - a}\\
\mathbf{else}:\\
\;\;\;\;\left(x - a \cdot \frac{y}{t}\right) + y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if t < -2.20000000000000007e102Initial program 44.6%
Taylor expanded in y around 0 44.6%
associate-*l/56.4%
Simplified56.4%
*-commutative56.4%
clear-num56.8%
un-div-inv56.4%
Applied egg-rr56.4%
Taylor expanded in t around -inf 67.9%
mul-1-neg67.9%
unsub-neg67.9%
*-commutative67.9%
cancel-sign-sub-inv67.9%
neg-mul-167.9%
distribute-rgt-in68.0%
associate-/l*90.7%
neg-mul-190.7%
sub-neg90.7%
Simplified90.7%
if -2.20000000000000007e102 < t < 3.89999999999999974e-33Initial program 93.0%
Taylor expanded in y around 0 93.0%
associate-*l/95.9%
Simplified95.9%
Taylor expanded in z around inf 92.6%
associate-*l/95.1%
*-commutative95.1%
Simplified95.1%
if 3.89999999999999974e-33 < t Initial program 56.1%
Taylor expanded in y around 0 56.1%
associate-*l/64.2%
Simplified64.2%
Taylor expanded in t around inf 73.9%
sub-neg73.9%
mul-1-neg73.9%
unsub-neg73.9%
associate-/l*78.0%
mul-1-neg78.0%
remove-double-neg78.0%
associate-/l*82.1%
Simplified82.1%
Final simplification91.0%
(FPCore (x y z t a)
:precision binary64
(if (<= t -6.8e+107)
(- x (* y (/ (- a z) t)))
(if (<= t 3.9e-33)
(+ (+ x y) (* z (/ y (- t a))))
(+ x (* y (- (/ z t) (/ a t)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -6.8e+107) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (z * (y / (t - a)));
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-6.8d+107)) then
tmp = x - (y * ((a - z) / t))
else if (t <= 3.9d-33) then
tmp = (x + y) + (z * (y / (t - a)))
else
tmp = x + (y * ((z / t) - (a / t)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -6.8e+107) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (z * (y / (t - a)));
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -6.8e+107: tmp = x - (y * ((a - z) / t)) elif t <= 3.9e-33: tmp = (x + y) + (z * (y / (t - a))) else: tmp = x + (y * ((z / t) - (a / t))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -6.8e+107) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= 3.9e-33) tmp = Float64(Float64(x + y) + Float64(z * Float64(y / Float64(t - a)))); else tmp = Float64(x + Float64(y * Float64(Float64(z / t) - Float64(a / t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -6.8e+107) tmp = x - (y * ((a - z) / t)); elseif (t <= 3.9e-33) tmp = (x + y) + (z * (y / (t - a))); else tmp = x + (y * ((z / t) - (a / t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -6.8e+107], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.9e-33], N[(N[(x + y), $MachinePrecision] + N[(z * N[(y / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(z / t), $MachinePrecision] - N[(a / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -6.8 \cdot 10^{+107}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq 3.9 \cdot 10^{-33}:\\
\;\;\;\;\left(x + y\right) + z \cdot \frac{y}{t - a}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \left(\frac{z}{t} - \frac{a}{t}\right)\\
\end{array}
\end{array}
if t < -6.7999999999999994e107Initial program 44.6%
Taylor expanded in y around 0 44.6%
associate-*l/56.4%
Simplified56.4%
*-commutative56.4%
clear-num56.8%
un-div-inv56.4%
Applied egg-rr56.4%
Taylor expanded in t around -inf 67.9%
mul-1-neg67.9%
unsub-neg67.9%
*-commutative67.9%
cancel-sign-sub-inv67.9%
neg-mul-167.9%
distribute-rgt-in68.0%
associate-/l*90.7%
neg-mul-190.7%
sub-neg90.7%
Simplified90.7%
if -6.7999999999999994e107 < t < 3.89999999999999974e-33Initial program 93.0%
Taylor expanded in y around 0 93.0%
associate-*l/95.9%
Simplified95.9%
Taylor expanded in z around inf 92.6%
associate-*l/95.1%
*-commutative95.1%
Simplified95.1%
if 3.89999999999999974e-33 < t Initial program 56.1%
Taylor expanded in y around 0 56.1%
associate-*l/64.2%
Simplified64.2%
*-commutative64.2%
clear-num64.1%
un-div-inv64.1%
Applied egg-rr64.1%
Taylor expanded in t around -inf 73.9%
mul-1-neg73.9%
unsub-neg73.9%
*-commutative73.9%
cancel-sign-sub-inv73.9%
neg-mul-173.9%
distribute-rgt-in74.0%
associate-/l*82.1%
neg-mul-182.1%
sub-neg82.1%
Simplified82.1%
div-sub82.1%
Applied egg-rr82.1%
Final simplification91.0%
(FPCore (x y z t a)
:precision binary64
(if (<= t -7.2e+100)
(- x (* y (/ (- a z) t)))
(if (<= t 3.9e-33)
(+ (+ x y) (* y (/ z (- t a))))
(+ x (* y (- (/ z t) (/ a t)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -7.2e+100) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (y * (z / (t - a)));
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-7.2d+100)) then
tmp = x - (y * ((a - z) / t))
else if (t <= 3.9d-33) then
tmp = (x + y) + (y * (z / (t - a)))
else
tmp = x + (y * ((z / t) - (a / t)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -7.2e+100) {
tmp = x - (y * ((a - z) / t));
} else if (t <= 3.9e-33) {
tmp = (x + y) + (y * (z / (t - a)));
} else {
tmp = x + (y * ((z / t) - (a / t)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -7.2e+100: tmp = x - (y * ((a - z) / t)) elif t <= 3.9e-33: tmp = (x + y) + (y * (z / (t - a))) else: tmp = x + (y * ((z / t) - (a / t))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -7.2e+100) tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); elseif (t <= 3.9e-33) tmp = Float64(Float64(x + y) + Float64(y * Float64(z / Float64(t - a)))); else tmp = Float64(x + Float64(y * Float64(Float64(z / t) - Float64(a / t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -7.2e+100) tmp = x - (y * ((a - z) / t)); elseif (t <= 3.9e-33) tmp = (x + y) + (y * (z / (t - a))); else tmp = x + (y * ((z / t) - (a / t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -7.2e+100], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.9e-33], N[(N[(x + y), $MachinePrecision] + N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(z / t), $MachinePrecision] - N[(a / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -7.2 \cdot 10^{+100}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\mathbf{elif}\;t \leq 3.9 \cdot 10^{-33}:\\
\;\;\;\;\left(x + y\right) + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \left(\frac{z}{t} - \frac{a}{t}\right)\\
\end{array}
\end{array}
if t < -7.2e100Initial program 44.6%
Taylor expanded in y around 0 44.6%
associate-*l/56.4%
Simplified56.4%
*-commutative56.4%
clear-num56.8%
un-div-inv56.4%
Applied egg-rr56.4%
Taylor expanded in t around -inf 67.9%
mul-1-neg67.9%
unsub-neg67.9%
*-commutative67.9%
cancel-sign-sub-inv67.9%
neg-mul-167.9%
distribute-rgt-in68.0%
associate-/l*90.7%
neg-mul-190.7%
sub-neg90.7%
Simplified90.7%
if -7.2e100 < t < 3.89999999999999974e-33Initial program 93.0%
Taylor expanded in z around inf 92.6%
associate-/l*91.3%
Simplified91.3%
if 3.89999999999999974e-33 < t Initial program 56.1%
Taylor expanded in y around 0 56.1%
associate-*l/64.2%
Simplified64.2%
*-commutative64.2%
clear-num64.1%
un-div-inv64.1%
Applied egg-rr64.1%
Taylor expanded in t around -inf 73.9%
mul-1-neg73.9%
unsub-neg73.9%
*-commutative73.9%
cancel-sign-sub-inv73.9%
neg-mul-173.9%
distribute-rgt-in74.0%
associate-/l*82.1%
neg-mul-182.1%
sub-neg82.1%
Simplified82.1%
div-sub82.1%
Applied egg-rr82.1%
Final simplification88.8%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -5.6e+51) (not (<= a 1.65e+21))) (+ x y) (- x (* y (/ (- a z) t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -5.6e+51) || !(a <= 1.65e+21)) {
tmp = x + y;
} else {
tmp = x - (y * ((a - z) / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a <= (-5.6d+51)) .or. (.not. (a <= 1.65d+21))) then
tmp = x + y
else
tmp = x - (y * ((a - z) / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -5.6e+51) || !(a <= 1.65e+21)) {
tmp = x + y;
} else {
tmp = x - (y * ((a - z) / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -5.6e+51) or not (a <= 1.65e+21): tmp = x + y else: tmp = x - (y * ((a - z) / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -5.6e+51) || !(a <= 1.65e+21)) tmp = Float64(x + y); else tmp = Float64(x - Float64(y * Float64(Float64(a - z) / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -5.6e+51) || ~((a <= 1.65e+21))) tmp = x + y; else tmp = x - (y * ((a - z) / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -5.6e+51], N[Not[LessEqual[a, 1.65e+21]], $MachinePrecision]], N[(x + y), $MachinePrecision], N[(x - N[(y * N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.6 \cdot 10^{+51} \lor \neg \left(a \leq 1.65 \cdot 10^{+21}\right):\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{a - z}{t}\\
\end{array}
\end{array}
if a < -5.60000000000000009e51 or 1.65e21 < a Initial program 75.5%
Taylor expanded in a around inf 75.0%
+-commutative75.0%
Simplified75.0%
if -5.60000000000000009e51 < a < 1.65e21Initial program 75.3%
Taylor expanded in y around 0 75.3%
associate-*l/76.8%
Simplified76.8%
*-commutative76.8%
clear-num76.7%
un-div-inv77.4%
Applied egg-rr77.4%
Taylor expanded in t around -inf 75.8%
mul-1-neg75.8%
unsub-neg75.8%
*-commutative75.8%
cancel-sign-sub-inv75.8%
neg-mul-175.8%
distribute-rgt-in75.8%
associate-/l*77.0%
neg-mul-177.0%
sub-neg77.0%
Simplified77.0%
Final simplification76.2%
(FPCore (x y z t a) :precision binary64 (if (<= z -2.4e+196) (* z (/ y t)) (if (or (<= z -5.6e-189) (not (<= z -1.02e-290))) (+ x y) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -2.4e+196) {
tmp = z * (y / t);
} else if ((z <= -5.6e-189) || !(z <= -1.02e-290)) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-2.4d+196)) then
tmp = z * (y / t)
else if ((z <= (-5.6d-189)) .or. (.not. (z <= (-1.02d-290)))) then
tmp = x + y
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -2.4e+196) {
tmp = z * (y / t);
} else if ((z <= -5.6e-189) || !(z <= -1.02e-290)) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -2.4e+196: tmp = z * (y / t) elif (z <= -5.6e-189) or not (z <= -1.02e-290): tmp = x + y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -2.4e+196) tmp = Float64(z * Float64(y / t)); elseif ((z <= -5.6e-189) || !(z <= -1.02e-290)) tmp = Float64(x + y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -2.4e+196) tmp = z * (y / t); elseif ((z <= -5.6e-189) || ~((z <= -1.02e-290))) tmp = x + y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -2.4e+196], N[(z * N[(y / t), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -5.6e-189], N[Not[LessEqual[z, -1.02e-290]], $MachinePrecision]], N[(x + y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.4 \cdot 10^{+196}:\\
\;\;\;\;z \cdot \frac{y}{t}\\
\mathbf{elif}\;z \leq -5.6 \cdot 10^{-189} \lor \neg \left(z \leq -1.02 \cdot 10^{-290}\right):\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -2.4e196Initial program 77.4%
sub-neg77.4%
+-commutative77.4%
distribute-frac-neg77.4%
distribute-rgt-neg-out77.4%
associate-/l*96.3%
fma-define96.5%
distribute-frac-neg96.5%
distribute-neg-frac296.5%
sub-neg96.5%
distribute-neg-in96.5%
remove-double-neg96.5%
+-commutative96.5%
sub-neg96.5%
Simplified96.5%
Taylor expanded in z around inf 60.7%
associate-/l*76.6%
Simplified76.6%
Taylor expanded in t around inf 44.2%
*-commutative44.2%
associate-/l*53.7%
un-div-inv53.7%
*-commutative53.7%
un-div-inv53.7%
Applied egg-rr53.7%
if -2.4e196 < z < -5.5999999999999999e-189 or -1.02e-290 < z Initial program 75.9%
Taylor expanded in a around inf 58.5%
+-commutative58.5%
Simplified58.5%
if -5.5999999999999999e-189 < z < -1.02e-290Initial program 67.7%
Taylor expanded in x around inf 93.1%
Final simplification60.7%
(FPCore (x y z t a) :precision binary64 (if (<= z -6.6e+188) (* y (/ z t)) (if (or (<= z -3.2e-189) (not (<= z -2.9e-292))) (+ x y) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -6.6e+188) {
tmp = y * (z / t);
} else if ((z <= -3.2e-189) || !(z <= -2.9e-292)) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-6.6d+188)) then
tmp = y * (z / t)
else if ((z <= (-3.2d-189)) .or. (.not. (z <= (-2.9d-292)))) then
tmp = x + y
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -6.6e+188) {
tmp = y * (z / t);
} else if ((z <= -3.2e-189) || !(z <= -2.9e-292)) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -6.6e+188: tmp = y * (z / t) elif (z <= -3.2e-189) or not (z <= -2.9e-292): tmp = x + y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -6.6e+188) tmp = Float64(y * Float64(z / t)); elseif ((z <= -3.2e-189) || !(z <= -2.9e-292)) tmp = Float64(x + y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -6.6e+188) tmp = y * (z / t); elseif ((z <= -3.2e-189) || ~((z <= -2.9e-292))) tmp = x + y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -6.6e+188], N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -3.2e-189], N[Not[LessEqual[z, -2.9e-292]], $MachinePrecision]], N[(x + y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6.6 \cdot 10^{+188}:\\
\;\;\;\;y \cdot \frac{z}{t}\\
\mathbf{elif}\;z \leq -3.2 \cdot 10^{-189} \lor \neg \left(z \leq -2.9 \cdot 10^{-292}\right):\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -6.59999999999999966e188Initial program 78.2%
sub-neg78.2%
+-commutative78.2%
distribute-frac-neg78.2%
distribute-rgt-neg-out78.2%
associate-/l*96.4%
fma-define96.6%
distribute-frac-neg96.6%
distribute-neg-frac296.6%
sub-neg96.6%
distribute-neg-in96.6%
remove-double-neg96.6%
+-commutative96.6%
sub-neg96.6%
Simplified96.6%
Taylor expanded in z around inf 62.0%
associate-/l*77.4%
Simplified77.4%
Taylor expanded in t around inf 46.0%
associate-/l*52.1%
Simplified52.1%
if -6.59999999999999966e188 < z < -3.2000000000000001e-189 or -2.89999999999999993e-292 < z Initial program 75.8%
Taylor expanded in a around inf 58.8%
+-commutative58.8%
Simplified58.8%
if -3.2000000000000001e-189 < z < -2.89999999999999993e-292Initial program 67.7%
Taylor expanded in x around inf 93.1%
Final simplification60.7%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -3.2e+59) (not (<= a 2e-6))) (+ x y) (+ x (* y (/ z t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.2e+59) || !(a <= 2e-6)) {
tmp = x + y;
} else {
tmp = x + (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a <= (-3.2d+59)) .or. (.not. (a <= 2d-6))) then
tmp = x + y
else
tmp = x + (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.2e+59) || !(a <= 2e-6)) {
tmp = x + y;
} else {
tmp = x + (y * (z / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -3.2e+59) or not (a <= 2e-6): tmp = x + y else: tmp = x + (y * (z / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -3.2e+59) || !(a <= 2e-6)) tmp = Float64(x + y); else tmp = Float64(x + Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -3.2e+59) || ~((a <= 2e-6))) tmp = x + y; else tmp = x + (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -3.2e+59], N[Not[LessEqual[a, 2e-6]], $MachinePrecision]], N[(x + y), $MachinePrecision], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.2 \cdot 10^{+59} \lor \neg \left(a \leq 2 \cdot 10^{-6}\right):\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if a < -3.19999999999999982e59 or 1.99999999999999991e-6 < a Initial program 74.8%
Taylor expanded in a around inf 74.3%
+-commutative74.3%
Simplified74.3%
if -3.19999999999999982e59 < a < 1.99999999999999991e-6Initial program 75.8%
Taylor expanded in y around 0 75.8%
associate-*l/77.3%
Simplified77.3%
*-commutative77.3%
clear-num77.2%
un-div-inv77.9%
Applied egg-rr77.9%
Taylor expanded in t around -inf 74.3%
mul-1-neg74.3%
unsub-neg74.3%
*-commutative74.3%
cancel-sign-sub-inv74.3%
neg-mul-174.3%
distribute-rgt-in74.3%
associate-/l*75.5%
neg-mul-175.5%
sub-neg75.5%
Simplified75.5%
Taylor expanded in a around 0 74.1%
neg-mul-174.1%
distribute-neg-frac274.1%
Simplified74.1%
sub-neg74.1%
+-commutative74.1%
distribute-lft-neg-in74.1%
add-sqr-sqrt38.4%
sqrt-unprod53.0%
sqr-neg53.0%
sqrt-unprod24.2%
add-sqr-sqrt50.3%
add-sqr-sqrt26.0%
sqrt-unprod54.6%
sqr-neg54.6%
sqrt-unprod35.7%
add-sqr-sqrt74.1%
Applied egg-rr74.1%
Final simplification74.2%
(FPCore (x y z t a) :precision binary64 (if (<= t -1.05e+225) x (if (<= t 1e+93) (+ x y) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.05e+225) {
tmp = x;
} else if (t <= 1e+93) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-1.05d+225)) then
tmp = x
else if (t <= 1d+93) then
tmp = x + y
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.05e+225) {
tmp = x;
} else if (t <= 1e+93) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.05e+225: tmp = x elif t <= 1e+93: tmp = x + y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.05e+225) tmp = x; elseif (t <= 1e+93) tmp = Float64(x + y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.05e+225) tmp = x; elseif (t <= 1e+93) tmp = x + y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.05e+225], x, If[LessEqual[t, 1e+93], N[(x + y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.05 \cdot 10^{+225}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 10^{+93}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if t < -1.05e225 or 1.00000000000000004e93 < t Initial program 45.4%
Taylor expanded in x around inf 61.6%
if -1.05e225 < t < 1.00000000000000004e93Initial program 84.4%
Taylor expanded in a around inf 58.2%
+-commutative58.2%
Simplified58.2%
Final simplification59.0%
(FPCore (x y z t a) :precision binary64 (if (<= a -3.2e+299) y x))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -3.2e+299) {
tmp = y;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-3.2d+299)) then
tmp = y
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -3.2e+299) {
tmp = y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -3.2e+299: tmp = y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -3.2e+299) tmp = y; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -3.2e+299) tmp = y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -3.2e+299], y, x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.2 \cdot 10^{+299}:\\
\;\;\;\;y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -3.1999999999999999e299Initial program 82.4%
Taylor expanded in t around 0 82.4%
Taylor expanded in x around 0 82.4%
Taylor expanded in z around 0 100.0%
if -3.1999999999999999e299 < a Initial program 75.3%
Taylor expanded in x around inf 49.5%
(FPCore (x y z t a) :precision binary64 x)
double code(double x, double y, double z, double t, double a) {
return x;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x
end function
public static double code(double x, double y, double z, double t, double a) {
return x;
}
def code(x, y, z, t, a): return x
function code(x, y, z, t, a) return x end
function tmp = code(x, y, z, t, a) tmp = x; end
code[x_, y_, z_, t_, a_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 75.4%
Taylor expanded in x around inf 48.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y)))
(t_2 (- (+ x y) (/ (* (- z t) y) (- a t)))))
(if (< t_2 -1.3664970889390727e-7)
t_1
(if (< t_2 1.4754293444577233e-239)
(/ (- (* y (- a z)) (* x t)) (- a t))
t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y);
double t_2 = (x + y) - (((z - t) * y) / (a - t));
double tmp;
if (t_2 < -1.3664970889390727e-7) {
tmp = t_1;
} else if (t_2 < 1.4754293444577233e-239) {
tmp = ((y * (a - z)) - (x * t)) / (a - t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (y + x) - (((z - t) * (1.0d0 / (a - t))) * y)
t_2 = (x + y) - (((z - t) * y) / (a - t))
if (t_2 < (-1.3664970889390727d-7)) then
tmp = t_1
else if (t_2 < 1.4754293444577233d-239) then
tmp = ((y * (a - z)) - (x * t)) / (a - t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y);
double t_2 = (x + y) - (((z - t) * y) / (a - t));
double tmp;
if (t_2 < -1.3664970889390727e-7) {
tmp = t_1;
} else if (t_2 < 1.4754293444577233e-239) {
tmp = ((y * (a - z)) - (x * t)) / (a - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y) t_2 = (x + y) - (((z - t) * y) / (a - t)) tmp = 0 if t_2 < -1.3664970889390727e-7: tmp = t_1 elif t_2 < 1.4754293444577233e-239: tmp = ((y * (a - z)) - (x * t)) / (a - t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y + x) - Float64(Float64(Float64(z - t) * Float64(1.0 / Float64(a - t))) * y)) t_2 = Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) tmp = 0.0 if (t_2 < -1.3664970889390727e-7) tmp = t_1; elseif (t_2 < 1.4754293444577233e-239) tmp = Float64(Float64(Float64(y * Float64(a - z)) - Float64(x * t)) / Float64(a - t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y); t_2 = (x + y) - (((z - t) * y) / (a - t)); tmp = 0.0; if (t_2 < -1.3664970889390727e-7) tmp = t_1; elseif (t_2 < 1.4754293444577233e-239) tmp = ((y * (a - z)) - (x * t)) / (a - t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y + x), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * N[(1.0 / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t$95$2, -1.3664970889390727e-7], t$95$1, If[Less[t$95$2, 1.4754293444577233e-239], N[(N[(N[(y * N[(a - z), $MachinePrecision]), $MachinePrecision] - N[(x * t), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y + x\right) - \left(\left(z - t\right) \cdot \frac{1}{a - t}\right) \cdot y\\
t_2 := \left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\\
\mathbf{if}\;t\_2 < -1.3664970889390727 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 < 1.4754293444577233 \cdot 10^{-239}:\\
\;\;\;\;\frac{y \cdot \left(a - z\right) - x \cdot t}{a - t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
herbie shell --seed 2024107
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
:precision binary64
:alt
(if (< (- (+ x y) (/ (* (- z t) y) (- a t))) -1.3664970889390727e-7) (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y)) (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) 1.4754293444577233e-239) (/ (- (* y (- a z)) (* x t)) (- a t)) (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y))))
(- (+ x y) (/ (* (- z t) y) (- a t))))