
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (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 - x) / (a - z)))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
def code(x, y, z, t, a): return x + ((y - z) * ((t - x) / (a - z)))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function tmp = code(x, y, z, t, a) tmp = x + ((y - z) * ((t - x) / (a - z))); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (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 - x) / (a - z)))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
def code(x, y, z, t, a): return x + ((y - z) * ((t - x) / (a - z)))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function tmp = code(x, y, z, t, a) tmp = x + ((y - z) * ((t - x) / (a - z))); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- y z) (- a z)))
(t_2 (- x (* (- y z) (/ (- x t) (- a z)))))
(t_3 (cbrt (- t x))))
(if (<= t_2 -5e-274)
(fma (- t x) t_1 x)
(if (<= t_2 3e-279)
(+ (+ (- t (/ (* y t) z)) (/ (- y a) (/ z x))) (/ a (/ z t)))
(if (<= t_2 2e+145)
(+
(/ (* (- y z) t) (- a z))
(* x (- (+ (/ z (- a z)) 1.0) (/ y (- a z)))))
(if (<= t_2 5e+303)
(fma (/ t (- a z)) (- y z) (* x (- (/ (- z y) (- a z)) -1.0)))
(fma (pow t_3 2.0) (* t_1 t_3) x)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y - z) / (a - z);
double t_2 = x - ((y - z) * ((x - t) / (a - z)));
double t_3 = cbrt((t - x));
double tmp;
if (t_2 <= -5e-274) {
tmp = fma((t - x), t_1, x);
} else if (t_2 <= 3e-279) {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
} else if (t_2 <= 2e+145) {
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z))));
} else if (t_2 <= 5e+303) {
tmp = fma((t / (a - z)), (y - z), (x * (((z - y) / (a - z)) - -1.0)));
} else {
tmp = fma(pow(t_3, 2.0), (t_1 * t_3), x);
}
return tmp;
}
function code(x, y, z, t, a) t_1 = Float64(Float64(y - z) / Float64(a - z)) t_2 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) t_3 = cbrt(Float64(t - x)) tmp = 0.0 if (t_2 <= -5e-274) tmp = fma(Float64(t - x), t_1, x); elseif (t_2 <= 3e-279) tmp = Float64(Float64(Float64(t - Float64(Float64(y * t) / z)) + Float64(Float64(y - a) / Float64(z / x))) + Float64(a / Float64(z / t))); elseif (t_2 <= 2e+145) tmp = Float64(Float64(Float64(Float64(y - z) * t) / Float64(a - z)) + Float64(x * Float64(Float64(Float64(z / Float64(a - z)) + 1.0) - Float64(y / Float64(a - z))))); elseif (t_2 <= 5e+303) tmp = fma(Float64(t / Float64(a - z)), Float64(y - z), Float64(x * Float64(Float64(Float64(z - y) / Float64(a - z)) - -1.0))); else tmp = fma((t_3 ^ 2.0), Float64(t_1 * t_3), x); end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Power[N[(t - x), $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[t$95$2, -5e-274], N[(N[(t - x), $MachinePrecision] * t$95$1 + x), $MachinePrecision], If[LessEqual[t$95$2, 3e-279], N[(N[(N[(t - N[(N[(y * t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - a), $MachinePrecision] / N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a / N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 2e+145], N[(N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(N[(z / N[(a - z), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5e+303], N[(N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision] * N[(y - z), $MachinePrecision] + N[(x * N[(N[(N[(z - y), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[t$95$3, 2.0], $MachinePrecision] * N[(t$95$1 * t$95$3), $MachinePrecision] + x), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y - z}{a - z}\\
t_2 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
t_3 := \sqrt[3]{t - x}\\
\mathbf{if}\;t_2 \leq -5 \cdot 10^{-274}:\\
\;\;\;\;\mathsf{fma}\left(t - x, t_1, x\right)\\
\mathbf{elif}\;t_2 \leq 3 \cdot 10^{-279}:\\
\;\;\;\;\left(\left(t - \frac{y \cdot t}{z}\right) + \frac{y - a}{\frac{z}{x}}\right) + \frac{a}{\frac{z}{t}}\\
\mathbf{elif}\;t_2 \leq 2 \cdot 10^{+145}:\\
\;\;\;\;\frac{\left(y - z\right) \cdot t}{a - z} + x \cdot \left(\left(\frac{z}{a - z} + 1\right) - \frac{y}{a - z}\right)\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{+303}:\\
\;\;\;\;\mathsf{fma}\left(\frac{t}{a - z}, y - z, x \cdot \left(\frac{z - y}{a - z} - -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left({t_3}^{2}, t_1 \cdot t_3, x\right)\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -5e-274Initial program 91.0%
+-commutative91.0%
associate-*r/76.0%
*-commutative76.0%
associate-*r/94.1%
fma-def94.1%
Simplified94.1%
if -5e-274 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 3e-279Initial program 3.2%
Taylor expanded in x around -inf 56.6%
Taylor expanded in z around inf 84.0%
sub-neg84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
*-commutative84.0%
associate-/l*96.5%
mul-1-neg96.5%
unsub-neg96.5%
mul-1-neg96.5%
remove-double-neg96.5%
associate-/l*96.5%
Simplified96.5%
if 3e-279 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 2e145Initial program 86.0%
Taylor expanded in x around -inf 96.8%
if 2e145 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 4.9999999999999997e303Initial program 92.2%
Taylor expanded in x around -inf 63.2%
*-un-lft-identity63.2%
fma-def63.2%
associate-/l*94.6%
mul-1-neg94.6%
*-commutative94.6%
associate--r+94.6%
div-sub94.6%
Applied egg-rr94.6%
fma-udef94.6%
*-lft-identity94.6%
associate-/r/94.6%
fma-def94.6%
*-commutative94.6%
distribute-rgt-neg-in94.6%
sub-neg94.6%
metadata-eval94.6%
Simplified94.6%
if 4.9999999999999997e303 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 78.6%
+-commutative78.6%
*-commutative78.6%
associate-*l/88.3%
associate-*r/94.2%
add-cube-cbrt94.6%
associate-*l*94.6%
fma-def94.6%
pow294.6%
Applied egg-rr94.6%
Final simplification95.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (- t x) (/ (- y z) (- a z)) x))
(t_2 (- x (* (- y z) (/ (- x t) (- a z))))))
(if (<= t_2 -5e-274)
t_1
(if (<= t_2 3e-279)
(+ (+ (- t (/ (* y t) z)) (/ (- y a) (/ z x))) (/ a (/ z t)))
(if (<= t_2 2e+145)
(+
(/ (* (- y z) t) (- a z))
(* x (- (+ (/ z (- a z)) 1.0) (/ y (- a z)))))
(if (<= t_2 1e+281)
(fma (/ t (- a z)) (- y z) (* x (- (/ (- z y) (- a z)) -1.0)))
t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((t - x), ((y - z) / (a - z)), x);
double t_2 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if (t_2 <= -5e-274) {
tmp = t_1;
} else if (t_2 <= 3e-279) {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
} else if (t_2 <= 2e+145) {
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z))));
} else if (t_2 <= 1e+281) {
tmp = fma((t / (a - z)), (y - z), (x * (((z - y) / (a - z)) - -1.0)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(t - x), Float64(Float64(y - z) / Float64(a - z)), x) t_2 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if (t_2 <= -5e-274) tmp = t_1; elseif (t_2 <= 3e-279) tmp = Float64(Float64(Float64(t - Float64(Float64(y * t) / z)) + Float64(Float64(y - a) / Float64(z / x))) + Float64(a / Float64(z / t))); elseif (t_2 <= 2e+145) tmp = Float64(Float64(Float64(Float64(y - z) * t) / Float64(a - z)) + Float64(x * Float64(Float64(Float64(z / Float64(a - z)) + 1.0) - Float64(y / Float64(a - z))))); elseif (t_2 <= 1e+281) tmp = fma(Float64(t / Float64(a - z)), Float64(y - z), Float64(x * Float64(Float64(Float64(z - y) / Float64(a - z)) - -1.0))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]}, Block[{t$95$2 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -5e-274], t$95$1, If[LessEqual[t$95$2, 3e-279], N[(N[(N[(t - N[(N[(y * t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - a), $MachinePrecision] / N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a / N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 2e+145], N[(N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(N[(z / N[(a - z), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 1e+281], N[(N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision] * N[(y - z), $MachinePrecision] + N[(x * N[(N[(N[(z - y), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
t_2 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_2 \leq -5 \cdot 10^{-274}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 3 \cdot 10^{-279}:\\
\;\;\;\;\left(\left(t - \frac{y \cdot t}{z}\right) + \frac{y - a}{\frac{z}{x}}\right) + \frac{a}{\frac{z}{t}}\\
\mathbf{elif}\;t_2 \leq 2 \cdot 10^{+145}:\\
\;\;\;\;\frac{\left(y - z\right) \cdot t}{a - z} + x \cdot \left(\left(\frac{z}{a - z} + 1\right) - \frac{y}{a - z}\right)\\
\mathbf{elif}\;t_2 \leq 10^{+281}:\\
\;\;\;\;\mathsf{fma}\left(\frac{t}{a - z}, y - z, x \cdot \left(\frac{z - y}{a - z} - -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -5e-274 or 1e281 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 89.7%
+-commutative89.7%
associate-*r/75.2%
*-commutative75.2%
associate-*r/94.4%
fma-def94.4%
Simplified94.4%
if -5e-274 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 3e-279Initial program 3.2%
Taylor expanded in x around -inf 56.6%
Taylor expanded in z around inf 84.0%
sub-neg84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
*-commutative84.0%
associate-/l*96.5%
mul-1-neg96.5%
unsub-neg96.5%
mul-1-neg96.5%
remove-double-neg96.5%
associate-/l*96.5%
Simplified96.5%
if 3e-279 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 2e145Initial program 86.0%
Taylor expanded in x around -inf 96.8%
if 2e145 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 1e281Initial program 90.7%
Taylor expanded in x around -inf 64.6%
*-un-lft-identity64.6%
fma-def64.6%
associate-/l*93.6%
mul-1-neg93.6%
*-commutative93.6%
associate--r+93.6%
div-sub93.6%
Applied egg-rr93.6%
fma-udef93.6%
*-lft-identity93.6%
associate-/r/93.6%
fma-def93.6%
*-commutative93.6%
distribute-rgt-neg-in93.6%
sub-neg93.6%
metadata-eval93.6%
Simplified93.6%
Final simplification95.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* (- y z) (/ (- x t) (- a z))))))
(if (or (<= t_1 -5e-274) (not (<= t_1 0.0)))
(fma (- t x) (/ (- y z) (- a z)) x)
(+ (+ (- t (/ (* y t) z)) (/ (- y a) (/ z x))) (/ a (/ z t))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if ((t_1 <= -5e-274) || !(t_1 <= 0.0)) {
tmp = fma((t - x), ((y - z) / (a - z)), x);
} else {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
}
return tmp;
}
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if ((t_1 <= -5e-274) || !(t_1 <= 0.0)) tmp = fma(Float64(t - x), Float64(Float64(y - z) / Float64(a - z)), x); else tmp = Float64(Float64(Float64(t - Float64(Float64(y * t) / z)) + Float64(Float64(y - a) / Float64(z / x))) + Float64(a / Float64(z / t))); end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -5e-274], N[Not[LessEqual[t$95$1, 0.0]], $MachinePrecision]], N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(N[(N[(t - N[(N[(y * t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - a), $MachinePrecision] / N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a / N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_1 \leq -5 \cdot 10^{-274} \lor \neg \left(t_1 \leq 0\right):\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(t - \frac{y \cdot t}{z}\right) + \frac{y - a}{\frac{z}{x}}\right) + \frac{a}{\frac{z}{t}}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -5e-274 or 0.0 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 88.5%
+-commutative88.5%
associate-*r/75.7%
*-commutative75.7%
associate-*r/93.1%
fma-def93.2%
Simplified93.2%
if -5e-274 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.0Initial program 3.1%
Taylor expanded in x around -inf 55.1%
Taylor expanded in z around inf 83.4%
sub-neg83.4%
+-commutative83.4%
mul-1-neg83.4%
unsub-neg83.4%
+-commutative83.4%
mul-1-neg83.4%
unsub-neg83.4%
*-commutative83.4%
associate-/l*96.4%
mul-1-neg96.4%
unsub-neg96.4%
mul-1-neg96.4%
remove-double-neg96.4%
associate-/l*96.4%
Simplified96.4%
Final simplification93.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* (- y z) (/ (- x t) (- a z))))))
(if (<= t_1 -1e-198)
t_1
(if (<= t_1 3e-279)
(+ (+ (- t (/ (* y t) z)) (/ (- y a) (/ z x))) (/ a (/ z t)))
(if (<= t_1 5e+169)
(+
(/ (* (- y z) t) (- a z))
(* x (- (+ (/ z (- a z)) 1.0) (/ y (- a z)))))
t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if (t_1 <= -1e-198) {
tmp = t_1;
} else if (t_1 <= 3e-279) {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
} else if (t_1 <= 5e+169) {
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y - z) * ((x - t) / (a - z)))
if (t_1 <= (-1d-198)) then
tmp = t_1
else if (t_1 <= 3d-279) then
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t))
else if (t_1 <= 5d+169) then
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0d0) - (y / (a - z))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if (t_1 <= -1e-198) {
tmp = t_1;
} else if (t_1 <= 3e-279) {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
} else if (t_1 <= 5e+169) {
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - ((y - z) * ((x - t) / (a - z))) tmp = 0 if t_1 <= -1e-198: tmp = t_1 elif t_1 <= 3e-279: tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t)) elif t_1 <= 5e+169: tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if (t_1 <= -1e-198) tmp = t_1; elseif (t_1 <= 3e-279) tmp = Float64(Float64(Float64(t - Float64(Float64(y * t) / z)) + Float64(Float64(y - a) / Float64(z / x))) + Float64(a / Float64(z / t))); elseif (t_1 <= 5e+169) tmp = Float64(Float64(Float64(Float64(y - z) * t) / Float64(a - z)) + Float64(x * Float64(Float64(Float64(z / Float64(a - z)) + 1.0) - Float64(y / Float64(a - z))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - ((y - z) * ((x - t) / (a - z))); tmp = 0.0; if (t_1 <= -1e-198) tmp = t_1; elseif (t_1 <= 3e-279) tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t)); elseif (t_1 <= 5e+169) tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-198], t$95$1, If[LessEqual[t$95$1, 3e-279], N[(N[(N[(t - N[(N[(y * t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - a), $MachinePrecision] / N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a / N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e+169], N[(N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(N[(z / N[(a - z), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_1 \leq -1 \cdot 10^{-198}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_1 \leq 3 \cdot 10^{-279}:\\
\;\;\;\;\left(\left(t - \frac{y \cdot t}{z}\right) + \frac{y - a}{\frac{z}{x}}\right) + \frac{a}{\frac{z}{t}}\\
\mathbf{elif}\;t_1 \leq 5 \cdot 10^{+169}:\\
\;\;\;\;\frac{\left(y - z\right) \cdot t}{a - z} + x \cdot \left(\left(\frac{z}{a - z} + 1\right) - \frac{y}{a - z}\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -9.9999999999999991e-199 or 5.00000000000000017e169 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 91.1%
if -9.9999999999999991e-199 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 3e-279Initial program 6.6%
Taylor expanded in x around -inf 61.7%
Taylor expanded in z around inf 83.0%
sub-neg83.0%
+-commutative83.0%
mul-1-neg83.0%
unsub-neg83.0%
+-commutative83.0%
mul-1-neg83.0%
unsub-neg83.0%
*-commutative83.0%
associate-/l*94.0%
mul-1-neg94.0%
unsub-neg94.0%
mul-1-neg94.0%
remove-double-neg94.0%
associate-/l*94.0%
Simplified94.0%
if 3e-279 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 5.00000000000000017e169Initial program 87.0%
Taylor expanded in x around -inf 97.0%
Final simplification93.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* (- y z) (/ (- x t) (- a z))))))
(if (or (<= t_1 -1e-198) (not (<= t_1 5e-271)))
t_1
(+ (+ (- t (/ (* y t) z)) (/ (- y a) (/ z x))) (/ a (/ z t))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) {
tmp = t_1;
} else {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - ((y - z) * ((x - t) / (a - z)))
if ((t_1 <= (-1d-198)) .or. (.not. (t_1 <= 5d-271))) then
tmp = t_1
else
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) {
tmp = t_1;
} else {
tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - ((y - z) * ((x - t) / (a - z))) tmp = 0 if (t_1 <= -1e-198) or not (t_1 <= 5e-271): tmp = t_1 else: tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t)) return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) tmp = t_1; else tmp = Float64(Float64(Float64(t - Float64(Float64(y * t) / z)) + Float64(Float64(y - a) / Float64(z / x))) + Float64(a / Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - ((y - z) * ((x - t) / (a - z))); tmp = 0.0; if ((t_1 <= -1e-198) || ~((t_1 <= 5e-271))) tmp = t_1; else tmp = ((t - ((y * t) / z)) + ((y - a) / (z / x))) + (a / (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -1e-198], N[Not[LessEqual[t$95$1, 5e-271]], $MachinePrecision]], t$95$1, N[(N[(N[(t - N[(N[(y * t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - a), $MachinePrecision] / N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a / N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_1 \leq -1 \cdot 10^{-198} \lor \neg \left(t_1 \leq 5 \cdot 10^{-271}\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(t - \frac{y \cdot t}{z}\right) + \frac{y - a}{\frac{z}{x}}\right) + \frac{a}{\frac{z}{t}}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -9.9999999999999991e-199 or 5.0000000000000002e-271 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 90.6%
if -9.9999999999999991e-199 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 5.0000000000000002e-271Initial program 6.5%
Taylor expanded in x around -inf 63.8%
Taylor expanded in z around inf 83.9%
sub-neg83.9%
+-commutative83.9%
mul-1-neg83.9%
unsub-neg83.9%
+-commutative83.9%
mul-1-neg83.9%
unsub-neg83.9%
*-commutative83.9%
associate-/l*94.4%
mul-1-neg94.4%
unsub-neg94.4%
mul-1-neg94.4%
remove-double-neg94.4%
associate-/l*94.4%
Simplified94.4%
Final simplification91.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* (- y z) (/ (- x t) (- a z))))))
(if (or (<= t_1 -1e-198) (not (<= t_1 5e-271)))
t_1
(+ t (/ (* (- y a) (- x t)) z)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((y - z) * ((x - t) / (a - z)));
double tmp;
if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) {
tmp = t_1;
} else {
tmp = t + (((y - a) * (x - 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 - ((y - z) * ((x - t) / (a - z)))
if ((t_1 <= (-1d-198)) .or. (.not. (t_1 <= 5d-271))) then
tmp = t_1
else
tmp = t + (((y - a) * (x - 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 - ((y - z) * ((x - t) / (a - z)));
double tmp;
if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) {
tmp = t_1;
} else {
tmp = t + (((y - a) * (x - t)) / z);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - ((y - z) * ((x - t) / (a - z))) tmp = 0 if (t_1 <= -1e-198) or not (t_1 <= 5e-271): tmp = t_1 else: tmp = t + (((y - a) * (x - t)) / z) return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(y - z) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if ((t_1 <= -1e-198) || !(t_1 <= 5e-271)) tmp = t_1; else tmp = Float64(t + Float64(Float64(Float64(y - a) * Float64(x - t)) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - ((y - z) * ((x - t) / (a - z))); tmp = 0.0; if ((t_1 <= -1e-198) || ~((t_1 <= 5e-271))) tmp = t_1; else tmp = t + (((y - a) * (x - t)) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(y - z), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -1e-198], N[Not[LessEqual[t$95$1, 5e-271]], $MachinePrecision]], t$95$1, N[(t + N[(N[(N[(y - a), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(y - z\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_1 \leq -1 \cdot 10^{-198} \lor \neg \left(t_1 \leq 5 \cdot 10^{-271}\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t + \frac{\left(y - a\right) \cdot \left(x - t\right)}{z}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -9.9999999999999991e-199 or 5.0000000000000002e-271 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 90.6%
if -9.9999999999999991e-199 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 5.0000000000000002e-271Initial program 6.5%
Taylor expanded in z around inf 84.1%
+-commutative84.1%
associate--l+84.1%
associate-*r/84.1%
associate-*r/84.1%
div-sub84.1%
distribute-lft-out--84.1%
mul-1-neg84.1%
distribute-neg-frac84.1%
unsub-neg84.1%
distribute-rgt-out--84.1%
Simplified84.1%
Final simplification89.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ (- y z) (- a z))))
(t_2 (* (- t x) (/ y (- a z))))
(t_3 (- x (/ y (/ a x)))))
(if (<= t -6500000000.0)
t_1
(if (<= t -4.5e-5)
t_2
(if (<= t -5e-12)
(/ (- t) (/ z (- y z)))
(if (<= t -2.7e-186)
x
(if (<= t 1.3e-287)
t_3
(if (<= t 8.5e-156) t_2 (if (<= t 3.2e-19) t_3 t_1)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * ((y - z) / (a - z));
double t_2 = (t - x) * (y / (a - z));
double t_3 = x - (y / (a / x));
double tmp;
if (t <= -6500000000.0) {
tmp = t_1;
} else if (t <= -4.5e-5) {
tmp = t_2;
} else if (t <= -5e-12) {
tmp = -t / (z / (y - z));
} else if (t <= -2.7e-186) {
tmp = x;
} else if (t <= 1.3e-287) {
tmp = t_3;
} else if (t <= 8.5e-156) {
tmp = t_2;
} else if (t <= 3.2e-19) {
tmp = t_3;
} 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) :: t_3
real(8) :: tmp
t_1 = t * ((y - z) / (a - z))
t_2 = (t - x) * (y / (a - z))
t_3 = x - (y / (a / x))
if (t <= (-6500000000.0d0)) then
tmp = t_1
else if (t <= (-4.5d-5)) then
tmp = t_2
else if (t <= (-5d-12)) then
tmp = -t / (z / (y - z))
else if (t <= (-2.7d-186)) then
tmp = x
else if (t <= 1.3d-287) then
tmp = t_3
else if (t <= 8.5d-156) then
tmp = t_2
else if (t <= 3.2d-19) then
tmp = t_3
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 = t * ((y - z) / (a - z));
double t_2 = (t - x) * (y / (a - z));
double t_3 = x - (y / (a / x));
double tmp;
if (t <= -6500000000.0) {
tmp = t_1;
} else if (t <= -4.5e-5) {
tmp = t_2;
} else if (t <= -5e-12) {
tmp = -t / (z / (y - z));
} else if (t <= -2.7e-186) {
tmp = x;
} else if (t <= 1.3e-287) {
tmp = t_3;
} else if (t <= 8.5e-156) {
tmp = t_2;
} else if (t <= 3.2e-19) {
tmp = t_3;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * ((y - z) / (a - z)) t_2 = (t - x) * (y / (a - z)) t_3 = x - (y / (a / x)) tmp = 0 if t <= -6500000000.0: tmp = t_1 elif t <= -4.5e-5: tmp = t_2 elif t <= -5e-12: tmp = -t / (z / (y - z)) elif t <= -2.7e-186: tmp = x elif t <= 1.3e-287: tmp = t_3 elif t <= 8.5e-156: tmp = t_2 elif t <= 3.2e-19: tmp = t_3 else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(Float64(y - z) / Float64(a - z))) t_2 = Float64(Float64(t - x) * Float64(y / Float64(a - z))) t_3 = Float64(x - Float64(y / Float64(a / x))) tmp = 0.0 if (t <= -6500000000.0) tmp = t_1; elseif (t <= -4.5e-5) tmp = t_2; elseif (t <= -5e-12) tmp = Float64(Float64(-t) / Float64(z / Float64(y - z))); elseif (t <= -2.7e-186) tmp = x; elseif (t <= 1.3e-287) tmp = t_3; elseif (t <= 8.5e-156) tmp = t_2; elseif (t <= 3.2e-19) tmp = t_3; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * ((y - z) / (a - z)); t_2 = (t - x) * (y / (a - z)); t_3 = x - (y / (a / x)); tmp = 0.0; if (t <= -6500000000.0) tmp = t_1; elseif (t <= -4.5e-5) tmp = t_2; elseif (t <= -5e-12) tmp = -t / (z / (y - z)); elseif (t <= -2.7e-186) tmp = x; elseif (t <= 1.3e-287) tmp = t_3; elseif (t <= 8.5e-156) tmp = t_2; elseif (t <= 3.2e-19) tmp = t_3; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t - x), $MachinePrecision] * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -6500000000.0], t$95$1, If[LessEqual[t, -4.5e-5], t$95$2, If[LessEqual[t, -5e-12], N[((-t) / N[(z / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -2.7e-186], x, If[LessEqual[t, 1.3e-287], t$95$3, If[LessEqual[t, 8.5e-156], t$95$2, If[LessEqual[t, 3.2e-19], t$95$3, t$95$1]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y - z}{a - z}\\
t_2 := \left(t - x\right) \cdot \frac{y}{a - z}\\
t_3 := x - \frac{y}{\frac{a}{x}}\\
\mathbf{if}\;t \leq -6500000000:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -4.5 \cdot 10^{-5}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq -5 \cdot 10^{-12}:\\
\;\;\;\;\frac{-t}{\frac{z}{y - z}}\\
\mathbf{elif}\;t \leq -2.7 \cdot 10^{-186}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 1.3 \cdot 10^{-287}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;t \leq 8.5 \cdot 10^{-156}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 3.2 \cdot 10^{-19}:\\
\;\;\;\;t_3\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -6.5e9 or 3.19999999999999982e-19 < t Initial program 86.3%
Taylor expanded in t around inf 77.7%
div-sub77.7%
Simplified77.7%
if -6.5e9 < t < -4.50000000000000028e-5 or 1.3e-287 < t < 8.5e-156Initial program 72.6%
add-cube-cbrt72.6%
fma-def72.6%
pow272.6%
*-commutative72.6%
associate-*l/63.0%
associate-*r/78.6%
Applied egg-rr78.6%
Taylor expanded in y around -inf 62.3%
associate-*l/73.7%
Simplified73.7%
if -4.50000000000000028e-5 < t < -4.9999999999999997e-12Initial program 99.2%
Taylor expanded in a around 0 99.2%
associate-*r/99.2%
neg-mul-199.2%
Simplified99.2%
Taylor expanded in x around 0 77.1%
mul-1-neg77.1%
associate-/l*77.1%
Simplified77.1%
if -4.9999999999999997e-12 < t < -2.6999999999999999e-186Initial program 68.1%
Taylor expanded in a around inf 44.1%
if -2.6999999999999999e-186 < t < 1.3e-287 or 8.5e-156 < t < 3.19999999999999982e-19Initial program 73.6%
Taylor expanded in z around 0 55.1%
Taylor expanded in t around 0 53.7%
+-commutative53.7%
mul-1-neg53.7%
unsub-neg53.7%
associate-/l*56.7%
Simplified56.7%
Final simplification67.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* y (/ (- x t) z))) (t_2 (* t (- 1.0 (/ y z)))))
(if (<= a -5.6e+119)
x
(if (<= a -2.5e-62)
(* (- y z) (/ t a))
(if (<= a -6e-122)
t_1
(if (<= a -4.2e-268)
t_2
(if (<= a 1.15e-49) t_1 (if (<= a 1.32e+38) t_2 x))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * ((x - t) / z);
double t_2 = t * (1.0 - (y / z));
double tmp;
if (a <= -5.6e+119) {
tmp = x;
} else if (a <= -2.5e-62) {
tmp = (y - z) * (t / a);
} else if (a <= -6e-122) {
tmp = t_1;
} else if (a <= -4.2e-268) {
tmp = t_2;
} else if (a <= 1.15e-49) {
tmp = t_1;
} else if (a <= 1.32e+38) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = y * ((x - t) / z)
t_2 = t * (1.0d0 - (y / z))
if (a <= (-5.6d+119)) then
tmp = x
else if (a <= (-2.5d-62)) then
tmp = (y - z) * (t / a)
else if (a <= (-6d-122)) then
tmp = t_1
else if (a <= (-4.2d-268)) then
tmp = t_2
else if (a <= 1.15d-49) then
tmp = t_1
else if (a <= 1.32d+38) then
tmp = t_2
else
tmp = 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 * ((x - t) / z);
double t_2 = t * (1.0 - (y / z));
double tmp;
if (a <= -5.6e+119) {
tmp = x;
} else if (a <= -2.5e-62) {
tmp = (y - z) * (t / a);
} else if (a <= -6e-122) {
tmp = t_1;
} else if (a <= -4.2e-268) {
tmp = t_2;
} else if (a <= 1.15e-49) {
tmp = t_1;
} else if (a <= 1.32e+38) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * ((x - t) / z) t_2 = t * (1.0 - (y / z)) tmp = 0 if a <= -5.6e+119: tmp = x elif a <= -2.5e-62: tmp = (y - z) * (t / a) elif a <= -6e-122: tmp = t_1 elif a <= -4.2e-268: tmp = t_2 elif a <= 1.15e-49: tmp = t_1 elif a <= 1.32e+38: tmp = t_2 else: tmp = x return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(Float64(x - t) / z)) t_2 = Float64(t * Float64(1.0 - Float64(y / z))) tmp = 0.0 if (a <= -5.6e+119) tmp = x; elseif (a <= -2.5e-62) tmp = Float64(Float64(y - z) * Float64(t / a)); elseif (a <= -6e-122) tmp = t_1; elseif (a <= -4.2e-268) tmp = t_2; elseif (a <= 1.15e-49) tmp = t_1; elseif (a <= 1.32e+38) tmp = t_2; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * ((x - t) / z); t_2 = t * (1.0 - (y / z)); tmp = 0.0; if (a <= -5.6e+119) tmp = x; elseif (a <= -2.5e-62) tmp = (y - z) * (t / a); elseif (a <= -6e-122) tmp = t_1; elseif (a <= -4.2e-268) tmp = t_2; elseif (a <= 1.15e-49) tmp = t_1; elseif (a <= 1.32e+38) tmp = t_2; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -5.6e+119], x, If[LessEqual[a, -2.5e-62], N[(N[(y - z), $MachinePrecision] * N[(t / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -6e-122], t$95$1, If[LessEqual[a, -4.2e-268], t$95$2, If[LessEqual[a, 1.15e-49], t$95$1, If[LessEqual[a, 1.32e+38], t$95$2, x]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{x - t}{z}\\
t_2 := t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{if}\;a \leq -5.6 \cdot 10^{+119}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -2.5 \cdot 10^{-62}:\\
\;\;\;\;\left(y - z\right) \cdot \frac{t}{a}\\
\mathbf{elif}\;a \leq -6 \cdot 10^{-122}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -4.2 \cdot 10^{-268}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;a \leq 1.15 \cdot 10^{-49}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq 1.32 \cdot 10^{+38}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -5.60000000000000026e119 or 1.32e38 < a Initial program 89.6%
Taylor expanded in a around inf 56.3%
if -5.60000000000000026e119 < a < -2.5000000000000001e-62Initial program 85.4%
add-cube-cbrt85.6%
fma-def85.6%
pow285.6%
*-commutative85.6%
associate-*l/70.8%
associate-*r/88.9%
Applied egg-rr88.9%
Taylor expanded in x around 0 50.8%
associate-/l*59.9%
associate-/r/56.7%
Simplified56.7%
Taylor expanded in a around inf 41.9%
if -2.5000000000000001e-62 < a < -6.00000000000000009e-122 or -4.19999999999999996e-268 < a < 1.15e-49Initial program 66.8%
Taylor expanded in a around 0 52.7%
associate-*r/52.7%
neg-mul-152.7%
Simplified52.7%
Taylor expanded in y around inf 52.0%
div-sub52.1%
Simplified52.1%
if -6.00000000000000009e-122 < a < -4.19999999999999996e-268 or 1.15e-49 < a < 1.32e38Initial program 76.2%
Taylor expanded in a around 0 49.3%
+-commutative49.3%
mul-1-neg49.3%
unsub-neg49.3%
associate-/l*50.7%
Simplified50.7%
Taylor expanded in t around inf 58.5%
Final simplification53.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- t x) (/ y (- a z)))))
(if (<= y -1.02e+155)
t_1
(if (<= y -1.9e-65)
(+ x (/ y (/ a (- t x))))
(if (<= y 3.3e-77)
(- x (/ t (/ (- a z) z)))
(if (or (<= y 5.6e+30) (not (<= y 4.8e+97)))
t_1
(* t (/ (- y z) (- a z)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -1.02e+155) {
tmp = t_1;
} else if (y <= -1.9e-65) {
tmp = x + (y / (a / (t - x)));
} else if (y <= 3.3e-77) {
tmp = x - (t / ((a - z) / z));
} else if ((y <= 5.6e+30) || !(y <= 4.8e+97)) {
tmp = t_1;
} else {
tmp = t * ((y - z) / (a - 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 = (t - x) * (y / (a - z))
if (y <= (-1.02d+155)) then
tmp = t_1
else if (y <= (-1.9d-65)) then
tmp = x + (y / (a / (t - x)))
else if (y <= 3.3d-77) then
tmp = x - (t / ((a - z) / z))
else if ((y <= 5.6d+30) .or. (.not. (y <= 4.8d+97))) then
tmp = t_1
else
tmp = t * ((y - z) / (a - z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -1.02e+155) {
tmp = t_1;
} else if (y <= -1.9e-65) {
tmp = x + (y / (a / (t - x)));
} else if (y <= 3.3e-77) {
tmp = x - (t / ((a - z) / z));
} else if ((y <= 5.6e+30) || !(y <= 4.8e+97)) {
tmp = t_1;
} else {
tmp = t * ((y - z) / (a - z));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (t - x) * (y / (a - z)) tmp = 0 if y <= -1.02e+155: tmp = t_1 elif y <= -1.9e-65: tmp = x + (y / (a / (t - x))) elif y <= 3.3e-77: tmp = x - (t / ((a - z) / z)) elif (y <= 5.6e+30) or not (y <= 4.8e+97): tmp = t_1 else: tmp = t * ((y - z) / (a - z)) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(t - x) * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -1.02e+155) tmp = t_1; elseif (y <= -1.9e-65) tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); elseif (y <= 3.3e-77) tmp = Float64(x - Float64(t / Float64(Float64(a - z) / z))); elseif ((y <= 5.6e+30) || !(y <= 4.8e+97)) tmp = t_1; else tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (t - x) * (y / (a - z)); tmp = 0.0; if (y <= -1.02e+155) tmp = t_1; elseif (y <= -1.9e-65) tmp = x + (y / (a / (t - x))); elseif (y <= 3.3e-77) tmp = x - (t / ((a - z) / z)); elseif ((y <= 5.6e+30) || ~((y <= 4.8e+97))) tmp = t_1; else tmp = t * ((y - z) / (a - z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.02e+155], t$95$1, If[LessEqual[y, -1.9e-65], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.3e-77], N[(x - N[(t / N[(N[(a - z), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[y, 5.6e+30], N[Not[LessEqual[y, 4.8e+97]], $MachinePrecision]], t$95$1, N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t - x\right) \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -1.02 \cdot 10^{+155}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq -1.9 \cdot 10^{-65}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\mathbf{elif}\;y \leq 3.3 \cdot 10^{-77}:\\
\;\;\;\;x - \frac{t}{\frac{a - z}{z}}\\
\mathbf{elif}\;y \leq 5.6 \cdot 10^{+30} \lor \neg \left(y \leq 4.8 \cdot 10^{+97}\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\end{array}
\end{array}
if y < -1.02e155 or 3.29999999999999991e-77 < y < 5.59999999999999966e30 or 4.8e97 < y Initial program 82.8%
add-cube-cbrt82.9%
fma-def82.9%
pow282.9%
*-commutative82.9%
associate-*l/64.6%
associate-*r/84.2%
Applied egg-rr84.2%
Taylor expanded in y around -inf 61.4%
associate-*l/76.5%
Simplified76.5%
if -1.02e155 < y < -1.9000000000000001e-65Initial program 84.6%
Taylor expanded in z around 0 56.5%
+-commutative56.5%
associate-/l*60.2%
Simplified60.2%
if -1.9000000000000001e-65 < y < 3.29999999999999991e-77Initial program 71.4%
Taylor expanded in x around -inf 75.9%
*-un-lft-identity75.9%
fma-def75.9%
associate-/l*87.6%
mul-1-neg87.6%
*-commutative87.6%
associate--r+84.7%
div-sub84.7%
Applied egg-rr84.7%
fma-udef84.7%
*-lft-identity84.7%
associate-/r/80.7%
fma-def80.7%
*-commutative80.7%
distribute-rgt-neg-in80.7%
sub-neg80.7%
metadata-eval80.7%
Simplified80.7%
Taylor expanded in z around 0 68.2%
associate-*r*68.2%
sub-neg68.2%
metadata-eval68.2%
distribute-lft-in68.2%
metadata-eval68.2%
+-commutative68.2%
*-commutative68.2%
mul-1-neg68.2%
unsub-neg68.2%
Simplified68.2%
Taylor expanded in y around 0 59.8%
+-commutative59.8%
mul-1-neg59.8%
unsub-neg59.8%
associate-/l*70.5%
Simplified70.5%
if 5.59999999999999966e30 < y < 4.8e97Initial program 81.7%
Taylor expanded in t around inf 75.8%
div-sub75.8%
Simplified75.8%
Final simplification70.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- t x) (/ y (- a z))))
(t_2 (* t (/ (- y z) (- a z))))
(t_3 (+ x (/ y (/ a (- t x))))))
(if (<= z -1.85e+93)
t_2
(if (<= z -7.5e-150)
t_1
(if (<= z 5e-14)
t_3
(if (<= z 2.45e+58) t_1 (if (<= z 6.5e+129) t_3 t_2)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double t_2 = t * ((y - z) / (a - z));
double t_3 = x + (y / (a / (t - x)));
double tmp;
if (z <= -1.85e+93) {
tmp = t_2;
} else if (z <= -7.5e-150) {
tmp = t_1;
} else if (z <= 5e-14) {
tmp = t_3;
} else if (z <= 2.45e+58) {
tmp = t_1;
} else if (z <= 6.5e+129) {
tmp = t_3;
} 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) :: t_3
real(8) :: tmp
t_1 = (t - x) * (y / (a - z))
t_2 = t * ((y - z) / (a - z))
t_3 = x + (y / (a / (t - x)))
if (z <= (-1.85d+93)) then
tmp = t_2
else if (z <= (-7.5d-150)) then
tmp = t_1
else if (z <= 5d-14) then
tmp = t_3
else if (z <= 2.45d+58) then
tmp = t_1
else if (z <= 6.5d+129) then
tmp = t_3
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 - x) * (y / (a - z));
double t_2 = t * ((y - z) / (a - z));
double t_3 = x + (y / (a / (t - x)));
double tmp;
if (z <= -1.85e+93) {
tmp = t_2;
} else if (z <= -7.5e-150) {
tmp = t_1;
} else if (z <= 5e-14) {
tmp = t_3;
} else if (z <= 2.45e+58) {
tmp = t_1;
} else if (z <= 6.5e+129) {
tmp = t_3;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (t - x) * (y / (a - z)) t_2 = t * ((y - z) / (a - z)) t_3 = x + (y / (a / (t - x))) tmp = 0 if z <= -1.85e+93: tmp = t_2 elif z <= -7.5e-150: tmp = t_1 elif z <= 5e-14: tmp = t_3 elif z <= 2.45e+58: tmp = t_1 elif z <= 6.5e+129: tmp = t_3 else: tmp = t_2 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(t - x) * Float64(y / Float64(a - z))) t_2 = Float64(t * Float64(Float64(y - z) / Float64(a - z))) t_3 = Float64(x + Float64(y / Float64(a / Float64(t - x)))) tmp = 0.0 if (z <= -1.85e+93) tmp = t_2; elseif (z <= -7.5e-150) tmp = t_1; elseif (z <= 5e-14) tmp = t_3; elseif (z <= 2.45e+58) tmp = t_1; elseif (z <= 6.5e+129) tmp = t_3; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (t - x) * (y / (a - z)); t_2 = t * ((y - z) / (a - z)); t_3 = x + (y / (a / (t - x))); tmp = 0.0; if (z <= -1.85e+93) tmp = t_2; elseif (z <= -7.5e-150) tmp = t_1; elseif (z <= 5e-14) tmp = t_3; elseif (z <= 2.45e+58) tmp = t_1; elseif (z <= 6.5e+129) tmp = t_3; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.85e+93], t$95$2, If[LessEqual[z, -7.5e-150], t$95$1, If[LessEqual[z, 5e-14], t$95$3, If[LessEqual[z, 2.45e+58], t$95$1, If[LessEqual[z, 6.5e+129], t$95$3, t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t - x\right) \cdot \frac{y}{a - z}\\
t_2 := t \cdot \frac{y - z}{a - z}\\
t_3 := x + \frac{y}{\frac{a}{t - x}}\\
\mathbf{if}\;z \leq -1.85 \cdot 10^{+93}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;z \leq -7.5 \cdot 10^{-150}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 5 \cdot 10^{-14}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;z \leq 2.45 \cdot 10^{+58}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 6.5 \cdot 10^{+129}:\\
\;\;\;\;t_3\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
if z < -1.84999999999999994e93 or 6.4999999999999995e129 < z Initial program 56.3%
Taylor expanded in t around inf 64.9%
div-sub64.9%
Simplified64.9%
if -1.84999999999999994e93 < z < -7.5000000000000004e-150 or 5.0000000000000002e-14 < z < 2.45000000000000009e58Initial program 84.3%
add-cube-cbrt84.1%
fma-def84.1%
pow284.1%
*-commutative84.1%
associate-*l/79.8%
associate-*r/86.6%
Applied egg-rr86.6%
Taylor expanded in y around -inf 57.0%
associate-*l/60.9%
Simplified60.9%
if -7.5000000000000004e-150 < z < 5.0000000000000002e-14 or 2.45000000000000009e58 < z < 6.4999999999999995e129Initial program 92.7%
Taylor expanded in z around 0 74.2%
+-commutative74.2%
associate-/l*77.7%
Simplified77.7%
Final simplification69.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- t x) (/ y (- a z)))))
(if (<= y -1e+155)
t_1
(if (<= y -4.8e-65)
(+ x (/ y (/ a (- t x))))
(if (<= y 3.5e-77)
(- x (/ t (/ (- a z) z)))
(if (<= y 8.8e+29)
t_1
(if (<= y 1.2e+94)
(* t (/ (- y z) (- a z)))
(/ y (/ (- a z) (- t x))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -1e+155) {
tmp = t_1;
} else if (y <= -4.8e-65) {
tmp = x + (y / (a / (t - x)));
} else if (y <= 3.5e-77) {
tmp = x - (t / ((a - z) / z));
} else if (y <= 8.8e+29) {
tmp = t_1;
} else if (y <= 1.2e+94) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = y / ((a - z) / (t - 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 = (t - x) * (y / (a - z))
if (y <= (-1d+155)) then
tmp = t_1
else if (y <= (-4.8d-65)) then
tmp = x + (y / (a / (t - x)))
else if (y <= 3.5d-77) then
tmp = x - (t / ((a - z) / z))
else if (y <= 8.8d+29) then
tmp = t_1
else if (y <= 1.2d+94) then
tmp = t * ((y - z) / (a - z))
else
tmp = y / ((a - z) / (t - x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -1e+155) {
tmp = t_1;
} else if (y <= -4.8e-65) {
tmp = x + (y / (a / (t - x)));
} else if (y <= 3.5e-77) {
tmp = x - (t / ((a - z) / z));
} else if (y <= 8.8e+29) {
tmp = t_1;
} else if (y <= 1.2e+94) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = y / ((a - z) / (t - x));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (t - x) * (y / (a - z)) tmp = 0 if y <= -1e+155: tmp = t_1 elif y <= -4.8e-65: tmp = x + (y / (a / (t - x))) elif y <= 3.5e-77: tmp = x - (t / ((a - z) / z)) elif y <= 8.8e+29: tmp = t_1 elif y <= 1.2e+94: tmp = t * ((y - z) / (a - z)) else: tmp = y / ((a - z) / (t - x)) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(t - x) * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -1e+155) tmp = t_1; elseif (y <= -4.8e-65) tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); elseif (y <= 3.5e-77) tmp = Float64(x - Float64(t / Float64(Float64(a - z) / z))); elseif (y <= 8.8e+29) tmp = t_1; elseif (y <= 1.2e+94) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(y / Float64(Float64(a - z) / Float64(t - x))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (t - x) * (y / (a - z)); tmp = 0.0; if (y <= -1e+155) tmp = t_1; elseif (y <= -4.8e-65) tmp = x + (y / (a / (t - x))); elseif (y <= 3.5e-77) tmp = x - (t / ((a - z) / z)); elseif (y <= 8.8e+29) tmp = t_1; elseif (y <= 1.2e+94) tmp = t * ((y - z) / (a - z)); else tmp = y / ((a - z) / (t - x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1e+155], t$95$1, If[LessEqual[y, -4.8e-65], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.5e-77], N[(x - N[(t / N[(N[(a - z), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.8e+29], t$95$1, If[LessEqual[y, 1.2e+94], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y / N[(N[(a - z), $MachinePrecision] / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t - x\right) \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -1 \cdot 10^{+155}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq -4.8 \cdot 10^{-65}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{-77}:\\
\;\;\;\;x - \frac{t}{\frac{a - z}{z}}\\
\mathbf{elif}\;y \leq 8.8 \cdot 10^{+29}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{+94}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{\frac{a - z}{t - x}}\\
\end{array}
\end{array}
if y < -1.00000000000000001e155 or 3.50000000000000013e-77 < y < 8.8000000000000005e29Initial program 79.0%
add-cube-cbrt79.0%
fma-def79.0%
pow279.0%
*-commutative79.0%
associate-*l/63.1%
associate-*r/79.7%
Applied egg-rr79.7%
Taylor expanded in y around -inf 63.0%
associate-*l/78.7%
Simplified78.7%
if -1.00000000000000001e155 < y < -4.8000000000000003e-65Initial program 84.6%
Taylor expanded in z around 0 56.5%
+-commutative56.5%
associate-/l*60.2%
Simplified60.2%
if -4.8000000000000003e-65 < y < 3.50000000000000013e-77Initial program 71.4%
Taylor expanded in x around -inf 75.9%
*-un-lft-identity75.9%
fma-def75.9%
associate-/l*87.6%
mul-1-neg87.6%
*-commutative87.6%
associate--r+84.7%
div-sub84.7%
Applied egg-rr84.7%
fma-udef84.7%
*-lft-identity84.7%
associate-/r/80.7%
fma-def80.7%
*-commutative80.7%
distribute-rgt-neg-in80.7%
sub-neg80.7%
metadata-eval80.7%
Simplified80.7%
Taylor expanded in z around 0 68.2%
associate-*r*68.2%
sub-neg68.2%
metadata-eval68.2%
distribute-lft-in68.2%
metadata-eval68.2%
+-commutative68.2%
*-commutative68.2%
mul-1-neg68.2%
unsub-neg68.2%
Simplified68.2%
Taylor expanded in y around 0 59.8%
+-commutative59.8%
mul-1-neg59.8%
unsub-neg59.8%
associate-/l*70.5%
Simplified70.5%
if 8.8000000000000005e29 < y < 1.19999999999999991e94Initial program 81.7%
Taylor expanded in t around inf 75.8%
div-sub75.8%
Simplified75.8%
if 1.19999999999999991e94 < y Initial program 86.0%
Taylor expanded in y around inf 76.5%
div-sub76.5%
*-commutative76.5%
associate-*r/60.1%
associate-/l*76.5%
Simplified76.5%
Final simplification71.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- t x) (/ y (- a z)))))
(if (<= y -2.7e+155)
t_1
(if (<= y -2.5e-65)
(+ x (* (- z y) (/ (- x t) a)))
(if (<= y 5e-77)
(- x (/ t (/ (- a z) z)))
(if (<= y 3.4e+30)
t_1
(if (<= y 2.75e+94)
(* t (/ (- y z) (- a z)))
(/ y (/ (- a z) (- t x))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -2.7e+155) {
tmp = t_1;
} else if (y <= -2.5e-65) {
tmp = x + ((z - y) * ((x - t) / a));
} else if (y <= 5e-77) {
tmp = x - (t / ((a - z) / z));
} else if (y <= 3.4e+30) {
tmp = t_1;
} else if (y <= 2.75e+94) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = y / ((a - z) / (t - 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 = (t - x) * (y / (a - z))
if (y <= (-2.7d+155)) then
tmp = t_1
else if (y <= (-2.5d-65)) then
tmp = x + ((z - y) * ((x - t) / a))
else if (y <= 5d-77) then
tmp = x - (t / ((a - z) / z))
else if (y <= 3.4d+30) then
tmp = t_1
else if (y <= 2.75d+94) then
tmp = t * ((y - z) / (a - z))
else
tmp = y / ((a - z) / (t - x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (t - x) * (y / (a - z));
double tmp;
if (y <= -2.7e+155) {
tmp = t_1;
} else if (y <= -2.5e-65) {
tmp = x + ((z - y) * ((x - t) / a));
} else if (y <= 5e-77) {
tmp = x - (t / ((a - z) / z));
} else if (y <= 3.4e+30) {
tmp = t_1;
} else if (y <= 2.75e+94) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = y / ((a - z) / (t - x));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (t - x) * (y / (a - z)) tmp = 0 if y <= -2.7e+155: tmp = t_1 elif y <= -2.5e-65: tmp = x + ((z - y) * ((x - t) / a)) elif y <= 5e-77: tmp = x - (t / ((a - z) / z)) elif y <= 3.4e+30: tmp = t_1 elif y <= 2.75e+94: tmp = t * ((y - z) / (a - z)) else: tmp = y / ((a - z) / (t - x)) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(t - x) * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -2.7e+155) tmp = t_1; elseif (y <= -2.5e-65) tmp = Float64(x + Float64(Float64(z - y) * Float64(Float64(x - t) / a))); elseif (y <= 5e-77) tmp = Float64(x - Float64(t / Float64(Float64(a - z) / z))); elseif (y <= 3.4e+30) tmp = t_1; elseif (y <= 2.75e+94) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(y / Float64(Float64(a - z) / Float64(t - x))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (t - x) * (y / (a - z)); tmp = 0.0; if (y <= -2.7e+155) tmp = t_1; elseif (y <= -2.5e-65) tmp = x + ((z - y) * ((x - t) / a)); elseif (y <= 5e-77) tmp = x - (t / ((a - z) / z)); elseif (y <= 3.4e+30) tmp = t_1; elseif (y <= 2.75e+94) tmp = t * ((y - z) / (a - z)); else tmp = y / ((a - z) / (t - x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.7e+155], t$95$1, If[LessEqual[y, -2.5e-65], N[(x + N[(N[(z - y), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e-77], N[(x - N[(t / N[(N[(a - z), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.4e+30], t$95$1, If[LessEqual[y, 2.75e+94], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y / N[(N[(a - z), $MachinePrecision] / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t - x\right) \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -2.7 \cdot 10^{+155}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq -2.5 \cdot 10^{-65}:\\
\;\;\;\;x + \left(z - y\right) \cdot \frac{x - t}{a}\\
\mathbf{elif}\;y \leq 5 \cdot 10^{-77}:\\
\;\;\;\;x - \frac{t}{\frac{a - z}{z}}\\
\mathbf{elif}\;y \leq 3.4 \cdot 10^{+30}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq 2.75 \cdot 10^{+94}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{\frac{a - z}{t - x}}\\
\end{array}
\end{array}
if y < -2.69999999999999994e155 or 4.99999999999999963e-77 < y < 3.4000000000000002e30Initial program 79.0%
add-cube-cbrt79.0%
fma-def79.0%
pow279.0%
*-commutative79.0%
associate-*l/63.1%
associate-*r/79.7%
Applied egg-rr79.7%
Taylor expanded in y around -inf 63.0%
associate-*l/78.7%
Simplified78.7%
if -2.69999999999999994e155 < y < -2.49999999999999991e-65Initial program 84.6%
Taylor expanded in a around inf 61.9%
if -2.49999999999999991e-65 < y < 4.99999999999999963e-77Initial program 71.4%
Taylor expanded in x around -inf 75.9%
*-un-lft-identity75.9%
fma-def75.9%
associate-/l*87.6%
mul-1-neg87.6%
*-commutative87.6%
associate--r+84.7%
div-sub84.7%
Applied egg-rr84.7%
fma-udef84.7%
*-lft-identity84.7%
associate-/r/80.7%
fma-def80.7%
*-commutative80.7%
distribute-rgt-neg-in80.7%
sub-neg80.7%
metadata-eval80.7%
Simplified80.7%
Taylor expanded in z around 0 68.2%
associate-*r*68.2%
sub-neg68.2%
metadata-eval68.2%
distribute-lft-in68.2%
metadata-eval68.2%
+-commutative68.2%
*-commutative68.2%
mul-1-neg68.2%
unsub-neg68.2%
Simplified68.2%
Taylor expanded in y around 0 59.8%
+-commutative59.8%
mul-1-neg59.8%
unsub-neg59.8%
associate-/l*70.5%
Simplified70.5%
if 3.4000000000000002e30 < y < 2.7499999999999999e94Initial program 81.7%
Taylor expanded in t around inf 75.8%
div-sub75.8%
Simplified75.8%
if 2.7499999999999999e94 < y Initial program 86.0%
Taylor expanded in y around inf 76.5%
div-sub76.5%
*-commutative76.5%
associate-*r/60.1%
associate-/l*76.5%
Simplified76.5%
Final simplification71.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ y (/ a x)))))
(if (<= a -4.7e+60)
t_1
(if (<= a -2.45e-231)
(* t (- 1.0 (/ y z)))
(if (<= a 1.06e-149)
(/ (* y (- x t)) z)
(if (<= a 1.1e-119)
t
(if (<= a 6.8e-63) (* y (/ (- x t) z)) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -4.7e+60) {
tmp = t_1;
} else if (a <= -2.45e-231) {
tmp = t * (1.0 - (y / z));
} else if (a <= 1.06e-149) {
tmp = (y * (x - t)) / z;
} else if (a <= 1.1e-119) {
tmp = t;
} else if (a <= 6.8e-63) {
tmp = y * ((x - t) / z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y / (a / x))
if (a <= (-4.7d+60)) then
tmp = t_1
else if (a <= (-2.45d-231)) then
tmp = t * (1.0d0 - (y / z))
else if (a <= 1.06d-149) then
tmp = (y * (x - t)) / z
else if (a <= 1.1d-119) then
tmp = t
else if (a <= 6.8d-63) then
tmp = y * ((x - t) / z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -4.7e+60) {
tmp = t_1;
} else if (a <= -2.45e-231) {
tmp = t * (1.0 - (y / z));
} else if (a <= 1.06e-149) {
tmp = (y * (x - t)) / z;
} else if (a <= 1.1e-119) {
tmp = t;
} else if (a <= 6.8e-63) {
tmp = y * ((x - t) / z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y / (a / x)) tmp = 0 if a <= -4.7e+60: tmp = t_1 elif a <= -2.45e-231: tmp = t * (1.0 - (y / z)) elif a <= 1.06e-149: tmp = (y * (x - t)) / z elif a <= 1.1e-119: tmp = t elif a <= 6.8e-63: tmp = y * ((x - t) / z) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y / Float64(a / x))) tmp = 0.0 if (a <= -4.7e+60) tmp = t_1; elseif (a <= -2.45e-231) tmp = Float64(t * Float64(1.0 - Float64(y / z))); elseif (a <= 1.06e-149) tmp = Float64(Float64(y * Float64(x - t)) / z); elseif (a <= 1.1e-119) tmp = t; elseif (a <= 6.8e-63) tmp = Float64(y * Float64(Float64(x - t) / z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y / (a / x)); tmp = 0.0; if (a <= -4.7e+60) tmp = t_1; elseif (a <= -2.45e-231) tmp = t * (1.0 - (y / z)); elseif (a <= 1.06e-149) tmp = (y * (x - t)) / z; elseif (a <= 1.1e-119) tmp = t; elseif (a <= 6.8e-63) tmp = y * ((x - t) / z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -4.7e+60], t$95$1, If[LessEqual[a, -2.45e-231], N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.06e-149], N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[a, 1.1e-119], t, If[LessEqual[a, 6.8e-63], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{y}{\frac{a}{x}}\\
\mathbf{if}\;a \leq -4.7 \cdot 10^{+60}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -2.45 \cdot 10^{-231}:\\
\;\;\;\;t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{elif}\;a \leq 1.06 \cdot 10^{-149}:\\
\;\;\;\;\frac{y \cdot \left(x - t\right)}{z}\\
\mathbf{elif}\;a \leq 1.1 \cdot 10^{-119}:\\
\;\;\;\;t\\
\mathbf{elif}\;a \leq 6.8 \cdot 10^{-63}:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -4.6999999999999998e60 or 6.79999999999999997e-63 < a Initial program 86.2%
Taylor expanded in z around 0 65.1%
Taylor expanded in t around 0 52.9%
+-commutative52.9%
mul-1-neg52.9%
unsub-neg52.9%
associate-/l*57.6%
Simplified57.6%
if -4.6999999999999998e60 < a < -2.45000000000000002e-231Initial program 76.0%
Taylor expanded in a around 0 50.5%
+-commutative50.5%
mul-1-neg50.5%
unsub-neg50.5%
associate-/l*56.8%
Simplified56.8%
Taylor expanded in t around inf 52.9%
if -2.45000000000000002e-231 < a < 1.05999999999999998e-149Initial program 72.5%
Taylor expanded in y around inf 61.9%
div-sub61.9%
*-commutative61.9%
associate-*r/66.8%
associate-/l*62.3%
Simplified62.3%
Taylor expanded in a around 0 63.7%
mul-1-neg63.7%
associate-/l*60.2%
associate-/r/58.9%
distribute-rgt-neg-in58.9%
Simplified58.9%
Taylor expanded in y around 0 63.7%
if 1.05999999999999998e-149 < a < 1.1e-119Initial program 63.9%
Taylor expanded in z around inf 63.4%
if 1.1e-119 < a < 6.79999999999999997e-63Initial program 59.4%
Taylor expanded in a around 0 34.9%
associate-*r/34.9%
neg-mul-134.9%
Simplified34.9%
Taylor expanded in y around inf 45.6%
div-sub45.6%
Simplified45.6%
Final simplification57.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ y (/ a x)))))
(if (<= a -4.7e+60)
t_1
(if (<= a -6.6e-233)
(* t (- 1.0 (/ y z)))
(if (<= a 8.5e-152)
(/ (* y (- x t)) z)
(if (<= a 9.2e-120)
t
(if (<= a 1.22e-63) (* (/ y z) (- x t)) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -4.7e+60) {
tmp = t_1;
} else if (a <= -6.6e-233) {
tmp = t * (1.0 - (y / z));
} else if (a <= 8.5e-152) {
tmp = (y * (x - t)) / z;
} else if (a <= 9.2e-120) {
tmp = t;
} else if (a <= 1.22e-63) {
tmp = (y / z) * (x - t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y / (a / x))
if (a <= (-4.7d+60)) then
tmp = t_1
else if (a <= (-6.6d-233)) then
tmp = t * (1.0d0 - (y / z))
else if (a <= 8.5d-152) then
tmp = (y * (x - t)) / z
else if (a <= 9.2d-120) then
tmp = t
else if (a <= 1.22d-63) then
tmp = (y / z) * (x - t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -4.7e+60) {
tmp = t_1;
} else if (a <= -6.6e-233) {
tmp = t * (1.0 - (y / z));
} else if (a <= 8.5e-152) {
tmp = (y * (x - t)) / z;
} else if (a <= 9.2e-120) {
tmp = t;
} else if (a <= 1.22e-63) {
tmp = (y / z) * (x - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y / (a / x)) tmp = 0 if a <= -4.7e+60: tmp = t_1 elif a <= -6.6e-233: tmp = t * (1.0 - (y / z)) elif a <= 8.5e-152: tmp = (y * (x - t)) / z elif a <= 9.2e-120: tmp = t elif a <= 1.22e-63: tmp = (y / z) * (x - t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y / Float64(a / x))) tmp = 0.0 if (a <= -4.7e+60) tmp = t_1; elseif (a <= -6.6e-233) tmp = Float64(t * Float64(1.0 - Float64(y / z))); elseif (a <= 8.5e-152) tmp = Float64(Float64(y * Float64(x - t)) / z); elseif (a <= 9.2e-120) tmp = t; elseif (a <= 1.22e-63) tmp = Float64(Float64(y / z) * Float64(x - t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y / (a / x)); tmp = 0.0; if (a <= -4.7e+60) tmp = t_1; elseif (a <= -6.6e-233) tmp = t * (1.0 - (y / z)); elseif (a <= 8.5e-152) tmp = (y * (x - t)) / z; elseif (a <= 9.2e-120) tmp = t; elseif (a <= 1.22e-63) tmp = (y / z) * (x - t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -4.7e+60], t$95$1, If[LessEqual[a, -6.6e-233], N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 8.5e-152], N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[a, 9.2e-120], t, If[LessEqual[a, 1.22e-63], N[(N[(y / z), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{y}{\frac{a}{x}}\\
\mathbf{if}\;a \leq -4.7 \cdot 10^{+60}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -6.6 \cdot 10^{-233}:\\
\;\;\;\;t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{elif}\;a \leq 8.5 \cdot 10^{-152}:\\
\;\;\;\;\frac{y \cdot \left(x - t\right)}{z}\\
\mathbf{elif}\;a \leq 9.2 \cdot 10^{-120}:\\
\;\;\;\;t\\
\mathbf{elif}\;a \leq 1.22 \cdot 10^{-63}:\\
\;\;\;\;\frac{y}{z} \cdot \left(x - t\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -4.6999999999999998e60 or 1.2199999999999999e-63 < a Initial program 86.2%
Taylor expanded in z around 0 65.1%
Taylor expanded in t around 0 52.9%
+-commutative52.9%
mul-1-neg52.9%
unsub-neg52.9%
associate-/l*57.6%
Simplified57.6%
if -4.6999999999999998e60 < a < -6.6000000000000001e-233Initial program 76.0%
Taylor expanded in a around 0 50.5%
+-commutative50.5%
mul-1-neg50.5%
unsub-neg50.5%
associate-/l*56.8%
Simplified56.8%
Taylor expanded in t around inf 52.9%
if -6.6000000000000001e-233 < a < 8.5000000000000007e-152Initial program 72.5%
Taylor expanded in y around inf 61.9%
div-sub61.9%
*-commutative61.9%
associate-*r/66.8%
associate-/l*62.3%
Simplified62.3%
Taylor expanded in a around 0 63.7%
mul-1-neg63.7%
associate-/l*60.2%
associate-/r/58.9%
distribute-rgt-neg-in58.9%
Simplified58.9%
Taylor expanded in y around 0 63.7%
if 8.5000000000000007e-152 < a < 9.19999999999999946e-120Initial program 63.9%
Taylor expanded in z around inf 63.4%
if 9.19999999999999946e-120 < a < 1.2199999999999999e-63Initial program 59.4%
Taylor expanded in y around inf 58.3%
div-sub58.3%
*-commutative58.3%
associate-*r/47.1%
associate-/l*58.3%
Simplified58.3%
Taylor expanded in a around 0 34.2%
mul-1-neg34.2%
associate-/l*45.6%
associate-/r/45.6%
distribute-rgt-neg-in45.6%
Simplified45.6%
Final simplification57.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ y (/ a x)))))
(if (<= a -3.15e+59)
t_1
(if (<= a -1.8e-233)
(* t (- 1.0 (/ y z)))
(if (<= a 7.6e-150)
(/ (* y (- x t)) z)
(if (<= a 8e-120)
(/ (- t) (/ (- a z) z))
(if (<= a 6.5e-64) (* (/ y z) (- x t)) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -3.15e+59) {
tmp = t_1;
} else if (a <= -1.8e-233) {
tmp = t * (1.0 - (y / z));
} else if (a <= 7.6e-150) {
tmp = (y * (x - t)) / z;
} else if (a <= 8e-120) {
tmp = -t / ((a - z) / z);
} else if (a <= 6.5e-64) {
tmp = (y / z) * (x - t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y / (a / x))
if (a <= (-3.15d+59)) then
tmp = t_1
else if (a <= (-1.8d-233)) then
tmp = t * (1.0d0 - (y / z))
else if (a <= 7.6d-150) then
tmp = (y * (x - t)) / z
else if (a <= 8d-120) then
tmp = -t / ((a - z) / z)
else if (a <= 6.5d-64) then
tmp = (y / z) * (x - t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -3.15e+59) {
tmp = t_1;
} else if (a <= -1.8e-233) {
tmp = t * (1.0 - (y / z));
} else if (a <= 7.6e-150) {
tmp = (y * (x - t)) / z;
} else if (a <= 8e-120) {
tmp = -t / ((a - z) / z);
} else if (a <= 6.5e-64) {
tmp = (y / z) * (x - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y / (a / x)) tmp = 0 if a <= -3.15e+59: tmp = t_1 elif a <= -1.8e-233: tmp = t * (1.0 - (y / z)) elif a <= 7.6e-150: tmp = (y * (x - t)) / z elif a <= 8e-120: tmp = -t / ((a - z) / z) elif a <= 6.5e-64: tmp = (y / z) * (x - t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y / Float64(a / x))) tmp = 0.0 if (a <= -3.15e+59) tmp = t_1; elseif (a <= -1.8e-233) tmp = Float64(t * Float64(1.0 - Float64(y / z))); elseif (a <= 7.6e-150) tmp = Float64(Float64(y * Float64(x - t)) / z); elseif (a <= 8e-120) tmp = Float64(Float64(-t) / Float64(Float64(a - z) / z)); elseif (a <= 6.5e-64) tmp = Float64(Float64(y / z) * Float64(x - t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y / (a / x)); tmp = 0.0; if (a <= -3.15e+59) tmp = t_1; elseif (a <= -1.8e-233) tmp = t * (1.0 - (y / z)); elseif (a <= 7.6e-150) tmp = (y * (x - t)) / z; elseif (a <= 8e-120) tmp = -t / ((a - z) / z); elseif (a <= 6.5e-64) tmp = (y / z) * (x - t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3.15e+59], t$95$1, If[LessEqual[a, -1.8e-233], N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 7.6e-150], N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[a, 8e-120], N[((-t) / N[(N[(a - z), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 6.5e-64], N[(N[(y / z), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{y}{\frac{a}{x}}\\
\mathbf{if}\;a \leq -3.15 \cdot 10^{+59}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -1.8 \cdot 10^{-233}:\\
\;\;\;\;t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{elif}\;a \leq 7.6 \cdot 10^{-150}:\\
\;\;\;\;\frac{y \cdot \left(x - t\right)}{z}\\
\mathbf{elif}\;a \leq 8 \cdot 10^{-120}:\\
\;\;\;\;\frac{-t}{\frac{a - z}{z}}\\
\mathbf{elif}\;a \leq 6.5 \cdot 10^{-64}:\\
\;\;\;\;\frac{y}{z} \cdot \left(x - t\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -3.15e59 or 6.5000000000000004e-64 < a Initial program 86.2%
Taylor expanded in z around 0 65.1%
Taylor expanded in t around 0 52.9%
+-commutative52.9%
mul-1-neg52.9%
unsub-neg52.9%
associate-/l*57.6%
Simplified57.6%
if -3.15e59 < a < -1.80000000000000004e-233Initial program 76.0%
Taylor expanded in a around 0 50.5%
+-commutative50.5%
mul-1-neg50.5%
unsub-neg50.5%
associate-/l*56.8%
Simplified56.8%
Taylor expanded in t around inf 52.9%
if -1.80000000000000004e-233 < a < 7.5999999999999997e-150Initial program 72.5%
Taylor expanded in y around inf 61.9%
div-sub61.9%
*-commutative61.9%
associate-*r/66.8%
associate-/l*62.3%
Simplified62.3%
Taylor expanded in a around 0 63.7%
mul-1-neg63.7%
associate-/l*60.2%
associate-/r/58.9%
distribute-rgt-neg-in58.9%
Simplified58.9%
Taylor expanded in y around 0 63.7%
if 7.5999999999999997e-150 < a < 7.99999999999999983e-120Initial program 63.9%
add-cube-cbrt64.2%
fma-def64.2%
pow264.2%
*-commutative64.2%
associate-*l/51.7%
associate-*r/76.2%
Applied egg-rr76.2%
Taylor expanded in x around 0 64.7%
associate-/l*88.3%
associate-/r/75.9%
Simplified75.9%
Taylor expanded in y around 0 40.5%
mul-1-neg40.5%
associate-/l*64.1%
distribute-neg-frac64.1%
Simplified64.1%
if 7.99999999999999983e-120 < a < 6.5000000000000004e-64Initial program 59.4%
Taylor expanded in y around inf 58.3%
div-sub58.3%
*-commutative58.3%
associate-*r/47.1%
associate-/l*58.3%
Simplified58.3%
Taylor expanded in a around 0 34.2%
mul-1-neg34.2%
associate-/l*45.6%
associate-/r/45.6%
distribute-rgt-neg-in45.6%
Simplified45.6%
Final simplification57.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ y (/ a x)))) (t_2 (* t (/ (- y z) (- a z)))))
(if (<= t -3600.0)
t_2
(if (<= t 2.9e-260)
t_1
(if (<= t 4.8e-162) (* (/ y z) (- x t)) (if (<= t 2.7e-19) t_1 t_2))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double t_2 = t * ((y - z) / (a - z));
double tmp;
if (t <= -3600.0) {
tmp = t_2;
} else if (t <= 2.9e-260) {
tmp = t_1;
} else if (t <= 4.8e-162) {
tmp = (y / z) * (x - t);
} else if (t <= 2.7e-19) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = x - (y / (a / x))
t_2 = t * ((y - z) / (a - z))
if (t <= (-3600.0d0)) then
tmp = t_2
else if (t <= 2.9d-260) then
tmp = t_1
else if (t <= 4.8d-162) then
tmp = (y / z) * (x - t)
else if (t <= 2.7d-19) then
tmp = t_1
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double t_2 = t * ((y - z) / (a - z));
double tmp;
if (t <= -3600.0) {
tmp = t_2;
} else if (t <= 2.9e-260) {
tmp = t_1;
} else if (t <= 4.8e-162) {
tmp = (y / z) * (x - t);
} else if (t <= 2.7e-19) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y / (a / x)) t_2 = t * ((y - z) / (a - z)) tmp = 0 if t <= -3600.0: tmp = t_2 elif t <= 2.9e-260: tmp = t_1 elif t <= 4.8e-162: tmp = (y / z) * (x - t) elif t <= 2.7e-19: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y / Float64(a / x))) t_2 = Float64(t * Float64(Float64(y - z) / Float64(a - z))) tmp = 0.0 if (t <= -3600.0) tmp = t_2; elseif (t <= 2.9e-260) tmp = t_1; elseif (t <= 4.8e-162) tmp = Float64(Float64(y / z) * Float64(x - t)); elseif (t <= 2.7e-19) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y / (a / x)); t_2 = t * ((y - z) / (a - z)); tmp = 0.0; if (t <= -3600.0) tmp = t_2; elseif (t <= 2.9e-260) tmp = t_1; elseif (t <= 4.8e-162) tmp = (y / z) * (x - t); elseif (t <= 2.7e-19) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3600.0], t$95$2, If[LessEqual[t, 2.9e-260], t$95$1, If[LessEqual[t, 4.8e-162], N[(N[(y / z), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.7e-19], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{y}{\frac{a}{x}}\\
t_2 := t \cdot \frac{y - z}{a - z}\\
\mathbf{if}\;t \leq -3600:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 2.9 \cdot 10^{-260}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 4.8 \cdot 10^{-162}:\\
\;\;\;\;\frac{y}{z} \cdot \left(x - t\right)\\
\mathbf{elif}\;t \leq 2.7 \cdot 10^{-19}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
if t < -3600 or 2.7000000000000001e-19 < t Initial program 86.3%
Taylor expanded in t around inf 77.7%
div-sub77.7%
Simplified77.7%
if -3600 < t < 2.8999999999999999e-260 or 4.8000000000000004e-162 < t < 2.7000000000000001e-19Initial program 73.4%
Taylor expanded in z around 0 51.0%
Taylor expanded in t around 0 48.6%
+-commutative48.6%
mul-1-neg48.6%
unsub-neg48.6%
associate-/l*51.2%
Simplified51.2%
if 2.8999999999999999e-260 < t < 4.8000000000000004e-162Initial program 68.1%
Taylor expanded in y around inf 65.9%
div-sub65.9%
*-commutative65.9%
associate-*r/68.4%
associate-/l*66.6%
Simplified66.6%
Taylor expanded in a around 0 56.5%
mul-1-neg56.5%
associate-/l*53.1%
associate-/r/60.3%
distribute-rgt-neg-in60.3%
Simplified60.3%
Final simplification64.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* (- z y) (/ (- x t) a)))))
(if (<= a -1.8e+47)
t_1
(if (<= a -7.2e-80)
(/ y (/ (- a z) (- t x)))
(if (<= a 1.55e-30) (+ t (/ (* (- y a) (- x t)) z)) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((z - y) * ((x - t) / a));
double tmp;
if (a <= -1.8e+47) {
tmp = t_1;
} else if (a <= -7.2e-80) {
tmp = y / ((a - z) / (t - x));
} else if (a <= 1.55e-30) {
tmp = t + (((y - a) * (x - t)) / z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z - y) * ((x - t) / a))
if (a <= (-1.8d+47)) then
tmp = t_1
else if (a <= (-7.2d-80)) then
tmp = y / ((a - z) / (t - x))
else if (a <= 1.55d-30) then
tmp = t + (((y - a) * (x - t)) / z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((z - y) * ((x - t) / a));
double tmp;
if (a <= -1.8e+47) {
tmp = t_1;
} else if (a <= -7.2e-80) {
tmp = y / ((a - z) / (t - x));
} else if (a <= 1.55e-30) {
tmp = t + (((y - a) * (x - t)) / z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + ((z - y) * ((x - t) / a)) tmp = 0 if a <= -1.8e+47: tmp = t_1 elif a <= -7.2e-80: tmp = y / ((a - z) / (t - x)) elif a <= 1.55e-30: tmp = t + (((y - a) * (x - t)) / z) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(z - y) * Float64(Float64(x - t) / a))) tmp = 0.0 if (a <= -1.8e+47) tmp = t_1; elseif (a <= -7.2e-80) tmp = Float64(y / Float64(Float64(a - z) / Float64(t - x))); elseif (a <= 1.55e-30) tmp = Float64(t + Float64(Float64(Float64(y - a) * Float64(x - t)) / z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + ((z - y) * ((x - t) / a)); tmp = 0.0; if (a <= -1.8e+47) tmp = t_1; elseif (a <= -7.2e-80) tmp = y / ((a - z) / (t - x)); elseif (a <= 1.55e-30) tmp = t + (((y - a) * (x - t)) / z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(z - y), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.8e+47], t$95$1, If[LessEqual[a, -7.2e-80], N[(y / N[(N[(a - z), $MachinePrecision] / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.55e-30], N[(t + N[(N[(N[(y - a), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(z - y\right) \cdot \frac{x - t}{a}\\
\mathbf{if}\;a \leq -1.8 \cdot 10^{+47}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -7.2 \cdot 10^{-80}:\\
\;\;\;\;\frac{y}{\frac{a - z}{t - x}}\\
\mathbf{elif}\;a \leq 1.55 \cdot 10^{-30}:\\
\;\;\;\;t + \frac{\left(y - a\right) \cdot \left(x - t\right)}{z}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -1.80000000000000004e47 or 1.54999999999999995e-30 < a Initial program 87.5%
Taylor expanded in a around inf 77.9%
if -1.80000000000000004e47 < a < -7.2e-80Initial program 79.5%
Taylor expanded in y around inf 67.2%
div-sub67.2%
*-commutative67.2%
associate-*r/51.8%
associate-/l*67.3%
Simplified67.3%
if -7.2e-80 < a < 1.54999999999999995e-30Initial program 69.6%
Taylor expanded in z around inf 78.4%
+-commutative78.4%
associate--l+78.4%
associate-*r/78.4%
associate-*r/78.4%
div-sub79.3%
distribute-lft-out--79.3%
mul-1-neg79.3%
distribute-neg-frac79.3%
unsub-neg79.3%
distribute-rgt-out--79.3%
Simplified79.3%
Final simplification77.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (- 1.0 (/ y z)))))
(if (<= a -4.5e+59)
x
(if (<= a -6.2e-268)
t_1
(if (<= a 3.8e-51) (* y (/ (- x t) z)) (if (<= a 2.6e+38) t_1 x))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (1.0 - (y / z));
double tmp;
if (a <= -4.5e+59) {
tmp = x;
} else if (a <= -6.2e-268) {
tmp = t_1;
} else if (a <= 3.8e-51) {
tmp = y * ((x - t) / z);
} else if (a <= 2.6e+38) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = t * (1.0d0 - (y / z))
if (a <= (-4.5d+59)) then
tmp = x
else if (a <= (-6.2d-268)) then
tmp = t_1
else if (a <= 3.8d-51) then
tmp = y * ((x - t) / z)
else if (a <= 2.6d+38) then
tmp = t_1
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = t * (1.0 - (y / z));
double tmp;
if (a <= -4.5e+59) {
tmp = x;
} else if (a <= -6.2e-268) {
tmp = t_1;
} else if (a <= 3.8e-51) {
tmp = y * ((x - t) / z);
} else if (a <= 2.6e+38) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (1.0 - (y / z)) tmp = 0 if a <= -4.5e+59: tmp = x elif a <= -6.2e-268: tmp = t_1 elif a <= 3.8e-51: tmp = y * ((x - t) / z) elif a <= 2.6e+38: tmp = t_1 else: tmp = x return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(1.0 - Float64(y / z))) tmp = 0.0 if (a <= -4.5e+59) tmp = x; elseif (a <= -6.2e-268) tmp = t_1; elseif (a <= 3.8e-51) tmp = Float64(y * Float64(Float64(x - t) / z)); elseif (a <= 2.6e+38) tmp = t_1; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (1.0 - (y / z)); tmp = 0.0; if (a <= -4.5e+59) tmp = x; elseif (a <= -6.2e-268) tmp = t_1; elseif (a <= 3.8e-51) tmp = y * ((x - t) / z); elseif (a <= 2.6e+38) tmp = t_1; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -4.5e+59], x, If[LessEqual[a, -6.2e-268], t$95$1, If[LessEqual[a, 3.8e-51], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.6e+38], t$95$1, x]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{if}\;a \leq -4.5 \cdot 10^{+59}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -6.2 \cdot 10^{-268}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq 3.8 \cdot 10^{-51}:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{elif}\;a \leq 2.6 \cdot 10^{+38}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -4.49999999999999959e59 or 2.5999999999999999e38 < a Initial program 89.6%
Taylor expanded in a around inf 52.5%
if -4.49999999999999959e59 < a < -6.1999999999999996e-268 or 3.80000000000000003e-51 < a < 2.5999999999999999e38Initial program 75.3%
Taylor expanded in a around 0 44.7%
+-commutative44.7%
mul-1-neg44.7%
unsub-neg44.7%
associate-/l*48.9%
Simplified48.9%
Taylor expanded in t around inf 50.0%
if -6.1999999999999996e-268 < a < 3.80000000000000003e-51Initial program 68.1%
Taylor expanded in a around 0 53.3%
associate-*r/53.3%
neg-mul-153.3%
Simplified53.3%
Taylor expanded in y around inf 52.2%
div-sub52.2%
Simplified52.2%
Final simplification51.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ y (/ a x)))))
(if (<= a -1.3e+60)
t_1
(if (<= a -5.2e-268)
(* t (- 1.0 (/ y z)))
(if (<= a 2.3e-62) (* y (/ (- x t) z)) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -1.3e+60) {
tmp = t_1;
} else if (a <= -5.2e-268) {
tmp = t * (1.0 - (y / z));
} else if (a <= 2.3e-62) {
tmp = y * ((x - t) / z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (y / (a / x))
if (a <= (-1.3d+60)) then
tmp = t_1
else if (a <= (-5.2d-268)) then
tmp = t * (1.0d0 - (y / z))
else if (a <= 2.3d-62) then
tmp = y * ((x - t) / z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y / (a / x));
double tmp;
if (a <= -1.3e+60) {
tmp = t_1;
} else if (a <= -5.2e-268) {
tmp = t * (1.0 - (y / z));
} else if (a <= 2.3e-62) {
tmp = y * ((x - t) / z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y / (a / x)) tmp = 0 if a <= -1.3e+60: tmp = t_1 elif a <= -5.2e-268: tmp = t * (1.0 - (y / z)) elif a <= 2.3e-62: tmp = y * ((x - t) / z) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y / Float64(a / x))) tmp = 0.0 if (a <= -1.3e+60) tmp = t_1; elseif (a <= -5.2e-268) tmp = Float64(t * Float64(1.0 - Float64(y / z))); elseif (a <= 2.3e-62) tmp = Float64(y * Float64(Float64(x - t) / z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y / (a / x)); tmp = 0.0; if (a <= -1.3e+60) tmp = t_1; elseif (a <= -5.2e-268) tmp = t * (1.0 - (y / z)); elseif (a <= 2.3e-62) tmp = y * ((x - t) / z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.3e+60], t$95$1, If[LessEqual[a, -5.2e-268], N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.3e-62], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{y}{\frac{a}{x}}\\
\mathbf{if}\;a \leq -1.3 \cdot 10^{+60}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a \leq -5.2 \cdot 10^{-268}:\\
\;\;\;\;t \cdot \left(1 - \frac{y}{z}\right)\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{-62}:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if a < -1.30000000000000004e60 or 2.3e-62 < a Initial program 86.2%
Taylor expanded in z around 0 65.1%
Taylor expanded in t around 0 52.9%
+-commutative52.9%
mul-1-neg52.9%
unsub-neg52.9%
associate-/l*57.6%
Simplified57.6%
if -1.30000000000000004e60 < a < -5.20000000000000005e-268Initial program 75.7%
Taylor expanded in a around 0 53.3%
+-commutative53.3%
mul-1-neg53.3%
unsub-neg53.3%
associate-/l*57.5%
Simplified57.5%
Taylor expanded in t around inf 52.8%
if -5.20000000000000005e-268 < a < 2.3e-62Initial program 68.0%
Taylor expanded in a around 0 55.5%
associate-*r/55.5%
neg-mul-155.5%
Simplified55.5%
Taylor expanded in y around inf 52.9%
div-sub52.9%
Simplified52.9%
Final simplification55.1%
(FPCore (x y z t a)
:precision binary64
(if (<= a -6.2e+49)
x
(if (<= a -1.3e-242)
t
(if (<= a 1.85e-51) (* y (/ x z)) (if (<= a 6.2e+37) t x)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.2e+49) {
tmp = x;
} else if (a <= -1.3e-242) {
tmp = t;
} else if (a <= 1.85e-51) {
tmp = y * (x / z);
} else if (a <= 6.2e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-6.2d+49)) then
tmp = x
else if (a <= (-1.3d-242)) then
tmp = t
else if (a <= 1.85d-51) then
tmp = y * (x / z)
else if (a <= 6.2d+37) then
tmp = t
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.2e+49) {
tmp = x;
} else if (a <= -1.3e-242) {
tmp = t;
} else if (a <= 1.85e-51) {
tmp = y * (x / z);
} else if (a <= 6.2e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6.2e+49: tmp = x elif a <= -1.3e-242: tmp = t elif a <= 1.85e-51: tmp = y * (x / z) elif a <= 6.2e+37: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6.2e+49) tmp = x; elseif (a <= -1.3e-242) tmp = t; elseif (a <= 1.85e-51) tmp = Float64(y * Float64(x / z)); elseif (a <= 6.2e+37) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -6.2e+49) tmp = x; elseif (a <= -1.3e-242) tmp = t; elseif (a <= 1.85e-51) tmp = y * (x / z); elseif (a <= 6.2e+37) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6.2e+49], x, If[LessEqual[a, -1.3e-242], t, If[LessEqual[a, 1.85e-51], N[(y * N[(x / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 6.2e+37], t, x]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.2 \cdot 10^{+49}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -1.3 \cdot 10^{-242}:\\
\;\;\;\;t\\
\mathbf{elif}\;a \leq 1.85 \cdot 10^{-51}:\\
\;\;\;\;y \cdot \frac{x}{z}\\
\mathbf{elif}\;a \leq 6.2 \cdot 10^{+37}:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -6.19999999999999985e49 or 6.2000000000000004e37 < a Initial program 89.8%
Taylor expanded in a around inf 51.4%
if -6.19999999999999985e49 < a < -1.30000000000000009e-242 or 1.84999999999999987e-51 < a < 6.2000000000000004e37Initial program 74.4%
Taylor expanded in z around inf 36.4%
if -1.30000000000000009e-242 < a < 1.84999999999999987e-51Initial program 69.0%
Taylor expanded in a around 0 47.9%
+-commutative47.9%
mul-1-neg47.9%
unsub-neg47.9%
associate-/l*54.5%
Simplified54.5%
Taylor expanded in x around -inf 42.1%
expm1-log1p-u27.5%
expm1-udef14.1%
associate-/l*14.1%
Applied egg-rr14.1%
expm1-def26.3%
expm1-log1p43.4%
associate-/r/40.7%
*-commutative40.7%
associate-*r/42.1%
*-commutative42.1%
associate-*r/43.4%
Simplified43.4%
Final simplification44.1%
(FPCore (x y z t a)
:precision binary64
(if (<= a -1.45e+50)
x
(if (<= a -5.2e-248)
t
(if (<= a 9.2e-51) (/ y (/ z x)) (if (<= a 1.02e+38) t x)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.45e+50) {
tmp = x;
} else if (a <= -5.2e-248) {
tmp = t;
} else if (a <= 9.2e-51) {
tmp = y / (z / x);
} else if (a <= 1.02e+38) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-1.45d+50)) then
tmp = x
else if (a <= (-5.2d-248)) then
tmp = t
else if (a <= 9.2d-51) then
tmp = y / (z / x)
else if (a <= 1.02d+38) then
tmp = t
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.45e+50) {
tmp = x;
} else if (a <= -5.2e-248) {
tmp = t;
} else if (a <= 9.2e-51) {
tmp = y / (z / x);
} else if (a <= 1.02e+38) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.45e+50: tmp = x elif a <= -5.2e-248: tmp = t elif a <= 9.2e-51: tmp = y / (z / x) elif a <= 1.02e+38: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.45e+50) tmp = x; elseif (a <= -5.2e-248) tmp = t; elseif (a <= 9.2e-51) tmp = Float64(y / Float64(z / x)); elseif (a <= 1.02e+38) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.45e+50) tmp = x; elseif (a <= -5.2e-248) tmp = t; elseif (a <= 9.2e-51) tmp = y / (z / x); elseif (a <= 1.02e+38) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.45e+50], x, If[LessEqual[a, -5.2e-248], t, If[LessEqual[a, 9.2e-51], N[(y / N[(z / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.02e+38], t, x]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.45 \cdot 10^{+50}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -5.2 \cdot 10^{-248}:\\
\;\;\;\;t\\
\mathbf{elif}\;a \leq 9.2 \cdot 10^{-51}:\\
\;\;\;\;\frac{y}{\frac{z}{x}}\\
\mathbf{elif}\;a \leq 1.02 \cdot 10^{+38}:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -1.45e50 or 1.02000000000000006e38 < a Initial program 89.8%
Taylor expanded in a around inf 51.4%
if -1.45e50 < a < -5.20000000000000013e-248 or 9.20000000000000007e-51 < a < 1.02000000000000006e38Initial program 74.4%
Taylor expanded in z around inf 36.4%
if -5.20000000000000013e-248 < a < 9.20000000000000007e-51Initial program 69.0%
Taylor expanded in a around 0 47.9%
+-commutative47.9%
mul-1-neg47.9%
unsub-neg47.9%
associate-/l*54.5%
Simplified54.5%
Taylor expanded in x around -inf 42.1%
associate-/l*43.4%
Simplified43.4%
Final simplification44.1%
(FPCore (x y z t a) :precision binary64 (if (<= a -7e+73) x (if (<= a 3.2e-48) (* y (/ (- x t) z)) (if (<= a 1.42e+37) t x))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -7e+73) {
tmp = x;
} else if (a <= 3.2e-48) {
tmp = y * ((x - t) / z);
} else if (a <= 1.42e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-7d+73)) then
tmp = x
else if (a <= 3.2d-48) then
tmp = y * ((x - t) / z)
else if (a <= 1.42d+37) then
tmp = t
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -7e+73) {
tmp = x;
} else if (a <= 3.2e-48) {
tmp = y * ((x - t) / z);
} else if (a <= 1.42e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -7e+73: tmp = x elif a <= 3.2e-48: tmp = y * ((x - t) / z) elif a <= 1.42e+37: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -7e+73) tmp = x; elseif (a <= 3.2e-48) tmp = Float64(y * Float64(Float64(x - t) / z)); elseif (a <= 1.42e+37) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -7e+73) tmp = x; elseif (a <= 3.2e-48) tmp = y * ((x - t) / z); elseif (a <= 1.42e+37) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -7e+73], x, If[LessEqual[a, 3.2e-48], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.42e+37], t, x]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -7 \cdot 10^{+73}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 3.2 \cdot 10^{-48}:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{elif}\;a \leq 1.42 \cdot 10^{+37}:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -7.00000000000000004e73 or 1.4199999999999999e37 < a Initial program 90.2%
Taylor expanded in a around inf 53.9%
if -7.00000000000000004e73 < a < 3.1999999999999998e-48Initial program 72.1%
Taylor expanded in a around 0 55.2%
associate-*r/55.2%
neg-mul-155.2%
Simplified55.2%
Taylor expanded in y around inf 46.0%
div-sub46.1%
Simplified46.1%
if 3.1999999999999998e-48 < a < 1.4199999999999999e37Initial program 73.9%
Taylor expanded in z around inf 40.6%
Final simplification48.4%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.65e+50) x (if (<= a 1.45e+37) t x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.65e+50) {
tmp = x;
} else if (a <= 1.45e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-1.65d+50)) then
tmp = x
else if (a <= 1.45d+37) then
tmp = t
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.65e+50) {
tmp = x;
} else if (a <= 1.45e+37) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.65e+50: tmp = x elif a <= 1.45e+37: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.65e+50) tmp = x; elseif (a <= 1.45e+37) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.65e+50) tmp = x; elseif (a <= 1.45e+37) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.65e+50], x, If[LessEqual[a, 1.45e+37], t, x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.65 \cdot 10^{+50}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 1.45 \cdot 10^{+37}:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -1.65e50 or 1.44999999999999989e37 < a Initial program 89.8%
Taylor expanded in a around inf 51.4%
if -1.65e50 < a < 1.44999999999999989e37Initial program 72.0%
Taylor expanded in z around inf 32.9%
Final simplification40.0%
(FPCore (x y z t a) :precision binary64 t)
double code(double x, double y, double z, double t, double a) {
return 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 = t
end function
public static double code(double x, double y, double z, double t, double a) {
return t;
}
def code(x, y, z, t, a): return t
function code(x, y, z, t, a) return t end
function tmp = code(x, y, z, t, a) tmp = t; end
code[x_, y_, z_, t_, a_] := t
\begin{array}{l}
\\
t
\end{array}
Initial program 78.8%
Taylor expanded in z around inf 24.6%
Final simplification24.6%
herbie shell --seed 2023240
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))