
(FPCore (x y z t) :precision binary64 (/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1.0)))
double code(double x, double y, double z, double t) {
return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0d0)
end function
public static double code(double x, double y, double z, double t) {
return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0);
}
def code(x, y, z, t): return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0)
function code(x, y, z, t) return Float64(Float64(x + Float64(Float64(Float64(y * z) - x) / Float64(Float64(t * z) - x))) / Float64(x + 1.0)) end
function tmp = code(x, y, z, t) tmp = (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0); end
code[x_, y_, z_, t_] := N[(N[(x + N[(N[(N[(y * z), $MachinePrecision] - x), $MachinePrecision] / N[(N[(t * z), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x + \frac{y \cdot z - x}{t \cdot z - x}}{x + 1}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1.0)))
double code(double x, double y, double z, double t) {
return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0d0)
end function
public static double code(double x, double y, double z, double t) {
return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0);
}
def code(x, y, z, t): return (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0)
function code(x, y, z, t) return Float64(Float64(x + Float64(Float64(Float64(y * z) - x) / Float64(Float64(t * z) - x))) / Float64(x + 1.0)) end
function tmp = code(x, y, z, t) tmp = (x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0); end
code[x_, y_, z_, t_] := N[(N[(x + N[(N[(N[(y * z), $MachinePrecision] - x), $MachinePrecision] / N[(N[(t * z), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x + \frac{y \cdot z - x}{t \cdot z - x}}{x + 1}
\end{array}
(FPCore (x y z t) :precision binary64 (let* ((t_1 (/ (+ x (/ (- (* y z) x) (- (* z t) x))) (+ x 1.0)))) (if (<= t_1 5e+260) t_1 (/ (+ x (/ y t)) (+ x 1.0)))))
double code(double x, double y, double z, double t) {
double t_1 = (x + (((y * z) - x) / ((z * t) - x))) / (x + 1.0);
double tmp;
if (t_1 <= 5e+260) {
tmp = t_1;
} else {
tmp = (x + (y / t)) / (x + 1.0);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x + (((y * z) - x) / ((z * t) - x))) / (x + 1.0d0)
if (t_1 <= 5d+260) then
tmp = t_1
else
tmp = (x + (y / t)) / (x + 1.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x + (((y * z) - x) / ((z * t) - x))) / (x + 1.0);
double tmp;
if (t_1 <= 5e+260) {
tmp = t_1;
} else {
tmp = (x + (y / t)) / (x + 1.0);
}
return tmp;
}
def code(x, y, z, t): t_1 = (x + (((y * z) - x) / ((z * t) - x))) / (x + 1.0) tmp = 0 if t_1 <= 5e+260: tmp = t_1 else: tmp = (x + (y / t)) / (x + 1.0) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x + Float64(Float64(Float64(y * z) - x) / Float64(Float64(z * t) - x))) / Float64(x + 1.0)) tmp = 0.0 if (t_1 <= 5e+260) tmp = t_1; else tmp = Float64(Float64(x + Float64(y / t)) / Float64(x + 1.0)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x + (((y * z) - x) / ((z * t) - x))) / (x + 1.0); tmp = 0.0; if (t_1 <= 5e+260) tmp = t_1; else tmp = (x + (y / t)) / (x + 1.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x + N[(N[(N[(y * z), $MachinePrecision] - x), $MachinePrecision] / N[(N[(z * t), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 5e+260], t$95$1, N[(N[(x + N[(y / t), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x + \frac{y \cdot z - x}{z \cdot t - x}}{x + 1}\\
\mathbf{if}\;t_1 \leq 5 \cdot 10^{+260}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1}\\
\end{array}
\end{array}
if (/.f64 (+.f64 x (/.f64 (-.f64 (*.f64 y z) x) (-.f64 (*.f64 t z) x))) (+.f64 x 1)) < 4.9999999999999996e260Initial program 96.9%
if 4.9999999999999996e260 < (/.f64 (+.f64 x (/.f64 (-.f64 (*.f64 y z) x) (-.f64 (*.f64 t z) x))) (+.f64 x 1)) Initial program 9.7%
*-commutative9.7%
Simplified9.7%
Taylor expanded in z around inf 81.6%
Final simplification95.3%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ x (* z t)))))
(if (<= x -2.6e-13)
(/ x (+ x 1.0))
(if (<= x -5.8e-43)
(/ y (* t (+ x 1.0)))
(if (<= x -3.1e-71)
t_1
(if (<= x -3e-93)
(/ (- (* y z)) x)
(if (<= x -3.7e-190)
t_1
(if (<= x 5e-102) (/ y t) (if (<= x 2.25e-35) t_1 1.0)))))))))
double code(double x, double y, double z, double t) {
double t_1 = x - (x / (z * t));
double tmp;
if (x <= -2.6e-13) {
tmp = x / (x + 1.0);
} else if (x <= -5.8e-43) {
tmp = y / (t * (x + 1.0));
} else if (x <= -3.1e-71) {
tmp = t_1;
} else if (x <= -3e-93) {
tmp = -(y * z) / x;
} else if (x <= -3.7e-190) {
tmp = t_1;
} else if (x <= 5e-102) {
tmp = y / t;
} else if (x <= 2.25e-35) {
tmp = t_1;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - (x / (z * t))
if (x <= (-2.6d-13)) then
tmp = x / (x + 1.0d0)
else if (x <= (-5.8d-43)) then
tmp = y / (t * (x + 1.0d0))
else if (x <= (-3.1d-71)) then
tmp = t_1
else if (x <= (-3d-93)) then
tmp = -(y * z) / x
else if (x <= (-3.7d-190)) then
tmp = t_1
else if (x <= 5d-102) then
tmp = y / t
else if (x <= 2.25d-35) then
tmp = t_1
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - (x / (z * t));
double tmp;
if (x <= -2.6e-13) {
tmp = x / (x + 1.0);
} else if (x <= -5.8e-43) {
tmp = y / (t * (x + 1.0));
} else if (x <= -3.1e-71) {
tmp = t_1;
} else if (x <= -3e-93) {
tmp = -(y * z) / x;
} else if (x <= -3.7e-190) {
tmp = t_1;
} else if (x <= 5e-102) {
tmp = y / t;
} else if (x <= 2.25e-35) {
tmp = t_1;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - (x / (z * t)) tmp = 0 if x <= -2.6e-13: tmp = x / (x + 1.0) elif x <= -5.8e-43: tmp = y / (t * (x + 1.0)) elif x <= -3.1e-71: tmp = t_1 elif x <= -3e-93: tmp = -(y * z) / x elif x <= -3.7e-190: tmp = t_1 elif x <= 5e-102: tmp = y / t elif x <= 2.25e-35: tmp = t_1 else: tmp = 1.0 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(x / Float64(z * t))) tmp = 0.0 if (x <= -2.6e-13) tmp = Float64(x / Float64(x + 1.0)); elseif (x <= -5.8e-43) tmp = Float64(y / Float64(t * Float64(x + 1.0))); elseif (x <= -3.1e-71) tmp = t_1; elseif (x <= -3e-93) tmp = Float64(Float64(-Float64(y * z)) / x); elseif (x <= -3.7e-190) tmp = t_1; elseif (x <= 5e-102) tmp = Float64(y / t); elseif (x <= 2.25e-35) tmp = t_1; else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - (x / (z * t)); tmp = 0.0; if (x <= -2.6e-13) tmp = x / (x + 1.0); elseif (x <= -5.8e-43) tmp = y / (t * (x + 1.0)); elseif (x <= -3.1e-71) tmp = t_1; elseif (x <= -3e-93) tmp = -(y * z) / x; elseif (x <= -3.7e-190) tmp = t_1; elseif (x <= 5e-102) tmp = y / t; elseif (x <= 2.25e-35) tmp = t_1; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(x / N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -2.6e-13], N[(x / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -5.8e-43], N[(y / N[(t * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.1e-71], t$95$1, If[LessEqual[x, -3e-93], N[((-N[(y * z), $MachinePrecision]) / x), $MachinePrecision], If[LessEqual[x, -3.7e-190], t$95$1, If[LessEqual[x, 5e-102], N[(y / t), $MachinePrecision], If[LessEqual[x, 2.25e-35], t$95$1, 1.0]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{x}{z \cdot t}\\
\mathbf{if}\;x \leq -2.6 \cdot 10^{-13}:\\
\;\;\;\;\frac{x}{x + 1}\\
\mathbf{elif}\;x \leq -5.8 \cdot 10^{-43}:\\
\;\;\;\;\frac{y}{t \cdot \left(x + 1\right)}\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-71}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq -3 \cdot 10^{-93}:\\
\;\;\;\;\frac{-y \cdot z}{x}\\
\mathbf{elif}\;x \leq -3.7 \cdot 10^{-190}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq 5 \cdot 10^{-102}:\\
\;\;\;\;\frac{y}{t}\\
\mathbf{elif}\;x \leq 2.25 \cdot 10^{-35}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -2.6e-13Initial program 85.4%
*-commutative85.4%
Simplified85.4%
Taylor expanded in t around inf 91.8%
+-commutative91.8%
Simplified91.8%
if -2.6e-13 < x < -5.8000000000000003e-43Initial program 68.0%
*-commutative68.0%
Simplified68.0%
Taylor expanded in y around inf 68.0%
times-frac99.7%
+-commutative99.7%
Simplified99.7%
Taylor expanded in t around inf 68.4%
if -5.8000000000000003e-43 < x < -3.10000000000000002e-71 or -3.0000000000000001e-93 < x < -3.7000000000000002e-190 or 5.00000000000000026e-102 < x < 2.25000000000000005e-35Initial program 87.2%
*-commutative87.2%
Simplified87.2%
Taylor expanded in y around 0 69.1%
+-commutative69.1%
Simplified69.1%
clear-num69.1%
inv-pow69.1%
Applied egg-rr69.1%
unpow-169.1%
div-sub69.1%
sub-neg69.1%
associate-/l*69.0%
*-inverses69.0%
metadata-eval69.0%
Simplified69.0%
Taylor expanded in x around 0 52.4%
*-commutative52.4%
sub-neg52.4%
distribute-neg-frac52.4%
metadata-eval52.4%
distribute-rgt-in52.4%
*-lft-identity52.4%
associate-*l/52.5%
associate-*r/52.5%
mul-1-neg52.5%
unsub-neg52.5%
Simplified52.5%
if -3.10000000000000002e-71 < x < -3.0000000000000001e-93Initial program 87.5%
*-commutative87.5%
Simplified87.5%
Taylor expanded in y around inf 63.2%
times-frac63.4%
+-commutative63.4%
Simplified63.4%
Taylor expanded in x around 0 63.4%
Taylor expanded in t around 0 63.6%
associate-*r/63.6%
mul-1-neg63.6%
distribute-rgt-neg-out63.6%
Simplified63.6%
if -3.7000000000000002e-190 < x < 5.00000000000000026e-102Initial program 93.0%
*-commutative93.0%
Simplified93.0%
Taylor expanded in y around inf 60.2%
times-frac48.9%
+-commutative48.9%
Simplified48.9%
Taylor expanded in x around 0 54.5%
if 2.25000000000000005e-35 < x Initial program 88.4%
*-commutative88.4%
Simplified88.4%
Taylor expanded in y around 0 92.6%
+-commutative92.6%
Simplified92.6%
Taylor expanded in x around inf 88.9%
Final simplification73.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (+ x (/ y t)) (+ x 1.0))))
(if (<= z -1.3e-44)
t_1
(if (<= z 1.75e-228)
(/ (+ 1.0 (- x (/ y (/ x z)))) (+ x 1.0))
(if (<= z 1.1e+46) (/ (- x (/ x (- (* z t) x))) (+ x 1.0)) t_1)))))
double code(double x, double y, double z, double t) {
double t_1 = (x + (y / t)) / (x + 1.0);
double tmp;
if (z <= -1.3e-44) {
tmp = t_1;
} else if (z <= 1.75e-228) {
tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0);
} else if (z <= 1.1e+46) {
tmp = (x - (x / ((z * t) - x))) / (x + 1.0);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x + (y / t)) / (x + 1.0d0)
if (z <= (-1.3d-44)) then
tmp = t_1
else if (z <= 1.75d-228) then
tmp = (1.0d0 + (x - (y / (x / z)))) / (x + 1.0d0)
else if (z <= 1.1d+46) then
tmp = (x - (x / ((z * t) - x))) / (x + 1.0d0)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x + (y / t)) / (x + 1.0);
double tmp;
if (z <= -1.3e-44) {
tmp = t_1;
} else if (z <= 1.75e-228) {
tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0);
} else if (z <= 1.1e+46) {
tmp = (x - (x / ((z * t) - x))) / (x + 1.0);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (x + (y / t)) / (x + 1.0) tmp = 0 if z <= -1.3e-44: tmp = t_1 elif z <= 1.75e-228: tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0) elif z <= 1.1e+46: tmp = (x - (x / ((z * t) - x))) / (x + 1.0) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x + Float64(y / t)) / Float64(x + 1.0)) tmp = 0.0 if (z <= -1.3e-44) tmp = t_1; elseif (z <= 1.75e-228) tmp = Float64(Float64(1.0 + Float64(x - Float64(y / Float64(x / z)))) / Float64(x + 1.0)); elseif (z <= 1.1e+46) tmp = Float64(Float64(x - Float64(x / Float64(Float64(z * t) - x))) / Float64(x + 1.0)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x + (y / t)) / (x + 1.0); tmp = 0.0; if (z <= -1.3e-44) tmp = t_1; elseif (z <= 1.75e-228) tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0); elseif (z <= 1.1e+46) tmp = (x - (x / ((z * t) - x))) / (x + 1.0); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x + N[(y / t), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.3e-44], t$95$1, If[LessEqual[z, 1.75e-228], N[(N[(1.0 + N[(x - N[(y / N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.1e+46], N[(N[(x - N[(x / N[(N[(z * t), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x + \frac{y}{t}}{x + 1}\\
\mathbf{if}\;z \leq -1.3 \cdot 10^{-44}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{-228}:\\
\;\;\;\;\frac{1 + \left(x - \frac{y}{\frac{x}{z}}\right)}{x + 1}\\
\mathbf{elif}\;z \leq 1.1 \cdot 10^{+46}:\\
\;\;\;\;\frac{x - \frac{x}{z \cdot t - x}}{x + 1}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -1.2999999999999999e-44 or 1.1e46 < z Initial program 77.3%
*-commutative77.3%
Simplified77.3%
Taylor expanded in z around inf 89.6%
if -1.2999999999999999e-44 < z < 1.74999999999999987e-228Initial program 99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in t around 0 91.3%
+-commutative91.3%
mul-1-neg91.3%
unsub-neg91.3%
associate-/l*91.3%
+-commutative91.3%
Simplified91.3%
if 1.74999999999999987e-228 < z < 1.1e46Initial program 98.0%
*-commutative98.0%
Simplified98.0%
Taylor expanded in y around 0 82.3%
+-commutative82.3%
Simplified82.3%
Final simplification88.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- x (/ x (* z t)))))
(if (<= x -9.5e-67)
(/ x (+ x 1.0))
(if (<= x -3.55e-97)
(/ (- (* y z)) x)
(if (<= x -2.2e-190)
t_1
(if (<= x 5.8e-102) (/ y t) (if (<= x 2.4e-34) t_1 1.0)))))))
double code(double x, double y, double z, double t) {
double t_1 = x - (x / (z * t));
double tmp;
if (x <= -9.5e-67) {
tmp = x / (x + 1.0);
} else if (x <= -3.55e-97) {
tmp = -(y * z) / x;
} else if (x <= -2.2e-190) {
tmp = t_1;
} else if (x <= 5.8e-102) {
tmp = y / t;
} else if (x <= 2.4e-34) {
tmp = t_1;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x - (x / (z * t))
if (x <= (-9.5d-67)) then
tmp = x / (x + 1.0d0)
else if (x <= (-3.55d-97)) then
tmp = -(y * z) / x
else if (x <= (-2.2d-190)) then
tmp = t_1
else if (x <= 5.8d-102) then
tmp = y / t
else if (x <= 2.4d-34) then
tmp = t_1
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x - (x / (z * t));
double tmp;
if (x <= -9.5e-67) {
tmp = x / (x + 1.0);
} else if (x <= -3.55e-97) {
tmp = -(y * z) / x;
} else if (x <= -2.2e-190) {
tmp = t_1;
} else if (x <= 5.8e-102) {
tmp = y / t;
} else if (x <= 2.4e-34) {
tmp = t_1;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - (x / (z * t)) tmp = 0 if x <= -9.5e-67: tmp = x / (x + 1.0) elif x <= -3.55e-97: tmp = -(y * z) / x elif x <= -2.2e-190: tmp = t_1 elif x <= 5.8e-102: tmp = y / t elif x <= 2.4e-34: tmp = t_1 else: tmp = 1.0 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(x / Float64(z * t))) tmp = 0.0 if (x <= -9.5e-67) tmp = Float64(x / Float64(x + 1.0)); elseif (x <= -3.55e-97) tmp = Float64(Float64(-Float64(y * z)) / x); elseif (x <= -2.2e-190) tmp = t_1; elseif (x <= 5.8e-102) tmp = Float64(y / t); elseif (x <= 2.4e-34) tmp = t_1; else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x - (x / (z * t)); tmp = 0.0; if (x <= -9.5e-67) tmp = x / (x + 1.0); elseif (x <= -3.55e-97) tmp = -(y * z) / x; elseif (x <= -2.2e-190) tmp = t_1; elseif (x <= 5.8e-102) tmp = y / t; elseif (x <= 2.4e-34) tmp = t_1; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(x / N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -9.5e-67], N[(x / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.55e-97], N[((-N[(y * z), $MachinePrecision]) / x), $MachinePrecision], If[LessEqual[x, -2.2e-190], t$95$1, If[LessEqual[x, 5.8e-102], N[(y / t), $MachinePrecision], If[LessEqual[x, 2.4e-34], t$95$1, 1.0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{x}{z \cdot t}\\
\mathbf{if}\;x \leq -9.5 \cdot 10^{-67}:\\
\;\;\;\;\frac{x}{x + 1}\\
\mathbf{elif}\;x \leq -3.55 \cdot 10^{-97}:\\
\;\;\;\;\frac{-y \cdot z}{x}\\
\mathbf{elif}\;x \leq -2.2 \cdot 10^{-190}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq 5.8 \cdot 10^{-102}:\\
\;\;\;\;\frac{y}{t}\\
\mathbf{elif}\;x \leq 2.4 \cdot 10^{-34}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -9.4999999999999994e-67Initial program 85.0%
*-commutative85.0%
Simplified85.0%
Taylor expanded in t around inf 82.2%
+-commutative82.2%
Simplified82.2%
if -9.4999999999999994e-67 < x < -3.55000000000000008e-97Initial program 87.5%
*-commutative87.5%
Simplified87.5%
Taylor expanded in y around inf 63.2%
times-frac63.4%
+-commutative63.4%
Simplified63.4%
Taylor expanded in x around 0 63.4%
Taylor expanded in t around 0 63.6%
associate-*r/63.6%
mul-1-neg63.6%
distribute-rgt-neg-out63.6%
Simplified63.6%
if -3.55000000000000008e-97 < x < -2.20000000000000004e-190 or 5.79999999999999973e-102 < x < 2.39999999999999991e-34Initial program 85.6%
*-commutative85.6%
Simplified85.6%
Taylor expanded in y around 0 67.6%
+-commutative67.6%
Simplified67.6%
clear-num67.6%
inv-pow67.6%
Applied egg-rr67.6%
unpow-167.6%
div-sub67.6%
sub-neg67.6%
associate-/l*67.6%
*-inverses67.6%
metadata-eval67.6%
Simplified67.6%
Taylor expanded in x around 0 51.3%
*-commutative51.3%
sub-neg51.3%
distribute-neg-frac51.3%
metadata-eval51.3%
distribute-rgt-in51.3%
*-lft-identity51.3%
associate-*l/51.4%
associate-*r/51.4%
mul-1-neg51.4%
unsub-neg51.4%
Simplified51.4%
if -2.20000000000000004e-190 < x < 5.79999999999999973e-102Initial program 93.0%
*-commutative93.0%
Simplified93.0%
Taylor expanded in y around inf 60.2%
times-frac48.9%
+-commutative48.9%
Simplified48.9%
Taylor expanded in x around 0 54.5%
if 2.39999999999999991e-34 < x Initial program 88.4%
*-commutative88.4%
Simplified88.4%
Taylor expanded in y around 0 92.6%
+-commutative92.6%
Simplified92.6%
Taylor expanded in x around inf 88.9%
Final simplification72.3%
(FPCore (x y z t) :precision binary64 (if (or (<= z -5.2e-46) (not (<= z 9.2e+33))) (/ (+ x (/ y t)) (+ x 1.0)) (/ (+ 1.0 (- x (/ y (/ x z)))) (+ x 1.0))))
double code(double x, double y, double z, double t) {
double tmp;
if ((z <= -5.2e-46) || !(z <= 9.2e+33)) {
tmp = (x + (y / t)) / (x + 1.0);
} else {
tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((z <= (-5.2d-46)) .or. (.not. (z <= 9.2d+33))) then
tmp = (x + (y / t)) / (x + 1.0d0)
else
tmp = (1.0d0 + (x - (y / (x / z)))) / (x + 1.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((z <= -5.2e-46) || !(z <= 9.2e+33)) {
tmp = (x + (y / t)) / (x + 1.0);
} else {
tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (z <= -5.2e-46) or not (z <= 9.2e+33): tmp = (x + (y / t)) / (x + 1.0) else: tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0) return tmp
function code(x, y, z, t) tmp = 0.0 if ((z <= -5.2e-46) || !(z <= 9.2e+33)) tmp = Float64(Float64(x + Float64(y / t)) / Float64(x + 1.0)); else tmp = Float64(Float64(1.0 + Float64(x - Float64(y / Float64(x / z)))) / Float64(x + 1.0)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((z <= -5.2e-46) || ~((z <= 9.2e+33))) tmp = (x + (y / t)) / (x + 1.0); else tmp = (1.0 + (x - (y / (x / z)))) / (x + 1.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[z, -5.2e-46], N[Not[LessEqual[z, 9.2e+33]], $MachinePrecision]], N[(N[(x + N[(y / t), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(x - N[(y / N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.2 \cdot 10^{-46} \lor \neg \left(z \leq 9.2 \cdot 10^{+33}\right):\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \left(x - \frac{y}{\frac{x}{z}}\right)}{x + 1}\\
\end{array}
\end{array}
if z < -5.2000000000000004e-46 or 9.20000000000000042e33 < z Initial program 77.6%
*-commutative77.6%
Simplified77.6%
Taylor expanded in z around inf 89.2%
if -5.2000000000000004e-46 < z < 9.20000000000000042e33Initial program 99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in t around 0 82.1%
+-commutative82.1%
mul-1-neg82.1%
unsub-neg82.1%
associate-/l*82.1%
+-commutative82.1%
Simplified82.1%
Final simplification85.9%
(FPCore (x y z t)
:precision binary64
(if (<= x -8.5e-10)
(/ x (+ x 1.0))
(if (<= x -1.76e-94)
(* z (/ y (- (* z t) x)))
(if (<= x 6.2e-102) (/ y t) (if (<= x 4e-33) (- x (/ x (* z t))) 1.0)))))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -8.5e-10) {
tmp = x / (x + 1.0);
} else if (x <= -1.76e-94) {
tmp = z * (y / ((z * t) - x));
} else if (x <= 6.2e-102) {
tmp = y / t;
} else if (x <= 4e-33) {
tmp = x - (x / (z * t));
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-8.5d-10)) then
tmp = x / (x + 1.0d0)
else if (x <= (-1.76d-94)) then
tmp = z * (y / ((z * t) - x))
else if (x <= 6.2d-102) then
tmp = y / t
else if (x <= 4d-33) then
tmp = x - (x / (z * t))
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -8.5e-10) {
tmp = x / (x + 1.0);
} else if (x <= -1.76e-94) {
tmp = z * (y / ((z * t) - x));
} else if (x <= 6.2e-102) {
tmp = y / t;
} else if (x <= 4e-33) {
tmp = x - (x / (z * t));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -8.5e-10: tmp = x / (x + 1.0) elif x <= -1.76e-94: tmp = z * (y / ((z * t) - x)) elif x <= 6.2e-102: tmp = y / t elif x <= 4e-33: tmp = x - (x / (z * t)) else: tmp = 1.0 return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -8.5e-10) tmp = Float64(x / Float64(x + 1.0)); elseif (x <= -1.76e-94) tmp = Float64(z * Float64(y / Float64(Float64(z * t) - x))); elseif (x <= 6.2e-102) tmp = Float64(y / t); elseif (x <= 4e-33) tmp = Float64(x - Float64(x / Float64(z * t))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -8.5e-10) tmp = x / (x + 1.0); elseif (x <= -1.76e-94) tmp = z * (y / ((z * t) - x)); elseif (x <= 6.2e-102) tmp = y / t; elseif (x <= 4e-33) tmp = x - (x / (z * t)); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -8.5e-10], N[(x / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -1.76e-94], N[(z * N[(y / N[(N[(z * t), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 6.2e-102], N[(y / t), $MachinePrecision], If[LessEqual[x, 4e-33], N[(x - N[(x / N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -8.5 \cdot 10^{-10}:\\
\;\;\;\;\frac{x}{x + 1}\\
\mathbf{elif}\;x \leq -1.76 \cdot 10^{-94}:\\
\;\;\;\;z \cdot \frac{y}{z \cdot t - x}\\
\mathbf{elif}\;x \leq 6.2 \cdot 10^{-102}:\\
\;\;\;\;\frac{y}{t}\\
\mathbf{elif}\;x \leq 4 \cdot 10^{-33}:\\
\;\;\;\;x - \frac{x}{z \cdot t}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -8.4999999999999996e-10Initial program 85.2%
*-commutative85.2%
Simplified85.2%
Taylor expanded in t around inf 93.2%
+-commutative93.2%
Simplified93.2%
if -8.4999999999999996e-10 < x < -1.7599999999999999e-94Initial program 85.3%
*-commutative85.3%
Simplified85.3%
Taylor expanded in y around inf 56.6%
times-frac66.2%
+-commutative66.2%
Simplified66.2%
Taylor expanded in x around 0 64.4%
if -1.7599999999999999e-94 < x < 6.20000000000000026e-102Initial program 90.2%
*-commutative90.2%
Simplified90.2%
Taylor expanded in y around inf 49.5%
times-frac42.7%
+-commutative42.7%
Simplified42.7%
Taylor expanded in x around 0 49.2%
if 6.20000000000000026e-102 < x < 4.0000000000000002e-33Initial program 88.5%
*-commutative88.5%
Simplified88.5%
Taylor expanded in y around 0 73.2%
+-commutative73.2%
Simplified73.2%
clear-num73.2%
inv-pow73.2%
Applied egg-rr73.2%
unpow-173.2%
div-sub73.2%
sub-neg73.2%
associate-/l*73.1%
*-inverses73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in x around 0 51.1%
*-commutative51.1%
sub-neg51.1%
distribute-neg-frac51.1%
metadata-eval51.1%
distribute-rgt-in51.1%
*-lft-identity51.1%
associate-*l/51.2%
associate-*r/51.2%
mul-1-neg51.2%
unsub-neg51.2%
Simplified51.2%
if 4.0000000000000002e-33 < x Initial program 88.4%
*-commutative88.4%
Simplified88.4%
Taylor expanded in y around 0 92.6%
+-commutative92.6%
Simplified92.6%
Taylor expanded in x around inf 88.9%
Final simplification72.6%
(FPCore (x y z t) :precision binary64 (if (<= x -3.65e+21) 1.0 (if (<= x 1.62e-35) (/ (+ x (/ y t)) (+ x 1.0)) 1.0)))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -3.65e+21) {
tmp = 1.0;
} else if (x <= 1.62e-35) {
tmp = (x + (y / t)) / (x + 1.0);
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-3.65d+21)) then
tmp = 1.0d0
else if (x <= 1.62d-35) then
tmp = (x + (y / t)) / (x + 1.0d0)
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -3.65e+21) {
tmp = 1.0;
} else if (x <= 1.62e-35) {
tmp = (x + (y / t)) / (x + 1.0);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -3.65e+21: tmp = 1.0 elif x <= 1.62e-35: tmp = (x + (y / t)) / (x + 1.0) else: tmp = 1.0 return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -3.65e+21) tmp = 1.0; elseif (x <= 1.62e-35) tmp = Float64(Float64(x + Float64(y / t)) / Float64(x + 1.0)); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -3.65e+21) tmp = 1.0; elseif (x <= 1.62e-35) tmp = (x + (y / t)) / (x + 1.0); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -3.65e+21], 1.0, If[LessEqual[x, 1.62e-35], N[(N[(x + N[(y / t), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.65 \cdot 10^{+21}:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 1.62 \cdot 10^{-35}:\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -3.65e21 or 1.62000000000000011e-35 < x Initial program 87.3%
*-commutative87.3%
Simplified87.3%
Taylor expanded in y around 0 93.7%
+-commutative93.7%
Simplified93.7%
Taylor expanded in x around inf 91.5%
if -3.65e21 < x < 1.62000000000000011e-35Initial program 88.8%
*-commutative88.8%
Simplified88.8%
Taylor expanded in z around inf 65.4%
Final simplification78.8%
(FPCore (x y z t) :precision binary64 (if (<= x -0.028) 1.0 (if (<= x 2.7e-80) (/ y t) 1.0)))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -0.028) {
tmp = 1.0;
} else if (x <= 2.7e-80) {
tmp = y / t;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-0.028d0)) then
tmp = 1.0d0
else if (x <= 2.7d-80) then
tmp = y / t
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -0.028) {
tmp = 1.0;
} else if (x <= 2.7e-80) {
tmp = y / t;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -0.028: tmp = 1.0 elif x <= 2.7e-80: tmp = y / t else: tmp = 1.0 return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -0.028) tmp = 1.0; elseif (x <= 2.7e-80) tmp = Float64(y / t); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -0.028) tmp = 1.0; elseif (x <= 2.7e-80) tmp = y / t; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -0.028], 1.0, If[LessEqual[x, 2.7e-80], N[(y / t), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.028:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 2.7 \cdot 10^{-80}:\\
\;\;\;\;\frac{y}{t}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -0.0280000000000000006 or 2.7000000000000002e-80 < x Initial program 86.7%
*-commutative86.7%
Simplified86.7%
Taylor expanded in y around 0 91.8%
+-commutative91.8%
Simplified91.8%
Taylor expanded in x around inf 85.2%
if -0.0280000000000000006 < x < 2.7000000000000002e-80Initial program 89.8%
*-commutative89.8%
Simplified89.8%
Taylor expanded in y around inf 50.0%
times-frac46.7%
+-commutative46.7%
Simplified46.7%
Taylor expanded in x around 0 45.3%
Final simplification68.3%
(FPCore (x y z t) :precision binary64 (if (<= x -8.8e-13) (/ x (+ x 1.0)) (if (<= x 2.75e-80) (/ y t) 1.0)))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -8.8e-13) {
tmp = x / (x + 1.0);
} else if (x <= 2.75e-80) {
tmp = y / t;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-8.8d-13)) then
tmp = x / (x + 1.0d0)
else if (x <= 2.75d-80) then
tmp = y / t
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -8.8e-13) {
tmp = x / (x + 1.0);
} else if (x <= 2.75e-80) {
tmp = y / t;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -8.8e-13: tmp = x / (x + 1.0) elif x <= 2.75e-80: tmp = y / t else: tmp = 1.0 return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -8.8e-13) tmp = Float64(x / Float64(x + 1.0)); elseif (x <= 2.75e-80) tmp = Float64(y / t); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -8.8e-13) tmp = x / (x + 1.0); elseif (x <= 2.75e-80) tmp = y / t; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -8.8e-13], N[(x / N[(x + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2.75e-80], N[(y / t), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -8.8 \cdot 10^{-13}:\\
\;\;\;\;\frac{x}{x + 1}\\
\mathbf{elif}\;x \leq 2.75 \cdot 10^{-80}:\\
\;\;\;\;\frac{y}{t}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < -8.79999999999999986e-13Initial program 85.4%
*-commutative85.4%
Simplified85.4%
Taylor expanded in t around inf 91.8%
+-commutative91.8%
Simplified91.8%
if -8.79999999999999986e-13 < x < 2.7499999999999998e-80Initial program 89.5%
*-commutative89.5%
Simplified89.5%
Taylor expanded in y around inf 50.4%
times-frac47.0%
+-commutative47.0%
Simplified47.0%
Taylor expanded in x around 0 46.5%
if 2.7499999999999998e-80 < x Initial program 88.0%
*-commutative88.0%
Simplified88.0%
Taylor expanded in y around 0 90.5%
+-commutative90.5%
Simplified90.5%
Taylor expanded in x around inf 79.6%
Final simplification68.9%
(FPCore (x y z t) :precision binary64 1.0)
double code(double x, double y, double z, double t) {
return 1.0;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 1.0d0
end function
public static double code(double x, double y, double z, double t) {
return 1.0;
}
def code(x, y, z, t): return 1.0
function code(x, y, z, t) return 1.0 end
function tmp = code(x, y, z, t) tmp = 1.0; end
code[x_, y_, z_, t_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 88.0%
*-commutative88.0%
Simplified88.0%
Taylor expanded in y around 0 71.3%
+-commutative71.3%
Simplified71.3%
Taylor expanded in x around inf 56.7%
Final simplification56.7%
(FPCore (x y z t) :precision binary64 (/ (+ x (- (/ y (- t (/ x z))) (/ x (- (* t z) x)))) (+ x 1.0)))
double code(double x, double y, double z, double t) {
return (x + ((y / (t - (x / z))) - (x / ((t * z) - x)))) / (x + 1.0);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x + ((y / (t - (x / z))) - (x / ((t * z) - x)))) / (x + 1.0d0)
end function
public static double code(double x, double y, double z, double t) {
return (x + ((y / (t - (x / z))) - (x / ((t * z) - x)))) / (x + 1.0);
}
def code(x, y, z, t): return (x + ((y / (t - (x / z))) - (x / ((t * z) - x)))) / (x + 1.0)
function code(x, y, z, t) return Float64(Float64(x + Float64(Float64(y / Float64(t - Float64(x / z))) - Float64(x / Float64(Float64(t * z) - x)))) / Float64(x + 1.0)) end
function tmp = code(x, y, z, t) tmp = (x + ((y / (t - (x / z))) - (x / ((t * z) - x)))) / (x + 1.0); end
code[x_, y_, z_, t_] := N[(N[(x + N[(N[(y / N[(t - N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x / N[(N[(t * z), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x + \left(\frac{y}{t - \frac{x}{z}} - \frac{x}{t \cdot z - x}\right)}{x + 1}
\end{array}
herbie shell --seed 2023257
(FPCore (x y z t)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, A"
:precision binary64
:herbie-target
(/ (+ x (- (/ y (- t (/ x z))) (/ x (- (* t z) x)))) (+ x 1.0))
(/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1.0)))