
(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(Float64(y - z) * 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[(N[(y - z), $MachinePrecision] * N[(t - x), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 27 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(Float64(y - z) * 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[(N[(y - z), $MachinePrecision] * N[(t - x), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\end{array}
(FPCore (x y z t a) :precision binary64 (if (or (<= z -7.2e+165) (not (<= z 3.5e+122))) (pow (cbrt (+ t (* (- t x) (/ (- a y) z)))) 3.0) (fma (- t x) (/ (- y z) (- a z)) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -7.2e+165) || !(z <= 3.5e+122)) {
tmp = pow(cbrt((t + ((t - x) * ((a - y) / z)))), 3.0);
} else {
tmp = fma((t - x), ((y - z) / (a - z)), x);
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -7.2e+165) || !(z <= 3.5e+122)) tmp = cbrt(Float64(t + Float64(Float64(t - x) * Float64(Float64(a - y) / z)))) ^ 3.0; else tmp = fma(Float64(t - x), Float64(Float64(y - z) / Float64(a - z)), x); end return tmp end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -7.2e+165], N[Not[LessEqual[z, 3.5e+122]], $MachinePrecision]], N[Power[N[Power[N[(t + N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision], N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -7.2 \cdot 10^{+165} \lor \neg \left(z \leq 3.5 \cdot 10^{+122}\right):\\
\;\;\;\;{\left(\sqrt[3]{t + \left(t - x\right) \cdot \frac{a - y}{z}}\right)}^{3}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
\end{array}
\end{array}
if z < -7.1999999999999996e165 or 3.50000000000000014e122 < z Initial program 36.5%
associate-/l*58.1%
Simplified58.1%
Taylor expanded in z around inf 72.7%
associate--l+72.7%
associate-*r/72.7%
associate-*r/72.7%
mul-1-neg72.7%
div-sub72.7%
mul-1-neg72.7%
distribute-lft-out--72.7%
associate-*r/72.7%
mul-1-neg72.7%
unsub-neg72.7%
distribute-rgt-out--72.7%
Simplified72.7%
add-cube-cbrt71.4%
pow371.4%
associate-/l*91.3%
Applied egg-rr91.3%
if -7.1999999999999996e165 < z < 3.50000000000000014e122Initial program 84.6%
+-commutative84.6%
*-commutative84.6%
associate-/l*91.8%
fma-define91.8%
Simplified91.8%
Final simplification91.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* (- y z) (/ (- t x) (- a z)))))
(t_2 (+ x (/ (* (- t x) (- y z)) (- a z)))))
(if (<= t_2 (- INFINITY))
t_1
(if (<= t_2 -1e-264)
t_2
(if (<= t_2 1e-267)
(+ t (/ (* (- t x) (- a y)) z))
(if (<= t_2 4e+266) t_2 t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((y - z) * ((t - x) / (a - z)));
double t_2 = x + (((t - x) * (y - z)) / (a - z));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = t_1;
} else if (t_2 <= -1e-264) {
tmp = t_2;
} else if (t_2 <= 1e-267) {
tmp = t + (((t - x) * (a - y)) / z);
} else if (t_2 <= 4e+266) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((y - z) * ((t - x) / (a - z)));
double t_2 = x + (((t - x) * (y - z)) / (a - z));
double tmp;
if (t_2 <= -Double.POSITIVE_INFINITY) {
tmp = t_1;
} else if (t_2 <= -1e-264) {
tmp = t_2;
} else if (t_2 <= 1e-267) {
tmp = t + (((t - x) * (a - y)) / z);
} else if (t_2 <= 4e+266) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + ((y - z) * ((t - x) / (a - z))) t_2 = x + (((t - x) * (y - z)) / (a - z)) tmp = 0 if t_2 <= -math.inf: tmp = t_1 elif t_2 <= -1e-264: tmp = t_2 elif t_2 <= 1e-267: tmp = t + (((t - x) * (a - y)) / z) elif t_2 <= 4e+266: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) t_2 = Float64(x + Float64(Float64(Float64(t - x) * Float64(y - z)) / Float64(a - z))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = t_1; elseif (t_2 <= -1e-264) tmp = t_2; elseif (t_2 <= 1e-267) tmp = Float64(t + Float64(Float64(Float64(t - x) * Float64(a - y)) / z)); elseif (t_2 <= 4e+266) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + ((y - z) * ((t - x) / (a - z))); t_2 = x + (((t - x) * (y - z)) / (a - z)); tmp = 0.0; if (t_2 <= -Inf) tmp = t_1; elseif (t_2 <= -1e-264) tmp = t_2; elseif (t_2 <= 1e-267) tmp = t + (((t - x) * (a - y)) / z); elseif (t_2 <= 4e+266) tmp = t_2; 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[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(N[(t - x), $MachinePrecision] * N[(y - z), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], t$95$1, If[LessEqual[t$95$2, -1e-264], t$95$2, If[LessEqual[t$95$2, 1e-267], N[(t + N[(N[(N[(t - x), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 4e+266], t$95$2, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
t_2 := x + \frac{\left(t - x\right) \cdot \left(y - z\right)}{a - z}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-264}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_2 \leq 10^{-267}:\\
\;\;\;\;t + \frac{\left(t - x\right) \cdot \left(a - y\right)}{z}\\
\mathbf{elif}\;t\_2 \leq 4 \cdot 10^{+266}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -inf.0 or 4.0000000000000001e266 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) Initial program 47.8%
associate-/l*81.7%
Simplified81.7%
if -inf.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -1e-264 or 9.9999999999999998e-268 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 4.0000000000000001e266Initial program 97.6%
if -1e-264 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 9.9999999999999998e-268Initial program 8.5%
associate-/l*4.6%
Simplified4.6%
Taylor expanded in z around inf 99.6%
associate--l+99.6%
associate-*r/99.6%
associate-*r/99.6%
mul-1-neg99.6%
div-sub99.6%
mul-1-neg99.6%
distribute-lft-out--99.6%
associate-*r/99.6%
mul-1-neg99.6%
unsub-neg99.6%
distribute-rgt-out--99.6%
Simplified99.6%
Final simplification90.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (/ (* (- t x) (- y z)) (- a z)))))
(if (or (<= t_1 -1e-264) (not (<= t_1 1e-267)))
(+ x (/ (- t x) (/ (- a z) (- y z))))
(+ t (/ (* (- t x) (- a y)) z)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((t - x) * (y - z)) / (a - z));
double tmp;
if ((t_1 <= -1e-264) || !(t_1 <= 1e-267)) {
tmp = x + ((t - x) / ((a - z) / (y - z)));
} else {
tmp = t + (((t - x) * (a - y)) / 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 + (((t - x) * (y - z)) / (a - z))
if ((t_1 <= (-1d-264)) .or. (.not. (t_1 <= 1d-267))) then
tmp = x + ((t - x) / ((a - z) / (y - z)))
else
tmp = t + (((t - x) * (a - y)) / 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 + (((t - x) * (y - z)) / (a - z));
double tmp;
if ((t_1 <= -1e-264) || !(t_1 <= 1e-267)) {
tmp = x + ((t - x) / ((a - z) / (y - z)));
} else {
tmp = t + (((t - x) * (a - y)) / z);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (((t - x) * (y - z)) / (a - z)) tmp = 0 if (t_1 <= -1e-264) or not (t_1 <= 1e-267): tmp = x + ((t - x) / ((a - z) / (y - z))) else: tmp = t + (((t - x) * (a - y)) / z) return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(Float64(t - x) * Float64(y - z)) / Float64(a - z))) tmp = 0.0 if ((t_1 <= -1e-264) || !(t_1 <= 1e-267)) tmp = Float64(x + Float64(Float64(t - x) / Float64(Float64(a - z) / Float64(y - z)))); else tmp = Float64(t + Float64(Float64(Float64(t - x) * Float64(a - y)) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (((t - x) * (y - z)) / (a - z)); tmp = 0.0; if ((t_1 <= -1e-264) || ~((t_1 <= 1e-267))) tmp = x + ((t - x) / ((a - z) / (y - z))); else tmp = t + (((t - x) * (a - y)) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(N[(t - x), $MachinePrecision] * N[(y - z), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -1e-264], N[Not[LessEqual[t$95$1, 1e-267]], $MachinePrecision]], N[(x + N[(N[(t - x), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t + N[(N[(N[(t - x), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{\left(t - x\right) \cdot \left(y - z\right)}{a - z}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-264} \lor \neg \left(t\_1 \leq 10^{-267}\right):\\
\;\;\;\;x + \frac{t - x}{\frac{a - z}{y - z}}\\
\mathbf{else}:\\
\;\;\;\;t + \frac{\left(t - x\right) \cdot \left(a - y\right)}{z}\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -1e-264 or 9.9999999999999998e-268 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) Initial program 74.5%
associate-/l*85.5%
Simplified85.5%
*-commutative85.5%
associate-*l/74.5%
associate-*r/89.4%
clear-num89.4%
un-div-inv89.4%
Applied egg-rr89.4%
if -1e-264 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 9.9999999999999998e-268Initial program 8.5%
associate-/l*4.6%
Simplified4.6%
Taylor expanded in z around inf 99.6%
associate--l+99.6%
associate-*r/99.6%
associate-*r/99.6%
mul-1-neg99.6%
div-sub99.6%
mul-1-neg99.6%
distribute-lft-out--99.6%
associate-*r/99.6%
mul-1-neg99.6%
unsub-neg99.6%
distribute-rgt-out--99.6%
Simplified99.6%
Final simplification90.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* t (/ y a)))))
(if (<= z -750000000.0)
(/ t (/ (- z a) z))
(if (<= z 3.8e-87)
t_1
(if (<= z 116000000000.0)
(* y (/ (- x t) z))
(if (<= z 4.4e+45)
t_1
(if (<= z 1.9e+91)
(/ (* t (- z y)) z)
(if (<= z 2.05e+227)
(* t (/ z (- z a)))
(if (<= z 7.5e+234) (* x (/ (- y a) z)) (+ t x))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (t * (y / a));
double tmp;
if (z <= -750000000.0) {
tmp = t / ((z - a) / z);
} else if (z <= 3.8e-87) {
tmp = t_1;
} else if (z <= 116000000000.0) {
tmp = y * ((x - t) / z);
} else if (z <= 4.4e+45) {
tmp = t_1;
} else if (z <= 1.9e+91) {
tmp = (t * (z - y)) / z;
} else if (z <= 2.05e+227) {
tmp = t * (z / (z - a));
} else if (z <= 7.5e+234) {
tmp = x * ((y - a) / z);
} else {
tmp = 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 = x + (t * (y / a))
if (z <= (-750000000.0d0)) then
tmp = t / ((z - a) / z)
else if (z <= 3.8d-87) then
tmp = t_1
else if (z <= 116000000000.0d0) then
tmp = y * ((x - t) / z)
else if (z <= 4.4d+45) then
tmp = t_1
else if (z <= 1.9d+91) then
tmp = (t * (z - y)) / z
else if (z <= 2.05d+227) then
tmp = t * (z / (z - a))
else if (z <= 7.5d+234) then
tmp = x * ((y - a) / z)
else
tmp = 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 = x + (t * (y / a));
double tmp;
if (z <= -750000000.0) {
tmp = t / ((z - a) / z);
} else if (z <= 3.8e-87) {
tmp = t_1;
} else if (z <= 116000000000.0) {
tmp = y * ((x - t) / z);
} else if (z <= 4.4e+45) {
tmp = t_1;
} else if (z <= 1.9e+91) {
tmp = (t * (z - y)) / z;
} else if (z <= 2.05e+227) {
tmp = t * (z / (z - a));
} else if (z <= 7.5e+234) {
tmp = x * ((y - a) / z);
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (t * (y / a)) tmp = 0 if z <= -750000000.0: tmp = t / ((z - a) / z) elif z <= 3.8e-87: tmp = t_1 elif z <= 116000000000.0: tmp = y * ((x - t) / z) elif z <= 4.4e+45: tmp = t_1 elif z <= 1.9e+91: tmp = (t * (z - y)) / z elif z <= 2.05e+227: tmp = t * (z / (z - a)) elif z <= 7.5e+234: tmp = x * ((y - a) / z) else: tmp = t + x return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(t * Float64(y / a))) tmp = 0.0 if (z <= -750000000.0) tmp = Float64(t / Float64(Float64(z - a) / z)); elseif (z <= 3.8e-87) tmp = t_1; elseif (z <= 116000000000.0) tmp = Float64(y * Float64(Float64(x - t) / z)); elseif (z <= 4.4e+45) tmp = t_1; elseif (z <= 1.9e+91) tmp = Float64(Float64(t * Float64(z - y)) / z); elseif (z <= 2.05e+227) tmp = Float64(t * Float64(z / Float64(z - a))); elseif (z <= 7.5e+234) tmp = Float64(x * Float64(Float64(y - a) / z)); else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (t * (y / a)); tmp = 0.0; if (z <= -750000000.0) tmp = t / ((z - a) / z); elseif (z <= 3.8e-87) tmp = t_1; elseif (z <= 116000000000.0) tmp = y * ((x - t) / z); elseif (z <= 4.4e+45) tmp = t_1; elseif (z <= 1.9e+91) tmp = (t * (z - y)) / z; elseif (z <= 2.05e+227) tmp = t * (z / (z - a)); elseif (z <= 7.5e+234) tmp = x * ((y - a) / z); else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -750000000.0], N[(t / N[(N[(z - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.8e-87], t$95$1, If[LessEqual[z, 116000000000.0], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.4e+45], t$95$1, If[LessEqual[z, 1.9e+91], N[(N[(t * N[(z - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[z, 2.05e+227], N[(t * N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 7.5e+234], N[(x * N[(N[(y - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(t + x), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + t \cdot \frac{y}{a}\\
\mathbf{if}\;z \leq -750000000:\\
\;\;\;\;\frac{t}{\frac{z - a}{z}}\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{-87}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 116000000000:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{elif}\;z \leq 4.4 \cdot 10^{+45}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.9 \cdot 10^{+91}:\\
\;\;\;\;\frac{t \cdot \left(z - y\right)}{z}\\
\mathbf{elif}\;z \leq 2.05 \cdot 10^{+227}:\\
\;\;\;\;t \cdot \frac{z}{z - a}\\
\mathbf{elif}\;z \leq 7.5 \cdot 10^{+234}:\\
\;\;\;\;x \cdot \frac{y - a}{z}\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if z < -7.5e8Initial program 47.2%
associate-/l*64.0%
Simplified64.0%
Taylor expanded in x around 0 45.9%
associate-/l*65.2%
Simplified65.2%
clear-num65.2%
un-div-inv65.2%
Applied egg-rr65.2%
Taylor expanded in y around 0 51.0%
associate-*r/51.0%
mul-1-neg51.0%
Simplified51.0%
if -7.5e8 < z < 3.8e-87 or 1.16e11 < z < 4.4000000000000001e45Initial program 94.8%
associate-/l*94.9%
Simplified94.9%
Taylor expanded in z around 0 72.8%
Taylor expanded in t around inf 63.1%
associate-/l*69.0%
Simplified69.0%
if 3.8e-87 < z < 1.16e11Initial program 77.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in z around inf 72.7%
associate--l+72.7%
associate-*r/72.7%
associate-*r/72.7%
mul-1-neg72.7%
div-sub77.5%
mul-1-neg77.5%
distribute-lft-out--77.5%
associate-*r/77.5%
mul-1-neg77.5%
unsub-neg77.5%
distribute-rgt-out--82.3%
Simplified82.3%
Taylor expanded in y around -inf 59.5%
mul-1-neg59.5%
associate-/l*63.6%
*-commutative63.6%
distribute-rgt-neg-in63.6%
Simplified63.6%
if 4.4000000000000001e45 < z < 1.8999999999999999e91Initial program 79.3%
associate-/l*88.8%
Simplified88.8%
Taylor expanded in x around 0 57.9%
associate-/l*67.5%
Simplified67.5%
Taylor expanded in a around 0 57.5%
associate-*r/57.5%
mul-1-neg57.5%
*-commutative57.5%
distribute-rgt-neg-out57.5%
Simplified57.5%
if 1.8999999999999999e91 < z < 2.04999999999999998e227Initial program 46.3%
associate-/l*64.1%
Simplified64.1%
Taylor expanded in x around 0 40.8%
associate-/l*57.4%
Simplified57.4%
Taylor expanded in y around 0 41.0%
mul-1-neg41.0%
associate-/l*54.9%
distribute-rgt-neg-in54.9%
distribute-frac-neg254.9%
Simplified54.9%
if 2.04999999999999998e227 < z < 7.5000000000000004e234Initial program 2.5%
associate-/l*3.7%
Simplified3.7%
Taylor expanded in z around inf 80.8%
associate--l+80.8%
associate-*r/80.8%
associate-*r/80.8%
mul-1-neg80.8%
div-sub80.8%
mul-1-neg80.8%
distribute-lft-out--80.8%
associate-*r/80.8%
mul-1-neg80.8%
unsub-neg80.8%
distribute-rgt-out--80.8%
Simplified80.8%
Taylor expanded in t around 0 62.7%
associate-/l*80.8%
Simplified80.8%
if 7.5000000000000004e234 < z Initial program 28.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in t around inf 70.8%
Taylor expanded in z around inf 79.3%
Final simplification63.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ z (- z a)))))
(if (<= z -125000000.0)
t_1
(if (<= z 1.8e-88)
(+ x (* t (/ y a)))
(if (<= z 410000000.0)
(* y (/ (- x t) z))
(if (or (<= z 2.6e+115) (not (<= z 1.38e+234))) (+ t x) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (z / (z - a));
double tmp;
if (z <= -125000000.0) {
tmp = t_1;
} else if (z <= 1.8e-88) {
tmp = x + (t * (y / a));
} else if (z <= 410000000.0) {
tmp = y * ((x - t) / z);
} else if ((z <= 2.6e+115) || !(z <= 1.38e+234)) {
tmp = t + x;
} 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 = t * (z / (z - a))
if (z <= (-125000000.0d0)) then
tmp = t_1
else if (z <= 1.8d-88) then
tmp = x + (t * (y / a))
else if (z <= 410000000.0d0) then
tmp = y * ((x - t) / z)
else if ((z <= 2.6d+115) .or. (.not. (z <= 1.38d+234))) then
tmp = t + x
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 * (z / (z - a));
double tmp;
if (z <= -125000000.0) {
tmp = t_1;
} else if (z <= 1.8e-88) {
tmp = x + (t * (y / a));
} else if (z <= 410000000.0) {
tmp = y * ((x - t) / z);
} else if ((z <= 2.6e+115) || !(z <= 1.38e+234)) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (z / (z - a)) tmp = 0 if z <= -125000000.0: tmp = t_1 elif z <= 1.8e-88: tmp = x + (t * (y / a)) elif z <= 410000000.0: tmp = y * ((x - t) / z) elif (z <= 2.6e+115) or not (z <= 1.38e+234): tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(z / Float64(z - a))) tmp = 0.0 if (z <= -125000000.0) tmp = t_1; elseif (z <= 1.8e-88) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (z <= 410000000.0) tmp = Float64(y * Float64(Float64(x - t) / z)); elseif ((z <= 2.6e+115) || !(z <= 1.38e+234)) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (z / (z - a)); tmp = 0.0; if (z <= -125000000.0) tmp = t_1; elseif (z <= 1.8e-88) tmp = x + (t * (y / a)); elseif (z <= 410000000.0) tmp = y * ((x - t) / z); elseif ((z <= 2.6e+115) || ~((z <= 1.38e+234))) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -125000000.0], t$95$1, If[LessEqual[z, 1.8e-88], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 410000000.0], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, 2.6e+115], N[Not[LessEqual[z, 1.38e+234]], $MachinePrecision]], N[(t + x), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{z}{z - a}\\
\mathbf{if}\;z \leq -125000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.8 \cdot 10^{-88}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq 410000000:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+115} \lor \neg \left(z \leq 1.38 \cdot 10^{+234}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.25e8 or 2.6e115 < z < 1.3799999999999999e234Initial program 42.9%
associate-/l*60.0%
Simplified60.0%
Taylor expanded in x around 0 44.3%
associate-/l*62.8%
Simplified62.8%
Taylor expanded in y around 0 38.7%
mul-1-neg38.7%
associate-/l*53.0%
distribute-rgt-neg-in53.0%
distribute-frac-neg253.0%
Simplified53.0%
if -1.25e8 < z < 1.8e-88Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
Taylor expanded in t around inf 61.6%
associate-/l*67.9%
Simplified67.9%
if 1.8e-88 < z < 4.1e8Initial program 77.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in z around inf 72.7%
associate--l+72.7%
associate-*r/72.7%
associate-*r/72.7%
mul-1-neg72.7%
div-sub77.5%
mul-1-neg77.5%
distribute-lft-out--77.5%
associate-*r/77.5%
mul-1-neg77.5%
unsub-neg77.5%
distribute-rgt-out--82.3%
Simplified82.3%
Taylor expanded in y around -inf 59.5%
mul-1-neg59.5%
associate-/l*63.6%
*-commutative63.6%
distribute-rgt-neg-in63.6%
Simplified63.6%
if 4.1e8 < z < 2.6e115 or 1.3799999999999999e234 < z Initial program 56.6%
associate-/l*80.6%
Simplified80.6%
Taylor expanded in t around inf 74.2%
Taylor expanded in z around inf 64.3%
Final simplification61.9%
(FPCore (x y z t a)
:precision binary64
(if (<= z -750000000.0)
(/ t (/ (- z a) z))
(if (<= z 3.8e-87)
(+ x (* t (/ y a)))
(if (<= z 58000000.0)
(* y (/ (- x t) z))
(if (or (<= z 1.2e+115) (not (<= z 1.4e+234)))
(+ t x)
(* t (/ z (- z a))))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -750000000.0) {
tmp = t / ((z - a) / z);
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 58000000.0) {
tmp = y * ((x - t) / z);
} else if ((z <= 1.2e+115) || !(z <= 1.4e+234)) {
tmp = t + x;
} else {
tmp = t * (z / (z - a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-750000000.0d0)) then
tmp = t / ((z - a) / z)
else if (z <= 3.8d-87) then
tmp = x + (t * (y / a))
else if (z <= 58000000.0d0) then
tmp = y * ((x - t) / z)
else if ((z <= 1.2d+115) .or. (.not. (z <= 1.4d+234))) then
tmp = t + x
else
tmp = t * (z / (z - a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -750000000.0) {
tmp = t / ((z - a) / z);
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 58000000.0) {
tmp = y * ((x - t) / z);
} else if ((z <= 1.2e+115) || !(z <= 1.4e+234)) {
tmp = t + x;
} else {
tmp = t * (z / (z - a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -750000000.0: tmp = t / ((z - a) / z) elif z <= 3.8e-87: tmp = x + (t * (y / a)) elif z <= 58000000.0: tmp = y * ((x - t) / z) elif (z <= 1.2e+115) or not (z <= 1.4e+234): tmp = t + x else: tmp = t * (z / (z - a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -750000000.0) tmp = Float64(t / Float64(Float64(z - a) / z)); elseif (z <= 3.8e-87) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (z <= 58000000.0) tmp = Float64(y * Float64(Float64(x - t) / z)); elseif ((z <= 1.2e+115) || !(z <= 1.4e+234)) tmp = Float64(t + x); else tmp = Float64(t * Float64(z / Float64(z - a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -750000000.0) tmp = t / ((z - a) / z); elseif (z <= 3.8e-87) tmp = x + (t * (y / a)); elseif (z <= 58000000.0) tmp = y * ((x - t) / z); elseif ((z <= 1.2e+115) || ~((z <= 1.4e+234))) tmp = t + x; else tmp = t * (z / (z - a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -750000000.0], N[(t / N[(N[(z - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.8e-87], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 58000000.0], N[(y * N[(N[(x - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, 1.2e+115], N[Not[LessEqual[z, 1.4e+234]], $MachinePrecision]], N[(t + x), $MachinePrecision], N[(t * N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -750000000:\\
\;\;\;\;\frac{t}{\frac{z - a}{z}}\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{-87}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq 58000000:\\
\;\;\;\;y \cdot \frac{x - t}{z}\\
\mathbf{elif}\;z \leq 1.2 \cdot 10^{+115} \lor \neg \left(z \leq 1.4 \cdot 10^{+234}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t \cdot \frac{z}{z - a}\\
\end{array}
\end{array}
if z < -7.5e8Initial program 47.2%
associate-/l*64.0%
Simplified64.0%
Taylor expanded in x around 0 45.9%
associate-/l*65.2%
Simplified65.2%
clear-num65.2%
un-div-inv65.2%
Applied egg-rr65.2%
Taylor expanded in y around 0 51.0%
associate-*r/51.0%
mul-1-neg51.0%
Simplified51.0%
if -7.5e8 < z < 3.8e-87Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
Taylor expanded in t around inf 61.6%
associate-/l*67.9%
Simplified67.9%
if 3.8e-87 < z < 5.8e7Initial program 77.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in z around inf 72.7%
associate--l+72.7%
associate-*r/72.7%
associate-*r/72.7%
mul-1-neg72.7%
div-sub77.5%
mul-1-neg77.5%
distribute-lft-out--77.5%
associate-*r/77.5%
mul-1-neg77.5%
unsub-neg77.5%
distribute-rgt-out--82.3%
Simplified82.3%
Taylor expanded in y around -inf 59.5%
mul-1-neg59.5%
associate-/l*63.6%
*-commutative63.6%
distribute-rgt-neg-in63.6%
Simplified63.6%
if 5.8e7 < z < 1.2e115 or 1.3999999999999999e234 < z Initial program 56.6%
associate-/l*80.6%
Simplified80.6%
Taylor expanded in t around inf 74.2%
Taylor expanded in z around inf 64.3%
if 1.2e115 < z < 1.3999999999999999e234Initial program 36.0%
associate-/l*53.5%
Simplified53.5%
Taylor expanded in x around 0 41.7%
associate-/l*58.8%
Simplified58.8%
Taylor expanded in y around 0 41.9%
mul-1-neg41.9%
associate-/l*56.3%
distribute-rgt-neg-in56.3%
distribute-frac-neg256.3%
Simplified56.3%
Final simplification61.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ y (- a z)))))
(if (<= y -2.55e+182)
t_1
(if (<= y -1.9e+85)
(* x (/ y z))
(if (<= y -2.1e+58)
(+ t x)
(if (<= y -6.5e-6)
(* t (/ (- y z) a))
(if (<= y 1.9e+35) (+ t x) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (y / (a - z));
double tmp;
if (y <= -2.55e+182) {
tmp = t_1;
} else if (y <= -1.9e+85) {
tmp = x * (y / z);
} else if (y <= -2.1e+58) {
tmp = t + x;
} else if (y <= -6.5e-6) {
tmp = t * ((y - z) / a);
} else if (y <= 1.9e+35) {
tmp = t + x;
} 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 = t * (y / (a - z))
if (y <= (-2.55d+182)) then
tmp = t_1
else if (y <= (-1.9d+85)) then
tmp = x * (y / z)
else if (y <= (-2.1d+58)) then
tmp = t + x
else if (y <= (-6.5d-6)) then
tmp = t * ((y - z) / a)
else if (y <= 1.9d+35) then
tmp = t + x
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 / (a - z));
double tmp;
if (y <= -2.55e+182) {
tmp = t_1;
} else if (y <= -1.9e+85) {
tmp = x * (y / z);
} else if (y <= -2.1e+58) {
tmp = t + x;
} else if (y <= -6.5e-6) {
tmp = t * ((y - z) / a);
} else if (y <= 1.9e+35) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (y / (a - z)) tmp = 0 if y <= -2.55e+182: tmp = t_1 elif y <= -1.9e+85: tmp = x * (y / z) elif y <= -2.1e+58: tmp = t + x elif y <= -6.5e-6: tmp = t * ((y - z) / a) elif y <= 1.9e+35: tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -2.55e+182) tmp = t_1; elseif (y <= -1.9e+85) tmp = Float64(x * Float64(y / z)); elseif (y <= -2.1e+58) tmp = Float64(t + x); elseif (y <= -6.5e-6) tmp = Float64(t * Float64(Float64(y - z) / a)); elseif (y <= 1.9e+35) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (y / (a - z)); tmp = 0.0; if (y <= -2.55e+182) tmp = t_1; elseif (y <= -1.9e+85) tmp = x * (y / z); elseif (y <= -2.1e+58) tmp = t + x; elseif (y <= -6.5e-6) tmp = t * ((y - z) / a); elseif (y <= 1.9e+35) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.55e+182], t$95$1, If[LessEqual[y, -1.9e+85], N[(x * N[(y / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.1e+58], N[(t + x), $MachinePrecision], If[LessEqual[y, -6.5e-6], N[(t * N[(N[(y - z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.9e+35], N[(t + x), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -2.55 \cdot 10^{+182}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.9 \cdot 10^{+85}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq -2.1 \cdot 10^{+58}:\\
\;\;\;\;t + x\\
\mathbf{elif}\;y \leq -6.5 \cdot 10^{-6}:\\
\;\;\;\;t \cdot \frac{y - z}{a}\\
\mathbf{elif}\;y \leq 1.9 \cdot 10^{+35}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.55000000000000005e182 or 1.9e35 < y Initial program 76.0%
associate-/l*85.2%
Simplified85.2%
Taylor expanded in x around 0 53.7%
associate-/l*66.5%
Simplified66.5%
Taylor expanded in y around inf 49.4%
associate-/l*58.6%
Simplified58.6%
if -2.55000000000000005e182 < y < -1.89999999999999996e85Initial program 68.2%
associate-/l*83.9%
Simplified83.9%
Taylor expanded in z around inf 63.9%
associate--l+63.9%
associate-*r/63.9%
associate-*r/63.9%
mul-1-neg63.9%
div-sub64.1%
mul-1-neg64.1%
distribute-lft-out--64.1%
associate-*r/64.1%
mul-1-neg64.1%
unsub-neg64.1%
distribute-rgt-out--64.1%
Simplified64.1%
Taylor expanded in t around 0 40.1%
associate-/l*51.6%
Simplified51.6%
Taylor expanded in y around inf 40.4%
associate-/l*48.4%
Simplified48.4%
if -1.89999999999999996e85 < y < -2.10000000000000012e58 or -6.4999999999999996e-6 < y < 1.9e35Initial program 67.0%
associate-/l*75.7%
Simplified75.7%
Taylor expanded in t around inf 68.6%
Taylor expanded in z around inf 56.8%
if -2.10000000000000012e58 < y < -6.4999999999999996e-6Initial program 60.0%
associate-/l*75.6%
Simplified75.6%
Taylor expanded in x around 0 44.7%
associate-/l*52.3%
Simplified52.3%
Taylor expanded in a around inf 37.9%
associate-/l*41.4%
Simplified41.4%
Final simplification55.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ y (- z)))))
(if (<= y -1.05e+183)
t_1
(if (<= y -1.65e+82)
(* x (/ y z))
(if (<= y 2.35e+35)
(+ t x)
(if (<= y 5.4e+153)
(* t (/ y a))
(if (<= y 4.8e+180) (+ t x) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (y / -z);
double tmp;
if (y <= -1.05e+183) {
tmp = t_1;
} else if (y <= -1.65e+82) {
tmp = x * (y / z);
} else if (y <= 2.35e+35) {
tmp = t + x;
} else if (y <= 5.4e+153) {
tmp = t * (y / a);
} else if (y <= 4.8e+180) {
tmp = t + x;
} 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 = t * (y / -z)
if (y <= (-1.05d+183)) then
tmp = t_1
else if (y <= (-1.65d+82)) then
tmp = x * (y / z)
else if (y <= 2.35d+35) then
tmp = t + x
else if (y <= 5.4d+153) then
tmp = t * (y / a)
else if (y <= 4.8d+180) then
tmp = t + x
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);
double tmp;
if (y <= -1.05e+183) {
tmp = t_1;
} else if (y <= -1.65e+82) {
tmp = x * (y / z);
} else if (y <= 2.35e+35) {
tmp = t + x;
} else if (y <= 5.4e+153) {
tmp = t * (y / a);
} else if (y <= 4.8e+180) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (y / -z) tmp = 0 if y <= -1.05e+183: tmp = t_1 elif y <= -1.65e+82: tmp = x * (y / z) elif y <= 2.35e+35: tmp = t + x elif y <= 5.4e+153: tmp = t * (y / a) elif y <= 4.8e+180: tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(y / Float64(-z))) tmp = 0.0 if (y <= -1.05e+183) tmp = t_1; elseif (y <= -1.65e+82) tmp = Float64(x * Float64(y / z)); elseif (y <= 2.35e+35) tmp = Float64(t + x); elseif (y <= 5.4e+153) tmp = Float64(t * Float64(y / a)); elseif (y <= 4.8e+180) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (y / -z); tmp = 0.0; if (y <= -1.05e+183) tmp = t_1; elseif (y <= -1.65e+82) tmp = x * (y / z); elseif (y <= 2.35e+35) tmp = t + x; elseif (y <= 5.4e+153) tmp = t * (y / a); elseif (y <= 4.8e+180) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / (-z)), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.05e+183], t$95$1, If[LessEqual[y, -1.65e+82], N[(x * N[(y / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.35e+35], N[(t + x), $MachinePrecision], If[LessEqual[y, 5.4e+153], N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+180], N[(t + x), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y}{-z}\\
\mathbf{if}\;y \leq -1.05 \cdot 10^{+183}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.65 \cdot 10^{+82}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq 2.35 \cdot 10^{+35}:\\
\;\;\;\;t + x\\
\mathbf{elif}\;y \leq 5.4 \cdot 10^{+153}:\\
\;\;\;\;t \cdot \frac{y}{a}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+180}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.05e183 or 4.7999999999999997e180 < y Initial program 78.5%
associate-/l*87.9%
Simplified87.9%
Taylor expanded in z around inf 56.3%
associate--l+56.3%
associate-*r/56.3%
associate-*r/56.3%
mul-1-neg56.3%
div-sub58.3%
mul-1-neg58.3%
distribute-lft-out--58.3%
associate-*r/58.3%
mul-1-neg58.3%
unsub-neg58.3%
distribute-rgt-out--60.4%
Simplified60.4%
Taylor expanded in y around -inf 60.3%
mul-1-neg60.3%
associate-/l*68.0%
*-commutative68.0%
distribute-rgt-neg-in68.0%
Simplified68.0%
Taylor expanded in t around inf 46.2%
mul-1-neg46.2%
associate-/l*55.5%
distribute-lft-neg-in55.5%
Simplified55.5%
if -1.05e183 < y < -1.6499999999999999e82Initial program 68.2%
associate-/l*83.9%
Simplified83.9%
Taylor expanded in z around inf 63.9%
associate--l+63.9%
associate-*r/63.9%
associate-*r/63.9%
mul-1-neg63.9%
div-sub64.1%
mul-1-neg64.1%
distribute-lft-out--64.1%
associate-*r/64.1%
mul-1-neg64.1%
unsub-neg64.1%
distribute-rgt-out--64.1%
Simplified64.1%
Taylor expanded in t around 0 40.1%
associate-/l*51.6%
Simplified51.6%
Taylor expanded in y around inf 40.4%
associate-/l*48.4%
Simplified48.4%
if -1.6499999999999999e82 < y < 2.35000000000000017e35 or 5.4000000000000001e153 < y < 4.7999999999999997e180Initial program 65.4%
associate-/l*75.1%
Simplified75.1%
Taylor expanded in t around inf 67.7%
Taylor expanded in z around inf 51.9%
if 2.35000000000000017e35 < y < 5.4000000000000001e153Initial program 84.8%
associate-/l*89.4%
Simplified89.4%
Taylor expanded in x around 0 54.4%
associate-/l*59.8%
Simplified59.8%
Taylor expanded in z around 0 44.2%
Final simplification51.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (/ y z))) (t_2 (* t (/ y a))))
(if (<= y -1.02e+183)
t_2
(if (<= y -1.05e+85)
t_1
(if (<= y 2.5e+35)
(+ t x)
(if (<= y 3.7e+154) t_2 (if (<= y 3.6e+180) (+ t x) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (y / z);
double t_2 = t * (y / a);
double tmp;
if (y <= -1.02e+183) {
tmp = t_2;
} else if (y <= -1.05e+85) {
tmp = t_1;
} else if (y <= 2.5e+35) {
tmp = t + x;
} else if (y <= 3.7e+154) {
tmp = t_2;
} else if (y <= 3.6e+180) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = x * (y / z)
t_2 = t * (y / a)
if (y <= (-1.02d+183)) then
tmp = t_2
else if (y <= (-1.05d+85)) then
tmp = t_1
else if (y <= 2.5d+35) then
tmp = t + x
else if (y <= 3.7d+154) then
tmp = t_2
else if (y <= 3.6d+180) then
tmp = t + x
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);
double t_2 = t * (y / a);
double tmp;
if (y <= -1.02e+183) {
tmp = t_2;
} else if (y <= -1.05e+85) {
tmp = t_1;
} else if (y <= 2.5e+35) {
tmp = t + x;
} else if (y <= 3.7e+154) {
tmp = t_2;
} else if (y <= 3.6e+180) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (y / z) t_2 = t * (y / a) tmp = 0 if y <= -1.02e+183: tmp = t_2 elif y <= -1.05e+85: tmp = t_1 elif y <= 2.5e+35: tmp = t + x elif y <= 3.7e+154: tmp = t_2 elif y <= 3.6e+180: tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(y / z)) t_2 = Float64(t * Float64(y / a)) tmp = 0.0 if (y <= -1.02e+183) tmp = t_2; elseif (y <= -1.05e+85) tmp = t_1; elseif (y <= 2.5e+35) tmp = Float64(t + x); elseif (y <= 3.7e+154) tmp = t_2; elseif (y <= 3.6e+180) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (y / z); t_2 = t * (y / a); tmp = 0.0; if (y <= -1.02e+183) tmp = t_2; elseif (y <= -1.05e+85) tmp = t_1; elseif (y <= 2.5e+35) tmp = t + x; elseif (y <= 3.7e+154) tmp = t_2; elseif (y <= 3.6e+180) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(y / z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.02e+183], t$95$2, If[LessEqual[y, -1.05e+85], t$95$1, If[LessEqual[y, 2.5e+35], N[(t + x), $MachinePrecision], If[LessEqual[y, 3.7e+154], t$95$2, If[LessEqual[y, 3.6e+180], N[(t + x), $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{y}{z}\\
t_2 := t \cdot \frac{y}{a}\\
\mathbf{if}\;y \leq -1.02 \cdot 10^{+183}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.05 \cdot 10^{+85}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{+35}:\\
\;\;\;\;t + x\\
\mathbf{elif}\;y \leq 3.7 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+180}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.02000000000000002e183 or 2.50000000000000011e35 < y < 3.69999999999999994e154Initial program 84.9%
associate-/l*93.2%
Simplified93.2%
Taylor expanded in x around 0 58.7%
associate-/l*71.6%
Simplified71.6%
Taylor expanded in z around 0 49.4%
if -1.02000000000000002e183 < y < -1.05000000000000005e85 or 3.6000000000000002e180 < y Initial program 70.0%
associate-/l*81.7%
Simplified81.7%
Taylor expanded in z around inf 57.6%
associate--l+57.6%
associate-*r/57.6%
associate-*r/57.6%
mul-1-neg57.6%
div-sub57.7%
mul-1-neg57.7%
distribute-lft-out--57.7%
associate-*r/57.7%
mul-1-neg57.7%
unsub-neg57.7%
distribute-rgt-out--57.7%
Simplified57.7%
Taylor expanded in t around 0 33.9%
associate-/l*47.4%
Simplified47.4%
Taylor expanded in y around inf 34.1%
associate-/l*45.8%
Simplified45.8%
if -1.05000000000000005e85 < y < 2.50000000000000011e35 or 3.69999999999999994e154 < y < 3.6000000000000002e180Initial program 65.4%
associate-/l*75.1%
Simplified75.1%
Taylor expanded in t around inf 67.7%
Taylor expanded in z around inf 51.9%
Final simplification50.3%
(FPCore (x y z t a)
:precision binary64
(if (<= x -3e+160)
(* x (- 1.0 (/ y a)))
(if (<= x -2.7e+90)
(/ x (/ z (- y a)))
(if (<= x -4.6e+49)
(+ x (* t (/ y a)))
(if (<= x 1.2e+141) (* t (/ (- y z) (- a z))) (* x (/ (- y a) z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -3e+160) {
tmp = x * (1.0 - (y / a));
} else if (x <= -2.7e+90) {
tmp = x / (z / (y - a));
} else if (x <= -4.6e+49) {
tmp = x + (t * (y / a));
} else if (x <= 1.2e+141) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x * ((y - 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) :: tmp
if (x <= (-3d+160)) then
tmp = x * (1.0d0 - (y / a))
else if (x <= (-2.7d+90)) then
tmp = x / (z / (y - a))
else if (x <= (-4.6d+49)) then
tmp = x + (t * (y / a))
else if (x <= 1.2d+141) then
tmp = t * ((y - z) / (a - z))
else
tmp = x * ((y - a) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -3e+160) {
tmp = x * (1.0 - (y / a));
} else if (x <= -2.7e+90) {
tmp = x / (z / (y - a));
} else if (x <= -4.6e+49) {
tmp = x + (t * (y / a));
} else if (x <= 1.2e+141) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x * ((y - a) / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -3e+160: tmp = x * (1.0 - (y / a)) elif x <= -2.7e+90: tmp = x / (z / (y - a)) elif x <= -4.6e+49: tmp = x + (t * (y / a)) elif x <= 1.2e+141: tmp = t * ((y - z) / (a - z)) else: tmp = x * ((y - a) / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -3e+160) tmp = Float64(x * Float64(1.0 - Float64(y / a))); elseif (x <= -2.7e+90) tmp = Float64(x / Float64(z / Float64(y - a))); elseif (x <= -4.6e+49) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (x <= 1.2e+141) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(x * Float64(Float64(y - a) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -3e+160) tmp = x * (1.0 - (y / a)); elseif (x <= -2.7e+90) tmp = x / (z / (y - a)); elseif (x <= -4.6e+49) tmp = x + (t * (y / a)); elseif (x <= 1.2e+141) tmp = t * ((y - z) / (a - z)); else tmp = x * ((y - a) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -3e+160], N[(x * N[(1.0 - N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -2.7e+90], N[(x / N[(z / N[(y - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -4.6e+49], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.2e+141], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(y - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3 \cdot 10^{+160}:\\
\;\;\;\;x \cdot \left(1 - \frac{y}{a}\right)\\
\mathbf{elif}\;x \leq -2.7 \cdot 10^{+90}:\\
\;\;\;\;\frac{x}{\frac{z}{y - a}}\\
\mathbf{elif}\;x \leq -4.6 \cdot 10^{+49}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;x \leq 1.2 \cdot 10^{+141}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{y - a}{z}\\
\end{array}
\end{array}
if x < -2.9999999999999999e160Initial program 51.5%
associate-/l*66.9%
Simplified66.9%
Taylor expanded in z around 0 49.7%
Taylor expanded in x around inf 55.1%
mul-1-neg55.1%
unsub-neg55.1%
Simplified55.1%
if -2.9999999999999999e160 < x < -2.7e90Initial program 51.5%
associate-/l*65.3%
Simplified65.3%
Taylor expanded in z around inf 57.9%
associate--l+57.9%
associate-*r/57.9%
associate-*r/57.9%
mul-1-neg57.9%
div-sub57.9%
mul-1-neg57.9%
distribute-lft-out--57.9%
associate-*r/57.9%
mul-1-neg57.9%
unsub-neg57.9%
distribute-rgt-out--65.1%
Simplified65.1%
Taylor expanded in t around 0 65.3%
associate-/l*71.9%
Simplified71.9%
clear-num71.7%
un-div-inv72.0%
Applied egg-rr72.0%
if -2.7e90 < x < -4.60000000000000004e49Initial program 92.0%
associate-/l*92.1%
Simplified92.1%
Taylor expanded in z around 0 77.7%
Taylor expanded in t around inf 69.6%
associate-/l*77.7%
Simplified77.7%
if -4.60000000000000004e49 < x < 1.19999999999999999e141Initial program 76.1%
associate-/l*84.9%
Simplified84.9%
Taylor expanded in x around 0 56.8%
associate-/l*70.0%
Simplified70.0%
if 1.19999999999999999e141 < x Initial program 54.9%
associate-/l*66.5%
Simplified66.5%
Taylor expanded in z around inf 39.4%
associate--l+39.4%
associate-*r/39.4%
associate-*r/39.4%
mul-1-neg39.4%
div-sub42.5%
mul-1-neg42.5%
distribute-lft-out--42.5%
associate-*r/42.5%
mul-1-neg42.5%
unsub-neg42.5%
distribute-rgt-out--42.5%
Simplified42.5%
Taylor expanded in t around 0 36.9%
associate-/l*57.4%
Simplified57.4%
Final simplification66.9%
(FPCore (x y z t a)
:precision binary64
(if (<= x -6.6e+160)
(* x (- 1.0 (/ y a)))
(if (<= x -1.1e+91)
(* y (/ (- t x) (- a z)))
(if (<= x -3.2e+49)
(+ x (* t (/ y a)))
(if (<= x 1.15e+144) (* t (/ (- y z) (- a z))) (* x (/ (- y a) z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -6.6e+160) {
tmp = x * (1.0 - (y / a));
} else if (x <= -1.1e+91) {
tmp = y * ((t - x) / (a - z));
} else if (x <= -3.2e+49) {
tmp = x + (t * (y / a));
} else if (x <= 1.15e+144) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x * ((y - 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) :: tmp
if (x <= (-6.6d+160)) then
tmp = x * (1.0d0 - (y / a))
else if (x <= (-1.1d+91)) then
tmp = y * ((t - x) / (a - z))
else if (x <= (-3.2d+49)) then
tmp = x + (t * (y / a))
else if (x <= 1.15d+144) then
tmp = t * ((y - z) / (a - z))
else
tmp = x * ((y - a) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -6.6e+160) {
tmp = x * (1.0 - (y / a));
} else if (x <= -1.1e+91) {
tmp = y * ((t - x) / (a - z));
} else if (x <= -3.2e+49) {
tmp = x + (t * (y / a));
} else if (x <= 1.15e+144) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x * ((y - a) / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -6.6e+160: tmp = x * (1.0 - (y / a)) elif x <= -1.1e+91: tmp = y * ((t - x) / (a - z)) elif x <= -3.2e+49: tmp = x + (t * (y / a)) elif x <= 1.15e+144: tmp = t * ((y - z) / (a - z)) else: tmp = x * ((y - a) / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -6.6e+160) tmp = Float64(x * Float64(1.0 - Float64(y / a))); elseif (x <= -1.1e+91) tmp = Float64(y * Float64(Float64(t - x) / Float64(a - z))); elseif (x <= -3.2e+49) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (x <= 1.15e+144) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(x * Float64(Float64(y - a) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -6.6e+160) tmp = x * (1.0 - (y / a)); elseif (x <= -1.1e+91) tmp = y * ((t - x) / (a - z)); elseif (x <= -3.2e+49) tmp = x + (t * (y / a)); elseif (x <= 1.15e+144) tmp = t * ((y - z) / (a - z)); else tmp = x * ((y - a) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -6.6e+160], N[(x * N[(1.0 - N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -1.1e+91], N[(y * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.2e+49], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.15e+144], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(y - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6.6 \cdot 10^{+160}:\\
\;\;\;\;x \cdot \left(1 - \frac{y}{a}\right)\\
\mathbf{elif}\;x \leq -1.1 \cdot 10^{+91}:\\
\;\;\;\;y \cdot \frac{t - x}{a - z}\\
\mathbf{elif}\;x \leq -3.2 \cdot 10^{+49}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;x \leq 1.15 \cdot 10^{+144}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{y - a}{z}\\
\end{array}
\end{array}
if x < -6.5999999999999994e160Initial program 52.8%
associate-/l*65.9%
Simplified65.9%
Taylor expanded in z around 0 51.0%
Taylor expanded in x around inf 56.8%
mul-1-neg56.8%
unsub-neg56.8%
Simplified56.8%
if -6.5999999999999994e160 < x < -1.1e91Initial program 48.6%
associate-/l*67.6%
Simplified67.6%
Taylor expanded in y around inf 68.2%
div-sub68.2%
Simplified68.2%
if -1.1e91 < x < -3.20000000000000014e49Initial program 92.0%
associate-/l*92.1%
Simplified92.1%
Taylor expanded in z around 0 77.7%
Taylor expanded in t around inf 69.6%
associate-/l*77.7%
Simplified77.7%
if -3.20000000000000014e49 < x < 1.1500000000000001e144Initial program 76.1%
associate-/l*84.9%
Simplified84.9%
Taylor expanded in x around 0 56.8%
associate-/l*70.0%
Simplified70.0%
if 1.1500000000000001e144 < x Initial program 54.9%
associate-/l*66.5%
Simplified66.5%
Taylor expanded in z around inf 39.4%
associate--l+39.4%
associate-*r/39.4%
associate-*r/39.4%
mul-1-neg39.4%
div-sub42.5%
mul-1-neg42.5%
distribute-lft-out--42.5%
associate-*r/42.5%
mul-1-neg42.5%
unsub-neg42.5%
distribute-rgt-out--42.5%
Simplified42.5%
Taylor expanded in t around 0 36.9%
associate-/l*57.4%
Simplified57.4%
Final simplification67.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ (- y z) (- a z)))))
(if (<= z -280000000.0)
t_1
(if (<= z 3.5e-87)
(- x (/ (* y (- x t)) a))
(if (<= z 800000000.0)
(* y (/ (- t x) (- a z)))
(if (<= z 2e+45) (+ x (* t (/ y a))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * ((y - z) / (a - z));
double tmp;
if (z <= -280000000.0) {
tmp = t_1;
} else if (z <= 3.5e-87) {
tmp = x - ((y * (x - t)) / a);
} else if (z <= 800000000.0) {
tmp = y * ((t - x) / (a - z));
} else if (z <= 2e+45) {
tmp = x + (t * (y / a));
} 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 = t * ((y - z) / (a - z))
if (z <= (-280000000.0d0)) then
tmp = t_1
else if (z <= 3.5d-87) then
tmp = x - ((y * (x - t)) / a)
else if (z <= 800000000.0d0) then
tmp = y * ((t - x) / (a - z))
else if (z <= 2d+45) then
tmp = x + (t * (y / a))
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 tmp;
if (z <= -280000000.0) {
tmp = t_1;
} else if (z <= 3.5e-87) {
tmp = x - ((y * (x - t)) / a);
} else if (z <= 800000000.0) {
tmp = y * ((t - x) / (a - z));
} else if (z <= 2e+45) {
tmp = x + (t * (y / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * ((y - z) / (a - z)) tmp = 0 if z <= -280000000.0: tmp = t_1 elif z <= 3.5e-87: tmp = x - ((y * (x - t)) / a) elif z <= 800000000.0: tmp = y * ((t - x) / (a - z)) elif z <= 2e+45: tmp = x + (t * (y / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(Float64(y - z) / Float64(a - z))) tmp = 0.0 if (z <= -280000000.0) tmp = t_1; elseif (z <= 3.5e-87) tmp = Float64(x - Float64(Float64(y * Float64(x - t)) / a)); elseif (z <= 800000000.0) tmp = Float64(y * Float64(Float64(t - x) / Float64(a - z))); elseif (z <= 2e+45) tmp = Float64(x + Float64(t * Float64(y / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * ((y - z) / (a - z)); tmp = 0.0; if (z <= -280000000.0) tmp = t_1; elseif (z <= 3.5e-87) tmp = x - ((y * (x - t)) / a); elseif (z <= 800000000.0) tmp = y * ((t - x) / (a - z)); elseif (z <= 2e+45) tmp = x + (t * (y / a)); 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]}, If[LessEqual[z, -280000000.0], t$95$1, If[LessEqual[z, 3.5e-87], N[(x - N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 800000000.0], N[(y * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2e+45], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y - z}{a - z}\\
\mathbf{if}\;z \leq -280000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{-87}:\\
\;\;\;\;x - \frac{y \cdot \left(x - t\right)}{a}\\
\mathbf{elif}\;z \leq 800000000:\\
\;\;\;\;y \cdot \frac{t - x}{a - z}\\
\mathbf{elif}\;z \leq 2 \cdot 10^{+45}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -2.8e8 or 1.9999999999999999e45 < z Initial program 44.5%
associate-/l*64.5%
Simplified64.5%
Taylor expanded in x around 0 40.7%
associate-/l*63.1%
Simplified63.1%
if -2.8e8 < z < 3.50000000000000012e-87Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
if 3.50000000000000012e-87 < z < 8e8Initial program 77.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in y around inf 72.4%
div-sub72.4%
Simplified72.4%
if 8e8 < z < 1.9999999999999999e45Initial program 88.9%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in z around 0 70.1%
Taylor expanded in t around inf 80.9%
associate-/l*81.2%
Simplified81.2%
Final simplification68.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ (- y z) (- a z)))))
(if (<= z -85000000.0)
t_1
(if (<= z 3.3e-87)
(+ x (/ (- t x) (/ a y)))
(if (<= z 46000000.0)
(* y (/ (- t x) (- a z)))
(if (<= z 2.15e+45) (+ x (* t (/ y a))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * ((y - z) / (a - z));
double tmp;
if (z <= -85000000.0) {
tmp = t_1;
} else if (z <= 3.3e-87) {
tmp = x + ((t - x) / (a / y));
} else if (z <= 46000000.0) {
tmp = y * ((t - x) / (a - z));
} else if (z <= 2.15e+45) {
tmp = x + (t * (y / a));
} 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 = t * ((y - z) / (a - z))
if (z <= (-85000000.0d0)) then
tmp = t_1
else if (z <= 3.3d-87) then
tmp = x + ((t - x) / (a / y))
else if (z <= 46000000.0d0) then
tmp = y * ((t - x) / (a - z))
else if (z <= 2.15d+45) then
tmp = x + (t * (y / a))
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 tmp;
if (z <= -85000000.0) {
tmp = t_1;
} else if (z <= 3.3e-87) {
tmp = x + ((t - x) / (a / y));
} else if (z <= 46000000.0) {
tmp = y * ((t - x) / (a - z));
} else if (z <= 2.15e+45) {
tmp = x + (t * (y / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * ((y - z) / (a - z)) tmp = 0 if z <= -85000000.0: tmp = t_1 elif z <= 3.3e-87: tmp = x + ((t - x) / (a / y)) elif z <= 46000000.0: tmp = y * ((t - x) / (a - z)) elif z <= 2.15e+45: tmp = x + (t * (y / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(Float64(y - z) / Float64(a - z))) tmp = 0.0 if (z <= -85000000.0) tmp = t_1; elseif (z <= 3.3e-87) tmp = Float64(x + Float64(Float64(t - x) / Float64(a / y))); elseif (z <= 46000000.0) tmp = Float64(y * Float64(Float64(t - x) / Float64(a - z))); elseif (z <= 2.15e+45) tmp = Float64(x + Float64(t * Float64(y / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * ((y - z) / (a - z)); tmp = 0.0; if (z <= -85000000.0) tmp = t_1; elseif (z <= 3.3e-87) tmp = x + ((t - x) / (a / y)); elseif (z <= 46000000.0) tmp = y * ((t - x) / (a - z)); elseif (z <= 2.15e+45) tmp = x + (t * (y / a)); 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]}, If[LessEqual[z, -85000000.0], t$95$1, If[LessEqual[z, 3.3e-87], N[(x + N[(N[(t - x), $MachinePrecision] / N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 46000000.0], N[(y * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.15e+45], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y - z}{a - z}\\
\mathbf{if}\;z \leq -85000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{-87}:\\
\;\;\;\;x + \frac{t - x}{\frac{a}{y}}\\
\mathbf{elif}\;z \leq 46000000:\\
\;\;\;\;y \cdot \frac{t - x}{a - z}\\
\mathbf{elif}\;z \leq 2.15 \cdot 10^{+45}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -8.5e7 or 2.1500000000000002e45 < z Initial program 44.5%
associate-/l*64.5%
Simplified64.5%
Taylor expanded in x around 0 40.7%
associate-/l*63.1%
Simplified63.1%
if -8.5e7 < z < 3.3e-87Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
*-commutative94.5%
associate-*l/95.3%
associate-*r/97.2%
clear-num97.2%
un-div-inv97.2%
Applied egg-rr97.2%
Taylor expanded in z around 0 76.8%
if 3.3e-87 < z < 4.6e7Initial program 77.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in y around inf 72.4%
div-sub72.4%
Simplified72.4%
if 4.6e7 < z < 2.1500000000000002e45Initial program 88.9%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in z around 0 70.1%
Taylor expanded in t around inf 80.9%
associate-/l*81.2%
Simplified81.2%
Final simplification70.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (/ (- t x) (/ a y)))))
(if (<= a -6e-34)
t_1
(if (<= a 2.9e-64)
(+ t (/ (* y (- x t)) z))
(if (<= a 1.7e-35)
(* y (/ (- t x) (- a z)))
(if (<= a 4.5e+59) (* t (/ (- y z) (- a z))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((t - x) / (a / y));
double tmp;
if (a <= -6e-34) {
tmp = t_1;
} else if (a <= 2.9e-64) {
tmp = t + ((y * (x - t)) / z);
} else if (a <= 1.7e-35) {
tmp = y * ((t - x) / (a - z));
} else if (a <= 4.5e+59) {
tmp = t * ((y - z) / (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 + ((t - x) / (a / y))
if (a <= (-6d-34)) then
tmp = t_1
else if (a <= 2.9d-64) then
tmp = t + ((y * (x - t)) / z)
else if (a <= 1.7d-35) then
tmp = y * ((t - x) / (a - z))
else if (a <= 4.5d+59) then
tmp = t * ((y - z) / (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 + ((t - x) / (a / y));
double tmp;
if (a <= -6e-34) {
tmp = t_1;
} else if (a <= 2.9e-64) {
tmp = t + ((y * (x - t)) / z);
} else if (a <= 1.7e-35) {
tmp = y * ((t - x) / (a - z));
} else if (a <= 4.5e+59) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + ((t - x) / (a / y)) tmp = 0 if a <= -6e-34: tmp = t_1 elif a <= 2.9e-64: tmp = t + ((y * (x - t)) / z) elif a <= 1.7e-35: tmp = y * ((t - x) / (a - z)) elif a <= 4.5e+59: tmp = t * ((y - z) / (a - z)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(t - x) / Float64(a / y))) tmp = 0.0 if (a <= -6e-34) tmp = t_1; elseif (a <= 2.9e-64) tmp = Float64(t + Float64(Float64(y * Float64(x - t)) / z)); elseif (a <= 1.7e-35) tmp = Float64(y * Float64(Float64(t - x) / Float64(a - z))); elseif (a <= 4.5e+59) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + ((t - x) / (a / y)); tmp = 0.0; if (a <= -6e-34) tmp = t_1; elseif (a <= 2.9e-64) tmp = t + ((y * (x - t)) / z); elseif (a <= 1.7e-35) tmp = y * ((t - x) / (a - z)); elseif (a <= 4.5e+59) tmp = t * ((y - z) / (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[(t - x), $MachinePrecision] / N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -6e-34], t$95$1, If[LessEqual[a, 2.9e-64], N[(t + N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.7e-35], N[(y * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 4.5e+59], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{t - x}{\frac{a}{y}}\\
\mathbf{if}\;a \leq -6 \cdot 10^{-34}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2.9 \cdot 10^{-64}:\\
\;\;\;\;t + \frac{y \cdot \left(x - t\right)}{z}\\
\mathbf{elif}\;a \leq 1.7 \cdot 10^{-35}:\\
\;\;\;\;y \cdot \frac{t - x}{a - z}\\
\mathbf{elif}\;a \leq 4.5 \cdot 10^{+59}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -6e-34 or 4.49999999999999959e59 < a Initial program 70.5%
associate-/l*84.6%
Simplified84.6%
*-commutative84.6%
associate-*l/70.5%
associate-*r/87.9%
clear-num87.9%
un-div-inv87.9%
Applied egg-rr87.9%
Taylor expanded in z around 0 65.7%
if -6e-34 < a < 2.8999999999999999e-64Initial program 67.7%
associate-/l*73.4%
Simplified73.4%
Taylor expanded in z around inf 84.5%
associate--l+84.5%
associate-*r/84.5%
associate-*r/84.5%
mul-1-neg84.5%
div-sub85.3%
mul-1-neg85.3%
distribute-lft-out--85.3%
associate-*r/85.3%
mul-1-neg85.3%
unsub-neg85.3%
distribute-rgt-out--85.3%
Simplified85.3%
Taylor expanded in y around inf 83.0%
if 2.8999999999999999e-64 < a < 1.7000000000000001e-35Initial program 78.8%
associate-/l*88.8%
Simplified88.8%
Taylor expanded in y around inf 78.4%
div-sub78.4%
Simplified78.4%
if 1.7000000000000001e-35 < a < 4.49999999999999959e59Initial program 72.4%
associate-/l*83.2%
Simplified83.2%
Taylor expanded in x around 0 56.6%
associate-/l*73.0%
Simplified73.0%
Final simplification74.7%
(FPCore (x y z t a)
:precision binary64
(if (<= z -200000000.0)
t
(if (<= z 3.8e-87)
(+ x (* t (/ y a)))
(if (<= z 4.2e-23) (* x (/ y z)) (if (<= z 3.5e+116) (+ t x) t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -200000000.0) {
tmp = t;
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = x * (y / z);
} else if (z <= 3.5e+116) {
tmp = t + x;
} else {
tmp = t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-200000000.0d0)) then
tmp = t
else if (z <= 3.8d-87) then
tmp = x + (t * (y / a))
else if (z <= 4.2d-23) then
tmp = x * (y / z)
else if (z <= 3.5d+116) then
tmp = t + x
else
tmp = t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -200000000.0) {
tmp = t;
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = x * (y / z);
} else if (z <= 3.5e+116) {
tmp = t + x;
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -200000000.0: tmp = t elif z <= 3.8e-87: tmp = x + (t * (y / a)) elif z <= 4.2e-23: tmp = x * (y / z) elif z <= 3.5e+116: tmp = t + x else: tmp = t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -200000000.0) tmp = t; elseif (z <= 3.8e-87) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (z <= 4.2e-23) tmp = Float64(x * Float64(y / z)); elseif (z <= 3.5e+116) tmp = Float64(t + x); else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -200000000.0) tmp = t; elseif (z <= 3.8e-87) tmp = x + (t * (y / a)); elseif (z <= 4.2e-23) tmp = x * (y / z); elseif (z <= 3.5e+116) tmp = t + x; else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -200000000.0], t, If[LessEqual[z, 3.8e-87], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.2e-23], N[(x * N[(y / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.5e+116], N[(t + x), $MachinePrecision], t]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -200000000:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{-87}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq 4.2 \cdot 10^{-23}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{+116}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -2e8 or 3.49999999999999997e116 < z Initial program 39.9%
associate-/l*61.5%
Simplified61.5%
Taylor expanded in z around inf 51.3%
if -2e8 < z < 3.8e-87Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
Taylor expanded in t around inf 61.6%
associate-/l*67.9%
Simplified67.9%
if 3.8e-87 < z < 4.2000000000000002e-23Initial program 82.1%
associate-/l*81.6%
Simplified81.6%
Taylor expanded in z around inf 75.9%
associate--l+75.9%
associate-*r/75.9%
associate-*r/75.9%
mul-1-neg75.9%
div-sub75.9%
mul-1-neg75.9%
distribute-lft-out--75.9%
associate-*r/75.9%
mul-1-neg75.9%
unsub-neg75.9%
distribute-rgt-out--82.2%
Simplified82.2%
Taylor expanded in t around 0 58.0%
associate-/l*58.1%
Simplified58.1%
Taylor expanded in y around inf 58.6%
associate-/l*58.7%
Simplified58.7%
if 4.2000000000000002e-23 < z < 3.49999999999999997e116Initial program 79.7%
associate-/l*90.5%
Simplified90.5%
Taylor expanded in t around inf 76.7%
Taylor expanded in z around inf 48.4%
Final simplification58.3%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -6.1e-65) (not (<= a 2.35e-64))) (+ x (* (- y z) (/ (- t x) (- a z)))) (+ t (/ (* (- t x) (- a y)) z))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -6.1e-65) || !(a <= 2.35e-64)) {
tmp = x + ((y - z) * ((t - x) / (a - z)));
} else {
tmp = t + (((t - x) * (a - y)) / z);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a <= (-6.1d-65)) .or. (.not. (a <= 2.35d-64))) then
tmp = x + ((y - z) * ((t - x) / (a - z)))
else
tmp = t + (((t - x) * (a - y)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -6.1e-65) || !(a <= 2.35e-64)) {
tmp = x + ((y - z) * ((t - x) / (a - z)));
} else {
tmp = t + (((t - x) * (a - y)) / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -6.1e-65) or not (a <= 2.35e-64): tmp = x + ((y - z) * ((t - x) / (a - z))) else: tmp = t + (((t - x) * (a - y)) / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -6.1e-65) || !(a <= 2.35e-64)) tmp = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))); else tmp = Float64(t + Float64(Float64(Float64(t - x) * Float64(a - y)) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -6.1e-65) || ~((a <= 2.35e-64))) tmp = x + ((y - z) * ((t - x) / (a - z))); else tmp = t + (((t - x) * (a - y)) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -6.1e-65], N[Not[LessEqual[a, 2.35e-64]], $MachinePrecision]], N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t + N[(N[(N[(t - x), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.1 \cdot 10^{-65} \lor \neg \left(a \leq 2.35 \cdot 10^{-64}\right):\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t + \frac{\left(t - x\right) \cdot \left(a - y\right)}{z}\\
\end{array}
\end{array}
if a < -6.10000000000000014e-65 or 2.3499999999999999e-64 < a Initial program 72.1%
associate-/l*85.1%
Simplified85.1%
if -6.10000000000000014e-65 < a < 2.3499999999999999e-64Initial program 66.6%
associate-/l*72.5%
Simplified72.5%
Taylor expanded in z around inf 85.6%
associate--l+85.6%
associate-*r/85.6%
associate-*r/85.6%
mul-1-neg85.6%
div-sub86.5%
mul-1-neg86.5%
distribute-lft-out--86.5%
associate-*r/86.5%
mul-1-neg86.5%
unsub-neg86.5%
distribute-rgt-out--86.5%
Simplified86.5%
Final simplification85.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ y (- a z)))))
(if (<= y -1.15e+181)
t_1
(if (<= y -8.5e+83) (* x (/ y z)) (if (<= y 1.95e+35) (+ t x) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (y / (a - z));
double tmp;
if (y <= -1.15e+181) {
tmp = t_1;
} else if (y <= -8.5e+83) {
tmp = x * (y / z);
} else if (y <= 1.95e+35) {
tmp = t + x;
} 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 = t * (y / (a - z))
if (y <= (-1.15d+181)) then
tmp = t_1
else if (y <= (-8.5d+83)) then
tmp = x * (y / z)
else if (y <= 1.95d+35) then
tmp = t + x
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 / (a - z));
double tmp;
if (y <= -1.15e+181) {
tmp = t_1;
} else if (y <= -8.5e+83) {
tmp = x * (y / z);
} else if (y <= 1.95e+35) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (y / (a - z)) tmp = 0 if y <= -1.15e+181: tmp = t_1 elif y <= -8.5e+83: tmp = x * (y / z) elif y <= 1.95e+35: tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -1.15e+181) tmp = t_1; elseif (y <= -8.5e+83) tmp = Float64(x * Float64(y / z)); elseif (y <= 1.95e+35) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (y / (a - z)); tmp = 0.0; if (y <= -1.15e+181) tmp = t_1; elseif (y <= -8.5e+83) tmp = x * (y / z); elseif (y <= 1.95e+35) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.15e+181], t$95$1, If[LessEqual[y, -8.5e+83], N[(x * N[(y / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.95e+35], N[(t + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -1.15 \cdot 10^{+181}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -8.5 \cdot 10^{+83}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq 1.95 \cdot 10^{+35}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.1499999999999999e181 or 1.95e35 < y Initial program 76.0%
associate-/l*85.2%
Simplified85.2%
Taylor expanded in x around 0 53.7%
associate-/l*66.5%
Simplified66.5%
Taylor expanded in y around inf 49.4%
associate-/l*58.6%
Simplified58.6%
if -1.1499999999999999e181 < y < -8.4999999999999995e83Initial program 68.2%
associate-/l*83.9%
Simplified83.9%
Taylor expanded in z around inf 63.9%
associate--l+63.9%
associate-*r/63.9%
associate-*r/63.9%
mul-1-neg63.9%
div-sub64.1%
mul-1-neg64.1%
distribute-lft-out--64.1%
associate-*r/64.1%
mul-1-neg64.1%
unsub-neg64.1%
distribute-rgt-out--64.1%
Simplified64.1%
Taylor expanded in t around 0 40.1%
associate-/l*51.6%
Simplified51.6%
Taylor expanded in y around inf 40.4%
associate-/l*48.4%
Simplified48.4%
if -8.4999999999999995e83 < y < 1.95e35Initial program 66.4%
associate-/l*75.7%
Simplified75.7%
Taylor expanded in t around inf 68.5%
Taylor expanded in z around inf 52.7%
Final simplification54.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ y (- a z)))))
(if (<= y -4e+182)
t_1
(if (<= y -8.5e+83)
(* x (/ (- y a) z))
(if (<= y 2.5e+35) (+ t x) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (y / (a - z));
double tmp;
if (y <= -4e+182) {
tmp = t_1;
} else if (y <= -8.5e+83) {
tmp = x * ((y - a) / z);
} else if (y <= 2.5e+35) {
tmp = t + x;
} 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 = t * (y / (a - z))
if (y <= (-4d+182)) then
tmp = t_1
else if (y <= (-8.5d+83)) then
tmp = x * ((y - a) / z)
else if (y <= 2.5d+35) then
tmp = t + x
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 / (a - z));
double tmp;
if (y <= -4e+182) {
tmp = t_1;
} else if (y <= -8.5e+83) {
tmp = x * ((y - a) / z);
} else if (y <= 2.5e+35) {
tmp = t + x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (y / (a - z)) tmp = 0 if y <= -4e+182: tmp = t_1 elif y <= -8.5e+83: tmp = x * ((y - a) / z) elif y <= 2.5e+35: tmp = t + x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(y / Float64(a - z))) tmp = 0.0 if (y <= -4e+182) tmp = t_1; elseif (y <= -8.5e+83) tmp = Float64(x * Float64(Float64(y - a) / z)); elseif (y <= 2.5e+35) tmp = Float64(t + x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (y / (a - z)); tmp = 0.0; if (y <= -4e+182) tmp = t_1; elseif (y <= -8.5e+83) tmp = x * ((y - a) / z); elseif (y <= 2.5e+35) tmp = t + x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+182], t$95$1, If[LessEqual[y, -8.5e+83], N[(x * N[(N[(y - a), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.5e+35], N[(t + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y}{a - z}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+182}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -8.5 \cdot 10^{+83}:\\
\;\;\;\;x \cdot \frac{y - a}{z}\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{+35}:\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -4.0000000000000003e182 or 2.50000000000000011e35 < y Initial program 76.0%
associate-/l*85.2%
Simplified85.2%
Taylor expanded in x around 0 53.7%
associate-/l*66.5%
Simplified66.5%
Taylor expanded in y around inf 49.4%
associate-/l*58.6%
Simplified58.6%
if -4.0000000000000003e182 < y < -8.4999999999999995e83Initial program 68.2%
associate-/l*83.9%
Simplified83.9%
Taylor expanded in z around inf 63.9%
associate--l+63.9%
associate-*r/63.9%
associate-*r/63.9%
mul-1-neg63.9%
div-sub64.1%
mul-1-neg64.1%
distribute-lft-out--64.1%
associate-*r/64.1%
mul-1-neg64.1%
unsub-neg64.1%
distribute-rgt-out--64.1%
Simplified64.1%
Taylor expanded in t around 0 40.1%
associate-/l*51.6%
Simplified51.6%
if -8.4999999999999995e83 < y < 2.50000000000000011e35Initial program 66.4%
associate-/l*75.7%
Simplified75.7%
Taylor expanded in t around inf 68.5%
Taylor expanded in z around inf 52.7%
Final simplification54.5%
(FPCore (x y z t a)
:precision binary64
(if (<= z -700000000.0)
t
(if (<= z 9.5e-88)
(+ x (* t (/ y a)))
(if (<= z 4.2e-23) (/ (* x y) (- z a)) (+ t x)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -700000000.0) {
tmp = t;
} else if (z <= 9.5e-88) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = (x * y) / (z - a);
} else {
tmp = 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) :: tmp
if (z <= (-700000000.0d0)) then
tmp = t
else if (z <= 9.5d-88) then
tmp = x + (t * (y / a))
else if (z <= 4.2d-23) then
tmp = (x * y) / (z - a)
else
tmp = t + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -700000000.0) {
tmp = t;
} else if (z <= 9.5e-88) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = (x * y) / (z - a);
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -700000000.0: tmp = t elif z <= 9.5e-88: tmp = x + (t * (y / a)) elif z <= 4.2e-23: tmp = (x * y) / (z - a) else: tmp = t + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -700000000.0) tmp = t; elseif (z <= 9.5e-88) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (z <= 4.2e-23) tmp = Float64(Float64(x * y) / Float64(z - a)); else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -700000000.0) tmp = t; elseif (z <= 9.5e-88) tmp = x + (t * (y / a)); elseif (z <= 4.2e-23) tmp = (x * y) / (z - a); else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -700000000.0], t, If[LessEqual[z, 9.5e-88], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.2e-23], N[(N[(x * y), $MachinePrecision] / N[(z - a), $MachinePrecision]), $MachinePrecision], N[(t + x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -700000000:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq 9.5 \cdot 10^{-88}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq 4.2 \cdot 10^{-23}:\\
\;\;\;\;\frac{x \cdot y}{z - a}\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if z < -7e8Initial program 47.2%
associate-/l*64.0%
Simplified64.0%
Taylor expanded in z around inf 44.0%
if -7e8 < z < 9.5e-88Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
Taylor expanded in t around inf 61.6%
associate-/l*67.9%
Simplified67.9%
if 9.5e-88 < z < 4.2000000000000002e-23Initial program 82.1%
+-commutative82.1%
associate-/l*81.6%
fma-define81.4%
Simplified81.4%
Taylor expanded in t around 0 57.3%
neg-mul-157.3%
distribute-neg-frac257.3%
Simplified57.3%
Taylor expanded in y around inf 69.3%
*-commutative69.3%
Simplified69.3%
if 4.2000000000000002e-23 < z Initial program 48.7%
associate-/l*69.7%
Simplified69.7%
Taylor expanded in t around inf 63.5%
Taylor expanded in z around inf 52.6%
Final simplification58.2%
(FPCore (x y z t a)
:precision binary64
(if (<= z -210000000.0)
(* t (/ z (- z a)))
(if (<= z 3.8e-87)
(+ x (* t (/ y a)))
(if (<= z 4.2e-23) (/ (* x y) (- z a)) (+ t x)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -210000000.0) {
tmp = t * (z / (z - a));
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = (x * y) / (z - a);
} else {
tmp = 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) :: tmp
if (z <= (-210000000.0d0)) then
tmp = t * (z / (z - a))
else if (z <= 3.8d-87) then
tmp = x + (t * (y / a))
else if (z <= 4.2d-23) then
tmp = (x * y) / (z - a)
else
tmp = t + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -210000000.0) {
tmp = t * (z / (z - a));
} else if (z <= 3.8e-87) {
tmp = x + (t * (y / a));
} else if (z <= 4.2e-23) {
tmp = (x * y) / (z - a);
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -210000000.0: tmp = t * (z / (z - a)) elif z <= 3.8e-87: tmp = x + (t * (y / a)) elif z <= 4.2e-23: tmp = (x * y) / (z - a) else: tmp = t + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -210000000.0) tmp = Float64(t * Float64(z / Float64(z - a))); elseif (z <= 3.8e-87) tmp = Float64(x + Float64(t * Float64(y / a))); elseif (z <= 4.2e-23) tmp = Float64(Float64(x * y) / Float64(z - a)); else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -210000000.0) tmp = t * (z / (z - a)); elseif (z <= 3.8e-87) tmp = x + (t * (y / a)); elseif (z <= 4.2e-23) tmp = (x * y) / (z - a); else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -210000000.0], N[(t * N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.8e-87], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.2e-23], N[(N[(x * y), $MachinePrecision] / N[(z - a), $MachinePrecision]), $MachinePrecision], N[(t + x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -210000000:\\
\;\;\;\;t \cdot \frac{z}{z - a}\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{-87}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq 4.2 \cdot 10^{-23}:\\
\;\;\;\;\frac{x \cdot y}{z - a}\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if z < -2.1e8Initial program 47.2%
associate-/l*64.0%
Simplified64.0%
Taylor expanded in x around 0 45.9%
associate-/l*65.2%
Simplified65.2%
Taylor expanded in y around 0 36.8%
mul-1-neg36.8%
associate-/l*51.0%
distribute-rgt-neg-in51.0%
distribute-frac-neg251.0%
Simplified51.0%
if -2.1e8 < z < 3.8e-87Initial program 95.3%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in z around 0 73.1%
Taylor expanded in t around inf 61.6%
associate-/l*67.9%
Simplified67.9%
if 3.8e-87 < z < 4.2000000000000002e-23Initial program 82.1%
+-commutative82.1%
associate-/l*81.6%
fma-define81.4%
Simplified81.4%
Taylor expanded in t around 0 57.3%
neg-mul-157.3%
distribute-neg-frac257.3%
Simplified57.3%
Taylor expanded in y around inf 69.3%
*-commutative69.3%
Simplified69.3%
if 4.2000000000000002e-23 < z Initial program 48.7%
associate-/l*69.7%
Simplified69.7%
Taylor expanded in t around inf 63.5%
Taylor expanded in z around inf 52.6%
Final simplification59.6%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -3.1e-40) (not (<= a 7.8e-63))) (+ x (* (- y z) (/ t (- a z)))) (+ t (/ (* y (- x t)) z))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.1e-40) || !(a <= 7.8e-63)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t + ((y * (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) :: tmp
if ((a <= (-3.1d-40)) .or. (.not. (a <= 7.8d-63))) then
tmp = x + ((y - z) * (t / (a - z)))
else
tmp = t + ((y * (x - t)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.1e-40) || !(a <= 7.8e-63)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t + ((y * (x - t)) / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -3.1e-40) or not (a <= 7.8e-63): tmp = x + ((y - z) * (t / (a - z))) else: tmp = t + ((y * (x - t)) / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -3.1e-40) || !(a <= 7.8e-63)) tmp = Float64(x + Float64(Float64(y - z) * Float64(t / Float64(a - z)))); else tmp = Float64(t + Float64(Float64(y * Float64(x - t)) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -3.1e-40) || ~((a <= 7.8e-63))) tmp = x + ((y - z) * (t / (a - z))); else tmp = t + ((y * (x - t)) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -3.1e-40], N[Not[LessEqual[a, 7.8e-63]], $MachinePrecision]], N[(x + N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t + N[(N[(y * N[(x - t), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.1 \cdot 10^{-40} \lor \neg \left(a \leq 7.8 \cdot 10^{-63}\right):\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t + \frac{y \cdot \left(x - t\right)}{z}\\
\end{array}
\end{array}
if a < -3.10000000000000011e-40 or 7.80000000000000044e-63 < a Initial program 71.5%
associate-/l*84.8%
Simplified84.8%
Taylor expanded in t around inf 76.0%
if -3.10000000000000011e-40 < a < 7.80000000000000044e-63Initial program 67.4%
associate-/l*73.2%
Simplified73.2%
Taylor expanded in z around inf 85.1%
associate--l+85.1%
associate-*r/85.1%
associate-*r/85.1%
mul-1-neg85.1%
div-sub86.0%
mul-1-neg86.0%
distribute-lft-out--86.0%
associate-*r/86.0%
mul-1-neg86.0%
unsub-neg86.0%
distribute-rgt-out--86.0%
Simplified86.0%
Taylor expanded in y around inf 83.6%
Final simplification79.5%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -3.35e-37) (not (<= a 1.08e-61))) (+ x (* (- y z) (/ t (- a z)))) (+ t (/ (* (- t x) (- a y)) z))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.35e-37) || !(a <= 1.08e-61)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t + (((t - x) * (a - y)) / z);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a <= (-3.35d-37)) .or. (.not. (a <= 1.08d-61))) then
tmp = x + ((y - z) * (t / (a - z)))
else
tmp = t + (((t - x) * (a - y)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -3.35e-37) || !(a <= 1.08e-61)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t + (((t - x) * (a - y)) / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -3.35e-37) or not (a <= 1.08e-61): tmp = x + ((y - z) * (t / (a - z))) else: tmp = t + (((t - x) * (a - y)) / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -3.35e-37) || !(a <= 1.08e-61)) tmp = Float64(x + Float64(Float64(y - z) * Float64(t / Float64(a - z)))); else tmp = Float64(t + Float64(Float64(Float64(t - x) * Float64(a - y)) / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -3.35e-37) || ~((a <= 1.08e-61))) tmp = x + ((y - z) * (t / (a - z))); else tmp = t + (((t - x) * (a - y)) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -3.35e-37], N[Not[LessEqual[a, 1.08e-61]], $MachinePrecision]], N[(x + N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t + N[(N[(N[(t - x), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.35 \cdot 10^{-37} \lor \neg \left(a \leq 1.08 \cdot 10^{-61}\right):\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t + \frac{\left(t - x\right) \cdot \left(a - y\right)}{z}\\
\end{array}
\end{array}
if a < -3.35e-37 or 1.08000000000000008e-61 < a Initial program 71.5%
associate-/l*84.8%
Simplified84.8%
Taylor expanded in t around inf 76.0%
if -3.35e-37 < a < 1.08000000000000008e-61Initial program 67.4%
associate-/l*73.2%
Simplified73.2%
Taylor expanded in z around inf 85.1%
associate--l+85.1%
associate-*r/85.1%
associate-*r/85.1%
mul-1-neg85.1%
div-sub86.0%
mul-1-neg86.0%
distribute-lft-out--86.0%
associate-*r/86.0%
mul-1-neg86.0%
unsub-neg86.0%
distribute-rgt-out--86.0%
Simplified86.0%
Final simplification80.6%
(FPCore (x y z t a) :precision binary64 (if (<= a -5.9e+171) x (if (or (<= a -7.5e-80) (not (<= a 3.6e-63))) (+ t x) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -5.9e+171) {
tmp = x;
} else if ((a <= -7.5e-80) || !(a <= 3.6e-63)) {
tmp = t + x;
} else {
tmp = t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-5.9d+171)) then
tmp = x
else if ((a <= (-7.5d-80)) .or. (.not. (a <= 3.6d-63))) then
tmp = t + x
else
tmp = t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -5.9e+171) {
tmp = x;
} else if ((a <= -7.5e-80) || !(a <= 3.6e-63)) {
tmp = t + x;
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -5.9e+171: tmp = x elif (a <= -7.5e-80) or not (a <= 3.6e-63): tmp = t + x else: tmp = t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -5.9e+171) tmp = x; elseif ((a <= -7.5e-80) || !(a <= 3.6e-63)) tmp = Float64(t + x); else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -5.9e+171) tmp = x; elseif ((a <= -7.5e-80) || ~((a <= 3.6e-63))) tmp = t + x; else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -5.9e+171], x, If[Or[LessEqual[a, -7.5e-80], N[Not[LessEqual[a, 3.6e-63]], $MachinePrecision]], N[(t + x), $MachinePrecision], t]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.9 \cdot 10^{+171}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq -7.5 \cdot 10^{-80} \lor \neg \left(a \leq 3.6 \cdot 10^{-63}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if a < -5.90000000000000035e171Initial program 67.7%
associate-/l*95.8%
Simplified95.8%
Taylor expanded in a around inf 67.8%
if -5.90000000000000035e171 < a < -7.49999999999999999e-80 or 3.60000000000000008e-63 < a Initial program 72.7%
associate-/l*82.2%
Simplified82.2%
Taylor expanded in t around inf 70.0%
Taylor expanded in z around inf 43.0%
if -7.49999999999999999e-80 < a < 3.60000000000000008e-63Initial program 66.6%
associate-/l*72.8%
Simplified72.8%
Taylor expanded in z around inf 42.1%
Final simplification44.9%
(FPCore (x y z t a) :precision binary64 (if (or (<= y -2.4e+100) (not (<= y 2.5e+35))) (* t (/ y a)) (+ t x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y <= -2.4e+100) || !(y <= 2.5e+35)) {
tmp = t * (y / a);
} else {
tmp = 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) :: tmp
if ((y <= (-2.4d+100)) .or. (.not. (y <= 2.5d+35))) then
tmp = t * (y / a)
else
tmp = t + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y <= -2.4e+100) || !(y <= 2.5e+35)) {
tmp = t * (y / a);
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (y <= -2.4e+100) or not (y <= 2.5e+35): tmp = t * (y / a) else: tmp = t + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((y <= -2.4e+100) || !(y <= 2.5e+35)) tmp = Float64(t * Float64(y / a)); else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((y <= -2.4e+100) || ~((y <= 2.5e+35))) tmp = t * (y / a); else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[y, -2.4e+100], N[Not[LessEqual[y, 2.5e+35]], $MachinePrecision]], N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(t + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.4 \cdot 10^{+100} \lor \neg \left(y \leq 2.5 \cdot 10^{+35}\right):\\
\;\;\;\;t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if y < -2.40000000000000012e100 or 2.50000000000000011e35 < y Initial program 76.2%
associate-/l*86.8%
Simplified86.8%
Taylor expanded in x around 0 48.9%
associate-/l*60.4%
Simplified60.4%
Taylor expanded in z around 0 35.9%
if -2.40000000000000012e100 < y < 2.50000000000000011e35Initial program 65.5%
associate-/l*75.0%
Simplified75.0%
Taylor expanded in t around inf 66.9%
Taylor expanded in z around inf 51.8%
Final simplification45.8%
(FPCore (x y z t a) :precision binary64 (if (<= a -4.8e-33) x (if (<= a 4.3e+85) t x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -4.8e-33) {
tmp = x;
} else if (a <= 4.3e+85) {
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 <= (-4.8d-33)) then
tmp = x
else if (a <= 4.3d+85) 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 <= -4.8e-33) {
tmp = x;
} else if (a <= 4.3e+85) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -4.8e-33: tmp = x elif a <= 4.3e+85: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -4.8e-33) tmp = x; elseif (a <= 4.3e+85) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -4.8e-33) tmp = x; elseif (a <= 4.3e+85) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -4.8e-33], x, If[LessEqual[a, 4.3e+85], t, x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4.8 \cdot 10^{-33}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 4.3 \cdot 10^{+85}:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -4.8e-33 or 4.2999999999999999e85 < a Initial program 70.5%
associate-/l*86.0%
Simplified86.0%
Taylor expanded in a around inf 45.3%
if -4.8e-33 < a < 4.2999999999999999e85Initial program 68.9%
associate-/l*74.7%
Simplified74.7%
Taylor expanded in z around inf 38.7%
Final simplification41.5%
(FPCore (x y z t a) :precision binary64 0.0)
double code(double x, double y, double z, double t, double a) {
return 0.0;
}
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 = 0.0d0
end function
public static double code(double x, double y, double z, double t, double a) {
return 0.0;
}
def code(x, y, z, t, a): return 0.0
function code(x, y, z, t, a) return 0.0 end
function tmp = code(x, y, z, t, a) tmp = 0.0; end
code[x_, y_, z_, t_, a_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 69.6%
+-commutative69.6%
associate-/l*79.5%
fma-define79.6%
Simplified79.6%
Taylor expanded in t around 0 39.2%
neg-mul-139.2%
distribute-neg-frac239.2%
Simplified39.2%
Taylor expanded in z around inf 2.7%
distribute-rgt1-in2.7%
metadata-eval2.7%
mul0-lft2.7%
Simplified2.7%
Final simplification2.7%
(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 69.6%
associate-/l*79.5%
Simplified79.5%
Taylor expanded in z around inf 28.7%
Final simplification28.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- t (* (/ y z) (- t x)))))
(if (< z -1.2536131056095036e+188)
t_1
(if (< z 4.446702369113811e+64)
(+ x (/ (- y z) (/ (- a z) (- t x))))
t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t - ((y / z) * (t - x));
double tmp;
if (z < -1.2536131056095036e+188) {
tmp = t_1;
} else if (z < 4.446702369113811e+64) {
tmp = x + ((y - z) / ((a - z) / (t - x)));
} 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 = t - ((y / z) * (t - x))
if (z < (-1.2536131056095036d+188)) then
tmp = t_1
else if (z < 4.446702369113811d+64) then
tmp = x + ((y - z) / ((a - z) / (t - x)))
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) * (t - x));
double tmp;
if (z < -1.2536131056095036e+188) {
tmp = t_1;
} else if (z < 4.446702369113811e+64) {
tmp = x + ((y - z) / ((a - z) / (t - x)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t - ((y / z) * (t - x)) tmp = 0 if z < -1.2536131056095036e+188: tmp = t_1 elif z < 4.446702369113811e+64: tmp = x + ((y - z) / ((a - z) / (t - x))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(t - Float64(Float64(y / z) * Float64(t - x))) tmp = 0.0 if (z < -1.2536131056095036e+188) tmp = t_1; elseif (z < 4.446702369113811e+64) tmp = Float64(x + Float64(Float64(y - z) / Float64(Float64(a - z) / Float64(t - x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t - ((y / z) * (t - x)); tmp = 0.0; if (z < -1.2536131056095036e+188) tmp = t_1; elseif (z < 4.446702369113811e+64) tmp = x + ((y - z) / ((a - z) / (t - x))); 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[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[z, -1.2536131056095036e+188], t$95$1, If[Less[z, 4.446702369113811e+64], N[(x + N[(N[(y - z), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t - \frac{y}{z} \cdot \left(t - x\right)\\
\mathbf{if}\;z < -1.2536131056095036 \cdot 10^{+188}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z < 4.446702369113811 \cdot 10^{+64}:\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t - x}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
herbie shell --seed 2024095
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:alt
(if (< z -1.2536131056095036e+188) (- t (* (/ y z) (- t x))) (if (< z 4.446702369113811e+64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))