
(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 22 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
(let* ((t_1 (+ x (/ (* (- y z) (- t x)) (- a z)))))
(if (<= t_1 -2e-294)
(+ x (* (- t x) (/ (- y z) (- a z))))
(if (<= t_1 0.0)
(- t (* x (/ (- a y) z)))
(if (<= t_1 1e+303) t_1 (- t (/ (- t x) (/ z (- y a)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((y - z) * (t - x)) / (a - z));
double tmp;
if (t_1 <= -2e-294) {
tmp = x + ((t - x) * ((y - z) / (a - z)));
} else if (t_1 <= 0.0) {
tmp = t - (x * ((a - y) / z));
} else if (t_1 <= 1e+303) {
tmp = t_1;
} else {
tmp = t - ((t - x) / (z / (y - a)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x + (((y - z) * (t - x)) / (a - z))
if (t_1 <= (-2d-294)) then
tmp = x + ((t - x) * ((y - z) / (a - z)))
else if (t_1 <= 0.0d0) then
tmp = t - (x * ((a - y) / z))
else if (t_1 <= 1d+303) then
tmp = t_1
else
tmp = t - ((t - x) / (z / (y - a)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((y - z) * (t - x)) / (a - z));
double tmp;
if (t_1 <= -2e-294) {
tmp = x + ((t - x) * ((y - z) / (a - z)));
} else if (t_1 <= 0.0) {
tmp = t - (x * ((a - y) / z));
} else if (t_1 <= 1e+303) {
tmp = t_1;
} else {
tmp = t - ((t - x) / (z / (y - a)));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (((y - z) * (t - x)) / (a - z)) tmp = 0 if t_1 <= -2e-294: tmp = x + ((t - x) * ((y - z) / (a - z))) elif t_1 <= 0.0: tmp = t - (x * ((a - y) / z)) elif t_1 <= 1e+303: tmp = t_1 else: tmp = t - ((t - x) / (z / (y - a))) return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(Float64(y - z) * Float64(t - x)) / Float64(a - z))) tmp = 0.0 if (t_1 <= -2e-294) tmp = Float64(x + Float64(Float64(t - x) * Float64(Float64(y - z) / Float64(a - z)))); elseif (t_1 <= 0.0) tmp = Float64(t - Float64(x * Float64(Float64(a - y) / z))); elseif (t_1 <= 1e+303) tmp = t_1; else tmp = Float64(t - Float64(Float64(t - x) / Float64(z / Float64(y - a)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (((y - z) * (t - x)) / (a - z)); tmp = 0.0; if (t_1 <= -2e-294) tmp = x + ((t - x) * ((y - z) / (a - z))); elseif (t_1 <= 0.0) tmp = t - (x * ((a - y) / z)); elseif (t_1 <= 1e+303) tmp = t_1; else tmp = t - ((t - x) / (z / (y - a))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(N[(y - z), $MachinePrecision] * N[(t - x), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e-294], N[(x + N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(t - N[(x * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+303], t$95$1, N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / N[(y - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\\
\mathbf{if}\;t_1 \leq -2 \cdot 10^{-294}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y - z}{a - z}\\
\mathbf{elif}\;t_1 \leq 0:\\
\;\;\;\;t - x \cdot \frac{a - y}{z}\\
\mathbf{elif}\;t_1 \leq 10^{+303}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y - a}}\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -2.00000000000000003e-294Initial program 63.9%
associate-*l/90.1%
Simplified90.1%
if -2.00000000000000003e-294 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 0.0Initial program 3.9%
associate-*l/3.9%
Simplified3.9%
Taylor expanded in z around inf 99.5%
associate--l+99.5%
associate-*r/99.5%
associate-*r/99.5%
div-sub99.5%
distribute-lft-out--99.5%
associate-*r/99.5%
distribute-rgt-out--99.5%
mul-1-neg99.5%
unsub-neg99.5%
associate-/l*99.6%
Simplified99.6%
Taylor expanded in t around 0 99.5%
mul-1-neg99.5%
associate-*r/100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
Simplified100.0%
if 0.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 1e303Initial program 97.6%
if 1e303 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) Initial program 35.1%
associate-*l/63.7%
Simplified63.7%
Taylor expanded in z around inf 64.2%
associate--l+64.2%
associate-*r/64.2%
associate-*r/64.2%
div-sub64.2%
distribute-lft-out--64.2%
associate-*r/64.2%
distribute-rgt-out--64.4%
mul-1-neg64.4%
unsub-neg64.4%
associate-/l*83.9%
Simplified83.9%
Final simplification92.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- t (/ y (/ z (- t x))))))
(if (<= z -9.2e+83)
t_1
(if (<= z -2.1e+21)
(+ x (* (- t x) (/ y a)))
(if (<= z -0.17)
(* t (/ (- y z) (- a z)))
(if (<= z -6.1e-6)
(+ x (* t (/ y a)))
(if (or (<= z -2.55e-76) (not (<= z 8.2e-18)))
t_1
(+ x (/ y (/ a (- t x)))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t - (y / (z / (t - x)));
double tmp;
if (z <= -9.2e+83) {
tmp = t_1;
} else if (z <= -2.1e+21) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.17) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -6.1e-6) {
tmp = x + (t * (y / a));
} else if ((z <= -2.55e-76) || !(z <= 8.2e-18)) {
tmp = t_1;
} else {
tmp = x + (y / (a / (t - x)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = t - (y / (z / (t - x)))
if (z <= (-9.2d+83)) then
tmp = t_1
else if (z <= (-2.1d+21)) then
tmp = x + ((t - x) * (y / a))
else if (z <= (-0.17d0)) then
tmp = t * ((y - z) / (a - z))
else if (z <= (-6.1d-6)) then
tmp = x + (t * (y / a))
else if ((z <= (-2.55d-76)) .or. (.not. (z <= 8.2d-18))) then
tmp = t_1
else
tmp = x + (y / (a / (t - x)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = t - (y / (z / (t - x)));
double tmp;
if (z <= -9.2e+83) {
tmp = t_1;
} else if (z <= -2.1e+21) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.17) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -6.1e-6) {
tmp = x + (t * (y / a));
} else if ((z <= -2.55e-76) || !(z <= 8.2e-18)) {
tmp = t_1;
} else {
tmp = x + (y / (a / (t - x)));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t - (y / (z / (t - x))) tmp = 0 if z <= -9.2e+83: tmp = t_1 elif z <= -2.1e+21: tmp = x + ((t - x) * (y / a)) elif z <= -0.17: tmp = t * ((y - z) / (a - z)) elif z <= -6.1e-6: tmp = x + (t * (y / a)) elif (z <= -2.55e-76) or not (z <= 8.2e-18): tmp = t_1 else: tmp = x + (y / (a / (t - x))) return tmp
function code(x, y, z, t, a) t_1 = Float64(t - Float64(y / Float64(z / Float64(t - x)))) tmp = 0.0 if (z <= -9.2e+83) tmp = t_1; elseif (z <= -2.1e+21) tmp = Float64(x + Float64(Float64(t - x) * Float64(y / a))); elseif (z <= -0.17) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); elseif (z <= -6.1e-6) tmp = Float64(x + Float64(t * Float64(y / a))); elseif ((z <= -2.55e-76) || !(z <= 8.2e-18)) tmp = t_1; else tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); 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 <= -9.2e+83) tmp = t_1; elseif (z <= -2.1e+21) tmp = x + ((t - x) * (y / a)); elseif (z <= -0.17) tmp = t * ((y - z) / (a - z)); elseif (z <= -6.1e-6) tmp = x + (t * (y / a)); elseif ((z <= -2.55e-76) || ~((z <= 8.2e-18))) tmp = t_1; else tmp = x + (y / (a / (t - x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t - N[(y / N[(z / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -9.2e+83], t$95$1, If[LessEqual[z, -2.1e+21], N[(x + N[(N[(t - x), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.17], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -6.1e-6], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -2.55e-76], N[Not[LessEqual[z, 8.2e-18]], $MachinePrecision]], t$95$1, N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t - \frac{y}{\frac{z}{t - x}}\\
\mathbf{if}\;z \leq -9.2 \cdot 10^{+83}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -2.1 \cdot 10^{+21}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -0.17:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{elif}\;z \leq -6.1 \cdot 10^{-6}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -2.55 \cdot 10^{-76} \lor \neg \left(z \leq 8.2 \cdot 10^{-18}\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\end{array}
\end{array}
if z < -9.1999999999999998e83 or -6.10000000000000004e-6 < z < -2.54999999999999993e-76 or 8.1999999999999995e-18 < z Initial program 40.1%
associate-*l/65.6%
Simplified65.6%
Taylor expanded in z around inf 68.1%
associate--l+68.1%
associate-*r/68.1%
associate-*r/68.1%
div-sub68.1%
distribute-lft-out--68.1%
associate-*r/68.1%
distribute-rgt-out--68.2%
mul-1-neg68.2%
unsub-neg68.2%
associate-/l*82.2%
Simplified82.2%
Taylor expanded in y around inf 62.4%
associate-/l*72.9%
Simplified72.9%
if -9.1999999999999998e83 < z < -2.1e21Initial program 81.2%
associate-*l/90.2%
Simplified90.2%
Taylor expanded in z around 0 71.7%
if -2.1e21 < z < -0.170000000000000012Initial program 70.0%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in x around 0 77.3%
associate-*r/97.3%
Simplified97.3%
if -0.170000000000000012 < z < -6.10000000000000004e-6Initial program 100.0%
associate-*l/100.0%
Simplified100.0%
Taylor expanded in z around 0 100.0%
Taylor expanded in t around inf 100.0%
associate-/l*99.2%
Simplified99.2%
div-inv99.2%
clear-num100.0%
Applied egg-rr100.0%
if -2.54999999999999993e-76 < z < 8.1999999999999995e-18Initial program 92.0%
associate-*l/94.9%
Simplified94.9%
Taylor expanded in z around 0 74.5%
associate-/l*77.8%
Simplified77.8%
Final simplification75.8%
(FPCore (x y z t a)
:precision binary64
(if (<= z -8.8e+83)
(- t (/ (- t x) (/ z y)))
(if (<= z -6.4e+18)
(+ x (* (- t x) (/ y a)))
(if (<= z -0.39)
(* t (/ (- y z) (- a z)))
(if (<= z -1.95e-6)
(+ x (* t (/ y a)))
(if (or (<= z -2.6e-73) (not (<= z 5e-17)))
(- t (/ y (/ z (- t x))))
(+ x (/ y (/ a (- t x))))))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -8.8e+83) {
tmp = t - ((t - x) / (z / y));
} else if (z <= -6.4e+18) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.39) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -1.95e-6) {
tmp = x + (t * (y / a));
} else if ((z <= -2.6e-73) || !(z <= 5e-17)) {
tmp = t - (y / (z / (t - x)));
} else {
tmp = x + (y / (a / (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 <= (-8.8d+83)) then
tmp = t - ((t - x) / (z / y))
else if (z <= (-6.4d+18)) then
tmp = x + ((t - x) * (y / a))
else if (z <= (-0.39d0)) then
tmp = t * ((y - z) / (a - z))
else if (z <= (-1.95d-6)) then
tmp = x + (t * (y / a))
else if ((z <= (-2.6d-73)) .or. (.not. (z <= 5d-17))) then
tmp = t - (y / (z / (t - x)))
else
tmp = x + (y / (a / (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 <= -8.8e+83) {
tmp = t - ((t - x) / (z / y));
} else if (z <= -6.4e+18) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.39) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -1.95e-6) {
tmp = x + (t * (y / a));
} else if ((z <= -2.6e-73) || !(z <= 5e-17)) {
tmp = t - (y / (z / (t - x)));
} else {
tmp = x + (y / (a / (t - x)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -8.8e+83: tmp = t - ((t - x) / (z / y)) elif z <= -6.4e+18: tmp = x + ((t - x) * (y / a)) elif z <= -0.39: tmp = t * ((y - z) / (a - z)) elif z <= -1.95e-6: tmp = x + (t * (y / a)) elif (z <= -2.6e-73) or not (z <= 5e-17): tmp = t - (y / (z / (t - x))) else: tmp = x + (y / (a / (t - x))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -8.8e+83) tmp = Float64(t - Float64(Float64(t - x) / Float64(z / y))); elseif (z <= -6.4e+18) tmp = Float64(x + Float64(Float64(t - x) * Float64(y / a))); elseif (z <= -0.39) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); elseif (z <= -1.95e-6) tmp = Float64(x + Float64(t * Float64(y / a))); elseif ((z <= -2.6e-73) || !(z <= 5e-17)) tmp = Float64(t - Float64(y / Float64(z / Float64(t - x)))); else tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -8.8e+83) tmp = t - ((t - x) / (z / y)); elseif (z <= -6.4e+18) tmp = x + ((t - x) * (y / a)); elseif (z <= -0.39) tmp = t * ((y - z) / (a - z)); elseif (z <= -1.95e-6) tmp = x + (t * (y / a)); elseif ((z <= -2.6e-73) || ~((z <= 5e-17))) tmp = t - (y / (z / (t - x))); else tmp = x + (y / (a / (t - x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -8.8e+83], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -6.4e+18], N[(x + N[(N[(t - x), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.39], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -1.95e-6], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -2.6e-73], N[Not[LessEqual[z, 5e-17]], $MachinePrecision]], N[(t - N[(y / N[(z / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -8.8 \cdot 10^{+83}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y}}\\
\mathbf{elif}\;z \leq -6.4 \cdot 10^{+18}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -0.39:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{elif}\;z \leq -1.95 \cdot 10^{-6}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -2.6 \cdot 10^{-73} \lor \neg \left(z \leq 5 \cdot 10^{-17}\right):\\
\;\;\;\;t - \frac{y}{\frac{z}{t - x}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\end{array}
\end{array}
if z < -8.79999999999999995e83Initial program 33.4%
associate-*l/70.1%
Simplified70.1%
Taylor expanded in z around inf 63.6%
associate--l+63.6%
associate-*r/63.6%
associate-*r/63.6%
div-sub63.6%
distribute-lft-out--63.6%
associate-*r/63.6%
distribute-rgt-out--63.7%
mul-1-neg63.7%
unsub-neg63.7%
associate-/l*82.6%
Simplified82.6%
Taylor expanded in y around inf 76.9%
if -8.79999999999999995e83 < z < -6.4e18Initial program 81.2%
associate-*l/90.2%
Simplified90.2%
Taylor expanded in z around 0 71.7%
if -6.4e18 < z < -0.39000000000000001Initial program 70.0%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in x around 0 77.3%
associate-*r/97.3%
Simplified97.3%
if -0.39000000000000001 < z < -1.95e-6Initial program 100.0%
associate-*l/100.0%
Simplified100.0%
Taylor expanded in z around 0 100.0%
Taylor expanded in t around inf 100.0%
associate-/l*99.2%
Simplified99.2%
div-inv99.2%
clear-num100.0%
Applied egg-rr100.0%
if -1.95e-6 < z < -2.6000000000000001e-73 or 4.9999999999999999e-17 < z Initial program 44.8%
associate-*l/62.5%
Simplified62.5%
Taylor expanded in z around inf 71.2%
associate--l+71.2%
associate-*r/71.2%
associate-*r/71.2%
div-sub71.2%
distribute-lft-out--71.2%
associate-*r/71.2%
distribute-rgt-out--71.3%
mul-1-neg71.3%
unsub-neg71.3%
associate-/l*81.9%
Simplified81.9%
Taylor expanded in y around inf 61.9%
associate-/l*71.2%
Simplified71.2%
if -2.6000000000000001e-73 < z < 4.9999999999999999e-17Initial program 92.0%
associate-*l/94.9%
Simplified94.9%
Taylor expanded in z around 0 74.5%
associate-/l*77.8%
Simplified77.8%
Final simplification76.2%
(FPCore (x y z t a)
:precision binary64
(if (<= z -9.2e+83)
(- t (/ (- t x) (/ z y)))
(if (<= z -2.1e+21)
(+ x (* (- t x) (/ y a)))
(if (<= z -0.031)
(/ t (/ (- a z) (- y z)))
(if (<= z -1.55e-12)
(+ x (* t (/ y a)))
(if (or (<= z -8e-73) (not (<= z 3.4e-18)))
(- t (/ y (/ z (- t x))))
(+ x (/ y (/ a (- t x))))))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.2e+83) {
tmp = t - ((t - x) / (z / y));
} else if (z <= -2.1e+21) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.031) {
tmp = t / ((a - z) / (y - z));
} else if (z <= -1.55e-12) {
tmp = x + (t * (y / a));
} else if ((z <= -8e-73) || !(z <= 3.4e-18)) {
tmp = t - (y / (z / (t - x)));
} else {
tmp = x + (y / (a / (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 <= (-9.2d+83)) then
tmp = t - ((t - x) / (z / y))
else if (z <= (-2.1d+21)) then
tmp = x + ((t - x) * (y / a))
else if (z <= (-0.031d0)) then
tmp = t / ((a - z) / (y - z))
else if (z <= (-1.55d-12)) then
tmp = x + (t * (y / a))
else if ((z <= (-8d-73)) .or. (.not. (z <= 3.4d-18))) then
tmp = t - (y / (z / (t - x)))
else
tmp = x + (y / (a / (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 <= -9.2e+83) {
tmp = t - ((t - x) / (z / y));
} else if (z <= -2.1e+21) {
tmp = x + ((t - x) * (y / a));
} else if (z <= -0.031) {
tmp = t / ((a - z) / (y - z));
} else if (z <= -1.55e-12) {
tmp = x + (t * (y / a));
} else if ((z <= -8e-73) || !(z <= 3.4e-18)) {
tmp = t - (y / (z / (t - x)));
} else {
tmp = x + (y / (a / (t - x)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -9.2e+83: tmp = t - ((t - x) / (z / y)) elif z <= -2.1e+21: tmp = x + ((t - x) * (y / a)) elif z <= -0.031: tmp = t / ((a - z) / (y - z)) elif z <= -1.55e-12: tmp = x + (t * (y / a)) elif (z <= -8e-73) or not (z <= 3.4e-18): tmp = t - (y / (z / (t - x))) else: tmp = x + (y / (a / (t - x))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -9.2e+83) tmp = Float64(t - Float64(Float64(t - x) / Float64(z / y))); elseif (z <= -2.1e+21) tmp = Float64(x + Float64(Float64(t - x) * Float64(y / a))); elseif (z <= -0.031) tmp = Float64(t / Float64(Float64(a - z) / Float64(y - z))); elseif (z <= -1.55e-12) tmp = Float64(x + Float64(t * Float64(y / a))); elseif ((z <= -8e-73) || !(z <= 3.4e-18)) tmp = Float64(t - Float64(y / Float64(z / Float64(t - x)))); else tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -9.2e+83) tmp = t - ((t - x) / (z / y)); elseif (z <= -2.1e+21) tmp = x + ((t - x) * (y / a)); elseif (z <= -0.031) tmp = t / ((a - z) / (y - z)); elseif (z <= -1.55e-12) tmp = x + (t * (y / a)); elseif ((z <= -8e-73) || ~((z <= 3.4e-18))) tmp = t - (y / (z / (t - x))); else tmp = x + (y / (a / (t - x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -9.2e+83], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -2.1e+21], N[(x + N[(N[(t - x), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.031], N[(t / N[(N[(a - z), $MachinePrecision] / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -1.55e-12], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -8e-73], N[Not[LessEqual[z, 3.4e-18]], $MachinePrecision]], N[(t - N[(y / N[(z / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.2 \cdot 10^{+83}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y}}\\
\mathbf{elif}\;z \leq -2.1 \cdot 10^{+21}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -0.031:\\
\;\;\;\;\frac{t}{\frac{a - z}{y - z}}\\
\mathbf{elif}\;z \leq -1.55 \cdot 10^{-12}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{elif}\;z \leq -8 \cdot 10^{-73} \lor \neg \left(z \leq 3.4 \cdot 10^{-18}\right):\\
\;\;\;\;t - \frac{y}{\frac{z}{t - x}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\end{array}
\end{array}
if z < -9.1999999999999998e83Initial program 33.4%
associate-*l/70.1%
Simplified70.1%
Taylor expanded in z around inf 63.6%
associate--l+63.6%
associate-*r/63.6%
associate-*r/63.6%
div-sub63.6%
distribute-lft-out--63.6%
associate-*r/63.6%
distribute-rgt-out--63.7%
mul-1-neg63.7%
unsub-neg63.7%
associate-/l*82.6%
Simplified82.6%
Taylor expanded in y around inf 76.9%
if -9.1999999999999998e83 < z < -2.1e21Initial program 81.2%
associate-*l/90.2%
Simplified90.2%
Taylor expanded in z around 0 71.7%
if -2.1e21 < z < -0.031Initial program 70.0%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in x around 0 77.3%
associate-/l*97.5%
Simplified97.5%
if -0.031 < z < -1.5500000000000001e-12Initial program 100.0%
associate-*l/100.0%
Simplified100.0%
Taylor expanded in z around 0 100.0%
Taylor expanded in t around inf 100.0%
associate-/l*99.2%
Simplified99.2%
div-inv99.2%
clear-num100.0%
Applied egg-rr100.0%
if -1.5500000000000001e-12 < z < -7.99999999999999998e-73 or 3.40000000000000001e-18 < z Initial program 44.8%
associate-*l/62.5%
Simplified62.5%
Taylor expanded in z around inf 71.2%
associate--l+71.2%
associate-*r/71.2%
associate-*r/71.2%
div-sub71.2%
distribute-lft-out--71.2%
associate-*r/71.2%
distribute-rgt-out--71.3%
mul-1-neg71.3%
unsub-neg71.3%
associate-/l*81.9%
Simplified81.9%
Taylor expanded in y around inf 61.9%
associate-/l*71.2%
Simplified71.2%
if -7.99999999999999998e-73 < z < 3.40000000000000001e-18Initial program 92.0%
associate-*l/94.9%
Simplified94.9%
Taylor expanded in z around 0 74.5%
associate-/l*77.8%
Simplified77.8%
Final simplification76.2%
(FPCore (x y z t a)
:precision binary64
(if (or (<= z -8.2e+167)
(not (or (<= z -9.8e-11) (and (not (<= z -8e-73)) (<= z 6.5e+98)))))
(- t (/ (- t x) (/ z (- y a))))
(+ x (* (- t x) (/ (- y z) (- a z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -8.2e+167) || !((z <= -9.8e-11) || (!(z <= -8e-73) && (z <= 6.5e+98)))) {
tmp = t - ((t - x) / (z / (y - a)));
} else {
tmp = x + ((t - x) * ((y - z) / (a - z)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((z <= (-8.2d+167)) .or. (.not. (z <= (-9.8d-11)) .or. (.not. (z <= (-8d-73))) .and. (z <= 6.5d+98))) then
tmp = t - ((t - x) / (z / (y - a)))
else
tmp = x + ((t - x) * ((y - z) / (a - z)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -8.2e+167) || !((z <= -9.8e-11) || (!(z <= -8e-73) && (z <= 6.5e+98)))) {
tmp = t - ((t - x) / (z / (y - a)));
} else {
tmp = x + ((t - x) * ((y - z) / (a - z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -8.2e+167) or not ((z <= -9.8e-11) or (not (z <= -8e-73) and (z <= 6.5e+98))): tmp = t - ((t - x) / (z / (y - a))) else: tmp = x + ((t - x) * ((y - z) / (a - z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -8.2e+167) || !((z <= -9.8e-11) || (!(z <= -8e-73) && (z <= 6.5e+98)))) tmp = Float64(t - Float64(Float64(t - x) / Float64(z / Float64(y - a)))); else tmp = Float64(x + Float64(Float64(t - x) * Float64(Float64(y - z) / Float64(a - z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -8.2e+167) || ~(((z <= -9.8e-11) || (~((z <= -8e-73)) && (z <= 6.5e+98))))) tmp = t - ((t - x) / (z / (y - a))); else tmp = x + ((t - x) * ((y - z) / (a - z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -8.2e+167], N[Not[Or[LessEqual[z, -9.8e-11], And[N[Not[LessEqual[z, -8e-73]], $MachinePrecision], LessEqual[z, 6.5e+98]]]], $MachinePrecision]], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / N[(y - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -8.2 \cdot 10^{+167} \lor \neg \left(z \leq -9.8 \cdot 10^{-11} \lor \neg \left(z \leq -8 \cdot 10^{-73}\right) \land z \leq 6.5 \cdot 10^{+98}\right):\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y - a}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y - z}{a - z}\\
\end{array}
\end{array}
if z < -8.2e167 or -9.7999999999999998e-11 < z < -7.99999999999999998e-73 or 6.4999999999999999e98 < z Initial program 29.7%
associate-*l/57.5%
Simplified57.5%
Taylor expanded in z around inf 69.4%
associate--l+69.4%
associate-*r/69.4%
associate-*r/69.4%
div-sub69.4%
distribute-lft-out--69.4%
associate-*r/69.4%
distribute-rgt-out--69.6%
mul-1-neg69.6%
unsub-neg69.6%
associate-/l*89.0%
Simplified89.0%
if -8.2e167 < z < -9.7999999999999998e-11 or -7.99999999999999998e-73 < z < 6.4999999999999999e98Initial program 80.9%
associate-*l/90.3%
Simplified90.3%
Final simplification89.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ (- y z) a))))
(if (<= z -9e+75)
t
(if (<= z -2.6e-19)
t_1
(if (<= z -1e-38)
(/ x (/ z (- a)))
(if (<= z -6.8e-193) x (if (<= z 255.0) t_1 t)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * ((y - z) / a);
double tmp;
if (z <= -9e+75) {
tmp = t;
} else if (z <= -2.6e-19) {
tmp = t_1;
} else if (z <= -1e-38) {
tmp = x / (z / -a);
} else if (z <= -6.8e-193) {
tmp = x;
} else if (z <= 255.0) {
tmp = t_1;
} 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) :: t_1
real(8) :: tmp
t_1 = t * ((y - z) / a)
if (z <= (-9d+75)) then
tmp = t
else if (z <= (-2.6d-19)) then
tmp = t_1
else if (z <= (-1d-38)) then
tmp = x / (z / -a)
else if (z <= (-6.8d-193)) then
tmp = x
else if (z <= 255.0d0) then
tmp = t_1
else
tmp = t
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);
double tmp;
if (z <= -9e+75) {
tmp = t;
} else if (z <= -2.6e-19) {
tmp = t_1;
} else if (z <= -1e-38) {
tmp = x / (z / -a);
} else if (z <= -6.8e-193) {
tmp = x;
} else if (z <= 255.0) {
tmp = t_1;
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * ((y - z) / a) tmp = 0 if z <= -9e+75: tmp = t elif z <= -2.6e-19: tmp = t_1 elif z <= -1e-38: tmp = x / (z / -a) elif z <= -6.8e-193: tmp = x elif z <= 255.0: tmp = t_1 else: tmp = t return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(Float64(y - z) / a)) tmp = 0.0 if (z <= -9e+75) tmp = t; elseif (z <= -2.6e-19) tmp = t_1; elseif (z <= -1e-38) tmp = Float64(x / Float64(z / Float64(-a))); elseif (z <= -6.8e-193) tmp = x; elseif (z <= 255.0) tmp = t_1; else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * ((y - z) / a); tmp = 0.0; if (z <= -9e+75) tmp = t; elseif (z <= -2.6e-19) tmp = t_1; elseif (z <= -1e-38) tmp = x / (z / -a); elseif (z <= -6.8e-193) tmp = x; elseif (z <= 255.0) tmp = t_1; else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(N[(y - z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -9e+75], t, If[LessEqual[z, -2.6e-19], t$95$1, If[LessEqual[z, -1e-38], N[(x / N[(z / (-a)), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -6.8e-193], x, If[LessEqual[z, 255.0], t$95$1, t]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y - z}{a}\\
\mathbf{if}\;z \leq -9 \cdot 10^{+75}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq -2.6 \cdot 10^{-19}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -1 \cdot 10^{-38}:\\
\;\;\;\;\frac{x}{\frac{z}{-a}}\\
\mathbf{elif}\;z \leq -6.8 \cdot 10^{-193}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 255:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -9.0000000000000007e75 or 255 < z Initial program 39.1%
associate-*l/67.8%
Simplified67.8%
Taylor expanded in z around inf 51.1%
if -9.0000000000000007e75 < z < -2.60000000000000013e-19 or -6.8000000000000004e-193 < z < 255Initial program 87.3%
associate-*l/92.6%
Simplified92.6%
Taylor expanded in x around 0 52.1%
associate-*r/54.3%
Simplified54.3%
Taylor expanded in a around inf 48.0%
if -2.60000000000000013e-19 < z < -9.9999999999999996e-39Initial program 3.1%
associate-*l/3.1%
Simplified3.1%
Taylor expanded in y around 0 3.1%
mul-1-neg3.1%
associate-*r/3.1%
unsub-neg3.1%
Simplified3.1%
Taylor expanded in t around 0 3.1%
sub-neg3.1%
mul-1-neg3.1%
remove-double-neg3.1%
associate-/l*3.1%
div-sub3.1%
sub-neg3.1%
*-inverses3.1%
metadata-eval3.1%
Simplified3.1%
Taylor expanded in a around 0 3.1%
neg-mul-13.1%
associate-+r+76.2%
neg-mul-176.2%
mul-1-neg76.2%
unsub-neg76.2%
distribute-rgt1-in76.2%
metadata-eval76.2%
mul0-lft76.2%
neg-sub076.2%
distribute-frac-neg76.2%
distribute-lft-neg-in76.2%
*-commutative76.2%
associate-/l*76.2%
Simplified76.2%
if -9.9999999999999996e-39 < z < -6.8000000000000004e-193Initial program 88.2%
associate-*l/88.3%
Simplified88.3%
Taylor expanded in a around inf 49.7%
Final simplification50.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- t (* x (/ (- a y) z)))))
(if (<= z -3.8e+70)
t_1
(if (<= z -8e-17)
(* t (/ (- y z) (- a z)))
(if (<= z -8e-73)
t_1
(if (<= z 7e-19)
(+ x (/ y (/ a (- t x))))
(- t (/ y (/ z (- t x))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t - (x * ((a - y) / z));
double tmp;
if (z <= -3.8e+70) {
tmp = t_1;
} else if (z <= -8e-17) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -8e-73) {
tmp = t_1;
} else if (z <= 7e-19) {
tmp = x + (y / (a / (t - x)));
} else {
tmp = t - (y / (z / (t - x)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = t - (x * ((a - y) / z))
if (z <= (-3.8d+70)) then
tmp = t_1
else if (z <= (-8d-17)) then
tmp = t * ((y - z) / (a - z))
else if (z <= (-8d-73)) then
tmp = t_1
else if (z <= 7d-19) then
tmp = x + (y / (a / (t - x)))
else
tmp = t - (y / (z / (t - x)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = t - (x * ((a - y) / z));
double tmp;
if (z <= -3.8e+70) {
tmp = t_1;
} else if (z <= -8e-17) {
tmp = t * ((y - z) / (a - z));
} else if (z <= -8e-73) {
tmp = t_1;
} else if (z <= 7e-19) {
tmp = x + (y / (a / (t - x)));
} else {
tmp = t - (y / (z / (t - x)));
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t - (x * ((a - y) / z)) tmp = 0 if z <= -3.8e+70: tmp = t_1 elif z <= -8e-17: tmp = t * ((y - z) / (a - z)) elif z <= -8e-73: tmp = t_1 elif z <= 7e-19: tmp = x + (y / (a / (t - x))) else: tmp = t - (y / (z / (t - x))) return tmp
function code(x, y, z, t, a) t_1 = Float64(t - Float64(x * Float64(Float64(a - y) / z))) tmp = 0.0 if (z <= -3.8e+70) tmp = t_1; elseif (z <= -8e-17) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); elseif (z <= -8e-73) tmp = t_1; elseif (z <= 7e-19) tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); else tmp = Float64(t - Float64(y / Float64(z / Float64(t - x)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t - (x * ((a - y) / z)); tmp = 0.0; if (z <= -3.8e+70) tmp = t_1; elseif (z <= -8e-17) tmp = t * ((y - z) / (a - z)); elseif (z <= -8e-73) tmp = t_1; elseif (z <= 7e-19) tmp = x + (y / (a / (t - x))); else tmp = t - (y / (z / (t - x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t - N[(x * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.8e+70], t$95$1, If[LessEqual[z, -8e-17], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -8e-73], t$95$1, If[LessEqual[z, 7e-19], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t - N[(y / N[(z / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t - x \cdot \frac{a - y}{z}\\
\mathbf{if}\;z \leq -3.8 \cdot 10^{+70}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -8 \cdot 10^{-17}:\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{elif}\;z \leq -8 \cdot 10^{-73}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 7 \cdot 10^{-19}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\mathbf{else}:\\
\;\;\;\;t - \frac{y}{\frac{z}{t - x}}\\
\end{array}
\end{array}
if z < -3.7999999999999998e70 or -8.00000000000000057e-17 < z < -7.99999999999999998e-73Initial program 36.6%
associate-*l/66.8%
Simplified66.8%
Taylor expanded in z around inf 65.6%
associate--l+65.6%
associate-*r/65.6%
associate-*r/65.6%
div-sub65.6%
distribute-lft-out--65.6%
associate-*r/65.6%
distribute-rgt-out--65.9%
mul-1-neg65.9%
unsub-neg65.9%
associate-/l*80.8%
Simplified80.8%
Taylor expanded in t around 0 67.5%
mul-1-neg67.5%
associate-*r/75.4%
*-commutative75.4%
distribute-rgt-neg-in75.4%
Simplified75.4%
if -3.7999999999999998e70 < z < -8.00000000000000057e-17Initial program 79.5%
associate-*l/89.6%
Simplified89.6%
Taylor expanded in x around 0 66.9%
associate-*r/71.8%
Simplified71.8%
if -7.99999999999999998e-73 < z < 7.00000000000000031e-19Initial program 92.0%
associate-*l/94.9%
Simplified94.9%
Taylor expanded in z around 0 74.5%
associate-/l*77.8%
Simplified77.8%
if 7.00000000000000031e-19 < z Initial program 45.1%
associate-*l/65.9%
Simplified65.9%
Taylor expanded in z around inf 67.7%
associate--l+67.7%
associate-*r/67.7%
associate-*r/67.7%
div-sub67.7%
distribute-lft-out--67.7%
associate-*r/67.7%
distribute-rgt-out--67.8%
mul-1-neg67.8%
unsub-neg67.8%
associate-/l*80.3%
Simplified80.3%
Taylor expanded in y around inf 61.0%
associate-/l*71.9%
Simplified71.9%
Final simplification75.2%
(FPCore (x y z t a)
:precision binary64
(if (<= z -9.5e+105)
(/ (- t) (+ (/ a z) -1.0))
(if (<= z -3.9e+20)
(- x (/ x (/ a y)))
(if (or (<= z -3.8) (not (<= z 2.8e-10)))
(* t (- (/ (- y) z) -1.0))
(+ x (* t (/ y a)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.5e+105) {
tmp = -t / ((a / z) + -1.0);
} else if (z <= -3.9e+20) {
tmp = x - (x / (a / y));
} else if ((z <= -3.8) || !(z <= 2.8e-10)) {
tmp = t * ((-y / z) - -1.0);
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-9.5d+105)) then
tmp = -t / ((a / z) + (-1.0d0))
else if (z <= (-3.9d+20)) then
tmp = x - (x / (a / y))
else if ((z <= (-3.8d0)) .or. (.not. (z <= 2.8d-10))) then
tmp = t * ((-y / z) - (-1.0d0))
else
tmp = x + (t * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.5e+105) {
tmp = -t / ((a / z) + -1.0);
} else if (z <= -3.9e+20) {
tmp = x - (x / (a / y));
} else if ((z <= -3.8) || !(z <= 2.8e-10)) {
tmp = t * ((-y / z) - -1.0);
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -9.5e+105: tmp = -t / ((a / z) + -1.0) elif z <= -3.9e+20: tmp = x - (x / (a / y)) elif (z <= -3.8) or not (z <= 2.8e-10): tmp = t * ((-y / z) - -1.0) else: tmp = x + (t * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -9.5e+105) tmp = Float64(Float64(-t) / Float64(Float64(a / z) + -1.0)); elseif (z <= -3.9e+20) tmp = Float64(x - Float64(x / Float64(a / y))); elseif ((z <= -3.8) || !(z <= 2.8e-10)) tmp = Float64(t * Float64(Float64(Float64(-y) / z) - -1.0)); else tmp = Float64(x + Float64(t * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -9.5e+105) tmp = -t / ((a / z) + -1.0); elseif (z <= -3.9e+20) tmp = x - (x / (a / y)); elseif ((z <= -3.8) || ~((z <= 2.8e-10))) tmp = t * ((-y / z) - -1.0); else tmp = x + (t * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -9.5e+105], N[((-t) / N[(N[(a / z), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -3.9e+20], N[(x - N[(x / N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[z, -3.8], N[Not[LessEqual[z, 2.8e-10]], $MachinePrecision]], N[(t * N[(N[((-y) / z), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.5 \cdot 10^{+105}:\\
\;\;\;\;\frac{-t}{\frac{a}{z} + -1}\\
\mathbf{elif}\;z \leq -3.9 \cdot 10^{+20}:\\
\;\;\;\;x - \frac{x}{\frac{a}{y}}\\
\mathbf{elif}\;z \leq -3.8 \lor \neg \left(z \leq 2.8 \cdot 10^{-10}\right):\\
\;\;\;\;t \cdot \left(\frac{-y}{z} - -1\right)\\
\mathbf{else}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\end{array}
\end{array}
if z < -9.4999999999999995e105Initial program 33.9%
associate-*l/69.7%
Simplified69.7%
Taylor expanded in y around 0 28.3%
mul-1-neg28.3%
associate-*r/61.0%
unsub-neg61.0%
Simplified61.0%
Taylor expanded in x around 0 31.0%
mul-1-neg31.0%
associate-/l*64.8%
div-sub64.8%
sub-neg64.8%
*-inverses64.8%
metadata-eval64.8%
Simplified64.8%
if -9.4999999999999995e105 < z < -3.9e20Initial program 65.8%
associate-*l/86.0%
Simplified86.0%
Taylor expanded in z around 0 59.4%
Taylor expanded in t around 0 45.8%
mul-1-neg45.8%
unsub-neg45.8%
associate-/l*52.4%
Simplified52.4%
if -3.9e20 < z < -3.7999999999999998 or 2.80000000000000015e-10 < z Initial program 45.9%
associate-*l/67.4%
Simplified67.4%
Taylor expanded in x around 0 44.8%
associate-*r/66.2%
Simplified66.2%
Taylor expanded in a around 0 55.6%
mul-1-neg55.6%
div-sub55.6%
sub-neg55.6%
*-inverses55.6%
metadata-eval55.6%
Simplified55.6%
if -3.7999999999999998 < z < 2.80000000000000015e-10Initial program 87.3%
associate-*l/89.8%
Simplified89.8%
Taylor expanded in z around 0 70.7%
Taylor expanded in t around inf 60.7%
associate-/l*64.5%
Simplified64.5%
div-inv65.4%
clear-num65.4%
Applied egg-rr65.4%
Final simplification61.7%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* t (/ y a))))
(if (<= z -2.7e+106)
t
(if (<= z -2.45e-191)
x
(if (<= z 2.55e-176)
t_1
(if (<= z 1.22e-37) x (if (<= z 1.25e-8) t_1 t)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = t * (y / a);
double tmp;
if (z <= -2.7e+106) {
tmp = t;
} else if (z <= -2.45e-191) {
tmp = x;
} else if (z <= 2.55e-176) {
tmp = t_1;
} else if (z <= 1.22e-37) {
tmp = x;
} else if (z <= 1.25e-8) {
tmp = t_1;
} 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) :: t_1
real(8) :: tmp
t_1 = t * (y / a)
if (z <= (-2.7d+106)) then
tmp = t
else if (z <= (-2.45d-191)) then
tmp = x
else if (z <= 2.55d-176) then
tmp = t_1
else if (z <= 1.22d-37) then
tmp = x
else if (z <= 1.25d-8) then
tmp = t_1
else
tmp = t
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);
double tmp;
if (z <= -2.7e+106) {
tmp = t;
} else if (z <= -2.45e-191) {
tmp = x;
} else if (z <= 2.55e-176) {
tmp = t_1;
} else if (z <= 1.22e-37) {
tmp = x;
} else if (z <= 1.25e-8) {
tmp = t_1;
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = t * (y / a) tmp = 0 if z <= -2.7e+106: tmp = t elif z <= -2.45e-191: tmp = x elif z <= 2.55e-176: tmp = t_1 elif z <= 1.22e-37: tmp = x elif z <= 1.25e-8: tmp = t_1 else: tmp = t return tmp
function code(x, y, z, t, a) t_1 = Float64(t * Float64(y / a)) tmp = 0.0 if (z <= -2.7e+106) tmp = t; elseif (z <= -2.45e-191) tmp = x; elseif (z <= 2.55e-176) tmp = t_1; elseif (z <= 1.22e-37) tmp = x; elseif (z <= 1.25e-8) tmp = t_1; else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t * (y / a); tmp = 0.0; if (z <= -2.7e+106) tmp = t; elseif (z <= -2.45e-191) tmp = x; elseif (z <= 2.55e-176) tmp = t_1; elseif (z <= 1.22e-37) tmp = x; elseif (z <= 1.25e-8) tmp = t_1; else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2.7e+106], t, If[LessEqual[z, -2.45e-191], x, If[LessEqual[z, 2.55e-176], t$95$1, If[LessEqual[z, 1.22e-37], x, If[LessEqual[z, 1.25e-8], t$95$1, t]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{y}{a}\\
\mathbf{if}\;z \leq -2.7 \cdot 10^{+106}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq -2.45 \cdot 10^{-191}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 2.55 \cdot 10^{-176}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq 1.22 \cdot 10^{-37}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 1.25 \cdot 10^{-8}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -2.70000000000000006e106 or 1.2499999999999999e-8 < z Initial program 39.2%
associate-*l/67.0%
Simplified67.0%
Taylor expanded in z around inf 51.6%
if -2.70000000000000006e106 < z < -2.45e-191 or 2.5500000000000001e-176 < z < 1.21999999999999994e-37Initial program 76.5%
associate-*l/85.3%
Simplified85.3%
Taylor expanded in a around inf 34.9%
if -2.45e-191 < z < 2.5500000000000001e-176 or 1.21999999999999994e-37 < z < 1.2499999999999999e-8Initial program 96.1%
associate-*l/96.2%
Simplified96.2%
Taylor expanded in x around 0 55.4%
associate-*r/58.3%
Simplified58.3%
Taylor expanded in z around 0 51.6%
associate-*r/56.4%
Simplified56.4%
Final simplification47.0%
(FPCore (x y z t a)
:precision binary64
(if (<= z -1.1e+106)
t
(if (<= z -1.68e-192)
x
(if (<= z 1.75e-177)
(/ t (/ a y))
(if (<= z 3.6e-37) x (if (<= z 4e-12) (* t (/ y a)) t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.1e+106) {
tmp = t;
} else if (z <= -1.68e-192) {
tmp = x;
} else if (z <= 1.75e-177) {
tmp = t / (a / y);
} else if (z <= 3.6e-37) {
tmp = x;
} else if (z <= 4e-12) {
tmp = t * (y / a);
} 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 <= (-1.1d+106)) then
tmp = t
else if (z <= (-1.68d-192)) then
tmp = x
else if (z <= 1.75d-177) then
tmp = t / (a / y)
else if (z <= 3.6d-37) then
tmp = x
else if (z <= 4d-12) then
tmp = t * (y / a)
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 <= -1.1e+106) {
tmp = t;
} else if (z <= -1.68e-192) {
tmp = x;
} else if (z <= 1.75e-177) {
tmp = t / (a / y);
} else if (z <= 3.6e-37) {
tmp = x;
} else if (z <= 4e-12) {
tmp = t * (y / a);
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.1e+106: tmp = t elif z <= -1.68e-192: tmp = x elif z <= 1.75e-177: tmp = t / (a / y) elif z <= 3.6e-37: tmp = x elif z <= 4e-12: tmp = t * (y / a) else: tmp = t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.1e+106) tmp = t; elseif (z <= -1.68e-192) tmp = x; elseif (z <= 1.75e-177) tmp = Float64(t / Float64(a / y)); elseif (z <= 3.6e-37) tmp = x; elseif (z <= 4e-12) tmp = Float64(t * Float64(y / a)); else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.1e+106) tmp = t; elseif (z <= -1.68e-192) tmp = x; elseif (z <= 1.75e-177) tmp = t / (a / y); elseif (z <= 3.6e-37) tmp = x; elseif (z <= 4e-12) tmp = t * (y / a); else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.1e+106], t, If[LessEqual[z, -1.68e-192], x, If[LessEqual[z, 1.75e-177], N[(t / N[(a / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.6e-37], x, If[LessEqual[z, 4e-12], N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision], t]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.1 \cdot 10^{+106}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq -1.68 \cdot 10^{-192}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{-177}:\\
\;\;\;\;\frac{t}{\frac{a}{y}}\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-37}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 4 \cdot 10^{-12}:\\
\;\;\;\;t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -1.09999999999999996e106 or 3.99999999999999992e-12 < z Initial program 39.2%
associate-*l/67.0%
Simplified67.0%
Taylor expanded in z around inf 51.6%
if -1.09999999999999996e106 < z < -1.6799999999999999e-192 or 1.7500000000000001e-177 < z < 3.60000000000000007e-37Initial program 76.5%
associate-*l/85.3%
Simplified85.3%
Taylor expanded in a around inf 34.9%
if -1.6799999999999999e-192 < z < 1.7500000000000001e-177Initial program 97.7%
associate-*l/95.7%
Simplified95.7%
Taylor expanded in y around -inf 71.5%
Taylor expanded in t around inf 51.6%
associate-/l*53.0%
Simplified53.0%
Taylor expanded in a around inf 50.8%
if 3.60000000000000007e-37 < z < 3.99999999999999992e-12Initial program 84.4%
associate-*l/99.7%
Simplified99.7%
Taylor expanded in x around 0 67.9%
associate-*r/83.2%
Simplified83.2%
Taylor expanded in z around 0 68.6%
associate-*r/99.7%
Simplified99.7%
Final simplification47.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -5e-53) (not (<= a 3.3e+20))) (+ x (/ (- t x) (/ a (- y z)))) (- t (/ (- t x) (/ z y)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -5e-53) || !(a <= 3.3e+20)) {
tmp = x + ((t - x) / (a / (y - z)));
} else {
tmp = t - ((t - x) / (z / y));
}
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 <= (-5d-53)) .or. (.not. (a <= 3.3d+20))) then
tmp = x + ((t - x) / (a / (y - z)))
else
tmp = t - ((t - x) / (z / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -5e-53) || !(a <= 3.3e+20)) {
tmp = x + ((t - x) / (a / (y - z)));
} else {
tmp = t - ((t - x) / (z / y));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -5e-53) or not (a <= 3.3e+20): tmp = x + ((t - x) / (a / (y - z))) else: tmp = t - ((t - x) / (z / y)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -5e-53) || !(a <= 3.3e+20)) tmp = Float64(x + Float64(Float64(t - x) / Float64(a / Float64(y - z)))); else tmp = Float64(t - Float64(Float64(t - x) / Float64(z / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -5e-53) || ~((a <= 3.3e+20))) tmp = x + ((t - x) / (a / (y - z))); else tmp = t - ((t - x) / (z / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -5e-53], N[Not[LessEqual[a, 3.3e+20]], $MachinePrecision]], N[(x + N[(N[(t - x), $MachinePrecision] / N[(a / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5 \cdot 10^{-53} \lor \neg \left(a \leq 3.3 \cdot 10^{+20}\right):\\
\;\;\;\;x + \frac{t - x}{\frac{a}{y - z}}\\
\mathbf{else}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y}}\\
\end{array}
\end{array}
if a < -5e-53 or 3.3e20 < a Initial program 69.3%
associate-*l/89.4%
Simplified89.4%
Taylor expanded in a around inf 63.4%
associate-/l*73.4%
Simplified73.4%
if -5e-53 < a < 3.3e20Initial program 57.3%
associate-*l/68.6%
Simplified68.6%
Taylor expanded in z around inf 77.8%
associate--l+77.8%
associate-*r/77.8%
associate-*r/77.8%
div-sub78.6%
distribute-lft-out--78.6%
associate-*r/78.6%
distribute-rgt-out--78.6%
mul-1-neg78.6%
unsub-neg78.6%
associate-/l*84.7%
Simplified84.7%
Taylor expanded in y around inf 76.2%
Final simplification74.8%
(FPCore (x y z t a)
:precision binary64
(if (<= a -1.06e-51)
(+ x (/ (- y z) (/ (- a z) t)))
(if (<= a 1.9e+15)
(- t (/ (- t x) (/ z y)))
(+ x (/ (- t x) (/ a (- y z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.06e-51) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 1.9e+15) {
tmp = t - ((t - x) / (z / y));
} else {
tmp = x + ((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 <= (-1.06d-51)) then
tmp = x + ((y - z) / ((a - z) / t))
else if (a <= 1.9d+15) then
tmp = t - ((t - x) / (z / y))
else
tmp = x + ((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 <= -1.06e-51) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 1.9e+15) {
tmp = t - ((t - x) / (z / y));
} else {
tmp = x + ((t - x) / (a / (y - z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.06e-51: tmp = x + ((y - z) / ((a - z) / t)) elif a <= 1.9e+15: tmp = t - ((t - x) / (z / y)) else: tmp = x + ((t - x) / (a / (y - z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.06e-51) tmp = Float64(x + Float64(Float64(y - z) / Float64(Float64(a - z) / t))); elseif (a <= 1.9e+15) tmp = Float64(t - Float64(Float64(t - x) / Float64(z / y))); else tmp = Float64(x + Float64(Float64(t - x) / Float64(a / Float64(y - z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.06e-51) tmp = x + ((y - z) / ((a - z) / t)); elseif (a <= 1.9e+15) tmp = t - ((t - x) / (z / y)); else tmp = x + ((t - x) / (a / (y - z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.06e-51], N[(x + N[(N[(y - z), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.9e+15], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] / N[(a / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.06 \cdot 10^{-51}:\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t}}\\
\mathbf{elif}\;a \leq 1.9 \cdot 10^{+15}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t - x}{\frac{a}{y - z}}\\
\end{array}
\end{array}
if a < -1.0600000000000001e-51Initial program 73.2%
associate-/l*87.9%
Simplified87.9%
Taylor expanded in t around inf 82.0%
if -1.0600000000000001e-51 < a < 1.9e15Initial program 57.3%
associate-*l/68.6%
Simplified68.6%
Taylor expanded in z around inf 77.8%
associate--l+77.8%
associate-*r/77.8%
associate-*r/77.8%
div-sub78.6%
distribute-lft-out--78.6%
associate-*r/78.6%
distribute-rgt-out--78.6%
mul-1-neg78.6%
unsub-neg78.6%
associate-/l*84.7%
Simplified84.7%
Taylor expanded in y around inf 76.2%
if 1.9e15 < a Initial program 61.9%
associate-*l/89.9%
Simplified89.9%
Taylor expanded in a around inf 61.5%
associate-/l*76.6%
Simplified76.6%
Final simplification78.2%
(FPCore (x y z t a)
:precision binary64
(if (<= a -5.8e-63)
(+ x (/ (- y z) (/ (- a z) t)))
(if (<= a 3.5e+20)
(+ t (/ (* (- t x) (- a y)) z))
(+ x (/ (- t x) (/ a (- y z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -5.8e-63) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 3.5e+20) {
tmp = t + (((t - x) * (a - y)) / z);
} else {
tmp = x + ((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 <= (-5.8d-63)) then
tmp = x + ((y - z) / ((a - z) / t))
else if (a <= 3.5d+20) then
tmp = t + (((t - x) * (a - y)) / z)
else
tmp = x + ((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 <= -5.8e-63) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 3.5e+20) {
tmp = t + (((t - x) * (a - y)) / z);
} else {
tmp = x + ((t - x) / (a / (y - z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -5.8e-63: tmp = x + ((y - z) / ((a - z) / t)) elif a <= 3.5e+20: tmp = t + (((t - x) * (a - y)) / z) else: tmp = x + ((t - x) / (a / (y - z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -5.8e-63) tmp = Float64(x + Float64(Float64(y - z) / Float64(Float64(a - z) / t))); elseif (a <= 3.5e+20) tmp = Float64(t + Float64(Float64(Float64(t - x) * Float64(a - y)) / z)); else tmp = Float64(x + Float64(Float64(t - x) / Float64(a / Float64(y - z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -5.8e-63) tmp = x + ((y - z) / ((a - z) / t)); elseif (a <= 3.5e+20) tmp = t + (((t - x) * (a - y)) / z); else tmp = x + ((t - x) / (a / (y - z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -5.8e-63], N[(x + N[(N[(y - z), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 3.5e+20], N[(t + N[(N[(N[(t - x), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] / N[(a / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.8 \cdot 10^{-63}:\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t}}\\
\mathbf{elif}\;a \leq 3.5 \cdot 10^{+20}:\\
\;\;\;\;t + \frac{\left(t - x\right) \cdot \left(a - y\right)}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t - x}{\frac{a}{y - z}}\\
\end{array}
\end{array}
if a < -5.7999999999999995e-63Initial program 73.1%
associate-/l*88.3%
Simplified88.3%
Taylor expanded in t around inf 82.6%
if -5.7999999999999995e-63 < a < 3.5e20Initial program 57.0%
associate-*l/67.8%
Simplified67.8%
associate-*l/57.0%
clear-num56.8%
associate-/r*67.8%
Applied egg-rr67.8%
Taylor expanded in z around inf 78.1%
associate--l+78.1%
associate-*r/78.1%
associate-*r/78.1%
div-sub78.9%
distribute-lft-out--78.9%
associate-*r/78.9%
mul-1-neg78.9%
unsub-neg78.9%
distribute-rgt-out--78.9%
Simplified78.9%
if 3.5e20 < a Initial program 61.9%
associate-*l/89.9%
Simplified89.9%
Taylor expanded in a around inf 61.5%
associate-/l*76.6%
Simplified76.6%
Final simplification79.8%
(FPCore (x y z t a)
:precision binary64
(if (<= a -2.15e-50)
(+ x (/ (- y z) (/ (- a z) t)))
(if (<= a 4.4e+18)
(- t (/ (- t x) (/ z (- y a))))
(+ x (/ (- t x) (/ a (- y z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.15e-50) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 4.4e+18) {
tmp = t - ((t - x) / (z / (y - a)));
} else {
tmp = x + ((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 <= (-2.15d-50)) then
tmp = x + ((y - z) / ((a - z) / t))
else if (a <= 4.4d+18) then
tmp = t - ((t - x) / (z / (y - a)))
else
tmp = x + ((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 <= -2.15e-50) {
tmp = x + ((y - z) / ((a - z) / t));
} else if (a <= 4.4e+18) {
tmp = t - ((t - x) / (z / (y - a)));
} else {
tmp = x + ((t - x) / (a / (y - z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.15e-50: tmp = x + ((y - z) / ((a - z) / t)) elif a <= 4.4e+18: tmp = t - ((t - x) / (z / (y - a))) else: tmp = x + ((t - x) / (a / (y - z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.15e-50) tmp = Float64(x + Float64(Float64(y - z) / Float64(Float64(a - z) / t))); elseif (a <= 4.4e+18) tmp = Float64(t - Float64(Float64(t - x) / Float64(z / Float64(y - a)))); else tmp = Float64(x + Float64(Float64(t - x) / Float64(a / Float64(y - z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -2.15e-50) tmp = x + ((y - z) / ((a - z) / t)); elseif (a <= 4.4e+18) tmp = t - ((t - x) / (z / (y - a))); else tmp = x + ((t - x) / (a / (y - z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.15e-50], N[(x + N[(N[(y - z), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 4.4e+18], N[(t - N[(N[(t - x), $MachinePrecision] / N[(z / N[(y - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] / N[(a / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.15 \cdot 10^{-50}:\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t}}\\
\mathbf{elif}\;a \leq 4.4 \cdot 10^{+18}:\\
\;\;\;\;t - \frac{t - x}{\frac{z}{y - a}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t - x}{\frac{a}{y - z}}\\
\end{array}
\end{array}
if a < -2.14999999999999999e-50Initial program 73.2%
associate-/l*87.9%
Simplified87.9%
Taylor expanded in t around inf 82.0%
if -2.14999999999999999e-50 < a < 4.4e18Initial program 57.3%
associate-*l/68.6%
Simplified68.6%
Taylor expanded in z around inf 77.8%
associate--l+77.8%
associate-*r/77.8%
associate-*r/77.8%
div-sub78.6%
distribute-lft-out--78.6%
associate-*r/78.6%
distribute-rgt-out--78.6%
mul-1-neg78.6%
unsub-neg78.6%
associate-/l*84.7%
Simplified84.7%
if 4.4e18 < a Initial program 61.9%
associate-*l/89.9%
Simplified89.9%
Taylor expanded in a around inf 61.5%
associate-/l*76.6%
Simplified76.6%
Final simplification82.4%
(FPCore (x y z t a) :precision binary64 (if (or (<= t -1.65e-71) (not (<= t 3.8e-64))) (* t (/ (- y z) (- a z))) (- x (/ x (/ a y)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((t <= -1.65e-71) || !(t <= 3.8e-64)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x - (x / (a / y));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((t <= (-1.65d-71)) .or. (.not. (t <= 3.8d-64))) then
tmp = t * ((y - z) / (a - z))
else
tmp = x - (x / (a / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((t <= -1.65e-71) || !(t <= 3.8e-64)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x - (x / (a / y));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (t <= -1.65e-71) or not (t <= 3.8e-64): tmp = t * ((y - z) / (a - z)) else: tmp = x - (x / (a / y)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((t <= -1.65e-71) || !(t <= 3.8e-64)) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(x - Float64(x / Float64(a / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((t <= -1.65e-71) || ~((t <= 3.8e-64))) tmp = t * ((y - z) / (a - z)); else tmp = x - (x / (a / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[t, -1.65e-71], N[Not[LessEqual[t, 3.8e-64]], $MachinePrecision]], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(x / N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.65 \cdot 10^{-71} \lor \neg \left(t \leq 3.8 \cdot 10^{-64}\right):\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{x}{\frac{a}{y}}\\
\end{array}
\end{array}
if t < -1.6500000000000001e-71 or 3.8000000000000002e-64 < t Initial program 64.2%
associate-*l/83.3%
Simplified83.3%
Taylor expanded in x around 0 54.1%
associate-*r/73.8%
Simplified73.8%
if -1.6500000000000001e-71 < t < 3.8000000000000002e-64Initial program 61.3%
associate-*l/69.9%
Simplified69.9%
Taylor expanded in z around 0 51.1%
Taylor expanded in t around 0 43.3%
mul-1-neg43.3%
unsub-neg43.3%
associate-/l*48.0%
Simplified48.0%
Final simplification65.6%
(FPCore (x y z t a) :precision binary64 (if (or (<= t -4.4e-70) (not (<= t 2.1e-37))) (* t (/ (- y z) (- a z))) (+ x (* (- t x) (/ y a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((t <= -4.4e-70) || !(t <= 2.1e-37)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x + ((t - x) * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((t <= (-4.4d-70)) .or. (.not. (t <= 2.1d-37))) then
tmp = t * ((y - z) / (a - z))
else
tmp = x + ((t - x) * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((t <= -4.4e-70) || !(t <= 2.1e-37)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x + ((t - x) * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (t <= -4.4e-70) or not (t <= 2.1e-37): tmp = t * ((y - z) / (a - z)) else: tmp = x + ((t - x) * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((t <= -4.4e-70) || !(t <= 2.1e-37)) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(x + Float64(Float64(t - x) * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((t <= -4.4e-70) || ~((t <= 2.1e-37))) tmp = t * ((y - z) / (a - z)); else tmp = x + ((t - x) * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[t, -4.4e-70], N[Not[LessEqual[t, 2.1e-37]], $MachinePrecision]], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.4 \cdot 10^{-70} \lor \neg \left(t \leq 2.1 \cdot 10^{-37}\right):\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;x + \left(t - x\right) \cdot \frac{y}{a}\\
\end{array}
\end{array}
if t < -4.3999999999999998e-70 or 2.1000000000000001e-37 < t Initial program 65.1%
associate-*l/84.3%
Simplified84.3%
Taylor expanded in x around 0 55.1%
associate-*r/75.6%
Simplified75.6%
if -4.3999999999999998e-70 < t < 2.1000000000000001e-37Initial program 59.9%
associate-*l/69.0%
Simplified69.0%
Taylor expanded in z around 0 50.6%
Final simplification67.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= t -8.8e-72) (not (<= t 1.8e-37))) (* t (/ (- y z) (- a z))) (+ x (/ y (/ a (- t x))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((t <= -8.8e-72) || !(t <= 1.8e-37)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x + (y / (a / (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 ((t <= (-8.8d-72)) .or. (.not. (t <= 1.8d-37))) then
tmp = t * ((y - z) / (a - z))
else
tmp = x + (y / (a / (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 ((t <= -8.8e-72) || !(t <= 1.8e-37)) {
tmp = t * ((y - z) / (a - z));
} else {
tmp = x + (y / (a / (t - x)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (t <= -8.8e-72) or not (t <= 1.8e-37): tmp = t * ((y - z) / (a - z)) else: tmp = x + (y / (a / (t - x))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((t <= -8.8e-72) || !(t <= 1.8e-37)) tmp = Float64(t * Float64(Float64(y - z) / Float64(a - z))); else tmp = Float64(x + Float64(y / Float64(a / Float64(t - x)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((t <= -8.8e-72) || ~((t <= 1.8e-37))) tmp = t * ((y - z) / (a - z)); else tmp = x + (y / (a / (t - x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[t, -8.8e-72], N[Not[LessEqual[t, 1.8e-37]], $MachinePrecision]], N[(t * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(a / N[(t - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8.8 \cdot 10^{-72} \lor \neg \left(t \leq 1.8 \cdot 10^{-37}\right):\\
\;\;\;\;t \cdot \frac{y - z}{a - z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - x}}\\
\end{array}
\end{array}
if t < -8.8000000000000001e-72 or 1.80000000000000004e-37 < t Initial program 65.1%
associate-*l/84.3%
Simplified84.3%
Taylor expanded in x around 0 55.1%
associate-*r/75.6%
Simplified75.6%
if -8.8000000000000001e-72 < t < 1.80000000000000004e-37Initial program 59.9%
associate-*l/69.0%
Simplified69.0%
Taylor expanded in z around 0 47.3%
associate-/l*50.7%
Simplified50.7%
Final simplification67.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -9.5e+105) (not (<= z 3e-6))) (/ (- t) (+ (/ a z) -1.0)) (+ x (* t (/ y a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -9.5e+105) || !(z <= 3e-6)) {
tmp = -t / ((a / z) + -1.0);
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((z <= (-9.5d+105)) .or. (.not. (z <= 3d-6))) then
tmp = -t / ((a / z) + (-1.0d0))
else
tmp = x + (t * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -9.5e+105) || !(z <= 3e-6)) {
tmp = -t / ((a / z) + -1.0);
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -9.5e+105) or not (z <= 3e-6): tmp = -t / ((a / z) + -1.0) else: tmp = x + (t * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -9.5e+105) || !(z <= 3e-6)) tmp = Float64(Float64(-t) / Float64(Float64(a / z) + -1.0)); else tmp = Float64(x + Float64(t * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -9.5e+105) || ~((z <= 3e-6))) tmp = -t / ((a / z) + -1.0); else tmp = x + (t * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -9.5e+105], N[Not[LessEqual[z, 3e-6]], $MachinePrecision]], N[((-t) / N[(N[(a / z), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.5 \cdot 10^{+105} \lor \neg \left(z \leq 3 \cdot 10^{-6}\right):\\
\;\;\;\;\frac{-t}{\frac{a}{z} + -1}\\
\mathbf{else}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\end{array}
\end{array}
if z < -9.4999999999999995e105 or 3.0000000000000001e-6 < z Initial program 39.2%
associate-*l/67.0%
Simplified67.0%
Taylor expanded in y around 0 29.7%
mul-1-neg29.7%
associate-*r/50.6%
unsub-neg50.6%
Simplified50.6%
Taylor expanded in x around 0 32.3%
mul-1-neg32.3%
associate-/l*56.0%
div-sub56.0%
sub-neg56.0%
*-inverses56.0%
metadata-eval56.0%
Simplified56.0%
if -9.4999999999999995e105 < z < 3.0000000000000001e-6Initial program 83.9%
associate-*l/89.4%
Simplified89.4%
Taylor expanded in z around 0 67.1%
Taylor expanded in t around inf 57.4%
associate-/l*60.6%
Simplified60.6%
div-inv61.3%
clear-num61.3%
Applied egg-rr61.3%
Final simplification58.9%
(FPCore (x y z t a) :precision binary64 (if (<= z -9.2e+106) t (if (<= z 2e-6) (+ x (* t (/ y a))) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.2e+106) {
tmp = t;
} else if (z <= 2e-6) {
tmp = x + (t * (y / a));
} 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 <= (-9.2d+106)) then
tmp = t
else if (z <= 2d-6) then
tmp = x + (t * (y / a))
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 <= -9.2e+106) {
tmp = t;
} else if (z <= 2e-6) {
tmp = x + (t * (y / a));
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -9.2e+106: tmp = t elif z <= 2e-6: tmp = x + (t * (y / a)) else: tmp = t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -9.2e+106) tmp = t; elseif (z <= 2e-6) tmp = Float64(x + Float64(t * Float64(y / a))); else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -9.2e+106) tmp = t; elseif (z <= 2e-6) tmp = x + (t * (y / a)); else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -9.2e+106], t, If[LessEqual[z, 2e-6], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.2 \cdot 10^{+106}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq 2 \cdot 10^{-6}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -9.2000000000000008e106 or 1.99999999999999991e-6 < z Initial program 39.2%
associate-*l/67.0%
Simplified67.0%
Taylor expanded in z around inf 51.6%
if -9.2000000000000008e106 < z < 1.99999999999999991e-6Initial program 83.9%
associate-*l/89.4%
Simplified89.4%
Taylor expanded in z around 0 67.1%
Taylor expanded in t around inf 57.4%
associate-/l*60.6%
Simplified60.6%
div-inv61.3%
clear-num61.3%
Applied egg-rr61.3%
Final simplification56.8%
(FPCore (x y z t a) :precision binary64 (if (<= z -5.3e+107) t (if (<= z 3.3e-27) x t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -5.3e+107) {
tmp = t;
} else if (z <= 3.3e-27) {
tmp = 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 <= (-5.3d+107)) then
tmp = t
else if (z <= 3.3d-27) then
tmp = 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 <= -5.3e+107) {
tmp = t;
} else if (z <= 3.3e-27) {
tmp = x;
} else {
tmp = t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -5.3e+107: tmp = t elif z <= 3.3e-27: tmp = x else: tmp = t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -5.3e+107) tmp = t; elseif (z <= 3.3e-27) tmp = x; else tmp = t; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -5.3e+107) tmp = t; elseif (z <= 3.3e-27) tmp = x; else tmp = t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -5.3e+107], t, If[LessEqual[z, 3.3e-27], x, t]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.3 \cdot 10^{+107}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{-27}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t\\
\end{array}
\end{array}
if z < -5.3e107 or 3.29999999999999998e-27 < z Initial program 41.2%
associate-*l/68.0%
Simplified68.0%
Taylor expanded in z around inf 50.0%
if -5.3e107 < z < 3.29999999999999998e-27Initial program 83.4%
associate-*l/89.1%
Simplified89.1%
Taylor expanded in a around inf 31.1%
Final simplification40.1%
(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 63.3%
associate-*l/79.1%
Simplified79.1%
Taylor expanded in y around 0 35.4%
mul-1-neg35.4%
associate-*r/47.8%
unsub-neg47.8%
Simplified47.8%
Taylor expanded in t around 0 20.1%
sub-neg20.1%
mul-1-neg20.1%
remove-double-neg20.1%
associate-/l*22.1%
div-sub22.1%
sub-neg22.1%
*-inverses22.1%
metadata-eval22.1%
Simplified22.1%
Taylor expanded in a around 0 2.8%
distribute-rgt1-in2.8%
metadata-eval2.8%
mul0-lft2.8%
Simplified2.8%
Final simplification2.8%
(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 63.3%
associate-*l/79.1%
Simplified79.1%
Taylor expanded in z around inf 27.9%
Final simplification27.9%
(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 2024010
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:herbie-target
(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))))