
(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 15 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 (+ (* y (- (+ (/ t (- a t)) 1.0) (/ z (- a t)))) x))
double code(double x, double y, double z, double t, double a) {
return (y * (((t / (a - t)) + 1.0) - (z / (a - t)))) + 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 = (y * (((t / (a - t)) + 1.0d0) - (z / (a - t)))) + x
end function
public static double code(double x, double y, double z, double t, double a) {
return (y * (((t / (a - t)) + 1.0) - (z / (a - t)))) + x;
}
def code(x, y, z, t, a): return (y * (((t / (a - t)) + 1.0) - (z / (a - t)))) + x
function code(x, y, z, t, a) return Float64(Float64(y * Float64(Float64(Float64(t / Float64(a - t)) + 1.0) - Float64(z / Float64(a - t)))) + x) end
function tmp = code(x, y, z, t, a) tmp = (y * (((t / (a - t)) + 1.0) - (z / (a - t)))) + x; end
code[x_, y_, z_, t_, a_] := N[(N[(y * N[(N[(N[(t / N[(a - t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[(z / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
y \cdot \left(\left(\frac{t}{a - t} + 1\right) - \frac{z}{a - t}\right) + x
\end{array}
Initial program 75.1%
associate--l+79.2%
sub-neg79.2%
+-commutative79.2%
associate-/l*86.6%
distribute-neg-frac86.6%
associate-/r/88.1%
fma-def88.1%
sub-neg88.1%
+-commutative88.1%
distribute-neg-in88.1%
unsub-neg88.1%
remove-double-neg88.1%
Simplified88.1%
Taylor expanded in y around 0 94.7%
Final simplification94.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (- x (/ y (/ a z))))))
(if (<= a -2.05e+49)
t_1
(if (<= a -1.15e-93)
(+ x (* y (/ z t)))
(if (<= a -2.6e-99)
(- y (/ (* y z) a))
(if (<= a 32.0) (+ x (* y (/ (- z a) t))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + (x - (y / (a / z)));
double tmp;
if (a <= -2.05e+49) {
tmp = t_1;
} else if (a <= -1.15e-93) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = y - ((y * z) / a);
} else if (a <= 32.0) {
tmp = x + (y * ((z - 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) :: tmp
t_1 = y + (x - (y / (a / z)))
if (a <= (-2.05d+49)) then
tmp = t_1
else if (a <= (-1.15d-93)) then
tmp = x + (y * (z / t))
else if (a <= (-2.6d-99)) then
tmp = y - ((y * z) / a)
else if (a <= 32.0d0) then
tmp = x + (y * ((z - 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 - (y / (a / z)));
double tmp;
if (a <= -2.05e+49) {
tmp = t_1;
} else if (a <= -1.15e-93) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = y - ((y * z) / a);
} else if (a <= 32.0) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + (x - (y / (a / z))) tmp = 0 if a <= -2.05e+49: tmp = t_1 elif a <= -1.15e-93: tmp = x + (y * (z / t)) elif a <= -2.6e-99: tmp = y - ((y * z) / a) elif a <= 32.0: tmp = x + (y * ((z - a) / t)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(x - Float64(y / Float64(a / z)))) tmp = 0.0 if (a <= -2.05e+49) tmp = t_1; elseif (a <= -1.15e-93) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= -2.6e-99) tmp = Float64(y - Float64(Float64(y * z) / a)); elseif (a <= 32.0) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + (x - (y / (a / z))); tmp = 0.0; if (a <= -2.05e+49) tmp = t_1; elseif (a <= -1.15e-93) tmp = x + (y * (z / t)); elseif (a <= -2.6e-99) tmp = y - ((y * z) / a); elseif (a <= 32.0) tmp = x + (y * ((z - a) / t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.05e+49], t$95$1, If[LessEqual[a, -1.15e-93], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -2.6e-99], N[(y - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 32.0], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{if}\;a \leq -2.05 \cdot 10^{+49}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -1.15 \cdot 10^{-93}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq -2.6 \cdot 10^{-99}:\\
\;\;\;\;y - \frac{y \cdot z}{a}\\
\mathbf{elif}\;a \leq 32:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -2.05e49 or 32 < a Initial program 79.8%
associate--l+80.7%
sub-neg80.7%
+-commutative80.7%
associate-/l*94.7%
distribute-neg-frac94.7%
associate-/r/95.6%
fma-def95.6%
sub-neg95.6%
+-commutative95.6%
distribute-neg-in95.6%
unsub-neg95.6%
remove-double-neg95.6%
Simplified95.6%
Taylor expanded in t around 0 79.2%
mul-1-neg79.2%
sub-neg79.2%
associate-/l*89.6%
Simplified89.6%
if -2.05e49 < a < -1.1499999999999999e-93Initial program 68.5%
associate--l+78.0%
sub-neg78.0%
+-commutative78.0%
associate-/l*77.2%
distribute-neg-frac77.2%
associate-/r/85.0%
fma-def85.0%
sub-neg85.0%
+-commutative85.0%
distribute-neg-in85.0%
unsub-neg85.0%
remove-double-neg85.0%
Simplified85.0%
Taylor expanded in y around 0 95.7%
Taylor expanded in a around 0 77.1%
if -1.1499999999999999e-93 < a < -2.60000000000000005e-99Initial program 99.6%
associate-*l/76.9%
Simplified76.9%
Taylor expanded in t around 0 76.9%
Taylor expanded in x around 0 97.0%
if -2.60000000000000005e-99 < a < 32Initial program 71.5%
associate--l+77.3%
sub-neg77.3%
+-commutative77.3%
associate-/l*80.8%
distribute-neg-frac80.8%
associate-/r/82.0%
fma-def82.0%
sub-neg82.0%
+-commutative82.0%
distribute-neg-in82.0%
unsub-neg82.0%
remove-double-neg82.0%
Simplified82.0%
Taylor expanded in y around 0 94.3%
Taylor expanded in t around inf 92.4%
associate-*r/92.4%
distribute-lft-in92.4%
neg-mul-192.4%
mul-1-neg92.4%
remove-double-neg92.4%
neg-mul-192.4%
sub-neg92.4%
Simplified92.4%
Final simplification89.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (- x (/ y (/ a z))))))
(if (<= a -2.05e+49)
t_1
(if (<= a -1.15e-93)
(+ x (* y (/ z t)))
(if (<= a -2.6e-99)
(- y (/ (* y z) a))
(if (<= a 4.8e+68) (+ x (/ (* y (- z a)) t)) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + (x - (y / (a / z)));
double tmp;
if (a <= -2.05e+49) {
tmp = t_1;
} else if (a <= -1.15e-93) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = y - ((y * z) / a);
} else if (a <= 4.8e+68) {
tmp = x + ((y * (z - 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) :: tmp
t_1 = y + (x - (y / (a / z)))
if (a <= (-2.05d+49)) then
tmp = t_1
else if (a <= (-1.15d-93)) then
tmp = x + (y * (z / t))
else if (a <= (-2.6d-99)) then
tmp = y - ((y * z) / a)
else if (a <= 4.8d+68) then
tmp = x + ((y * (z - 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 - (y / (a / z)));
double tmp;
if (a <= -2.05e+49) {
tmp = t_1;
} else if (a <= -1.15e-93) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = y - ((y * z) / a);
} else if (a <= 4.8e+68) {
tmp = x + ((y * (z - a)) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + (x - (y / (a / z))) tmp = 0 if a <= -2.05e+49: tmp = t_1 elif a <= -1.15e-93: tmp = x + (y * (z / t)) elif a <= -2.6e-99: tmp = y - ((y * z) / a) elif a <= 4.8e+68: tmp = x + ((y * (z - a)) / t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(x - Float64(y / Float64(a / z)))) tmp = 0.0 if (a <= -2.05e+49) tmp = t_1; elseif (a <= -1.15e-93) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= -2.6e-99) tmp = Float64(y - Float64(Float64(y * z) / a)); elseif (a <= 4.8e+68) tmp = Float64(x + Float64(Float64(y * Float64(z - a)) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + (x - (y / (a / z))); tmp = 0.0; if (a <= -2.05e+49) tmp = t_1; elseif (a <= -1.15e-93) tmp = x + (y * (z / t)); elseif (a <= -2.6e-99) tmp = y - ((y * z) / a); elseif (a <= 4.8e+68) tmp = x + ((y * (z - a)) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.05e+49], t$95$1, If[LessEqual[a, -1.15e-93], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -2.6e-99], N[(y - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 4.8e+68], N[(x + N[(N[(y * N[(z - a), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{if}\;a \leq -2.05 \cdot 10^{+49}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -1.15 \cdot 10^{-93}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq -2.6 \cdot 10^{-99}:\\
\;\;\;\;y - \frac{y \cdot z}{a}\\
\mathbf{elif}\;a \leq 4.8 \cdot 10^{+68}:\\
\;\;\;\;x + \frac{y \cdot \left(z - a\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -2.05e49 or 4.80000000000000016e68 < a Initial program 78.7%
associate--l+79.7%
sub-neg79.7%
+-commutative79.7%
associate-/l*95.1%
distribute-neg-frac95.1%
associate-/r/95.1%
fma-def95.1%
sub-neg95.1%
+-commutative95.1%
distribute-neg-in95.1%
unsub-neg95.1%
remove-double-neg95.1%
Simplified95.1%
Taylor expanded in t around 0 79.9%
mul-1-neg79.9%
sub-neg79.9%
associate-/l*91.3%
Simplified91.3%
if -2.05e49 < a < -1.1499999999999999e-93Initial program 68.5%
associate--l+78.0%
sub-neg78.0%
+-commutative78.0%
associate-/l*77.2%
distribute-neg-frac77.2%
associate-/r/85.0%
fma-def85.0%
sub-neg85.0%
+-commutative85.0%
distribute-neg-in85.0%
unsub-neg85.0%
remove-double-neg85.0%
Simplified85.0%
Taylor expanded in y around 0 95.7%
Taylor expanded in a around 0 77.1%
if -1.1499999999999999e-93 < a < -2.60000000000000005e-99Initial program 99.6%
associate-*l/76.9%
Simplified76.9%
Taylor expanded in t around 0 76.9%
Taylor expanded in x around 0 97.0%
if -2.60000000000000005e-99 < a < 4.80000000000000016e68Initial program 73.1%
associate--l+78.4%
associate-/l*81.7%
Simplified81.7%
Taylor expanded in t around -inf 91.7%
+-commutative91.7%
sub-neg91.7%
mul-1-neg91.7%
+-commutative91.7%
*-commutative91.7%
+-commutative91.7%
*-commutative91.7%
mul-1-neg91.7%
unsub-neg91.7%
mul-1-neg91.7%
sub-neg91.7%
distribute-lft-out--91.7%
Simplified91.7%
Final simplification89.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (- x (/ y (/ a z))))))
(if (<= a -3e+50)
t_1
(if (<= a -1.8e-84)
(+ x (* y (/ z t)))
(if (<= a -2.6e-99)
(- (+ y x) (/ (* y z) a))
(if (<= a 4.8e+68) (+ x (/ (* y (- z a)) t)) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + (x - (y / (a / z)));
double tmp;
if (a <= -3e+50) {
tmp = t_1;
} else if (a <= -1.8e-84) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = (y + x) - ((y * z) / a);
} else if (a <= 4.8e+68) {
tmp = x + ((y * (z - 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) :: tmp
t_1 = y + (x - (y / (a / z)))
if (a <= (-3d+50)) then
tmp = t_1
else if (a <= (-1.8d-84)) then
tmp = x + (y * (z / t))
else if (a <= (-2.6d-99)) then
tmp = (y + x) - ((y * z) / a)
else if (a <= 4.8d+68) then
tmp = x + ((y * (z - 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 - (y / (a / z)));
double tmp;
if (a <= -3e+50) {
tmp = t_1;
} else if (a <= -1.8e-84) {
tmp = x + (y * (z / t));
} else if (a <= -2.6e-99) {
tmp = (y + x) - ((y * z) / a);
} else if (a <= 4.8e+68) {
tmp = x + ((y * (z - a)) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + (x - (y / (a / z))) tmp = 0 if a <= -3e+50: tmp = t_1 elif a <= -1.8e-84: tmp = x + (y * (z / t)) elif a <= -2.6e-99: tmp = (y + x) - ((y * z) / a) elif a <= 4.8e+68: tmp = x + ((y * (z - a)) / t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(x - Float64(y / Float64(a / z)))) tmp = 0.0 if (a <= -3e+50) tmp = t_1; elseif (a <= -1.8e-84) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= -2.6e-99) tmp = Float64(Float64(y + x) - Float64(Float64(y * z) / a)); elseif (a <= 4.8e+68) tmp = Float64(x + Float64(Float64(y * Float64(z - a)) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + (x - (y / (a / z))); tmp = 0.0; if (a <= -3e+50) tmp = t_1; elseif (a <= -1.8e-84) tmp = x + (y * (z / t)); elseif (a <= -2.6e-99) tmp = (y + x) - ((y * z) / a); elseif (a <= 4.8e+68) tmp = x + ((y * (z - a)) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3e+50], t$95$1, If[LessEqual[a, -1.8e-84], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -2.6e-99], N[(N[(y + x), $MachinePrecision] - N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 4.8e+68], N[(x + N[(N[(y * N[(z - a), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{if}\;a \leq -3 \cdot 10^{+50}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -1.8 \cdot 10^{-84}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq -2.6 \cdot 10^{-99}:\\
\;\;\;\;\left(y + x\right) - \frac{y \cdot z}{a}\\
\mathbf{elif}\;a \leq 4.8 \cdot 10^{+68}:\\
\;\;\;\;x + \frac{y \cdot \left(z - a\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -2.9999999999999998e50 or 4.80000000000000016e68 < a Initial program 78.7%
associate--l+79.7%
sub-neg79.7%
+-commutative79.7%
associate-/l*95.1%
distribute-neg-frac95.1%
associate-/r/95.1%
fma-def95.1%
sub-neg95.1%
+-commutative95.1%
distribute-neg-in95.1%
unsub-neg95.1%
remove-double-neg95.1%
Simplified95.1%
Taylor expanded in t around 0 79.9%
mul-1-neg79.9%
sub-neg79.9%
associate-/l*91.3%
Simplified91.3%
if -2.9999999999999998e50 < a < -1.80000000000000002e-84Initial program 65.4%
associate--l+75.9%
sub-neg75.9%
+-commutative75.9%
associate-/l*75.0%
distribute-neg-frac75.0%
associate-/r/83.5%
fma-def83.5%
sub-neg83.5%
+-commutative83.5%
distribute-neg-in83.5%
unsub-neg83.5%
remove-double-neg83.5%
Simplified83.5%
Taylor expanded in y around 0 95.2%
Taylor expanded in a around 0 74.9%
if -1.80000000000000002e-84 < a < -2.60000000000000005e-99Initial program 99.8%
associate--l+99.8%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 98.3%
if -2.60000000000000005e-99 < a < 4.80000000000000016e68Initial program 73.1%
associate--l+78.4%
associate-/l*81.7%
Simplified81.7%
Taylor expanded in t around -inf 91.7%
+-commutative91.7%
sub-neg91.7%
mul-1-neg91.7%
+-commutative91.7%
*-commutative91.7%
+-commutative91.7%
*-commutative91.7%
mul-1-neg91.7%
unsub-neg91.7%
mul-1-neg91.7%
sub-neg91.7%
distribute-lft-out--91.7%
Simplified91.7%
Final simplification89.7%
(FPCore (x y z t a)
:precision binary64
(if (<= t -1.1e+181)
(- x (/ y (/ t (- a z))))
(if (<= t 2.7e+55)
(+ x (+ y (/ (- t z) (/ (- a t) y))))
(+ x (* y (/ (- z a) t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.1e+181) {
tmp = x - (y / (t / (a - z)));
} else if (t <= 2.7e+55) {
tmp = x + (y + ((t - z) / ((a - t) / y)));
} else {
tmp = x + (y * ((z - 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 <= (-1.1d+181)) then
tmp = x - (y / (t / (a - z)))
else if (t <= 2.7d+55) then
tmp = x + (y + ((t - z) / ((a - t) / y)))
else
tmp = x + (y * ((z - 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 <= -1.1e+181) {
tmp = x - (y / (t / (a - z)));
} else if (t <= 2.7e+55) {
tmp = x + (y + ((t - z) / ((a - t) / y)));
} else {
tmp = x + (y * ((z - a) / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.1e+181: tmp = x - (y / (t / (a - z))) elif t <= 2.7e+55: tmp = x + (y + ((t - z) / ((a - t) / y))) else: tmp = x + (y * ((z - a) / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.1e+181) tmp = Float64(x - Float64(y / Float64(t / Float64(a - z)))); elseif (t <= 2.7e+55) tmp = Float64(x + Float64(y + Float64(Float64(t - z) / Float64(Float64(a - t) / y)))); else tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.1e+181) tmp = x - (y / (t / (a - z))); elseif (t <= 2.7e+55) tmp = x + (y + ((t - z) / ((a - t) / y))); else tmp = x + (y * ((z - a) / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.1e+181], N[(x - N[(y / N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.7e+55], N[(x + N[(y + N[(N[(t - z), $MachinePrecision] / N[(N[(a - t), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.1 \cdot 10^{+181}:\\
\;\;\;\;x - \frac{y}{\frac{t}{a - z}}\\
\mathbf{elif}\;t \leq 2.7 \cdot 10^{+55}:\\
\;\;\;\;x + \left(y + \frac{t - z}{\frac{a - t}{y}}\right)\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\end{array}
\end{array}
if t < -1.1000000000000001e181Initial program 47.7%
associate--l+54.8%
associate-/l*57.8%
Simplified57.8%
Taylor expanded in t around -inf 85.6%
Taylor expanded in y around 0 85.6%
associate-/l*96.4%
Simplified96.4%
if -1.1000000000000001e181 < t < 2.69999999999999977e55Initial program 85.5%
associate--l+89.0%
associate-/l*94.3%
Simplified94.3%
if 2.69999999999999977e55 < t Initial program 52.4%
associate--l+57.1%
sub-neg57.1%
+-commutative57.1%
associate-/l*76.1%
distribute-neg-frac76.1%
associate-/r/80.0%
fma-def80.0%
sub-neg80.0%
+-commutative80.0%
distribute-neg-in80.0%
unsub-neg80.0%
remove-double-neg80.0%
Simplified80.0%
Taylor expanded in y around 0 88.8%
Taylor expanded in t around inf 93.2%
associate-*r/93.2%
distribute-lft-in93.2%
neg-mul-193.2%
mul-1-neg93.2%
remove-double-neg93.2%
neg-mul-193.2%
sub-neg93.2%
Simplified93.2%
Final simplification94.4%
(FPCore (x y z t a)
:precision binary64
(if (<= a -4e+50)
(+ y (- x (/ y (/ a z))))
(if (<= a 8e+107)
(- x (/ (* y z) (- a t)))
(+ x (- y (/ y (/ (- a t) z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -4e+50) {
tmp = y + (x - (y / (a / z)));
} else if (a <= 8e+107) {
tmp = x - ((y * z) / (a - t));
} else {
tmp = x + (y - (y / ((a - t) / z)));
}
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 <= (-4d+50)) then
tmp = y + (x - (y / (a / z)))
else if (a <= 8d+107) then
tmp = x - ((y * z) / (a - t))
else
tmp = x + (y - (y / ((a - t) / z)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -4e+50) {
tmp = y + (x - (y / (a / z)));
} else if (a <= 8e+107) {
tmp = x - ((y * z) / (a - t));
} else {
tmp = x + (y - (y / ((a - t) / z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -4e+50: tmp = y + (x - (y / (a / z))) elif a <= 8e+107: tmp = x - ((y * z) / (a - t)) else: tmp = x + (y - (y / ((a - t) / z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -4e+50) tmp = Float64(y + Float64(x - Float64(y / Float64(a / z)))); elseif (a <= 8e+107) tmp = Float64(x - Float64(Float64(y * z) / Float64(a - t))); else tmp = Float64(x + Float64(y - Float64(y / Float64(Float64(a - t) / z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -4e+50) tmp = y + (x - (y / (a / z))); elseif (a <= 8e+107) tmp = x - ((y * z) / (a - t)); else tmp = x + (y - (y / ((a - t) / z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -4e+50], N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 8e+107], N[(x - N[(N[(y * z), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y - N[(y / N[(N[(a - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4 \cdot 10^{+50}:\\
\;\;\;\;y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{elif}\;a \leq 8 \cdot 10^{+107}:\\
\;\;\;\;x - \frac{y \cdot z}{a - t}\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - \frac{y}{\frac{a - t}{z}}\right)\\
\end{array}
\end{array}
if a < -4.0000000000000003e50Initial program 78.0%
associate--l+79.9%
sub-neg79.9%
+-commutative79.9%
associate-/l*96.3%
distribute-neg-frac96.3%
associate-/r/96.3%
fma-def96.3%
sub-neg96.3%
+-commutative96.3%
distribute-neg-in96.3%
unsub-neg96.3%
remove-double-neg96.3%
Simplified96.3%
Taylor expanded in t around 0 80.0%
mul-1-neg80.0%
sub-neg80.0%
associate-/l*92.6%
Simplified92.6%
if -4.0000000000000003e50 < a < 7.9999999999999998e107Initial program 72.9%
associate--l+78.7%
sub-neg78.7%
+-commutative78.7%
associate-/l*81.5%
distribute-neg-frac81.5%
associate-/r/83.9%
fma-def83.9%
sub-neg83.9%
+-commutative83.9%
distribute-neg-in83.9%
unsub-neg83.9%
remove-double-neg83.9%
Simplified83.9%
Taylor expanded in z around inf 92.7%
associate-*r/92.7%
associate-*r*92.7%
neg-mul-192.7%
Simplified92.7%
if 7.9999999999999998e107 < a Initial program 80.5%
associate--l+80.4%
associate-/l*95.0%
Simplified95.0%
Taylor expanded in z around inf 80.4%
associate-/l*94.0%
Simplified94.0%
Final simplification92.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ y (/ t z))))
(if (<= a -4.7e-176)
(+ y x)
(if (<= a 1.5e-286)
t_1
(if (<= a 1.8e-229)
x
(if (<= a 3.7e-175) t_1 (if (<= a 8e+107) x (+ y x))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y / (t / z);
double tmp;
if (a <= -4.7e-176) {
tmp = y + x;
} else if (a <= 1.5e-286) {
tmp = t_1;
} else if (a <= 1.8e-229) {
tmp = x;
} else if (a <= 3.7e-175) {
tmp = t_1;
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + 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) :: t_1
real(8) :: tmp
t_1 = y / (t / z)
if (a <= (-4.7d-176)) then
tmp = y + x
else if (a <= 1.5d-286) then
tmp = t_1
else if (a <= 1.8d-229) then
tmp = x
else if (a <= 3.7d-175) then
tmp = t_1
else if (a <= 8d+107) then
tmp = x
else
tmp = y + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = y / (t / z);
double tmp;
if (a <= -4.7e-176) {
tmp = y + x;
} else if (a <= 1.5e-286) {
tmp = t_1;
} else if (a <= 1.8e-229) {
tmp = x;
} else if (a <= 3.7e-175) {
tmp = t_1;
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y / (t / z) tmp = 0 if a <= -4.7e-176: tmp = y + x elif a <= 1.5e-286: tmp = t_1 elif a <= 1.8e-229: tmp = x elif a <= 3.7e-175: tmp = t_1 elif a <= 8e+107: tmp = x else: tmp = y + x return tmp
function code(x, y, z, t, a) t_1 = Float64(y / Float64(t / z)) tmp = 0.0 if (a <= -4.7e-176) tmp = Float64(y + x); elseif (a <= 1.5e-286) tmp = t_1; elseif (a <= 1.8e-229) tmp = x; elseif (a <= 3.7e-175) tmp = t_1; elseif (a <= 8e+107) tmp = x; else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y / (t / z); tmp = 0.0; if (a <= -4.7e-176) tmp = y + x; elseif (a <= 1.5e-286) tmp = t_1; elseif (a <= 1.8e-229) tmp = x; elseif (a <= 3.7e-175) tmp = t_1; elseif (a <= 8e+107) tmp = x; else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y / N[(t / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -4.7e-176], N[(y + x), $MachinePrecision], If[LessEqual[a, 1.5e-286], t$95$1, If[LessEqual[a, 1.8e-229], x, If[LessEqual[a, 3.7e-175], t$95$1, If[LessEqual[a, 8e+107], x, N[(y + x), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{\frac{t}{z}}\\
\mathbf{if}\;a \leq -4.7 \cdot 10^{-176}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 1.5 \cdot 10^{-286}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq 1.8 \cdot 10^{-229}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 3.7 \cdot 10^{-175}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq 8 \cdot 10^{+107}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -4.69999999999999984e-176 or 7.9999999999999998e107 < a Initial program 77.1%
associate--l+80.3%
sub-neg80.3%
+-commutative80.3%
associate-/l*89.3%
distribute-neg-frac89.3%
associate-/r/91.2%
fma-def91.2%
sub-neg91.2%
+-commutative91.2%
distribute-neg-in91.2%
unsub-neg91.2%
remove-double-neg91.2%
Simplified91.2%
Taylor expanded in a around inf 72.0%
if -4.69999999999999984e-176 < a < 1.5e-286 or 1.80000000000000001e-229 < a < 3.69999999999999998e-175Initial program 70.1%
associate--l+74.7%
associate-/l*76.6%
Simplified76.6%
Taylor expanded in t around -inf 97.0%
Taylor expanded in y around 0 96.9%
associate-/l*94.4%
Simplified94.4%
Taylor expanded in z around inf 66.1%
associate-/l*64.1%
Simplified64.1%
if 1.5e-286 < a < 1.80000000000000001e-229 or 3.69999999999999998e-175 < a < 7.9999999999999998e107Initial program 73.9%
associate--l+79.7%
sub-neg79.7%
+-commutative79.7%
associate-/l*87.1%
distribute-neg-frac87.1%
associate-/r/87.1%
fma-def87.1%
sub-neg87.1%
+-commutative87.1%
distribute-neg-in87.1%
unsub-neg87.1%
remove-double-neg87.1%
Simplified87.1%
Taylor expanded in x around inf 59.8%
Final simplification67.4%
(FPCore (x y z t a)
:precision binary64
(if (<= a -2.5e-181)
(+ y x)
(if (<= a 7.5e-287)
(/ (* y z) t)
(if (<= a 1.9e-230)
x
(if (<= a 4e-176) (/ y (/ t z)) (if (<= a 8e+107) x (+ y x)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.5e-181) {
tmp = y + x;
} else if (a <= 7.5e-287) {
tmp = (y * z) / t;
} else if (a <= 1.9e-230) {
tmp = x;
} else if (a <= 4e-176) {
tmp = y / (t / z);
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + 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 <= (-2.5d-181)) then
tmp = y + x
else if (a <= 7.5d-287) then
tmp = (y * z) / t
else if (a <= 1.9d-230) then
tmp = x
else if (a <= 4d-176) then
tmp = y / (t / z)
else if (a <= 8d+107) then
tmp = x
else
tmp = y + 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 <= -2.5e-181) {
tmp = y + x;
} else if (a <= 7.5e-287) {
tmp = (y * z) / t;
} else if (a <= 1.9e-230) {
tmp = x;
} else if (a <= 4e-176) {
tmp = y / (t / z);
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.5e-181: tmp = y + x elif a <= 7.5e-287: tmp = (y * z) / t elif a <= 1.9e-230: tmp = x elif a <= 4e-176: tmp = y / (t / z) elif a <= 8e+107: tmp = x else: tmp = y + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.5e-181) tmp = Float64(y + x); elseif (a <= 7.5e-287) tmp = Float64(Float64(y * z) / t); elseif (a <= 1.9e-230) tmp = x; elseif (a <= 4e-176) tmp = Float64(y / Float64(t / z)); elseif (a <= 8e+107) tmp = x; else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -2.5e-181) tmp = y + x; elseif (a <= 7.5e-287) tmp = (y * z) / t; elseif (a <= 1.9e-230) tmp = x; elseif (a <= 4e-176) tmp = y / (t / z); elseif (a <= 8e+107) tmp = x; else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.5e-181], N[(y + x), $MachinePrecision], If[LessEqual[a, 7.5e-287], N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[a, 1.9e-230], x, If[LessEqual[a, 4e-176], N[(y / N[(t / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 8e+107], x, N[(y + x), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.5 \cdot 10^{-181}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 7.5 \cdot 10^{-287}:\\
\;\;\;\;\frac{y \cdot z}{t}\\
\mathbf{elif}\;a \leq 1.9 \cdot 10^{-230}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 4 \cdot 10^{-176}:\\
\;\;\;\;\frac{y}{\frac{t}{z}}\\
\mathbf{elif}\;a \leq 8 \cdot 10^{+107}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -2.5000000000000001e-181 or 7.9999999999999998e107 < a Initial program 77.1%
associate--l+80.3%
sub-neg80.3%
+-commutative80.3%
associate-/l*89.3%
distribute-neg-frac89.3%
associate-/r/91.2%
fma-def91.2%
sub-neg91.2%
+-commutative91.2%
distribute-neg-in91.2%
unsub-neg91.2%
remove-double-neg91.2%
Simplified91.2%
Taylor expanded in a around inf 72.0%
if -2.5000000000000001e-181 < a < 7.5000000000000001e-287Initial program 67.4%
associate--l+70.6%
associate-/l*73.1%
Simplified73.1%
Taylor expanded in t around -inf 96.0%
Taylor expanded in y around 0 95.9%
associate-/l*92.3%
Simplified92.3%
Taylor expanded in z around inf 62.6%
if 7.5000000000000001e-287 < a < 1.8999999999999999e-230 or 4e-176 < a < 7.9999999999999998e107Initial program 73.9%
associate--l+79.7%
sub-neg79.7%
+-commutative79.7%
associate-/l*87.1%
distribute-neg-frac87.1%
associate-/r/87.1%
fma-def87.1%
sub-neg87.1%
+-commutative87.1%
distribute-neg-in87.1%
unsub-neg87.1%
remove-double-neg87.1%
Simplified87.1%
Taylor expanded in x around inf 59.8%
if 1.8999999999999999e-230 < a < 4e-176Initial program 77.2%
associate--l+85.8%
associate-/l*85.9%
Simplified85.9%
Taylor expanded in t around -inf 99.6%
Taylor expanded in y around 0 99.6%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in z around inf 75.3%
associate-/l*75.7%
Simplified75.7%
Final simplification67.8%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -5.5e-31) (not (<= a 4.8e+68))) (+ y (- x (/ y (/ a z)))) (+ x (/ (* y z) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -5.5e-31) || !(a <= 4.8e+68)) {
tmp = y + (x - (y / (a / z)));
} 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 <= (-5.5d-31)) .or. (.not. (a <= 4.8d+68))) then
tmp = y + (x - (y / (a / z)))
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 <= -5.5e-31) || !(a <= 4.8e+68)) {
tmp = y + (x - (y / (a / z)));
} else {
tmp = x + ((y * z) / t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -5.5e-31) or not (a <= 4.8e+68): tmp = y + (x - (y / (a / z))) else: tmp = x + ((y * z) / t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -5.5e-31) || !(a <= 4.8e+68)) tmp = Float64(y + Float64(x - Float64(y / Float64(a / z)))); else tmp = Float64(x + Float64(Float64(y * z) / t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -5.5e-31) || ~((a <= 4.8e+68))) tmp = y + (x - (y / (a / z))); else tmp = x + ((y * z) / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -5.5e-31], N[Not[LessEqual[a, 4.8e+68]], $MachinePrecision]], N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.5 \cdot 10^{-31} \lor \neg \left(a \leq 4.8 \cdot 10^{+68}\right):\\
\;\;\;\;y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot z}{t}\\
\end{array}
\end{array}
if a < -5.49999999999999958e-31 or 4.80000000000000016e68 < a Initial program 77.0%
associate--l+78.2%
sub-neg78.2%
+-commutative78.2%
associate-/l*93.1%
distribute-neg-frac93.1%
associate-/r/93.1%
fma-def93.1%
sub-neg93.1%
+-commutative93.1%
distribute-neg-in93.1%
unsub-neg93.1%
remove-double-neg93.1%
Simplified93.1%
Taylor expanded in t around 0 78.5%
mul-1-neg78.5%
sub-neg78.5%
associate-/l*88.6%
Simplified88.6%
if -5.49999999999999958e-31 < a < 4.80000000000000016e68Initial program 73.6%
associate--l+80.0%
sub-neg80.0%
+-commutative80.0%
associate-/l*81.5%
distribute-neg-frac81.5%
associate-/r/84.2%
fma-def84.2%
sub-neg84.2%
+-commutative84.2%
distribute-neg-in84.2%
unsub-neg84.2%
remove-double-neg84.2%
Simplified84.2%
Taylor expanded in y around 0 94.7%
Taylor expanded in a around 0 85.5%
Final simplification86.8%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -1e+54) (not (<= a 8e+107))) (+ y (- x (/ y (/ a z)))) (- x (/ (* y z) (- a t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -1e+54) || !(a <= 8e+107)) {
tmp = y + (x - (y / (a / z)));
} else {
tmp = x - ((y * z) / (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 ((a <= (-1d+54)) .or. (.not. (a <= 8d+107))) then
tmp = y + (x - (y / (a / z)))
else
tmp = x - ((y * z) / (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 ((a <= -1e+54) || !(a <= 8e+107)) {
tmp = y + (x - (y / (a / z)));
} else {
tmp = x - ((y * z) / (a - t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -1e+54) or not (a <= 8e+107): tmp = y + (x - (y / (a / z))) else: tmp = x - ((y * z) / (a - t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -1e+54) || !(a <= 8e+107)) tmp = Float64(y + Float64(x - Float64(y / Float64(a / z)))); else tmp = Float64(x - Float64(Float64(y * z) / Float64(a - t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -1e+54) || ~((a <= 8e+107))) tmp = y + (x - (y / (a / z))); else tmp = x - ((y * z) / (a - t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -1e+54], N[Not[LessEqual[a, 8e+107]], $MachinePrecision]], N[(y + N[(x - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[(y * z), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1 \cdot 10^{+54} \lor \neg \left(a \leq 8 \cdot 10^{+107}\right):\\
\;\;\;\;y + \left(x - \frac{y}{\frac{a}{z}}\right)\\
\mathbf{else}:\\
\;\;\;\;x - \frac{y \cdot z}{a - t}\\
\end{array}
\end{array}
if a < -1.0000000000000001e54 or 7.9999999999999998e107 < a Initial program 79.0%
associate--l+80.1%
sub-neg80.1%
+-commutative80.1%
associate-/l*95.7%
distribute-neg-frac95.7%
associate-/r/95.7%
fma-def95.7%
sub-neg95.7%
+-commutative95.7%
distribute-neg-in95.7%
unsub-neg95.7%
remove-double-neg95.7%
Simplified95.7%
Taylor expanded in t around 0 80.2%
mul-1-neg80.2%
sub-neg80.2%
associate-/l*92.6%
Simplified92.6%
if -1.0000000000000001e54 < a < 7.9999999999999998e107Initial program 72.9%
associate--l+78.7%
sub-neg78.7%
+-commutative78.7%
associate-/l*81.5%
distribute-neg-frac81.5%
associate-/r/83.9%
fma-def83.9%
sub-neg83.9%
+-commutative83.9%
distribute-neg-in83.9%
unsub-neg83.9%
remove-double-neg83.9%
Simplified83.9%
Taylor expanded in z around inf 92.7%
associate-*r/92.7%
associate-*r*92.7%
neg-mul-192.7%
Simplified92.7%
Final simplification92.7%
(FPCore (x y z t a) :precision binary64 (if (<= a -8.5e+49) (+ y x) (if (<= a 3.8e+27) (+ x (* y (/ z t))) (+ y x))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -8.5e+49) {
tmp = y + x;
} else if (a <= 3.8e+27) {
tmp = x + (y * (z / t));
} else {
tmp = y + 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 <= (-8.5d+49)) then
tmp = y + x
else if (a <= 3.8d+27) then
tmp = x + (y * (z / t))
else
tmp = y + 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 <= -8.5e+49) {
tmp = y + x;
} else if (a <= 3.8e+27) {
tmp = x + (y * (z / t));
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -8.5e+49: tmp = y + x elif a <= 3.8e+27: tmp = x + (y * (z / t)) else: tmp = y + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -8.5e+49) tmp = Float64(y + x); elseif (a <= 3.8e+27) tmp = Float64(x + Float64(y * Float64(z / t))); else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -8.5e+49) tmp = y + x; elseif (a <= 3.8e+27) tmp = x + (y * (z / t)); else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -8.5e+49], N[(y + x), $MachinePrecision], If[LessEqual[a, 3.8e+27], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -8.5 \cdot 10^{+49}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 3.8 \cdot 10^{+27}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -8.4999999999999996e49 or 3.80000000000000022e27 < a Initial program 79.6%
associate--l+80.5%
sub-neg80.5%
+-commutative80.5%
associate-/l*95.3%
distribute-neg-frac95.3%
associate-/r/95.3%
fma-def95.3%
sub-neg95.3%
+-commutative95.3%
distribute-neg-in95.3%
unsub-neg95.3%
remove-double-neg95.3%
Simplified95.3%
Taylor expanded in a around inf 82.5%
if -8.4999999999999996e49 < a < 3.80000000000000022e27Initial program 72.1%
associate--l+78.3%
sub-neg78.3%
+-commutative78.3%
associate-/l*80.7%
distribute-neg-frac80.7%
associate-/r/83.2%
fma-def83.2%
sub-neg83.2%
+-commutative83.2%
distribute-neg-in83.2%
unsub-neg83.2%
remove-double-neg83.2%
Simplified83.2%
Taylor expanded in y around 0 94.3%
Taylor expanded in a around 0 83.1%
Final simplification82.9%
(FPCore (x y z t a) :precision binary64 (if (<= a -2.05e+49) (+ y x) (if (<= a 3.7e+27) (+ x (/ y (/ t z))) (+ y x))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.05e+49) {
tmp = y + x;
} else if (a <= 3.7e+27) {
tmp = x + (y / (t / z));
} else {
tmp = y + 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 <= (-2.05d+49)) then
tmp = y + x
else if (a <= 3.7d+27) then
tmp = x + (y / (t / z))
else
tmp = y + 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 <= -2.05e+49) {
tmp = y + x;
} else if (a <= 3.7e+27) {
tmp = x + (y / (t / z));
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.05e+49: tmp = y + x elif a <= 3.7e+27: tmp = x + (y / (t / z)) else: tmp = y + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.05e+49) tmp = Float64(y + x); elseif (a <= 3.7e+27) tmp = Float64(x + Float64(y / Float64(t / z))); else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -2.05e+49) tmp = y + x; elseif (a <= 3.7e+27) tmp = x + (y / (t / z)); else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.05e+49], N[(y + x), $MachinePrecision], If[LessEqual[a, 3.7e+27], N[(x + N[(y / N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.05 \cdot 10^{+49}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 3.7 \cdot 10^{+27}:\\
\;\;\;\;x + \frac{y}{\frac{t}{z}}\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -2.05e49 or 3.70000000000000002e27 < a Initial program 79.6%
associate--l+80.5%
sub-neg80.5%
+-commutative80.5%
associate-/l*95.3%
distribute-neg-frac95.3%
associate-/r/95.3%
fma-def95.3%
sub-neg95.3%
+-commutative95.3%
distribute-neg-in95.3%
unsub-neg95.3%
remove-double-neg95.3%
Simplified95.3%
Taylor expanded in a around inf 82.5%
if -2.05e49 < a < 3.70000000000000002e27Initial program 72.1%
associate--l+78.3%
sub-neg78.3%
+-commutative78.3%
associate-/l*80.7%
distribute-neg-frac80.7%
associate-/r/83.2%
fma-def83.2%
sub-neg83.2%
+-commutative83.2%
distribute-neg-in83.2%
unsub-neg83.2%
remove-double-neg83.2%
Simplified83.2%
Taylor expanded in y around 0 94.3%
Taylor expanded in a around 0 83.2%
associate-/l*83.1%
Simplified83.1%
Final simplification82.9%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.55e+50) (+ y x) (if (<= a 8e+107) (+ x (/ (* y z) t)) (+ y x))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.55e+50) {
tmp = y + x;
} else if (a <= 8e+107) {
tmp = x + ((y * z) / t);
} else {
tmp = y + 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 <= (-1.55d+50)) then
tmp = y + x
else if (a <= 8d+107) then
tmp = x + ((y * z) / t)
else
tmp = y + 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 <= -1.55e+50) {
tmp = y + x;
} else if (a <= 8e+107) {
tmp = x + ((y * z) / t);
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.55e+50: tmp = y + x elif a <= 8e+107: tmp = x + ((y * z) / t) else: tmp = y + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.55e+50) tmp = Float64(y + x); elseif (a <= 8e+107) tmp = Float64(x + Float64(Float64(y * z) / t)); else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.55e+50) tmp = y + x; elseif (a <= 8e+107) tmp = x + ((y * z) / t); else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.55e+50], N[(y + x), $MachinePrecision], If[LessEqual[a, 8e+107], N[(x + N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(y + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.55 \cdot 10^{+50}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 8 \cdot 10^{+107}:\\
\;\;\;\;x + \frac{y \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -1.55000000000000001e50 or 7.9999999999999998e107 < a Initial program 79.0%
associate--l+80.1%
sub-neg80.1%
+-commutative80.1%
associate-/l*95.7%
distribute-neg-frac95.7%
associate-/r/95.7%
fma-def95.7%
sub-neg95.7%
+-commutative95.7%
distribute-neg-in95.7%
unsub-neg95.7%
remove-double-neg95.7%
Simplified95.7%
Taylor expanded in a around inf 85.3%
if -1.55000000000000001e50 < a < 7.9999999999999998e107Initial program 72.9%
associate--l+78.7%
sub-neg78.7%
+-commutative78.7%
associate-/l*81.5%
distribute-neg-frac81.5%
associate-/r/83.9%
fma-def83.9%
sub-neg83.9%
+-commutative83.9%
distribute-neg-in83.9%
unsub-neg83.9%
remove-double-neg83.9%
Simplified83.9%
Taylor expanded in y around 0 94.2%
Taylor expanded in a around 0 82.1%
Final simplification83.2%
(FPCore (x y z t a) :precision binary64 (if (<= a -2.4e-23) (+ y x) (if (<= a 8e+107) x (+ y x))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.4e-23) {
tmp = y + x;
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + 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 <= (-2.4d-23)) then
tmp = y + x
else if (a <= 8d+107) then
tmp = x
else
tmp = y + 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 <= -2.4e-23) {
tmp = y + x;
} else if (a <= 8e+107) {
tmp = x;
} else {
tmp = y + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.4e-23: tmp = y + x elif a <= 8e+107: tmp = x else: tmp = y + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.4e-23) tmp = Float64(y + x); elseif (a <= 8e+107) tmp = x; else tmp = Float64(y + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -2.4e-23) tmp = y + x; elseif (a <= 8e+107) tmp = x; else tmp = y + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.4e-23], N[(y + x), $MachinePrecision], If[LessEqual[a, 8e+107], x, N[(y + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.4 \cdot 10^{-23}:\\
\;\;\;\;y + x\\
\mathbf{elif}\;a \leq 8 \cdot 10^{+107}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y + x\\
\end{array}
\end{array}
if a < -2.39999999999999996e-23 or 7.9999999999999998e107 < a Initial program 77.6%
associate--l+78.6%
sub-neg78.6%
+-commutative78.6%
associate-/l*93.9%
distribute-neg-frac93.9%
associate-/r/93.9%
fma-def93.9%
sub-neg93.9%
+-commutative93.9%
distribute-neg-in93.9%
unsub-neg93.9%
remove-double-neg93.9%
Simplified93.9%
Taylor expanded in a around inf 81.2%
if -2.39999999999999996e-23 < a < 7.9999999999999998e107Initial program 73.4%
associate--l+79.6%
sub-neg79.6%
+-commutative79.6%
associate-/l*81.7%
distribute-neg-frac81.7%
associate-/r/84.2%
fma-def84.2%
sub-neg84.2%
+-commutative84.2%
distribute-neg-in84.2%
unsub-neg84.2%
remove-double-neg84.2%
Simplified84.2%
Taylor expanded in x around inf 50.9%
Final simplification63.0%
(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.1%
associate--l+79.2%
sub-neg79.2%
+-commutative79.2%
associate-/l*86.6%
distribute-neg-frac86.6%
associate-/r/88.1%
fma-def88.1%
sub-neg88.1%
+-commutative88.1%
distribute-neg-in88.1%
unsub-neg88.1%
remove-double-neg88.1%
Simplified88.1%
Taylor expanded in x around inf 52.1%
Final simplification52.1%
(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 2023200
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(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))))