
(FPCore (x y z t a) :precision binary64 (/ (- x (* y z)) (- t (* a z))))
double code(double x, double y, double z, double t, double a) {
return (x - (y * z)) / (t - (a * z));
}
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 - (a * z))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x - (y * z)) / (t - (a * z));
}
def code(x, y, z, t, a): return (x - (y * z)) / (t - (a * z))
function code(x, y, z, t, a) return Float64(Float64(x - Float64(y * z)) / Float64(t - Float64(a * z))) end
function tmp = code(x, y, z, t, a) tmp = (x - (y * z)) / (t - (a * z)); end
code[x_, y_, z_, t_, a_] := N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / N[(t - N[(a * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y \cdot z}{t - a \cdot z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (/ (- x (* y z)) (- t (* a z))))
double code(double x, double y, double z, double t, double a) {
return (x - (y * z)) / (t - (a * z));
}
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 - (a * z))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x - (y * z)) / (t - (a * z));
}
def code(x, y, z, t, a): return (x - (y * z)) / (t - (a * z))
function code(x, y, z, t, a) return Float64(Float64(x - Float64(y * z)) / Float64(t - Float64(a * z))) end
function tmp = code(x, y, z, t, a) tmp = (x - (y * z)) / (t - (a * z)); end
code[x_, y_, z_, t_, a_] := N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / N[(t - N[(a * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y \cdot z}{t - a \cdot z}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- x (* y z)) (- t (* a z)))))
(if (<= t_1 -5e-324)
t_1
(if (<= t_1 0.0)
(/ y (- a (/ t z)))
(if (<= t_1 1e+278) t_1 (/ (- y (/ x z)) a))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (x - (y * z)) / (t - (a * z));
double tmp;
if (t_1 <= -5e-324) {
tmp = t_1;
} else if (t_1 <= 0.0) {
tmp = y / (a - (t / z));
} else if (t_1 <= 1e+278) {
tmp = t_1;
} else {
tmp = (y - (x / z)) / a;
}
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 * z)) / (t - (a * z))
if (t_1 <= (-5d-324)) then
tmp = t_1
else if (t_1 <= 0.0d0) then
tmp = y / (a - (t / z))
else if (t_1 <= 1d+278) then
tmp = t_1
else
tmp = (y - (x / z)) / a
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 * z)) / (t - (a * z));
double tmp;
if (t_1 <= -5e-324) {
tmp = t_1;
} else if (t_1 <= 0.0) {
tmp = y / (a - (t / z));
} else if (t_1 <= 1e+278) {
tmp = t_1;
} else {
tmp = (y - (x / z)) / a;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x - (y * z)) / (t - (a * z)) tmp = 0 if t_1 <= -5e-324: tmp = t_1 elif t_1 <= 0.0: tmp = y / (a - (t / z)) elif t_1 <= 1e+278: tmp = t_1 else: tmp = (y - (x / z)) / a return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(x - Float64(y * z)) / Float64(t - Float64(a * z))) tmp = 0.0 if (t_1 <= -5e-324) tmp = t_1; elseif (t_1 <= 0.0) tmp = Float64(y / Float64(a - Float64(t / z))); elseif (t_1 <= 1e+278) tmp = t_1; else tmp = Float64(Float64(y - Float64(x / z)) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (x - (y * z)) / (t - (a * z)); tmp = 0.0; if (t_1 <= -5e-324) tmp = t_1; elseif (t_1 <= 0.0) tmp = y / (a - (t / z)); elseif (t_1 <= 1e+278) tmp = t_1; else tmp = (y - (x / z)) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / N[(t - N[(a * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e-324], t$95$1, If[LessEqual[t$95$1, 0.0], N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+278], t$95$1, N[(N[(y - N[(x / z), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x - y \cdot z}{t - a \cdot z}\\
\mathbf{if}\;t_1 \leq -5 \cdot 10^{-324}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_1 \leq 0:\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\mathbf{elif}\;t_1 \leq 10^{+278}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y - \frac{x}{z}}{a}\\
\end{array}
\end{array}
if (/.f64 (-.f64 x (*.f64 y z)) (-.f64 t (*.f64 a z))) < -4.94066e-324 or 0.0 < (/.f64 (-.f64 x (*.f64 y z)) (-.f64 t (*.f64 a z))) < 9.99999999999999964e277Initial program 97.9%
if -4.94066e-324 < (/.f64 (-.f64 x (*.f64 y z)) (-.f64 t (*.f64 a z))) < 0.0Initial program 57.9%
sub-neg57.9%
+-commutative57.9%
neg-sub057.9%
associate-+l-57.9%
sub0-neg57.9%
neg-mul-157.9%
sub-neg57.9%
+-commutative57.9%
neg-sub057.9%
associate-+l-57.9%
sub0-neg57.9%
neg-mul-157.9%
times-frac57.9%
metadata-eval57.9%
*-lft-identity57.9%
*-commutative57.9%
Simplified57.9%
div-sub57.9%
associate-/l*57.9%
Applied egg-rr57.9%
Taylor expanded in z around 0 79.0%
neg-mul-179.0%
unsub-neg79.0%
Simplified79.0%
Taylor expanded in y around inf 79.0%
if 9.99999999999999964e277 < (/.f64 (-.f64 x (*.f64 y z)) (-.f64 t (*.f64 a z))) Initial program 46.1%
sub-neg46.1%
+-commutative46.1%
neg-sub046.1%
associate-+l-46.1%
sub0-neg46.1%
neg-mul-146.1%
sub-neg46.1%
+-commutative46.1%
neg-sub046.1%
associate-+l-46.1%
sub0-neg46.1%
neg-mul-146.1%
times-frac46.1%
metadata-eval46.1%
*-lft-identity46.1%
*-commutative46.1%
Simplified46.1%
div-sub46.1%
associate-/l*69.9%
Applied egg-rr69.9%
Taylor expanded in a around inf 89.7%
Final simplification95.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (* a z) t)) (t_2 (/ (- x) t_1)))
(if (<= z -1.35e+20)
(/ (- y (/ x z)) a)
(if (<= z -2.45e-188)
t_2
(if (<= z 1e-235)
(/ (- x (* y z)) t)
(if (<= z 5e-82)
t_2
(if (<= z 9.5e-25)
(* z (/ y t_1))
(if (<= z 60000000000.0) t_2 (/ y (- a (/ t z)))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (a * z) - t;
double t_2 = -x / t_1;
double tmp;
if (z <= -1.35e+20) {
tmp = (y - (x / z)) / a;
} else if (z <= -2.45e-188) {
tmp = t_2;
} else if (z <= 1e-235) {
tmp = (x - (y * z)) / t;
} else if (z <= 5e-82) {
tmp = t_2;
} else if (z <= 9.5e-25) {
tmp = z * (y / t_1);
} else if (z <= 60000000000.0) {
tmp = t_2;
} else {
tmp = 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (a * z) - t
t_2 = -x / t_1
if (z <= (-1.35d+20)) then
tmp = (y - (x / z)) / a
else if (z <= (-2.45d-188)) then
tmp = t_2
else if (z <= 1d-235) then
tmp = (x - (y * z)) / t
else if (z <= 5d-82) then
tmp = t_2
else if (z <= 9.5d-25) then
tmp = z * (y / t_1)
else if (z <= 60000000000.0d0) then
tmp = t_2
else
tmp = 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 t_1 = (a * z) - t;
double t_2 = -x / t_1;
double tmp;
if (z <= -1.35e+20) {
tmp = (y - (x / z)) / a;
} else if (z <= -2.45e-188) {
tmp = t_2;
} else if (z <= 1e-235) {
tmp = (x - (y * z)) / t;
} else if (z <= 5e-82) {
tmp = t_2;
} else if (z <= 9.5e-25) {
tmp = z * (y / t_1);
} else if (z <= 60000000000.0) {
tmp = t_2;
} else {
tmp = y / (a - (t / z));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (a * z) - t t_2 = -x / t_1 tmp = 0 if z <= -1.35e+20: tmp = (y - (x / z)) / a elif z <= -2.45e-188: tmp = t_2 elif z <= 1e-235: tmp = (x - (y * z)) / t elif z <= 5e-82: tmp = t_2 elif z <= 9.5e-25: tmp = z * (y / t_1) elif z <= 60000000000.0: tmp = t_2 else: tmp = y / (a - (t / z)) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(a * z) - t) t_2 = Float64(Float64(-x) / t_1) tmp = 0.0 if (z <= -1.35e+20) tmp = Float64(Float64(y - Float64(x / z)) / a); elseif (z <= -2.45e-188) tmp = t_2; elseif (z <= 1e-235) tmp = Float64(Float64(x - Float64(y * z)) / t); elseif (z <= 5e-82) tmp = t_2; elseif (z <= 9.5e-25) tmp = Float64(z * Float64(y / t_1)); elseif (z <= 60000000000.0) tmp = t_2; else tmp = Float64(y / Float64(a - Float64(t / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (a * z) - t; t_2 = -x / t_1; tmp = 0.0; if (z <= -1.35e+20) tmp = (y - (x / z)) / a; elseif (z <= -2.45e-188) tmp = t_2; elseif (z <= 1e-235) tmp = (x - (y * z)) / t; elseif (z <= 5e-82) tmp = t_2; elseif (z <= 9.5e-25) tmp = z * (y / t_1); elseif (z <= 60000000000.0) tmp = t_2; else tmp = y / (a - (t / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(a * z), $MachinePrecision] - t), $MachinePrecision]}, Block[{t$95$2 = N[((-x) / t$95$1), $MachinePrecision]}, If[LessEqual[z, -1.35e+20], N[(N[(y - N[(x / z), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[z, -2.45e-188], t$95$2, If[LessEqual[z, 1e-235], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[z, 5e-82], t$95$2, If[LessEqual[z, 9.5e-25], N[(z * N[(y / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 60000000000.0], t$95$2, N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot z - t\\
t_2 := \frac{-x}{t_1}\\
\mathbf{if}\;z \leq -1.35 \cdot 10^{+20}:\\
\;\;\;\;\frac{y - \frac{x}{z}}{a}\\
\mathbf{elif}\;z \leq -2.45 \cdot 10^{-188}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;z \leq 10^{-235}:\\
\;\;\;\;\frac{x - y \cdot z}{t}\\
\mathbf{elif}\;z \leq 5 \cdot 10^{-82}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;z \leq 9.5 \cdot 10^{-25}:\\
\;\;\;\;z \cdot \frac{y}{t_1}\\
\mathbf{elif}\;z \leq 60000000000:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\end{array}
\end{array}
if z < -1.35e20Initial program 71.7%
sub-neg71.7%
+-commutative71.7%
neg-sub071.7%
associate-+l-71.7%
sub0-neg71.7%
neg-mul-171.7%
sub-neg71.7%
+-commutative71.7%
neg-sub071.7%
associate-+l-71.7%
sub0-neg71.7%
neg-mul-171.7%
times-frac71.7%
metadata-eval71.7%
*-lft-identity71.7%
*-commutative71.7%
Simplified71.7%
div-sub71.7%
associate-/l*80.7%
Applied egg-rr80.7%
Taylor expanded in a around inf 84.0%
if -1.35e20 < z < -2.45000000000000002e-188 or 9.9999999999999996e-236 < z < 4.9999999999999998e-82 or 9.50000000000000065e-25 < z < 6e10Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in y around 0 90.0%
neg-mul-190.0%
Simplified90.0%
if -2.45000000000000002e-188 < z < 9.9999999999999996e-236Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in a around 0 97.7%
associate-*r/97.7%
neg-mul-197.7%
neg-sub097.7%
sub-neg97.7%
+-commutative97.7%
associate--r+97.7%
neg-sub097.7%
remove-double-neg97.7%
*-commutative97.7%
Simplified97.7%
if 4.9999999999999998e-82 < z < 9.50000000000000065e-25Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in y around inf 79.1%
associate-/l*72.6%
Simplified72.6%
Taylor expanded in y around 0 79.1%
fma-neg79.1%
associate-*l/79.1%
*-commutative79.1%
fma-neg79.1%
*-commutative79.1%
Simplified79.1%
if 6e10 < z Initial program 73.1%
sub-neg73.1%
+-commutative73.1%
neg-sub073.1%
associate-+l-73.1%
sub0-neg73.1%
neg-mul-173.1%
sub-neg73.1%
+-commutative73.1%
neg-sub073.1%
associate-+l-73.1%
sub0-neg73.1%
neg-mul-173.1%
times-frac73.1%
metadata-eval73.1%
*-lft-identity73.1%
*-commutative73.1%
Simplified73.1%
div-sub73.1%
associate-/l*83.6%
Applied egg-rr83.6%
Taylor expanded in z around 0 95.9%
neg-mul-195.9%
unsub-neg95.9%
Simplified95.9%
Taylor expanded in y around inf 71.1%
Final simplification85.8%
(FPCore (x y z t a)
:precision binary64
(if (<= z -2.75e+20)
(/ y a)
(if (<= z 3.5e-72)
(/ x t)
(if (<= z 5.2e-24)
(/ y a)
(if (<= z 130000000000.0)
(/ x t)
(if (<= z 5.5e+32) (* z (- (/ y t))) (/ y a)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -2.75e+20) {
tmp = y / a;
} else if (z <= 3.5e-72) {
tmp = x / t;
} else if (z <= 5.2e-24) {
tmp = y / a;
} else if (z <= 130000000000.0) {
tmp = x / t;
} else if (z <= 5.5e+32) {
tmp = z * -(y / t);
} else {
tmp = y / a;
}
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.75d+20)) then
tmp = y / a
else if (z <= 3.5d-72) then
tmp = x / t
else if (z <= 5.2d-24) then
tmp = y / a
else if (z <= 130000000000.0d0) then
tmp = x / t
else if (z <= 5.5d+32) then
tmp = z * -(y / t)
else
tmp = y / a
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.75e+20) {
tmp = y / a;
} else if (z <= 3.5e-72) {
tmp = x / t;
} else if (z <= 5.2e-24) {
tmp = y / a;
} else if (z <= 130000000000.0) {
tmp = x / t;
} else if (z <= 5.5e+32) {
tmp = z * -(y / t);
} else {
tmp = y / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -2.75e+20: tmp = y / a elif z <= 3.5e-72: tmp = x / t elif z <= 5.2e-24: tmp = y / a elif z <= 130000000000.0: tmp = x / t elif z <= 5.5e+32: tmp = z * -(y / t) else: tmp = y / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -2.75e+20) tmp = Float64(y / a); elseif (z <= 3.5e-72) tmp = Float64(x / t); elseif (z <= 5.2e-24) tmp = Float64(y / a); elseif (z <= 130000000000.0) tmp = Float64(x / t); elseif (z <= 5.5e+32) tmp = Float64(z * Float64(-Float64(y / t))); else tmp = Float64(y / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -2.75e+20) tmp = y / a; elseif (z <= 3.5e-72) tmp = x / t; elseif (z <= 5.2e-24) tmp = y / a; elseif (z <= 130000000000.0) tmp = x / t; elseif (z <= 5.5e+32) tmp = z * -(y / t); else tmp = y / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -2.75e+20], N[(y / a), $MachinePrecision], If[LessEqual[z, 3.5e-72], N[(x / t), $MachinePrecision], If[LessEqual[z, 5.2e-24], N[(y / a), $MachinePrecision], If[LessEqual[z, 130000000000.0], N[(x / t), $MachinePrecision], If[LessEqual[z, 5.5e+32], N[(z * (-N[(y / t), $MachinePrecision])), $MachinePrecision], N[(y / a), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.75 \cdot 10^{+20}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{-72}:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{elif}\;z \leq 5.2 \cdot 10^{-24}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 130000000000:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{elif}\;z \leq 5.5 \cdot 10^{+32}:\\
\;\;\;\;z \cdot \left(-\frac{y}{t}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a}\\
\end{array}
\end{array}
if z < -2.75e20 or 3.5e-72 < z < 5.2e-24 or 5.49999999999999984e32 < z Initial program 73.7%
sub-neg73.7%
+-commutative73.7%
neg-sub073.7%
associate-+l-73.7%
sub0-neg73.7%
neg-mul-173.7%
sub-neg73.7%
+-commutative73.7%
neg-sub073.7%
associate-+l-73.7%
sub0-neg73.7%
neg-mul-173.7%
times-frac73.7%
metadata-eval73.7%
*-lft-identity73.7%
*-commutative73.7%
Simplified73.7%
Taylor expanded in z around inf 59.3%
if -2.75e20 < z < 3.5e-72 or 5.2e-24 < z < 1.3e11Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in z around 0 71.8%
if 1.3e11 < z < 5.49999999999999984e32Initial program 99.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
times-frac99.4%
metadata-eval99.4%
*-lft-identity99.4%
*-commutative99.4%
Simplified99.4%
Taylor expanded in y around inf 99.4%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in y around 0 99.4%
fma-neg99.4%
associate-*l/99.4%
*-commutative99.4%
fma-neg99.4%
*-commutative99.4%
Simplified99.4%
Taylor expanded in z around 0 79.4%
associate-*r/79.4%
neg-mul-179.4%
Simplified79.4%
Final simplification66.2%
(FPCore (x y z t a)
:precision binary64
(if (<= z -1.7e+20)
(/ y a)
(if (<= z 3.6e-72)
(/ x t)
(if (<= z 1.05e-24)
(/ y a)
(if (<= z 4400000000.0)
(/ x t)
(if (<= z 6.8e+29) (/ y (/ (- t) z)) (/ y a)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.7e+20) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else if (z <= 1.05e-24) {
tmp = y / a;
} else if (z <= 4400000000.0) {
tmp = x / t;
} else if (z <= 6.8e+29) {
tmp = y / (-t / z);
} else {
tmp = y / a;
}
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 <= (-1.7d+20)) then
tmp = y / a
else if (z <= 3.6d-72) then
tmp = x / t
else if (z <= 1.05d-24) then
tmp = y / a
else if (z <= 4400000000.0d0) then
tmp = x / t
else if (z <= 6.8d+29) then
tmp = y / (-t / z)
else
tmp = y / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.7e+20) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else if (z <= 1.05e-24) {
tmp = y / a;
} else if (z <= 4400000000.0) {
tmp = x / t;
} else if (z <= 6.8e+29) {
tmp = y / (-t / z);
} else {
tmp = y / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.7e+20: tmp = y / a elif z <= 3.6e-72: tmp = x / t elif z <= 1.05e-24: tmp = y / a elif z <= 4400000000.0: tmp = x / t elif z <= 6.8e+29: tmp = y / (-t / z) else: tmp = y / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.7e+20) tmp = Float64(y / a); elseif (z <= 3.6e-72) tmp = Float64(x / t); elseif (z <= 1.05e-24) tmp = Float64(y / a); elseif (z <= 4400000000.0) tmp = Float64(x / t); elseif (z <= 6.8e+29) tmp = Float64(y / Float64(Float64(-t) / z)); else tmp = Float64(y / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.7e+20) tmp = y / a; elseif (z <= 3.6e-72) tmp = x / t; elseif (z <= 1.05e-24) tmp = y / a; elseif (z <= 4400000000.0) tmp = x / t; elseif (z <= 6.8e+29) tmp = y / (-t / z); else tmp = y / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.7e+20], N[(y / a), $MachinePrecision], If[LessEqual[z, 3.6e-72], N[(x / t), $MachinePrecision], If[LessEqual[z, 1.05e-24], N[(y / a), $MachinePrecision], If[LessEqual[z, 4400000000.0], N[(x / t), $MachinePrecision], If[LessEqual[z, 6.8e+29], N[(y / N[((-t) / z), $MachinePrecision]), $MachinePrecision], N[(y / a), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.7 \cdot 10^{+20}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-72}:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{elif}\;z \leq 1.05 \cdot 10^{-24}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 4400000000:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+29}:\\
\;\;\;\;\frac{y}{\frac{-t}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a}\\
\end{array}
\end{array}
if z < -1.7e20 or 3.6e-72 < z < 1.05e-24 or 6.79999999999999963e29 < z Initial program 73.7%
sub-neg73.7%
+-commutative73.7%
neg-sub073.7%
associate-+l-73.7%
sub0-neg73.7%
neg-mul-173.7%
sub-neg73.7%
+-commutative73.7%
neg-sub073.7%
associate-+l-73.7%
sub0-neg73.7%
neg-mul-173.7%
times-frac73.7%
metadata-eval73.7%
*-lft-identity73.7%
*-commutative73.7%
Simplified73.7%
Taylor expanded in z around inf 59.3%
if -1.7e20 < z < 3.6e-72 or 1.05e-24 < z < 4.4e9Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in z around 0 71.8%
if 4.4e9 < z < 6.79999999999999963e29Initial program 99.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
times-frac99.4%
metadata-eval99.4%
*-lft-identity99.4%
*-commutative99.4%
Simplified99.4%
Taylor expanded in y around inf 99.4%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in a around 0 80.0%
neg-mul-180.0%
distribute-neg-frac80.0%
Simplified80.0%
Final simplification66.2%
(FPCore (x y z t a)
:precision binary64
(if (<= z -1.1e+17)
(/ y a)
(if (<= z 3.6e-72)
(/ x t)
(if (<= z 8.1e-50)
(/ y a)
(if (<= z 820000000000.0)
(/ (- x) (* a z))
(if (<= z 8e+36) (/ y (/ (- t) z)) (/ y a)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.1e+17) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else if (z <= 8.1e-50) {
tmp = y / a;
} else if (z <= 820000000000.0) {
tmp = -x / (a * z);
} else if (z <= 8e+36) {
tmp = y / (-t / z);
} else {
tmp = y / a;
}
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 <= (-1.1d+17)) then
tmp = y / a
else if (z <= 3.6d-72) then
tmp = x / t
else if (z <= 8.1d-50) then
tmp = y / a
else if (z <= 820000000000.0d0) then
tmp = -x / (a * z)
else if (z <= 8d+36) then
tmp = y / (-t / z)
else
tmp = y / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.1e+17) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else if (z <= 8.1e-50) {
tmp = y / a;
} else if (z <= 820000000000.0) {
tmp = -x / (a * z);
} else if (z <= 8e+36) {
tmp = y / (-t / z);
} else {
tmp = y / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.1e+17: tmp = y / a elif z <= 3.6e-72: tmp = x / t elif z <= 8.1e-50: tmp = y / a elif z <= 820000000000.0: tmp = -x / (a * z) elif z <= 8e+36: tmp = y / (-t / z) else: tmp = y / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.1e+17) tmp = Float64(y / a); elseif (z <= 3.6e-72) tmp = Float64(x / t); elseif (z <= 8.1e-50) tmp = Float64(y / a); elseif (z <= 820000000000.0) tmp = Float64(Float64(-x) / Float64(a * z)); elseif (z <= 8e+36) tmp = Float64(y / Float64(Float64(-t) / z)); else tmp = Float64(y / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.1e+17) tmp = y / a; elseif (z <= 3.6e-72) tmp = x / t; elseif (z <= 8.1e-50) tmp = y / a; elseif (z <= 820000000000.0) tmp = -x / (a * z); elseif (z <= 8e+36) tmp = y / (-t / z); else tmp = y / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.1e+17], N[(y / a), $MachinePrecision], If[LessEqual[z, 3.6e-72], N[(x / t), $MachinePrecision], If[LessEqual[z, 8.1e-50], N[(y / a), $MachinePrecision], If[LessEqual[z, 820000000000.0], N[((-x) / N[(a * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 8e+36], N[(y / N[((-t) / z), $MachinePrecision]), $MachinePrecision], N[(y / a), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.1 \cdot 10^{+17}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-72}:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{elif}\;z \leq 8.1 \cdot 10^{-50}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 820000000000:\\
\;\;\;\;\frac{-x}{a \cdot z}\\
\mathbf{elif}\;z \leq 8 \cdot 10^{+36}:\\
\;\;\;\;\frac{y}{\frac{-t}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a}\\
\end{array}
\end{array}
if z < -1.1e17 or 3.6e-72 < z < 8.0999999999999999e-50 or 8.00000000000000034e36 < z Initial program 73.3%
sub-neg73.3%
+-commutative73.3%
neg-sub073.3%
associate-+l-73.3%
sub0-neg73.3%
neg-mul-173.3%
sub-neg73.3%
+-commutative73.3%
neg-sub073.3%
associate-+l-73.3%
sub0-neg73.3%
neg-mul-173.3%
times-frac73.3%
metadata-eval73.3%
*-lft-identity73.3%
*-commutative73.3%
Simplified73.3%
Taylor expanded in z around inf 60.3%
if -1.1e17 < z < 3.6e-72Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in z around 0 72.2%
if 8.0999999999999999e-50 < z < 8.2e11Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
*-commutative100.0%
Simplified100.0%
Taylor expanded in y around 0 85.9%
neg-mul-185.9%
Simplified85.9%
Taylor expanded in z around inf 56.4%
associate-*r/56.4%
mul-1-neg56.4%
*-commutative56.4%
Simplified56.4%
if 8.2e11 < z < 8.00000000000000034e36Initial program 99.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
sub-neg99.4%
+-commutative99.4%
neg-sub099.4%
associate-+l-99.4%
sub0-neg99.4%
neg-mul-199.4%
times-frac99.4%
metadata-eval99.4%
*-lft-identity99.4%
*-commutative99.4%
Simplified99.4%
Taylor expanded in y around inf 99.4%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in a around 0 80.0%
neg-mul-180.0%
distribute-neg-frac80.0%
Simplified80.0%
Final simplification66.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- x) (- (* a z) t))))
(if (<= z -8.5e+18)
(/ (- y (/ x z)) a)
(if (<= z -8.5e-193)
t_1
(if (<= z 6.8e-236)
(/ (- x (* y z)) t)
(if (<= z 4.1e-82) t_1 (/ y (- a (/ t z)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = -x / ((a * z) - t);
double tmp;
if (z <= -8.5e+18) {
tmp = (y - (x / z)) / a;
} else if (z <= -8.5e-193) {
tmp = t_1;
} else if (z <= 6.8e-236) {
tmp = (x - (y * z)) / t;
} else if (z <= 4.1e-82) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = -x / ((a * z) - t)
if (z <= (-8.5d+18)) then
tmp = (y - (x / z)) / a
else if (z <= (-8.5d-193)) then
tmp = t_1
else if (z <= 6.8d-236) then
tmp = (x - (y * z)) / t
else if (z <= 4.1d-82) then
tmp = t_1
else
tmp = 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 t_1 = -x / ((a * z) - t);
double tmp;
if (z <= -8.5e+18) {
tmp = (y - (x / z)) / a;
} else if (z <= -8.5e-193) {
tmp = t_1;
} else if (z <= 6.8e-236) {
tmp = (x - (y * z)) / t;
} else if (z <= 4.1e-82) {
tmp = t_1;
} else {
tmp = y / (a - (t / z));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = -x / ((a * z) - t) tmp = 0 if z <= -8.5e+18: tmp = (y - (x / z)) / a elif z <= -8.5e-193: tmp = t_1 elif z <= 6.8e-236: tmp = (x - (y * z)) / t elif z <= 4.1e-82: tmp = t_1 else: tmp = y / (a - (t / z)) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(-x) / Float64(Float64(a * z) - t)) tmp = 0.0 if (z <= -8.5e+18) tmp = Float64(Float64(y - Float64(x / z)) / a); elseif (z <= -8.5e-193) tmp = t_1; elseif (z <= 6.8e-236) tmp = Float64(Float64(x - Float64(y * z)) / t); elseif (z <= 4.1e-82) tmp = t_1; else tmp = Float64(y / Float64(a - Float64(t / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = -x / ((a * z) - t); tmp = 0.0; if (z <= -8.5e+18) tmp = (y - (x / z)) / a; elseif (z <= -8.5e-193) tmp = t_1; elseif (z <= 6.8e-236) tmp = (x - (y * z)) / t; elseif (z <= 4.1e-82) tmp = t_1; else tmp = y / (a - (t / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[((-x) / N[(N[(a * z), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.5e+18], N[(N[(y - N[(x / z), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[z, -8.5e-193], t$95$1, If[LessEqual[z, 6.8e-236], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[z, 4.1e-82], t$95$1, N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-x}{a \cdot z - t}\\
\mathbf{if}\;z \leq -8.5 \cdot 10^{+18}:\\
\;\;\;\;\frac{y - \frac{x}{z}}{a}\\
\mathbf{elif}\;z \leq -8.5 \cdot 10^{-193}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{-236}:\\
\;\;\;\;\frac{x - y \cdot z}{t}\\
\mathbf{elif}\;z \leq 4.1 \cdot 10^{-82}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\end{array}
\end{array}
if z < -8.5e18Initial program 71.7%
sub-neg71.7%
+-commutative71.7%
neg-sub071.7%
associate-+l-71.7%
sub0-neg71.7%
neg-mul-171.7%
sub-neg71.7%
+-commutative71.7%
neg-sub071.7%
associate-+l-71.7%
sub0-neg71.7%
neg-mul-171.7%
times-frac71.7%
metadata-eval71.7%
*-lft-identity71.7%
*-commutative71.7%
Simplified71.7%
div-sub71.7%
associate-/l*80.7%
Applied egg-rr80.7%
Taylor expanded in a around inf 84.0%
if -8.5e18 < z < -8.50000000000000004e-193 or 6.79999999999999961e-236 < z < 4.09999999999999996e-82Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in y around 0 89.3%
neg-mul-189.3%
Simplified89.3%
if -8.50000000000000004e-193 < z < 6.79999999999999961e-236Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in a around 0 97.7%
associate-*r/97.7%
neg-mul-197.7%
neg-sub097.7%
sub-neg97.7%
+-commutative97.7%
associate--r+97.7%
neg-sub097.7%
remove-double-neg97.7%
*-commutative97.7%
Simplified97.7%
if 4.09999999999999996e-82 < z Initial program 80.8%
sub-neg80.8%
+-commutative80.8%
neg-sub080.8%
associate-+l-80.8%
sub0-neg80.8%
neg-mul-180.8%
sub-neg80.8%
+-commutative80.8%
neg-sub080.8%
associate-+l-80.8%
sub0-neg80.8%
neg-mul-180.8%
times-frac80.8%
metadata-eval80.8%
*-lft-identity80.8%
*-commutative80.8%
Simplified80.8%
div-sub80.8%
associate-/l*86.9%
Applied egg-rr86.9%
Taylor expanded in z around 0 95.7%
neg-mul-195.7%
unsub-neg95.7%
Simplified95.7%
Taylor expanded in y around inf 67.8%
Final simplification83.9%
(FPCore (x y z t a)
:precision binary64
(if (<= z -2.3e+113)
(/ (- y (/ x z)) a)
(if (<= z 3.9e+133)
(/ (- x (* y z)) (- t (* a z)))
(+ (/ y (- a (/ t z))) (/ x t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -2.3e+113) {
tmp = (y - (x / z)) / a;
} else if (z <= 3.9e+133) {
tmp = (x - (y * z)) / (t - (a * z));
} else {
tmp = (y / (a - (t / z))) + (x / 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 (z <= (-2.3d+113)) then
tmp = (y - (x / z)) / a
else if (z <= 3.9d+133) then
tmp = (x - (y * z)) / (t - (a * z))
else
tmp = (y / (a - (t / z))) + (x / t)
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.3e+113) {
tmp = (y - (x / z)) / a;
} else if (z <= 3.9e+133) {
tmp = (x - (y * z)) / (t - (a * z));
} else {
tmp = (y / (a - (t / z))) + (x / t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -2.3e+113: tmp = (y - (x / z)) / a elif z <= 3.9e+133: tmp = (x - (y * z)) / (t - (a * z)) else: tmp = (y / (a - (t / z))) + (x / t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -2.3e+113) tmp = Float64(Float64(y - Float64(x / z)) / a); elseif (z <= 3.9e+133) tmp = Float64(Float64(x - Float64(y * z)) / Float64(t - Float64(a * z))); else tmp = Float64(Float64(y / Float64(a - Float64(t / z))) + Float64(x / t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -2.3e+113) tmp = (y - (x / z)) / a; elseif (z <= 3.9e+133) tmp = (x - (y * z)) / (t - (a * z)); else tmp = (y / (a - (t / z))) + (x / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -2.3e+113], N[(N[(y - N[(x / z), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[z, 3.9e+133], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / N[(t - N[(a * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.3 \cdot 10^{+113}:\\
\;\;\;\;\frac{y - \frac{x}{z}}{a}\\
\mathbf{elif}\;z \leq 3.9 \cdot 10^{+133}:\\
\;\;\;\;\frac{x - y \cdot z}{t - a \cdot z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a - \frac{t}{z}} + \frac{x}{t}\\
\end{array}
\end{array}
if z < -2.29999999999999997e113Initial program 60.6%
sub-neg60.6%
+-commutative60.6%
neg-sub060.6%
associate-+l-60.6%
sub0-neg60.6%
neg-mul-160.6%
sub-neg60.6%
+-commutative60.6%
neg-sub060.6%
associate-+l-60.6%
sub0-neg60.6%
neg-mul-160.6%
times-frac60.6%
metadata-eval60.6%
*-lft-identity60.6%
*-commutative60.6%
Simplified60.6%
div-sub60.6%
associate-/l*73.0%
Applied egg-rr73.0%
Taylor expanded in a around inf 91.7%
if -2.29999999999999997e113 < z < 3.90000000000000014e133Initial program 97.3%
if 3.90000000000000014e133 < z Initial program 61.8%
sub-neg61.8%
+-commutative61.8%
neg-sub061.8%
associate-+l-61.8%
sub0-neg61.8%
neg-mul-161.8%
sub-neg61.8%
+-commutative61.8%
neg-sub061.8%
associate-+l-61.8%
sub0-neg61.8%
neg-mul-161.8%
times-frac61.8%
metadata-eval61.8%
*-lft-identity61.8%
*-commutative61.8%
Simplified61.8%
div-sub61.8%
associate-/l*75.9%
Applied egg-rr75.9%
Taylor expanded in z around 0 96.6%
neg-mul-196.6%
unsub-neg96.6%
Simplified96.6%
Taylor expanded in z around 0 82.5%
associate-*r/82.5%
mul-1-neg82.5%
Simplified82.5%
Final simplification94.8%
(FPCore (x y z t a) :precision binary64 (- (/ y (- a (/ t z))) (/ x (- (* a z) t))))
double code(double x, double y, double z, double t, double a) {
return (y / (a - (t / z))) - (x / ((a * z) - 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 = (y / (a - (t / z))) - (x / ((a * z) - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (y / (a - (t / z))) - (x / ((a * z) - t));
}
def code(x, y, z, t, a): return (y / (a - (t / z))) - (x / ((a * z) - t))
function code(x, y, z, t, a) return Float64(Float64(y / Float64(a - Float64(t / z))) - Float64(x / Float64(Float64(a * z) - t))) end
function tmp = code(x, y, z, t, a) tmp = (y / (a - (t / z))) - (x / ((a * z) - t)); end
code[x_, y_, z_, t_, a_] := N[(N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x / N[(N[(a * z), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{y}{a - \frac{t}{z}} - \frac{x}{a \cdot z - t}
\end{array}
Initial program 87.8%
sub-neg87.8%
+-commutative87.8%
neg-sub087.8%
associate-+l-87.8%
sub0-neg87.8%
neg-mul-187.8%
sub-neg87.8%
+-commutative87.8%
neg-sub087.8%
associate-+l-87.8%
sub0-neg87.8%
neg-mul-187.8%
times-frac87.8%
metadata-eval87.8%
*-lft-identity87.8%
*-commutative87.8%
Simplified87.8%
div-sub87.0%
associate-/l*89.0%
Applied egg-rr89.0%
Taylor expanded in z around 0 94.7%
neg-mul-194.7%
unsub-neg94.7%
Simplified94.7%
Final simplification94.7%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -62000000000000.0) (not (<= z 1.65e-95))) (/ y (- a (/ t z))) (/ x t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -62000000000000.0) || !(z <= 1.65e-95)) {
tmp = y / (a - (t / z));
} else {
tmp = x / 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 ((z <= (-62000000000000.0d0)) .or. (.not. (z <= 1.65d-95))) then
tmp = y / (a - (t / z))
else
tmp = x / t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -62000000000000.0) || !(z <= 1.65e-95)) {
tmp = y / (a - (t / z));
} else {
tmp = x / t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -62000000000000.0) or not (z <= 1.65e-95): tmp = y / (a - (t / z)) else: tmp = x / t return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -62000000000000.0) || !(z <= 1.65e-95)) tmp = Float64(y / Float64(a - Float64(t / z))); else tmp = Float64(x / t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -62000000000000.0) || ~((z <= 1.65e-95))) tmp = y / (a - (t / z)); else tmp = x / t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -62000000000000.0], N[Not[LessEqual[z, 1.65e-95]], $MachinePrecision]], N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -62000000000000 \lor \neg \left(z \leq 1.65 \cdot 10^{-95}\right):\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t}\\
\end{array}
\end{array}
if z < -6.2e13 or 1.65e-95 < z Initial program 76.8%
sub-neg76.8%
+-commutative76.8%
neg-sub076.8%
associate-+l-76.8%
sub0-neg76.8%
neg-mul-176.8%
sub-neg76.8%
+-commutative76.8%
neg-sub076.8%
associate-+l-76.8%
sub0-neg76.8%
neg-mul-176.8%
times-frac76.8%
metadata-eval76.8%
*-lft-identity76.8%
*-commutative76.8%
Simplified76.8%
div-sub76.8%
associate-/l*84.2%
Applied egg-rr84.2%
Taylor expanded in z around 0 95.0%
neg-mul-195.0%
unsub-neg95.0%
Simplified95.0%
Taylor expanded in y around inf 67.7%
if -6.2e13 < z < 1.65e-95Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in z around 0 74.8%
Final simplification71.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -1.3e+22) (not (<= z 3.6e-72))) (/ y (- a (/ t z))) (/ (- x (* y z)) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -1.3e+22) || !(z <= 3.6e-72)) {
tmp = y / (a - (t / 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 ((z <= (-1.3d+22)) .or. (.not. (z <= 3.6d-72))) then
tmp = y / (a - (t / 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 ((z <= -1.3e+22) || !(z <= 3.6e-72)) {
tmp = y / (a - (t / z));
} else {
tmp = (x - (y * z)) / t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -1.3e+22) or not (z <= 3.6e-72): tmp = y / (a - (t / z)) else: tmp = (x - (y * z)) / t return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -1.3e+22) || !(z <= 3.6e-72)) tmp = Float64(y / Float64(a - Float64(t / z))); else tmp = Float64(Float64(x - Float64(y * z)) / t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -1.3e+22) || ~((z <= 3.6e-72))) tmp = y / (a - (t / z)); else tmp = (x - (y * z)) / t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -1.3e+22], N[Not[LessEqual[z, 3.6e-72]], $MachinePrecision]], N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.3 \cdot 10^{+22} \lor \neg \left(z \leq 3.6 \cdot 10^{-72}\right):\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x - y \cdot z}{t}\\
\end{array}
\end{array}
if z < -1.3e22 or 3.6e-72 < z Initial program 75.4%
sub-neg75.4%
+-commutative75.4%
neg-sub075.4%
associate-+l-75.4%
sub0-neg75.4%
neg-mul-175.4%
sub-neg75.4%
+-commutative75.4%
neg-sub075.4%
associate-+l-75.4%
sub0-neg75.4%
neg-mul-175.4%
times-frac75.4%
metadata-eval75.4%
*-lft-identity75.4%
*-commutative75.4%
Simplified75.4%
div-sub75.4%
associate-/l*83.9%
Applied egg-rr83.9%
Taylor expanded in z around 0 95.4%
neg-mul-195.4%
unsub-neg95.4%
Simplified95.4%
Taylor expanded in y around inf 68.6%
if -1.3e22 < z < 3.6e-72Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in a around 0 82.4%
associate-*r/82.4%
neg-mul-182.4%
neg-sub082.4%
sub-neg82.4%
+-commutative82.4%
associate--r+82.4%
neg-sub082.4%
remove-double-neg82.4%
*-commutative82.4%
Simplified82.4%
Final simplification75.7%
(FPCore (x y z t a) :precision binary64 (if (<= z -1.65e+34) (/ (- y (/ x z)) a) (if (<= z 3.6e-72) (/ (- x (* y z)) t) (/ y (- a (/ t z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.65e+34) {
tmp = (y - (x / z)) / a;
} else if (z <= 3.6e-72) {
tmp = (x - (y * z)) / t;
} else {
tmp = 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 (z <= (-1.65d+34)) then
tmp = (y - (x / z)) / a
else if (z <= 3.6d-72) then
tmp = (x - (y * z)) / t
else
tmp = 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 (z <= -1.65e+34) {
tmp = (y - (x / z)) / a;
} else if (z <= 3.6e-72) {
tmp = (x - (y * z)) / t;
} else {
tmp = y / (a - (t / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.65e+34: tmp = (y - (x / z)) / a elif z <= 3.6e-72: tmp = (x - (y * z)) / t else: tmp = y / (a - (t / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.65e+34) tmp = Float64(Float64(y - Float64(x / z)) / a); elseif (z <= 3.6e-72) tmp = Float64(Float64(x - Float64(y * z)) / t); else tmp = Float64(y / Float64(a - Float64(t / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.65e+34) tmp = (y - (x / z)) / a; elseif (z <= 3.6e-72) tmp = (x - (y * z)) / t; else tmp = y / (a - (t / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.65e+34], N[(N[(y - N[(x / z), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[z, 3.6e-72], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision], N[(y / N[(a - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.65 \cdot 10^{+34}:\\
\;\;\;\;\frac{y - \frac{x}{z}}{a}\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-72}:\\
\;\;\;\;\frac{x - y \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a - \frac{t}{z}}\\
\end{array}
\end{array}
if z < -1.64999999999999994e34Initial program 69.9%
sub-neg69.9%
+-commutative69.9%
neg-sub069.9%
associate-+l-69.9%
sub0-neg69.9%
neg-mul-169.9%
sub-neg69.9%
+-commutative69.9%
neg-sub069.9%
associate-+l-69.9%
sub0-neg69.9%
neg-mul-169.9%
times-frac69.9%
metadata-eval69.9%
*-lft-identity69.9%
*-commutative69.9%
Simplified69.9%
div-sub69.9%
associate-/l*79.5%
Applied egg-rr79.5%
Taylor expanded in a around inf 86.0%
if -1.64999999999999994e34 < z < 3.6e-72Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in a around 0 82.0%
associate-*r/82.0%
neg-mul-182.0%
neg-sub082.0%
sub-neg82.0%
+-commutative82.0%
associate--r+82.0%
neg-sub082.0%
remove-double-neg82.0%
*-commutative82.0%
Simplified82.0%
if 3.6e-72 < z Initial program 79.9%
sub-neg79.9%
+-commutative79.9%
neg-sub079.9%
associate-+l-79.9%
sub0-neg79.9%
neg-mul-179.9%
sub-neg79.9%
+-commutative79.9%
neg-sub079.9%
associate-+l-79.9%
sub0-neg79.9%
neg-mul-179.9%
times-frac79.9%
metadata-eval79.9%
*-lft-identity79.9%
*-commutative79.9%
Simplified79.9%
div-sub79.9%
associate-/l*87.8%
Applied egg-rr87.8%
Taylor expanded in z around 0 97.0%
neg-mul-197.0%
unsub-neg97.0%
Simplified97.0%
Taylor expanded in y around inf 69.3%
Final simplification79.8%
(FPCore (x y z t a) :precision binary64 (if (<= z -1.9e+19) (/ y a) (if (<= z 3.6e-72) (/ x t) (/ y a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.9e+19) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else {
tmp = y / a;
}
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 <= (-1.9d+19)) then
tmp = y / a
else if (z <= 3.6d-72) then
tmp = x / t
else
tmp = y / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.9e+19) {
tmp = y / a;
} else if (z <= 3.6e-72) {
tmp = x / t;
} else {
tmp = y / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.9e+19: tmp = y / a elif z <= 3.6e-72: tmp = x / t else: tmp = y / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.9e+19) tmp = Float64(y / a); elseif (z <= 3.6e-72) tmp = Float64(x / t); else tmp = Float64(y / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.9e+19) tmp = y / a; elseif (z <= 3.6e-72) tmp = x / t; else tmp = y / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.9e+19], N[(y / a), $MachinePrecision], If[LessEqual[z, 3.6e-72], N[(x / t), $MachinePrecision], N[(y / a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.9 \cdot 10^{+19}:\\
\;\;\;\;\frac{y}{a}\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-72}:\\
\;\;\;\;\frac{x}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{a}\\
\end{array}
\end{array}
if z < -1.9e19 or 3.6e-72 < z Initial program 75.7%
sub-neg75.7%
+-commutative75.7%
neg-sub075.7%
associate-+l-75.7%
sub0-neg75.7%
neg-mul-175.7%
sub-neg75.7%
+-commutative75.7%
neg-sub075.7%
associate-+l-75.7%
sub0-neg75.7%
neg-mul-175.7%
times-frac75.7%
metadata-eval75.7%
*-lft-identity75.7%
*-commutative75.7%
Simplified75.7%
Taylor expanded in z around inf 56.4%
if -1.9e19 < z < 3.6e-72Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in z around 0 72.2%
Final simplification64.3%
(FPCore (x y z t a) :precision binary64 (/ x t))
double code(double x, double y, double z, double t, double a) {
return x / 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 / t
end function
public static double code(double x, double y, double z, double t, double a) {
return x / t;
}
def code(x, y, z, t, a): return x / t
function code(x, y, z, t, a) return Float64(x / t) end
function tmp = code(x, y, z, t, a) tmp = x / t; end
code[x_, y_, z_, t_, a_] := N[(x / t), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{t}
\end{array}
Initial program 87.8%
sub-neg87.8%
+-commutative87.8%
neg-sub087.8%
associate-+l-87.8%
sub0-neg87.8%
neg-mul-187.8%
sub-neg87.8%
+-commutative87.8%
neg-sub087.8%
associate-+l-87.8%
sub0-neg87.8%
neg-mul-187.8%
times-frac87.8%
metadata-eval87.8%
*-lft-identity87.8%
*-commutative87.8%
Simplified87.8%
Taylor expanded in z around 0 44.7%
Final simplification44.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- t (* a z))) (t_2 (- (/ x t_1) (/ y (- (/ t z) a)))))
(if (< z -32113435955957344.0)
t_2
(if (< z 3.5139522372978296e-86) (* (- x (* y z)) (/ 1.0 t_1)) t_2))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t - (a * z);
double t_2 = (x / t_1) - (y / ((t / z) - a));
double tmp;
if (z < -32113435955957344.0) {
tmp = t_2;
} else if (z < 3.5139522372978296e-86) {
tmp = (x - (y * z)) * (1.0 / 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 = t - (a * z)
t_2 = (x / t_1) - (y / ((t / z) - a))
if (z < (-32113435955957344.0d0)) then
tmp = t_2
else if (z < 3.5139522372978296d-86) then
tmp = (x - (y * z)) * (1.0d0 / 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 = t - (a * z);
double t_2 = (x / t_1) - (y / ((t / z) - a));
double tmp;
if (z < -32113435955957344.0) {
tmp = t_2;
} else if (z < 3.5139522372978296e-86) {
tmp = (x - (y * z)) * (1.0 / t_1);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t - (a * z) t_2 = (x / t_1) - (y / ((t / z) - a)) tmp = 0 if z < -32113435955957344.0: tmp = t_2 elif z < 3.5139522372978296e-86: tmp = (x - (y * z)) * (1.0 / t_1) else: tmp = t_2 return tmp
function code(x, y, z, t, a) t_1 = Float64(t - Float64(a * z)) t_2 = Float64(Float64(x / t_1) - Float64(y / Float64(Float64(t / z) - a))) tmp = 0.0 if (z < -32113435955957344.0) tmp = t_2; elseif (z < 3.5139522372978296e-86) tmp = Float64(Float64(x - Float64(y * z)) * Float64(1.0 / t_1)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t - (a * z); t_2 = (x / t_1) - (y / ((t / z) - a)); tmp = 0.0; if (z < -32113435955957344.0) tmp = t_2; elseif (z < 3.5139522372978296e-86) tmp = (x - (y * z)) * (1.0 / t_1); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t - N[(a * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x / t$95$1), $MachinePrecision] - N[(y / N[(N[(t / z), $MachinePrecision] - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[z, -32113435955957344.0], t$95$2, If[Less[z, 3.5139522372978296e-86], N[(N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision] * N[(1.0 / t$95$1), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t - a \cdot z\\
t_2 := \frac{x}{t_1} - \frac{y}{\frac{t}{z} - a}\\
\mathbf{if}\;z < -32113435955957344:\\
\;\;\;\;t_2\\
\mathbf{elif}\;z < 3.5139522372978296 \cdot 10^{-86}:\\
\;\;\;\;\left(x - y \cdot z\right) \cdot \frac{1}{t_1}\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
herbie shell --seed 2023185
(FPCore (x y z t a)
:name "Diagrams.Solve.Tridiagonal:solveTriDiagonal from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< z -32113435955957344.0) (- (/ x (- t (* a z))) (/ y (- (/ t z) a))) (if (< z 3.5139522372978296e-86) (* (- x (* y z)) (/ 1.0 (- t (* a z)))) (- (/ x (- t (* a z))) (/ y (- (/ t z) a)))))
(/ (- x (* y z)) (- t (* a z))))