
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y x) (- z t)) (- a t))))
double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (((y - x) * (z - t)) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
def code(x, y, z, t, a): return x + (((y - x) * (z - t)) / (a - t))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - x) * Float64(z - t)) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = x + (((y - x) * (z - t)) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 28 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y x) (- z t)) (- a t))))
double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (((y - x) * (z - t)) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
def code(x, y, z, t, a): return x + (((y - x) * (z - t)) / (a - t))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - x) * Float64(z - t)) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = x + (((y - x) * (z - t)) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* (- z t) (/ (- x y) (- a t)))))
(t_2 (- x (/ (* (- y x) (- t z)) (- a t)))))
(if (<= t_2 -2e+287)
t_1
(if (<= t_2 -2e-292)
t_2
(if (<= t_2 5e-291)
(+ y (/ (* (- y x) (- a z)) t))
(if (<= t_2 2e+272) t_2 t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((z - t) * ((x - y) / (a - t)));
double t_2 = x - (((y - x) * (t - z)) / (a - t));
double tmp;
if (t_2 <= -2e+287) {
tmp = t_1;
} else if (t_2 <= -2e-292) {
tmp = t_2;
} else if (t_2 <= 5e-291) {
tmp = y + (((y - x) * (a - z)) / t);
} else if (t_2 <= 2e+272) {
tmp = t_2;
} 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 - ((z - t) * ((x - y) / (a - t)))
t_2 = x - (((y - x) * (t - z)) / (a - t))
if (t_2 <= (-2d+287)) then
tmp = t_1
else if (t_2 <= (-2d-292)) then
tmp = t_2
else if (t_2 <= 5d-291) then
tmp = y + (((y - x) * (a - z)) / t)
else if (t_2 <= 2d+272) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - ((z - t) * ((x - y) / (a - t)));
double t_2 = x - (((y - x) * (t - z)) / (a - t));
double tmp;
if (t_2 <= -2e+287) {
tmp = t_1;
} else if (t_2 <= -2e-292) {
tmp = t_2;
} else if (t_2 <= 5e-291) {
tmp = y + (((y - x) * (a - z)) / t);
} else if (t_2 <= 2e+272) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - ((z - t) * ((x - y) / (a - t))) t_2 = x - (((y - x) * (t - z)) / (a - t)) tmp = 0 if t_2 <= -2e+287: tmp = t_1 elif t_2 <= -2e-292: tmp = t_2 elif t_2 <= 5e-291: tmp = y + (((y - x) * (a - z)) / t) elif t_2 <= 2e+272: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(z - t) * Float64(Float64(x - y) / Float64(a - t)))) t_2 = Float64(x - Float64(Float64(Float64(y - x) * Float64(t - z)) / Float64(a - t))) tmp = 0.0 if (t_2 <= -2e+287) tmp = t_1; elseif (t_2 <= -2e-292) tmp = t_2; elseif (t_2 <= 5e-291) tmp = Float64(y + Float64(Float64(Float64(y - x) * Float64(a - z)) / t)); elseif (t_2 <= 2e+272) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - ((z - t) * ((x - y) / (a - t))); t_2 = x - (((y - x) * (t - z)) / (a - t)); tmp = 0.0; if (t_2 <= -2e+287) tmp = t_1; elseif (t_2 <= -2e-292) tmp = t_2; elseif (t_2 <= 5e-291) tmp = y + (((y - x) * (a - z)) / t); elseif (t_2 <= 2e+272) 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[(z - t), $MachinePrecision] * N[(N[(x - y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x - N[(N[(N[(y - x), $MachinePrecision] * N[(t - z), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -2e+287], t$95$1, If[LessEqual[t$95$2, -2e-292], t$95$2, If[LessEqual[t$95$2, 5e-291], N[(y + N[(N[(N[(y - x), $MachinePrecision] * N[(a - z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 2e+272], t$95$2, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(z - t\right) \cdot \frac{x - y}{a - t}\\
t_2 := x - \frac{\left(y - x\right) \cdot \left(t - z\right)}{a - t}\\
\mathbf{if}\;t_2 \leq -2 \cdot 10^{+287}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq -2 \cdot 10^{-292}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-291}:\\
\;\;\;\;y + \frac{\left(y - x\right) \cdot \left(a - z\right)}{t}\\
\mathbf{elif}\;t_2 \leq 2 \cdot 10^{+272}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < -2.0000000000000002e287 or 2.0000000000000001e272 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) Initial program 52.3%
associate-*l/85.1%
Simplified85.1%
if -2.0000000000000002e287 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < -2.0000000000000001e-292 or 5.0000000000000003e-291 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < 2.0000000000000001e272Initial program 97.5%
if -2.0000000000000001e-292 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < 5.0000000000000003e-291Initial program 7.2%
+-commutative7.2%
associate-*l/3.7%
fma-def3.7%
Simplified3.7%
fma-udef3.7%
associate-/r/7.2%
div-inv7.2%
clear-num7.2%
Applied egg-rr7.2%
Taylor expanded in t around inf 100.0%
associate--l+100.0%
associate-*r/100.0%
associate-*r/100.0%
div-sub99.9%
distribute-lft-out--99.9%
associate-*r/99.9%
mul-1-neg99.9%
unsub-neg99.9%
distribute-rgt-out--100.0%
Simplified100.0%
Final simplification93.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- z t) (- a t))) (t_2 (- x (/ (* (- y x) (- t z)) (- a t)))))
(if (<= t_2 -2e-292)
(fma t_1 (- y x) x)
(if (<= t_2 5e-291)
(+ y (/ (* (- y x) (- a z)) t))
(+ x (* (- y x) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (z - t) / (a - t);
double t_2 = x - (((y - x) * (t - z)) / (a - t));
double tmp;
if (t_2 <= -2e-292) {
tmp = fma(t_1, (y - x), x);
} else if (t_2 <= 5e-291) {
tmp = y + (((y - x) * (a - z)) / t);
} else {
tmp = x + ((y - x) * t_1);
}
return tmp;
}
function code(x, y, z, t, a) t_1 = Float64(Float64(z - t) / Float64(a - t)) t_2 = Float64(x - Float64(Float64(Float64(y - x) * Float64(t - z)) / Float64(a - t))) tmp = 0.0 if (t_2 <= -2e-292) tmp = fma(t_1, Float64(y - x), x); elseif (t_2 <= 5e-291) tmp = Float64(y + Float64(Float64(Float64(y - x) * Float64(a - z)) / t)); else tmp = Float64(x + Float64(Float64(y - x) * t_1)); end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x - N[(N[(N[(y - x), $MachinePrecision] * N[(t - z), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -2e-292], N[(t$95$1 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[t$95$2, 5e-291], N[(y + N[(N[(N[(y - x), $MachinePrecision] * N[(a - z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y - x), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z - t}{a - t}\\
t_2 := x - \frac{\left(y - x\right) \cdot \left(t - z\right)}{a - t}\\
\mathbf{if}\;t_2 \leq -2 \cdot 10^{-292}:\\
\;\;\;\;\mathsf{fma}\left(t_1, y - x, x\right)\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-291}:\\
\;\;\;\;y + \frac{\left(y - x\right) \cdot \left(a - z\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - x\right) \cdot t_1\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < -2.0000000000000001e-292Initial program 80.1%
+-commutative80.1%
*-commutative80.1%
associate-/l*85.9%
associate-/r/92.8%
fma-def92.8%
Simplified92.8%
if -2.0000000000000001e-292 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < 5.0000000000000003e-291Initial program 7.2%
+-commutative7.2%
associate-*l/3.7%
fma-def3.7%
Simplified3.7%
fma-udef3.7%
associate-/r/7.2%
div-inv7.2%
clear-num7.2%
Applied egg-rr7.2%
Taylor expanded in t around inf 100.0%
associate--l+100.0%
associate-*r/100.0%
associate-*r/100.0%
div-sub99.9%
distribute-lft-out--99.9%
associate-*r/99.9%
mul-1-neg99.9%
unsub-neg99.9%
distribute-rgt-out--100.0%
Simplified100.0%
if 5.0000000000000003e-291 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) Initial program 79.5%
+-commutative79.5%
associate-*l/87.4%
fma-def87.5%
Simplified87.5%
fma-udef87.4%
associate-/r/91.1%
div-inv91.0%
clear-num91.0%
Applied egg-rr91.0%
Final simplification92.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (- 1.0 (/ z a))))
(t_2 (* (- z a) (/ x t)))
(t_3 (* z (/ y (- a t)))))
(if (<= t -1.9e+161)
y
(if (<= t -7.2e+71)
t_2
(if (<= t -9.2e-32)
t_1
(if (<= t -9.4e-153)
t_3
(if (<= t 1.1e-183)
t_1
(if (<= t 3e-161)
t_3
(if (<= t 3e+52)
t_1
(if (<= t 6.2e+126) t_3 (if (<= t 1.7e+162) t_2 y)))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double t_2 = (z - a) * (x / t);
double t_3 = z * (y / (a - t));
double tmp;
if (t <= -1.9e+161) {
tmp = y;
} else if (t <= -7.2e+71) {
tmp = t_2;
} else if (t <= -9.2e-32) {
tmp = t_1;
} else if (t <= -9.4e-153) {
tmp = t_3;
} else if (t <= 1.1e-183) {
tmp = t_1;
} else if (t <= 3e-161) {
tmp = t_3;
} else if (t <= 3e+52) {
tmp = t_1;
} else if (t <= 6.2e+126) {
tmp = t_3;
} else if (t <= 1.7e+162) {
tmp = t_2;
} else {
tmp = 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) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = x * (1.0d0 - (z / a))
t_2 = (z - a) * (x / t)
t_3 = z * (y / (a - t))
if (t <= (-1.9d+161)) then
tmp = y
else if (t <= (-7.2d+71)) then
tmp = t_2
else if (t <= (-9.2d-32)) then
tmp = t_1
else if (t <= (-9.4d-153)) then
tmp = t_3
else if (t <= 1.1d-183) then
tmp = t_1
else if (t <= 3d-161) then
tmp = t_3
else if (t <= 3d+52) then
tmp = t_1
else if (t <= 6.2d+126) then
tmp = t_3
else if (t <= 1.7d+162) then
tmp = t_2
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double t_2 = (z - a) * (x / t);
double t_3 = z * (y / (a - t));
double tmp;
if (t <= -1.9e+161) {
tmp = y;
} else if (t <= -7.2e+71) {
tmp = t_2;
} else if (t <= -9.2e-32) {
tmp = t_1;
} else if (t <= -9.4e-153) {
tmp = t_3;
} else if (t <= 1.1e-183) {
tmp = t_1;
} else if (t <= 3e-161) {
tmp = t_3;
} else if (t <= 3e+52) {
tmp = t_1;
} else if (t <= 6.2e+126) {
tmp = t_3;
} else if (t <= 1.7e+162) {
tmp = t_2;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (1.0 - (z / a)) t_2 = (z - a) * (x / t) t_3 = z * (y / (a - t)) tmp = 0 if t <= -1.9e+161: tmp = y elif t <= -7.2e+71: tmp = t_2 elif t <= -9.2e-32: tmp = t_1 elif t <= -9.4e-153: tmp = t_3 elif t <= 1.1e-183: tmp = t_1 elif t <= 3e-161: tmp = t_3 elif t <= 3e+52: tmp = t_1 elif t <= 6.2e+126: tmp = t_3 elif t <= 1.7e+162: tmp = t_2 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(1.0 - Float64(z / a))) t_2 = Float64(Float64(z - a) * Float64(x / t)) t_3 = Float64(z * Float64(y / Float64(a - t))) tmp = 0.0 if (t <= -1.9e+161) tmp = y; elseif (t <= -7.2e+71) tmp = t_2; elseif (t <= -9.2e-32) tmp = t_1; elseif (t <= -9.4e-153) tmp = t_3; elseif (t <= 1.1e-183) tmp = t_1; elseif (t <= 3e-161) tmp = t_3; elseif (t <= 3e+52) tmp = t_1; elseif (t <= 6.2e+126) tmp = t_3; elseif (t <= 1.7e+162) tmp = t_2; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (1.0 - (z / a)); t_2 = (z - a) * (x / t); t_3 = z * (y / (a - t)); tmp = 0.0; if (t <= -1.9e+161) tmp = y; elseif (t <= -7.2e+71) tmp = t_2; elseif (t <= -9.2e-32) tmp = t_1; elseif (t <= -9.4e-153) tmp = t_3; elseif (t <= 1.1e-183) tmp = t_1; elseif (t <= 3e-161) tmp = t_3; elseif (t <= 3e+52) tmp = t_1; elseif (t <= 6.2e+126) tmp = t_3; elseif (t <= 1.7e+162) tmp = t_2; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(z * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.9e+161], y, If[LessEqual[t, -7.2e+71], t$95$2, If[LessEqual[t, -9.2e-32], t$95$1, If[LessEqual[t, -9.4e-153], t$95$3, If[LessEqual[t, 1.1e-183], t$95$1, If[LessEqual[t, 3e-161], t$95$3, If[LessEqual[t, 3e+52], t$95$1, If[LessEqual[t, 6.2e+126], t$95$3, If[LessEqual[t, 1.7e+162], t$95$2, y]]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(1 - \frac{z}{a}\right)\\
t_2 := \left(z - a\right) \cdot \frac{x}{t}\\
t_3 := z \cdot \frac{y}{a - t}\\
\mathbf{if}\;t \leq -1.9 \cdot 10^{+161}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -7.2 \cdot 10^{+71}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq -9.2 \cdot 10^{-32}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -9.4 \cdot 10^{-153}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;t \leq 1.1 \cdot 10^{-183}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 3 \cdot 10^{-161}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;t \leq 3 \cdot 10^{+52}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{+126}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;t \leq 1.7 \cdot 10^{+162}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -1.9000000000000001e161 or 1.70000000000000001e162 < t Initial program 32.0%
associate-*l/56.3%
Simplified56.3%
Taylor expanded in t around inf 57.9%
if -1.9000000000000001e161 < t < -7.1999999999999999e71 or 6.2e126 < t < 1.70000000000000001e162Initial program 76.5%
associate-*l/76.1%
Simplified76.1%
Taylor expanded in t around inf 66.5%
associate--l+66.5%
associate-*r/66.5%
associate-*r/66.5%
div-sub66.5%
distribute-lft-out--66.5%
associate-*r/66.5%
mul-1-neg66.5%
unsub-neg66.5%
distribute-rgt-out--66.5%
associate-/l*73.1%
Simplified73.1%
Taylor expanded in y around 0 47.0%
associate-*l/53.2%
Simplified53.2%
if -7.1999999999999999e71 < t < -9.2000000000000002e-32 or -9.3999999999999998e-153 < t < 1.1e-183 or 2.99999999999999989e-161 < t < 3e52Initial program 88.9%
associate-*l/89.4%
Simplified89.4%
Taylor expanded in t around 0 71.0%
associate-/l*72.7%
Simplified72.7%
Taylor expanded in x around inf 62.8%
mul-1-neg62.8%
unsub-neg62.8%
Simplified62.8%
if -9.2000000000000002e-32 < t < -9.3999999999999998e-153 or 1.1e-183 < t < 2.99999999999999989e-161 or 3e52 < t < 6.2e126Initial program 77.3%
+-commutative77.3%
associate-*l/85.9%
fma-def86.1%
Simplified86.1%
fma-udef85.9%
associate-/r/85.2%
div-inv83.7%
clear-num83.8%
Applied egg-rr83.8%
Taylor expanded in y around inf 73.9%
div-sub73.9%
Simplified73.9%
Taylor expanded in z around inf 49.6%
associate-/l*57.5%
associate-/r/58.3%
Simplified58.3%
Final simplification60.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ (* (- y x) (- t z)) (- a t)))))
(if (or (<= t_1 -2e-292) (not (<= t_1 5e-291)))
(+ x (* (- y x) (/ (- z t) (- a t))))
(+ y (/ (* (- y x) (- a z)) t)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (((y - x) * (t - z)) / (a - t));
double tmp;
if ((t_1 <= -2e-292) || !(t_1 <= 5e-291)) {
tmp = x + ((y - x) * ((z - t) / (a - t)));
} else {
tmp = y + (((y - x) * (a - z)) / t);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x - (((y - x) * (t - z)) / (a - t))
if ((t_1 <= (-2d-292)) .or. (.not. (t_1 <= 5d-291))) then
tmp = x + ((y - x) * ((z - t) / (a - t)))
else
tmp = y + (((y - x) * (a - z)) / t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (((y - x) * (t - z)) / (a - t));
double tmp;
if ((t_1 <= -2e-292) || !(t_1 <= 5e-291)) {
tmp = x + ((y - x) * ((z - t) / (a - t)));
} else {
tmp = y + (((y - x) * (a - z)) / t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (((y - x) * (t - z)) / (a - t)) tmp = 0 if (t_1 <= -2e-292) or not (t_1 <= 5e-291): tmp = x + ((y - x) * ((z - t) / (a - t))) else: tmp = y + (((y - x) * (a - z)) / t) return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(Float64(y - x) * Float64(t - z)) / Float64(a - t))) tmp = 0.0 if ((t_1 <= -2e-292) || !(t_1 <= 5e-291)) tmp = Float64(x + Float64(Float64(y - x) * Float64(Float64(z - t) / Float64(a - t)))); else tmp = Float64(y + Float64(Float64(Float64(y - x) * Float64(a - z)) / t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (((y - x) * (t - z)) / (a - t)); tmp = 0.0; if ((t_1 <= -2e-292) || ~((t_1 <= 5e-291))) tmp = x + ((y - x) * ((z - t) / (a - t))); else tmp = y + (((y - x) * (a - z)) / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(N[(y - x), $MachinePrecision] * N[(t - z), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -2e-292], N[Not[LessEqual[t$95$1, 5e-291]], $MachinePrecision]], N[(x + N[(N[(y - x), $MachinePrecision] * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + N[(N[(N[(y - x), $MachinePrecision] * N[(a - z), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\left(y - x\right) \cdot \left(t - z\right)}{a - t}\\
\mathbf{if}\;t_1 \leq -2 \cdot 10^{-292} \lor \neg \left(t_1 \leq 5 \cdot 10^{-291}\right):\\
\;\;\;\;x + \left(y - x\right) \cdot \frac{z - t}{a - t}\\
\mathbf{else}:\\
\;\;\;\;y + \frac{\left(y - x\right) \cdot \left(a - z\right)}{t}\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < -2.0000000000000001e-292 or 5.0000000000000003e-291 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) Initial program 79.8%
+-commutative79.8%
associate-*l/87.3%
fma-def87.4%
Simplified87.4%
fma-udef87.3%
associate-/r/92.2%
div-inv91.8%
clear-num91.9%
Applied egg-rr91.9%
if -2.0000000000000001e-292 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) (-.f64 z t)) (-.f64 a t))) < 5.0000000000000003e-291Initial program 7.2%
+-commutative7.2%
associate-*l/3.7%
fma-def3.7%
Simplified3.7%
fma-udef3.7%
associate-/r/7.2%
div-inv7.2%
clear-num7.2%
Applied egg-rr7.2%
Taylor expanded in t around inf 100.0%
associate--l+100.0%
associate-*r/100.0%
associate-*r/100.0%
div-sub99.9%
distribute-lft-out--99.9%
associate-*r/99.9%
mul-1-neg99.9%
unsub-neg99.9%
distribute-rgt-out--100.0%
Simplified100.0%
Final simplification92.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ y (/ a z))))
(if (<= t -1.45e+107)
y
(if (<= t -2.35e-22)
x
(if (<= t -1.9e-271)
t_1
(if (<= t 1.4e-268)
x
(if (<= t 3e-233)
t_1
(if (<= t 1.12e-184)
x
(if (<= t 1.3e-152) t_1 (if (<= t 6e+64) x y))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y / (a / z);
double tmp;
if (t <= -1.45e+107) {
tmp = y;
} else if (t <= -2.35e-22) {
tmp = x;
} else if (t <= -1.9e-271) {
tmp = t_1;
} else if (t <= 1.4e-268) {
tmp = x;
} else if (t <= 3e-233) {
tmp = t_1;
} else if (t <= 1.12e-184) {
tmp = x;
} else if (t <= 1.3e-152) {
tmp = t_1;
} else if (t <= 6e+64) {
tmp = x;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = y / (a / z)
if (t <= (-1.45d+107)) then
tmp = y
else if (t <= (-2.35d-22)) then
tmp = x
else if (t <= (-1.9d-271)) then
tmp = t_1
else if (t <= 1.4d-268) then
tmp = x
else if (t <= 3d-233) then
tmp = t_1
else if (t <= 1.12d-184) then
tmp = x
else if (t <= 1.3d-152) then
tmp = t_1
else if (t <= 6d+64) then
tmp = x
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = y / (a / z);
double tmp;
if (t <= -1.45e+107) {
tmp = y;
} else if (t <= -2.35e-22) {
tmp = x;
} else if (t <= -1.9e-271) {
tmp = t_1;
} else if (t <= 1.4e-268) {
tmp = x;
} else if (t <= 3e-233) {
tmp = t_1;
} else if (t <= 1.12e-184) {
tmp = x;
} else if (t <= 1.3e-152) {
tmp = t_1;
} else if (t <= 6e+64) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y / (a / z) tmp = 0 if t <= -1.45e+107: tmp = y elif t <= -2.35e-22: tmp = x elif t <= -1.9e-271: tmp = t_1 elif t <= 1.4e-268: tmp = x elif t <= 3e-233: tmp = t_1 elif t <= 1.12e-184: tmp = x elif t <= 1.3e-152: tmp = t_1 elif t <= 6e+64: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(y / Float64(a / z)) tmp = 0.0 if (t <= -1.45e+107) tmp = y; elseif (t <= -2.35e-22) tmp = x; elseif (t <= -1.9e-271) tmp = t_1; elseif (t <= 1.4e-268) tmp = x; elseif (t <= 3e-233) tmp = t_1; elseif (t <= 1.12e-184) tmp = x; elseif (t <= 1.3e-152) tmp = t_1; elseif (t <= 6e+64) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y / (a / z); tmp = 0.0; if (t <= -1.45e+107) tmp = y; elseif (t <= -2.35e-22) tmp = x; elseif (t <= -1.9e-271) tmp = t_1; elseif (t <= 1.4e-268) tmp = x; elseif (t <= 3e-233) tmp = t_1; elseif (t <= 1.12e-184) tmp = x; elseif (t <= 1.3e-152) tmp = t_1; elseif (t <= 6e+64) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.45e+107], y, If[LessEqual[t, -2.35e-22], x, If[LessEqual[t, -1.9e-271], t$95$1, If[LessEqual[t, 1.4e-268], x, If[LessEqual[t, 3e-233], t$95$1, If[LessEqual[t, 1.12e-184], x, If[LessEqual[t, 1.3e-152], t$95$1, If[LessEqual[t, 6e+64], x, y]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{\frac{a}{z}}\\
\mathbf{if}\;t \leq -1.45 \cdot 10^{+107}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -2.35 \cdot 10^{-22}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq -1.9 \cdot 10^{-271}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 1.4 \cdot 10^{-268}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 3 \cdot 10^{-233}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 1.12 \cdot 10^{-184}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 1.3 \cdot 10^{-152}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 6 \cdot 10^{+64}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -1.44999999999999994e107 or 6.0000000000000004e64 < t Initial program 47.0%
associate-*l/64.0%
Simplified64.0%
Taylor expanded in t around inf 45.2%
if -1.44999999999999994e107 < t < -2.3500000000000001e-22 or -1.90000000000000005e-271 < t < 1.40000000000000008e-268 or 2.99999999999999999e-233 < t < 1.11999999999999997e-184 or 1.30000000000000006e-152 < t < 6.0000000000000004e64Initial program 86.3%
associate-*l/89.3%
Simplified89.3%
Taylor expanded in a around inf 43.3%
if -2.3500000000000001e-22 < t < -1.90000000000000005e-271 or 1.40000000000000008e-268 < t < 2.99999999999999999e-233 or 1.11999999999999997e-184 < t < 1.30000000000000006e-152Initial program 89.6%
associate-*l/87.8%
Simplified87.8%
Taylor expanded in z around -inf 71.8%
Taylor expanded in y around inf 55.8%
*-commutative55.8%
Simplified55.8%
Taylor expanded in a around inf 49.0%
associate-/l*48.4%
Simplified48.4%
Final simplification45.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* y z) a)))
(if (<= t -2.9e+107)
y
(if (<= t -1.15e-21)
x
(if (<= t -2.55e-273)
(/ y (/ a z))
(if (<= t 2.9e-267)
x
(if (<= t 8e-234)
t_1
(if (<= t 1.05e-183)
x
(if (<= t 4e-154) t_1 (if (<= t 4e+68) x y))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y * z) / a;
double tmp;
if (t <= -2.9e+107) {
tmp = y;
} else if (t <= -1.15e-21) {
tmp = x;
} else if (t <= -2.55e-273) {
tmp = y / (a / z);
} else if (t <= 2.9e-267) {
tmp = x;
} else if (t <= 8e-234) {
tmp = t_1;
} else if (t <= 1.05e-183) {
tmp = x;
} else if (t <= 4e-154) {
tmp = t_1;
} else if (t <= 4e+68) {
tmp = x;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = (y * z) / a
if (t <= (-2.9d+107)) then
tmp = y
else if (t <= (-1.15d-21)) then
tmp = x
else if (t <= (-2.55d-273)) then
tmp = y / (a / z)
else if (t <= 2.9d-267) then
tmp = x
else if (t <= 8d-234) then
tmp = t_1
else if (t <= 1.05d-183) then
tmp = x
else if (t <= 4d-154) then
tmp = t_1
else if (t <= 4d+68) then
tmp = x
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (y * z) / a;
double tmp;
if (t <= -2.9e+107) {
tmp = y;
} else if (t <= -1.15e-21) {
tmp = x;
} else if (t <= -2.55e-273) {
tmp = y / (a / z);
} else if (t <= 2.9e-267) {
tmp = x;
} else if (t <= 8e-234) {
tmp = t_1;
} else if (t <= 1.05e-183) {
tmp = x;
} else if (t <= 4e-154) {
tmp = t_1;
} else if (t <= 4e+68) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y * z) / a tmp = 0 if t <= -2.9e+107: tmp = y elif t <= -1.15e-21: tmp = x elif t <= -2.55e-273: tmp = y / (a / z) elif t <= 2.9e-267: tmp = x elif t <= 8e-234: tmp = t_1 elif t <= 1.05e-183: tmp = x elif t <= 4e-154: tmp = t_1 elif t <= 4e+68: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y * z) / a) tmp = 0.0 if (t <= -2.9e+107) tmp = y; elseif (t <= -1.15e-21) tmp = x; elseif (t <= -2.55e-273) tmp = Float64(y / Float64(a / z)); elseif (t <= 2.9e-267) tmp = x; elseif (t <= 8e-234) tmp = t_1; elseif (t <= 1.05e-183) tmp = x; elseif (t <= 4e-154) tmp = t_1; elseif (t <= 4e+68) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y * z) / a; tmp = 0.0; if (t <= -2.9e+107) tmp = y; elseif (t <= -1.15e-21) tmp = x; elseif (t <= -2.55e-273) tmp = y / (a / z); elseif (t <= 2.9e-267) tmp = x; elseif (t <= 8e-234) tmp = t_1; elseif (t <= 1.05e-183) tmp = x; elseif (t <= 4e-154) tmp = t_1; elseif (t <= 4e+68) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]}, If[LessEqual[t, -2.9e+107], y, If[LessEqual[t, -1.15e-21], x, If[LessEqual[t, -2.55e-273], N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.9e-267], x, If[LessEqual[t, 8e-234], t$95$1, If[LessEqual[t, 1.05e-183], x, If[LessEqual[t, 4e-154], t$95$1, If[LessEqual[t, 4e+68], x, y]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y \cdot z}{a}\\
\mathbf{if}\;t \leq -2.9 \cdot 10^{+107}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -1.15 \cdot 10^{-21}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq -2.55 \cdot 10^{-273}:\\
\;\;\;\;\frac{y}{\frac{a}{z}}\\
\mathbf{elif}\;t \leq 2.9 \cdot 10^{-267}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 8 \cdot 10^{-234}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 1.05 \cdot 10^{-183}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 4 \cdot 10^{-154}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 4 \cdot 10^{+68}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -2.89999999999999988e107 or 3.99999999999999981e68 < t Initial program 47.0%
associate-*l/64.0%
Simplified64.0%
Taylor expanded in t around inf 45.2%
if -2.89999999999999988e107 < t < -1.15e-21 or -2.54999999999999983e-273 < t < 2.90000000000000021e-267 or 7.9999999999999997e-234 < t < 1.0500000000000001e-183 or 3.9999999999999999e-154 < t < 3.99999999999999981e68Initial program 86.3%
associate-*l/89.3%
Simplified89.3%
Taylor expanded in a around inf 43.3%
if -1.15e-21 < t < -2.54999999999999983e-273Initial program 87.0%
associate-*l/89.3%
Simplified89.3%
Taylor expanded in z around -inf 64.9%
Taylor expanded in y around inf 48.5%
*-commutative48.5%
Simplified48.5%
Taylor expanded in a around inf 38.6%
associate-/l*42.8%
Simplified42.8%
if 2.90000000000000021e-267 < t < 7.9999999999999997e-234 or 1.0500000000000001e-183 < t < 3.9999999999999999e-154Initial program 94.8%
associate-*l/84.8%
Simplified84.8%
Taylor expanded in z around -inf 85.1%
Taylor expanded in y around inf 70.0%
*-commutative70.0%
Simplified70.0%
Taylor expanded in a around inf 69.2%
Final simplification45.8%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ y (/ a z))))
(if (<= z -7.6e+143)
(/ (- y) (/ t z))
(if (<= z -2.7e-21)
t_1
(if (<= z -4.5e-157)
y
(if (<= z 2.6e-119)
x
(if (<= z 3.9e-41)
y
(if (<= z 1.2e+63)
x
(if (<= z 3.7e+193) t_1 (/ (- x) (/ a z)))))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y / (a / z);
double tmp;
if (z <= -7.6e+143) {
tmp = -y / (t / z);
} else if (z <= -2.7e-21) {
tmp = t_1;
} else if (z <= -4.5e-157) {
tmp = y;
} else if (z <= 2.6e-119) {
tmp = x;
} else if (z <= 3.9e-41) {
tmp = y;
} else if (z <= 1.2e+63) {
tmp = x;
} else if (z <= 3.7e+193) {
tmp = t_1;
} else {
tmp = -x / (a / z);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = y / (a / z)
if (z <= (-7.6d+143)) then
tmp = -y / (t / z)
else if (z <= (-2.7d-21)) then
tmp = t_1
else if (z <= (-4.5d-157)) then
tmp = y
else if (z <= 2.6d-119) then
tmp = x
else if (z <= 3.9d-41) then
tmp = y
else if (z <= 1.2d+63) then
tmp = x
else if (z <= 3.7d+193) then
tmp = t_1
else
tmp = -x / (a / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = y / (a / z);
double tmp;
if (z <= -7.6e+143) {
tmp = -y / (t / z);
} else if (z <= -2.7e-21) {
tmp = t_1;
} else if (z <= -4.5e-157) {
tmp = y;
} else if (z <= 2.6e-119) {
tmp = x;
} else if (z <= 3.9e-41) {
tmp = y;
} else if (z <= 1.2e+63) {
tmp = x;
} else if (z <= 3.7e+193) {
tmp = t_1;
} else {
tmp = -x / (a / z);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y / (a / z) tmp = 0 if z <= -7.6e+143: tmp = -y / (t / z) elif z <= -2.7e-21: tmp = t_1 elif z <= -4.5e-157: tmp = y elif z <= 2.6e-119: tmp = x elif z <= 3.9e-41: tmp = y elif z <= 1.2e+63: tmp = x elif z <= 3.7e+193: tmp = t_1 else: tmp = -x / (a / z) return tmp
function code(x, y, z, t, a) t_1 = Float64(y / Float64(a / z)) tmp = 0.0 if (z <= -7.6e+143) tmp = Float64(Float64(-y) / Float64(t / z)); elseif (z <= -2.7e-21) tmp = t_1; elseif (z <= -4.5e-157) tmp = y; elseif (z <= 2.6e-119) tmp = x; elseif (z <= 3.9e-41) tmp = y; elseif (z <= 1.2e+63) tmp = x; elseif (z <= 3.7e+193) tmp = t_1; else tmp = Float64(Float64(-x) / Float64(a / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y / (a / z); tmp = 0.0; if (z <= -7.6e+143) tmp = -y / (t / z); elseif (z <= -2.7e-21) tmp = t_1; elseif (z <= -4.5e-157) tmp = y; elseif (z <= 2.6e-119) tmp = x; elseif (z <= 3.9e-41) tmp = y; elseif (z <= 1.2e+63) tmp = x; elseif (z <= 3.7e+193) tmp = t_1; else tmp = -x / (a / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -7.6e+143], N[((-y) / N[(t / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -2.7e-21], t$95$1, If[LessEqual[z, -4.5e-157], y, If[LessEqual[z, 2.6e-119], x, If[LessEqual[z, 3.9e-41], y, If[LessEqual[z, 1.2e+63], x, If[LessEqual[z, 3.7e+193], t$95$1, N[((-x) / N[(a / z), $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{\frac{a}{z}}\\
\mathbf{if}\;z \leq -7.6 \cdot 10^{+143}:\\
\;\;\;\;\frac{-y}{\frac{t}{z}}\\
\mathbf{elif}\;z \leq -2.7 \cdot 10^{-21}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -4.5 \cdot 10^{-157}:\\
\;\;\;\;y\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{-119}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 3.9 \cdot 10^{-41}:\\
\;\;\;\;y\\
\mathbf{elif}\;z \leq 1.2 \cdot 10^{+63}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 3.7 \cdot 10^{+193}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{-x}{\frac{a}{z}}\\
\end{array}
\end{array}
if z < -7.60000000000000001e143Initial program 77.2%
associate-*l/89.4%
Simplified89.4%
Taylor expanded in z around -inf 64.0%
Taylor expanded in y around inf 40.5%
*-commutative40.5%
Simplified40.5%
Taylor expanded in a around 0 43.7%
mul-1-neg43.7%
associate-/l*50.0%
distribute-neg-frac50.0%
Simplified50.0%
if -7.60000000000000001e143 < z < -2.7000000000000001e-21 or 1.2e63 < z < 3.7000000000000003e193Initial program 75.4%
associate-*l/83.5%
Simplified83.5%
Taylor expanded in z around -inf 60.2%
Taylor expanded in y around inf 37.7%
*-commutative37.7%
Simplified37.7%
Taylor expanded in a around inf 32.8%
associate-/l*37.9%
Simplified37.9%
if -2.7000000000000001e-21 < z < -4.49999999999999999e-157 or 2.60000000000000012e-119 < z < 3.89999999999999991e-41Initial program 59.7%
associate-*l/71.4%
Simplified71.4%
Taylor expanded in t around inf 46.7%
if -4.49999999999999999e-157 < z < 2.60000000000000012e-119 or 3.89999999999999991e-41 < z < 1.2e63Initial program 74.9%
associate-*l/74.2%
Simplified74.2%
Taylor expanded in a around inf 48.3%
if 3.7000000000000003e193 < z Initial program 87.9%
associate-*l/99.4%
Simplified99.4%
Taylor expanded in t around 0 62.1%
associate-/l*70.7%
Simplified70.7%
Taylor expanded in x around inf 55.8%
mul-1-neg55.8%
unsub-neg55.8%
Simplified55.8%
Taylor expanded in z around inf 40.7%
mul-1-neg40.7%
associate-/l*52.8%
distribute-neg-frac52.8%
Simplified52.8%
Final simplification46.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* y (/ (- z t) (- a t)))) (t_2 (* (- z t) (/ y (- a t)))))
(if (<= t -2.25e+138)
t_1
(if (<= t -3.2e+66)
(* z (/ (- y x) (- a t)))
(if (<= t -1.52e-97)
t_2
(if (<= t 1.65e-68)
(+ x (/ z (/ a (- y x))))
(if (<= t 1.22e-8)
t_2
(if (<= t 4.6e+49) (+ x (/ z (- (/ a x)))) t_1))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * ((z - t) / (a - t));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -2.25e+138) {
tmp = t_1;
} else if (t <= -3.2e+66) {
tmp = z * ((y - x) / (a - t));
} else if (t <= -1.52e-97) {
tmp = t_2;
} else if (t <= 1.65e-68) {
tmp = x + (z / (a / (y - x)));
} else if (t <= 1.22e-8) {
tmp = t_2;
} else if (t <= 4.6e+49) {
tmp = x + (z / -(a / 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 = y * ((z - t) / (a - t))
t_2 = (z - t) * (y / (a - t))
if (t <= (-2.25d+138)) then
tmp = t_1
else if (t <= (-3.2d+66)) then
tmp = z * ((y - x) / (a - t))
else if (t <= (-1.52d-97)) then
tmp = t_2
else if (t <= 1.65d-68) then
tmp = x + (z / (a / (y - x)))
else if (t <= 1.22d-8) then
tmp = t_2
else if (t <= 4.6d+49) then
tmp = x + (z / -(a / 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 = y * ((z - t) / (a - t));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -2.25e+138) {
tmp = t_1;
} else if (t <= -3.2e+66) {
tmp = z * ((y - x) / (a - t));
} else if (t <= -1.52e-97) {
tmp = t_2;
} else if (t <= 1.65e-68) {
tmp = x + (z / (a / (y - x)));
} else if (t <= 1.22e-8) {
tmp = t_2;
} else if (t <= 4.6e+49) {
tmp = x + (z / -(a / x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * ((z - t) / (a - t)) t_2 = (z - t) * (y / (a - t)) tmp = 0 if t <= -2.25e+138: tmp = t_1 elif t <= -3.2e+66: tmp = z * ((y - x) / (a - t)) elif t <= -1.52e-97: tmp = t_2 elif t <= 1.65e-68: tmp = x + (z / (a / (y - x))) elif t <= 1.22e-8: tmp = t_2 elif t <= 4.6e+49: tmp = x + (z / -(a / x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(Float64(z - t) / Float64(a - t))) t_2 = Float64(Float64(z - t) * Float64(y / Float64(a - t))) tmp = 0.0 if (t <= -2.25e+138) tmp = t_1; elseif (t <= -3.2e+66) tmp = Float64(z * Float64(Float64(y - x) / Float64(a - t))); elseif (t <= -1.52e-97) tmp = t_2; elseif (t <= 1.65e-68) tmp = Float64(x + Float64(z / Float64(a / Float64(y - x)))); elseif (t <= 1.22e-8) tmp = t_2; elseif (t <= 4.6e+49) tmp = Float64(x + Float64(z / Float64(-Float64(a / x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * ((z - t) / (a - t)); t_2 = (z - t) * (y / (a - t)); tmp = 0.0; if (t <= -2.25e+138) tmp = t_1; elseif (t <= -3.2e+66) tmp = z * ((y - x) / (a - t)); elseif (t <= -1.52e-97) tmp = t_2; elseif (t <= 1.65e-68) tmp = x + (z / (a / (y - x))); elseif (t <= 1.22e-8) tmp = t_2; elseif (t <= 4.6e+49) tmp = x + (z / -(a / x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.25e+138], t$95$1, If[LessEqual[t, -3.2e+66], N[(z * N[(N[(y - x), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -1.52e-97], t$95$2, If[LessEqual[t, 1.65e-68], N[(x + N[(z / N[(a / N[(y - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.22e-8], t$95$2, If[LessEqual[t, 4.6e+49], N[(x + N[(z / (-N[(a / x), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z - t}{a - t}\\
t_2 := \left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{if}\;t \leq -2.25 \cdot 10^{+138}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -3.2 \cdot 10^{+66}:\\
\;\;\;\;z \cdot \frac{y - x}{a - t}\\
\mathbf{elif}\;t \leq -1.52 \cdot 10^{-97}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 1.65 \cdot 10^{-68}:\\
\;\;\;\;x + \frac{z}{\frac{a}{y - x}}\\
\mathbf{elif}\;t \leq 1.22 \cdot 10^{-8}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 4.6 \cdot 10^{+49}:\\
\;\;\;\;x + \frac{z}{-\frac{a}{x}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -2.24999999999999991e138 or 4.60000000000000004e49 < t Initial program 45.6%
+-commutative45.6%
associate-*l/64.6%
fma-def64.7%
Simplified64.7%
fma-udef64.6%
associate-/r/71.2%
div-inv71.2%
clear-num71.2%
Applied egg-rr71.2%
Taylor expanded in y around inf 67.9%
div-sub67.9%
Simplified67.9%
if -2.24999999999999991e138 < t < -3.2e66Initial program 83.7%
associate-*l/83.1%
Simplified83.1%
Taylor expanded in z around inf 61.7%
div-sub61.7%
Simplified61.7%
if -3.2e66 < t < -1.5200000000000001e-97 or 1.6499999999999999e-68 < t < 1.22e-8Initial program 78.3%
associate-*l/82.6%
Simplified82.6%
Taylor expanded in x around 0 61.4%
associate-/l*61.5%
associate-/r/63.6%
Simplified63.6%
if -1.5200000000000001e-97 < t < 1.6499999999999999e-68Initial program 93.8%
associate-*l/90.5%
Simplified90.5%
Taylor expanded in t around 0 86.3%
associate-/l*85.8%
Simplified85.8%
if 1.22e-8 < t < 4.60000000000000004e49Initial program 74.9%
associate-*l/93.3%
Simplified93.3%
Taylor expanded in t around 0 61.2%
associate-/l*73.8%
Simplified73.8%
Taylor expanded in y around 0 80.5%
associate-*r/80.5%
neg-mul-180.5%
Simplified80.5%
Final simplification74.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* y (/ (- z t) (- a t)))) (t_2 (* (- z t) (/ y (- a t)))))
(if (<= t -2.25e+138)
t_1
(if (<= t -3.3e+66)
(* z (/ (- y x) (- a t)))
(if (<= t -2.45e-98)
t_2
(if (<= t 1.42e-67)
(- x (/ (* z (- x y)) a))
(if (<= t 0.000235)
t_2
(if (<= t 1.45e+50) (+ x (/ z (- (/ a x)))) t_1))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * ((z - t) / (a - t));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -2.25e+138) {
tmp = t_1;
} else if (t <= -3.3e+66) {
tmp = z * ((y - x) / (a - t));
} else if (t <= -2.45e-98) {
tmp = t_2;
} else if (t <= 1.42e-67) {
tmp = x - ((z * (x - y)) / a);
} else if (t <= 0.000235) {
tmp = t_2;
} else if (t <= 1.45e+50) {
tmp = x + (z / -(a / 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 = y * ((z - t) / (a - t))
t_2 = (z - t) * (y / (a - t))
if (t <= (-2.25d+138)) then
tmp = t_1
else if (t <= (-3.3d+66)) then
tmp = z * ((y - x) / (a - t))
else if (t <= (-2.45d-98)) then
tmp = t_2
else if (t <= 1.42d-67) then
tmp = x - ((z * (x - y)) / a)
else if (t <= 0.000235d0) then
tmp = t_2
else if (t <= 1.45d+50) then
tmp = x + (z / -(a / 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 = y * ((z - t) / (a - t));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -2.25e+138) {
tmp = t_1;
} else if (t <= -3.3e+66) {
tmp = z * ((y - x) / (a - t));
} else if (t <= -2.45e-98) {
tmp = t_2;
} else if (t <= 1.42e-67) {
tmp = x - ((z * (x - y)) / a);
} else if (t <= 0.000235) {
tmp = t_2;
} else if (t <= 1.45e+50) {
tmp = x + (z / -(a / x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * ((z - t) / (a - t)) t_2 = (z - t) * (y / (a - t)) tmp = 0 if t <= -2.25e+138: tmp = t_1 elif t <= -3.3e+66: tmp = z * ((y - x) / (a - t)) elif t <= -2.45e-98: tmp = t_2 elif t <= 1.42e-67: tmp = x - ((z * (x - y)) / a) elif t <= 0.000235: tmp = t_2 elif t <= 1.45e+50: tmp = x + (z / -(a / x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(Float64(z - t) / Float64(a - t))) t_2 = Float64(Float64(z - t) * Float64(y / Float64(a - t))) tmp = 0.0 if (t <= -2.25e+138) tmp = t_1; elseif (t <= -3.3e+66) tmp = Float64(z * Float64(Float64(y - x) / Float64(a - t))); elseif (t <= -2.45e-98) tmp = t_2; elseif (t <= 1.42e-67) tmp = Float64(x - Float64(Float64(z * Float64(x - y)) / a)); elseif (t <= 0.000235) tmp = t_2; elseif (t <= 1.45e+50) tmp = Float64(x + Float64(z / Float64(-Float64(a / x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * ((z - t) / (a - t)); t_2 = (z - t) * (y / (a - t)); tmp = 0.0; if (t <= -2.25e+138) tmp = t_1; elseif (t <= -3.3e+66) tmp = z * ((y - x) / (a - t)); elseif (t <= -2.45e-98) tmp = t_2; elseif (t <= 1.42e-67) tmp = x - ((z * (x - y)) / a); elseif (t <= 0.000235) tmp = t_2; elseif (t <= 1.45e+50) tmp = x + (z / -(a / x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.25e+138], t$95$1, If[LessEqual[t, -3.3e+66], N[(z * N[(N[(y - x), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -2.45e-98], t$95$2, If[LessEqual[t, 1.42e-67], N[(x - N[(N[(z * N[(x - y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.000235], t$95$2, If[LessEqual[t, 1.45e+50], N[(x + N[(z / (-N[(a / x), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z - t}{a - t}\\
t_2 := \left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{if}\;t \leq -2.25 \cdot 10^{+138}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -3.3 \cdot 10^{+66}:\\
\;\;\;\;z \cdot \frac{y - x}{a - t}\\
\mathbf{elif}\;t \leq -2.45 \cdot 10^{-98}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 1.42 \cdot 10^{-67}:\\
\;\;\;\;x - \frac{z \cdot \left(x - y\right)}{a}\\
\mathbf{elif}\;t \leq 0.000235:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 1.45 \cdot 10^{+50}:\\
\;\;\;\;x + \frac{z}{-\frac{a}{x}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -2.24999999999999991e138 or 1.45e50 < t Initial program 45.6%
+-commutative45.6%
associate-*l/64.6%
fma-def64.7%
Simplified64.7%
fma-udef64.6%
associate-/r/71.2%
div-inv71.2%
clear-num71.2%
Applied egg-rr71.2%
Taylor expanded in y around inf 67.9%
div-sub67.9%
Simplified67.9%
if -2.24999999999999991e138 < t < -3.3000000000000001e66Initial program 83.7%
associate-*l/83.1%
Simplified83.1%
Taylor expanded in z around inf 61.7%
div-sub61.7%
Simplified61.7%
if -3.3000000000000001e66 < t < -2.45000000000000007e-98 or 1.42000000000000004e-67 < t < 2.34999999999999993e-4Initial program 78.3%
associate-*l/82.6%
Simplified82.6%
Taylor expanded in x around 0 61.4%
associate-/l*61.5%
associate-/r/63.6%
Simplified63.6%
if -2.45000000000000007e-98 < t < 1.42000000000000004e-67Initial program 93.8%
associate-*l/90.5%
Simplified90.5%
Taylor expanded in t around 0 86.3%
if 2.34999999999999993e-4 < t < 1.45e50Initial program 74.9%
associate-*l/93.3%
Simplified93.3%
Taylor expanded in t around 0 61.2%
associate-/l*73.8%
Simplified73.8%
Taylor expanded in y around 0 80.5%
associate-*r/80.5%
neg-mul-180.5%
Simplified80.5%
Final simplification74.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (/ (- x y) (/ t z)))) (t_2 (* (- z t) (/ y (- a t)))))
(if (<= t -3.5e+44)
t_1
(if (<= t -32000.0)
(+ x (/ (* y z) a))
(if (<= t -2.1e-65)
t_2
(if (<= t 8.2e-69)
(- x (/ (* z (- x y)) a))
(if (<= t 4.5e-11)
t_2
(if (<= t 5.5e+49) (+ x (/ z (- (/ a x)))) t_1))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + ((x - y) / (t / z));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -3.5e+44) {
tmp = t_1;
} else if (t <= -32000.0) {
tmp = x + ((y * z) / a);
} else if (t <= -2.1e-65) {
tmp = t_2;
} else if (t <= 8.2e-69) {
tmp = x - ((z * (x - y)) / a);
} else if (t <= 4.5e-11) {
tmp = t_2;
} else if (t <= 5.5e+49) {
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / z))
t_2 = (z - t) * (y / (a - t))
if (t <= (-3.5d+44)) then
tmp = t_1
else if (t <= (-32000.0d0)) then
tmp = x + ((y * z) / a)
else if (t <= (-2.1d-65)) then
tmp = t_2
else if (t <= 8.2d-69) then
tmp = x - ((z * (x - y)) / a)
else if (t <= 4.5d-11) then
tmp = t_2
else if (t <= 5.5d+49) then
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / z));
double t_2 = (z - t) * (y / (a - t));
double tmp;
if (t <= -3.5e+44) {
tmp = t_1;
} else if (t <= -32000.0) {
tmp = x + ((y * z) / a);
} else if (t <= -2.1e-65) {
tmp = t_2;
} else if (t <= 8.2e-69) {
tmp = x - ((z * (x - y)) / a);
} else if (t <= 4.5e-11) {
tmp = t_2;
} else if (t <= 5.5e+49) {
tmp = x + (z / -(a / x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + ((x - y) / (t / z)) t_2 = (z - t) * (y / (a - t)) tmp = 0 if t <= -3.5e+44: tmp = t_1 elif t <= -32000.0: tmp = x + ((y * z) / a) elif t <= -2.1e-65: tmp = t_2 elif t <= 8.2e-69: tmp = x - ((z * (x - y)) / a) elif t <= 4.5e-11: tmp = t_2 elif t <= 5.5e+49: tmp = x + (z / -(a / x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(Float64(x - y) / Float64(t / z))) t_2 = Float64(Float64(z - t) * Float64(y / Float64(a - t))) tmp = 0.0 if (t <= -3.5e+44) tmp = t_1; elseif (t <= -32000.0) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= -2.1e-65) tmp = t_2; elseif (t <= 8.2e-69) tmp = Float64(x - Float64(Float64(z * Float64(x - y)) / a)); elseif (t <= 4.5e-11) tmp = t_2; elseif (t <= 5.5e+49) tmp = Float64(x + Float64(z / Float64(-Float64(a / x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + ((x - y) / (t / z)); t_2 = (z - t) * (y / (a - t)); tmp = 0.0; if (t <= -3.5e+44) tmp = t_1; elseif (t <= -32000.0) tmp = x + ((y * z) / a); elseif (t <= -2.1e-65) tmp = t_2; elseif (t <= 8.2e-69) tmp = x - ((z * (x - y)) / a); elseif (t <= 4.5e-11) tmp = t_2; elseif (t <= 5.5e+49) tmp = x + (z / -(a / x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(N[(x - y), $MachinePrecision] / N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.5e+44], t$95$1, If[LessEqual[t, -32000.0], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -2.1e-65], t$95$2, If[LessEqual[t, 8.2e-69], N[(x - N[(N[(z * N[(x - y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.5e-11], t$95$2, If[LessEqual[t, 5.5e+49], N[(x + N[(z / (-N[(a / x), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \frac{x - y}{\frac{t}{z}}\\
t_2 := \left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{if}\;t \leq -3.5 \cdot 10^{+44}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -32000:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq -2.1 \cdot 10^{-65}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 8.2 \cdot 10^{-69}:\\
\;\;\;\;x - \frac{z \cdot \left(x - y\right)}{a}\\
\mathbf{elif}\;t \leq 4.5 \cdot 10^{-11}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq 5.5 \cdot 10^{+49}:\\
\;\;\;\;x + \frac{z}{-\frac{a}{x}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -3.4999999999999999e44 or 5.50000000000000042e49 < t Initial program 53.8%
associate-*l/69.5%
Simplified69.5%
Taylor expanded in t around inf 68.5%
associate--l+68.5%
associate-*r/68.5%
associate-*r/68.5%
div-sub68.5%
distribute-lft-out--68.5%
associate-*r/68.5%
mul-1-neg68.5%
unsub-neg68.5%
distribute-rgt-out--68.6%
associate-/l*81.9%
Simplified81.9%
Taylor expanded in z around inf 73.8%
if -3.4999999999999999e44 < t < -32000Initial program 56.4%
associate-*l/67.7%
Simplified67.7%
Taylor expanded in t around 0 57.1%
associate-/l*68.0%
Simplified68.0%
Taylor expanded in y around inf 68.2%
if -32000 < t < -2.10000000000000003e-65 or 8.1999999999999998e-69 < t < 4.5e-11Initial program 80.4%
associate-*l/80.3%
Simplified80.3%
Taylor expanded in x around 0 65.0%
associate-/l*65.1%
associate-/r/68.9%
Simplified68.9%
if -2.10000000000000003e-65 < t < 8.1999999999999998e-69Initial program 94.2%
associate-*l/91.1%
Simplified91.1%
Taylor expanded in t around 0 85.2%
if 4.5e-11 < t < 5.50000000000000042e49Initial program 74.9%
associate-*l/93.3%
Simplified93.3%
Taylor expanded in t around 0 61.2%
associate-/l*73.8%
Simplified73.8%
Taylor expanded in y around 0 80.5%
associate-*r/80.5%
neg-mul-180.5%
Simplified80.5%
Final simplification78.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- z a) (/ x t))))
(if (<= t -1.05e+160)
y
(if (<= t -9e+65)
t_1
(if (<= t 4.6e-117)
(+ x (/ z (/ a y)))
(if (<= t 4e+54)
(* x (- 1.0 (/ z a)))
(if (<= t 1.05e+127)
(* z (/ y (- a t)))
(if (<= t 7e+159) t_1 y))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -1.05e+160) {
tmp = y;
} else if (t <= -9e+65) {
tmp = t_1;
} else if (t <= 4.6e-117) {
tmp = x + (z / (a / y));
} else if (t <= 4e+54) {
tmp = x * (1.0 - (z / a));
} else if (t <= 1.05e+127) {
tmp = z * (y / (a - t));
} else if (t <= 7e+159) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = (z - a) * (x / t)
if (t <= (-1.05d+160)) then
tmp = y
else if (t <= (-9d+65)) then
tmp = t_1
else if (t <= 4.6d-117) then
tmp = x + (z / (a / y))
else if (t <= 4d+54) then
tmp = x * (1.0d0 - (z / a))
else if (t <= 1.05d+127) then
tmp = z * (y / (a - t))
else if (t <= 7d+159) then
tmp = t_1
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -1.05e+160) {
tmp = y;
} else if (t <= -9e+65) {
tmp = t_1;
} else if (t <= 4.6e-117) {
tmp = x + (z / (a / y));
} else if (t <= 4e+54) {
tmp = x * (1.0 - (z / a));
} else if (t <= 1.05e+127) {
tmp = z * (y / (a - t));
} else if (t <= 7e+159) {
tmp = t_1;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (z - a) * (x / t) tmp = 0 if t <= -1.05e+160: tmp = y elif t <= -9e+65: tmp = t_1 elif t <= 4.6e-117: tmp = x + (z / (a / y)) elif t <= 4e+54: tmp = x * (1.0 - (z / a)) elif t <= 1.05e+127: tmp = z * (y / (a - t)) elif t <= 7e+159: tmp = t_1 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(z - a) * Float64(x / t)) tmp = 0.0 if (t <= -1.05e+160) tmp = y; elseif (t <= -9e+65) tmp = t_1; elseif (t <= 4.6e-117) tmp = Float64(x + Float64(z / Float64(a / y))); elseif (t <= 4e+54) tmp = Float64(x * Float64(1.0 - Float64(z / a))); elseif (t <= 1.05e+127) tmp = Float64(z * Float64(y / Float64(a - t))); elseif (t <= 7e+159) tmp = t_1; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (z - a) * (x / t); tmp = 0.0; if (t <= -1.05e+160) tmp = y; elseif (t <= -9e+65) tmp = t_1; elseif (t <= 4.6e-117) tmp = x + (z / (a / y)); elseif (t <= 4e+54) tmp = x * (1.0 - (z / a)); elseif (t <= 1.05e+127) tmp = z * (y / (a - t)); elseif (t <= 7e+159) tmp = t_1; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.05e+160], y, If[LessEqual[t, -9e+65], t$95$1, If[LessEqual[t, 4.6e-117], N[(x + N[(z / N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4e+54], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.05e+127], N[(z * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 7e+159], t$95$1, y]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{if}\;t \leq -1.05 \cdot 10^{+160}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -9 \cdot 10^{+65}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 4.6 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{z}{\frac{a}{y}}\\
\mathbf{elif}\;t \leq 4 \cdot 10^{+54}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{elif}\;t \leq 1.05 \cdot 10^{+127}:\\
\;\;\;\;z \cdot \frac{y}{a - t}\\
\mathbf{elif}\;t \leq 7 \cdot 10^{+159}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -1.04999999999999998e160 or 6.9999999999999999e159 < t Initial program 32.0%
associate-*l/56.3%
Simplified56.3%
Taylor expanded in t around inf 57.9%
if -1.04999999999999998e160 < t < -9e65 or 1.04999999999999996e127 < t < 6.9999999999999999e159Initial program 79.3%
associate-*l/79.0%
Simplified79.0%
Taylor expanded in t around inf 67.4%
associate--l+67.4%
associate-*r/67.4%
associate-*r/67.4%
div-sub67.4%
distribute-lft-out--67.4%
associate-*r/67.4%
mul-1-neg67.4%
unsub-neg67.4%
distribute-rgt-out--67.4%
associate-/l*73.2%
Simplified73.2%
Taylor expanded in y around 0 47.6%
associate-*l/53.0%
Simplified53.0%
if -9e65 < t < 4.59999999999999989e-117Initial program 90.2%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in t around 0 76.8%
associate-/l*78.0%
Simplified78.0%
Taylor expanded in y around inf 66.1%
if 4.59999999999999989e-117 < t < 4.0000000000000003e54Initial program 78.4%
associate-*l/83.2%
Simplified83.2%
Taylor expanded in t around 0 49.3%
associate-/l*54.1%
Simplified54.1%
Taylor expanded in x around inf 59.9%
mul-1-neg59.9%
unsub-neg59.9%
Simplified59.9%
if 4.0000000000000003e54 < t < 1.04999999999999996e127Initial program 68.6%
+-commutative68.6%
associate-*l/86.8%
fma-def86.9%
Simplified86.9%
fma-udef86.8%
associate-/r/87.0%
div-inv87.0%
clear-num87.0%
Applied egg-rr87.0%
Taylor expanded in y around inf 80.9%
div-sub80.9%
Simplified80.9%
Taylor expanded in z around inf 29.7%
associate-/l*48.1%
associate-/r/48.0%
Simplified48.0%
Final simplification60.8%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- z a) (/ x t))))
(if (<= t -2.55e+160)
y
(if (<= t -3.4e+64)
t_1
(if (<= t 5.5e-117)
(+ x (/ (* y z) a))
(if (<= t 2.4e+57)
(* x (- 1.0 (/ z a)))
(if (<= t 6.2e+126)
(* z (/ y (- a t)))
(if (<= t 2.45e+160) t_1 y))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -2.55e+160) {
tmp = y;
} else if (t <= -3.4e+64) {
tmp = t_1;
} else if (t <= 5.5e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 2.4e+57) {
tmp = x * (1.0 - (z / a));
} else if (t <= 6.2e+126) {
tmp = z * (y / (a - t));
} else if (t <= 2.45e+160) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = (z - a) * (x / t)
if (t <= (-2.55d+160)) then
tmp = y
else if (t <= (-3.4d+64)) then
tmp = t_1
else if (t <= 5.5d-117) then
tmp = x + ((y * z) / a)
else if (t <= 2.4d+57) then
tmp = x * (1.0d0 - (z / a))
else if (t <= 6.2d+126) then
tmp = z * (y / (a - t))
else if (t <= 2.45d+160) then
tmp = t_1
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -2.55e+160) {
tmp = y;
} else if (t <= -3.4e+64) {
tmp = t_1;
} else if (t <= 5.5e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 2.4e+57) {
tmp = x * (1.0 - (z / a));
} else if (t <= 6.2e+126) {
tmp = z * (y / (a - t));
} else if (t <= 2.45e+160) {
tmp = t_1;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (z - a) * (x / t) tmp = 0 if t <= -2.55e+160: tmp = y elif t <= -3.4e+64: tmp = t_1 elif t <= 5.5e-117: tmp = x + ((y * z) / a) elif t <= 2.4e+57: tmp = x * (1.0 - (z / a)) elif t <= 6.2e+126: tmp = z * (y / (a - t)) elif t <= 2.45e+160: tmp = t_1 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(z - a) * Float64(x / t)) tmp = 0.0 if (t <= -2.55e+160) tmp = y; elseif (t <= -3.4e+64) tmp = t_1; elseif (t <= 5.5e-117) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= 2.4e+57) tmp = Float64(x * Float64(1.0 - Float64(z / a))); elseif (t <= 6.2e+126) tmp = Float64(z * Float64(y / Float64(a - t))); elseif (t <= 2.45e+160) tmp = t_1; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (z - a) * (x / t); tmp = 0.0; if (t <= -2.55e+160) tmp = y; elseif (t <= -3.4e+64) tmp = t_1; elseif (t <= 5.5e-117) tmp = x + ((y * z) / a); elseif (t <= 2.4e+57) tmp = x * (1.0 - (z / a)); elseif (t <= 6.2e+126) tmp = z * (y / (a - t)); elseif (t <= 2.45e+160) tmp = t_1; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.55e+160], y, If[LessEqual[t, -3.4e+64], t$95$1, If[LessEqual[t, 5.5e-117], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.4e+57], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.2e+126], N[(z * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.45e+160], t$95$1, y]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{if}\;t \leq -2.55 \cdot 10^{+160}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -3.4 \cdot 10^{+64}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 5.5 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{+57}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{+126}:\\
\;\;\;\;z \cdot \frac{y}{a - t}\\
\mathbf{elif}\;t \leq 2.45 \cdot 10^{+160}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -2.5500000000000001e160 or 2.4500000000000001e160 < t Initial program 32.0%
associate-*l/56.3%
Simplified56.3%
Taylor expanded in t around inf 57.9%
if -2.5500000000000001e160 < t < -3.4000000000000002e64 or 6.2e126 < t < 2.4500000000000001e160Initial program 79.3%
associate-*l/79.0%
Simplified79.0%
Taylor expanded in t around inf 67.4%
associate--l+67.4%
associate-*r/67.4%
associate-*r/67.4%
div-sub67.4%
distribute-lft-out--67.4%
associate-*r/67.4%
mul-1-neg67.4%
unsub-neg67.4%
distribute-rgt-out--67.4%
associate-/l*73.2%
Simplified73.2%
Taylor expanded in y around 0 47.6%
associate-*l/53.0%
Simplified53.0%
if -3.4000000000000002e64 < t < 5.50000000000000025e-117Initial program 90.2%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in t around 0 76.8%
associate-/l*78.0%
Simplified78.0%
Taylor expanded in y around inf 67.1%
if 5.50000000000000025e-117 < t < 2.40000000000000005e57Initial program 78.4%
associate-*l/83.2%
Simplified83.2%
Taylor expanded in t around 0 49.3%
associate-/l*54.1%
Simplified54.1%
Taylor expanded in x around inf 59.9%
mul-1-neg59.9%
unsub-neg59.9%
Simplified59.9%
if 2.40000000000000005e57 < t < 6.2e126Initial program 68.6%
+-commutative68.6%
associate-*l/86.8%
fma-def86.9%
Simplified86.9%
fma-udef86.8%
associate-/r/87.0%
div-inv87.0%
clear-num87.0%
Applied egg-rr87.0%
Taylor expanded in y around inf 80.9%
div-sub80.9%
Simplified80.9%
Taylor expanded in z around inf 29.7%
associate-/l*48.1%
associate-/r/48.0%
Simplified48.0%
Final simplification61.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (- z a) (/ x t))))
(if (<= t -3.1e+159)
y
(if (<= t -1e+65)
t_1
(if (<= t 7.2e-117)
(+ x (/ (* y z) a))
(if (<= t 5.8e+56)
(* x (- 1.0 (/ z a)))
(if (<= t 1.1e+127)
(/ y (/ (- a t) z))
(if (<= t 6.1e+159) t_1 y))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -3.1e+159) {
tmp = y;
} else if (t <= -1e+65) {
tmp = t_1;
} else if (t <= 7.2e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 5.8e+56) {
tmp = x * (1.0 - (z / a));
} else if (t <= 1.1e+127) {
tmp = y / ((a - t) / z);
} else if (t <= 6.1e+159) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = (z - a) * (x / t)
if (t <= (-3.1d+159)) then
tmp = y
else if (t <= (-1d+65)) then
tmp = t_1
else if (t <= 7.2d-117) then
tmp = x + ((y * z) / a)
else if (t <= 5.8d+56) then
tmp = x * (1.0d0 - (z / a))
else if (t <= 1.1d+127) then
tmp = y / ((a - t) / z)
else if (t <= 6.1d+159) then
tmp = t_1
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (z - a) * (x / t);
double tmp;
if (t <= -3.1e+159) {
tmp = y;
} else if (t <= -1e+65) {
tmp = t_1;
} else if (t <= 7.2e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 5.8e+56) {
tmp = x * (1.0 - (z / a));
} else if (t <= 1.1e+127) {
tmp = y / ((a - t) / z);
} else if (t <= 6.1e+159) {
tmp = t_1;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (z - a) * (x / t) tmp = 0 if t <= -3.1e+159: tmp = y elif t <= -1e+65: tmp = t_1 elif t <= 7.2e-117: tmp = x + ((y * z) / a) elif t <= 5.8e+56: tmp = x * (1.0 - (z / a)) elif t <= 1.1e+127: tmp = y / ((a - t) / z) elif t <= 6.1e+159: tmp = t_1 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(z - a) * Float64(x / t)) tmp = 0.0 if (t <= -3.1e+159) tmp = y; elseif (t <= -1e+65) tmp = t_1; elseif (t <= 7.2e-117) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= 5.8e+56) tmp = Float64(x * Float64(1.0 - Float64(z / a))); elseif (t <= 1.1e+127) tmp = Float64(y / Float64(Float64(a - t) / z)); elseif (t <= 6.1e+159) tmp = t_1; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (z - a) * (x / t); tmp = 0.0; if (t <= -3.1e+159) tmp = y; elseif (t <= -1e+65) tmp = t_1; elseif (t <= 7.2e-117) tmp = x + ((y * z) / a); elseif (t <= 5.8e+56) tmp = x * (1.0 - (z / a)); elseif (t <= 1.1e+127) tmp = y / ((a - t) / z); elseif (t <= 6.1e+159) tmp = t_1; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.1e+159], y, If[LessEqual[t, -1e+65], t$95$1, If[LessEqual[t, 7.2e-117], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5.8e+56], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.1e+127], N[(y / N[(N[(a - t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.1e+159], t$95$1, y]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{if}\;t \leq -3.1 \cdot 10^{+159}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -1 \cdot 10^{+65}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 7.2 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 5.8 \cdot 10^{+56}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{elif}\;t \leq 1.1 \cdot 10^{+127}:\\
\;\;\;\;\frac{y}{\frac{a - t}{z}}\\
\mathbf{elif}\;t \leq 6.1 \cdot 10^{+159}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -3.0999999999999998e159 or 6.1e159 < t Initial program 32.0%
associate-*l/56.3%
Simplified56.3%
Taylor expanded in t around inf 57.9%
if -3.0999999999999998e159 < t < -9.9999999999999999e64 or 1.1000000000000001e127 < t < 6.1e159Initial program 79.3%
associate-*l/79.0%
Simplified79.0%
Taylor expanded in t around inf 67.4%
associate--l+67.4%
associate-*r/67.4%
associate-*r/67.4%
div-sub67.4%
distribute-lft-out--67.4%
associate-*r/67.4%
mul-1-neg67.4%
unsub-neg67.4%
distribute-rgt-out--67.4%
associate-/l*73.2%
Simplified73.2%
Taylor expanded in y around 0 47.6%
associate-*l/53.0%
Simplified53.0%
if -9.9999999999999999e64 < t < 7.2000000000000001e-117Initial program 90.2%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in t around 0 76.8%
associate-/l*78.0%
Simplified78.0%
Taylor expanded in y around inf 67.1%
if 7.2000000000000001e-117 < t < 5.80000000000000014e56Initial program 78.4%
associate-*l/83.2%
Simplified83.2%
Taylor expanded in t around 0 49.3%
associate-/l*54.1%
Simplified54.1%
Taylor expanded in x around inf 59.9%
mul-1-neg59.9%
unsub-neg59.9%
Simplified59.9%
if 5.80000000000000014e56 < t < 1.1000000000000001e127Initial program 68.6%
associate-*l/86.8%
Simplified86.8%
Taylor expanded in z around -inf 29.8%
Taylor expanded in y around inf 29.7%
associate-/l*48.1%
Simplified48.1%
Final simplification61.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- x) (/ a z))))
(if (<= z -0.104)
t_1
(if (<= z -1.45e-155)
y
(if (<= z 1.36e-117)
x
(if (<= z 2.5e-40)
y
(if (<= z 2.85e+57) x (if (<= z 8e+189) (/ y (/ a z)) t_1))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = -x / (a / z);
double tmp;
if (z <= -0.104) {
tmp = t_1;
} else if (z <= -1.45e-155) {
tmp = y;
} else if (z <= 1.36e-117) {
tmp = x;
} else if (z <= 2.5e-40) {
tmp = y;
} else if (z <= 2.85e+57) {
tmp = x;
} else if (z <= 8e+189) {
tmp = y / (a / z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = -x / (a / z)
if (z <= (-0.104d0)) then
tmp = t_1
else if (z <= (-1.45d-155)) then
tmp = y
else if (z <= 1.36d-117) then
tmp = x
else if (z <= 2.5d-40) then
tmp = y
else if (z <= 2.85d+57) then
tmp = x
else if (z <= 8d+189) then
tmp = y / (a / z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = -x / (a / z);
double tmp;
if (z <= -0.104) {
tmp = t_1;
} else if (z <= -1.45e-155) {
tmp = y;
} else if (z <= 1.36e-117) {
tmp = x;
} else if (z <= 2.5e-40) {
tmp = y;
} else if (z <= 2.85e+57) {
tmp = x;
} else if (z <= 8e+189) {
tmp = y / (a / z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = -x / (a / z) tmp = 0 if z <= -0.104: tmp = t_1 elif z <= -1.45e-155: tmp = y elif z <= 1.36e-117: tmp = x elif z <= 2.5e-40: tmp = y elif z <= 2.85e+57: tmp = x elif z <= 8e+189: tmp = y / (a / z) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(-x) / Float64(a / z)) tmp = 0.0 if (z <= -0.104) tmp = t_1; elseif (z <= -1.45e-155) tmp = y; elseif (z <= 1.36e-117) tmp = x; elseif (z <= 2.5e-40) tmp = y; elseif (z <= 2.85e+57) tmp = x; elseif (z <= 8e+189) tmp = Float64(y / Float64(a / z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = -x / (a / z); tmp = 0.0; if (z <= -0.104) tmp = t_1; elseif (z <= -1.45e-155) tmp = y; elseif (z <= 1.36e-117) tmp = x; elseif (z <= 2.5e-40) tmp = y; elseif (z <= 2.85e+57) tmp = x; elseif (z <= 8e+189) tmp = y / (a / z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[((-x) / N[(a / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.104], t$95$1, If[LessEqual[z, -1.45e-155], y, If[LessEqual[z, 1.36e-117], x, If[LessEqual[z, 2.5e-40], y, If[LessEqual[z, 2.85e+57], x, If[LessEqual[z, 8e+189], N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-x}{\frac{a}{z}}\\
\mathbf{if}\;z \leq -0.104:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -1.45 \cdot 10^{-155}:\\
\;\;\;\;y\\
\mathbf{elif}\;z \leq 1.36 \cdot 10^{-117}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 2.5 \cdot 10^{-40}:\\
\;\;\;\;y\\
\mathbf{elif}\;z \leq 2.85 \cdot 10^{+57}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 8 \cdot 10^{+189}:\\
\;\;\;\;\frac{y}{\frac{a}{z}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -0.103999999999999995 or 8.0000000000000002e189 < z Initial program 82.1%
associate-*l/90.3%
Simplified90.3%
Taylor expanded in t around 0 52.1%
associate-/l*56.2%
Simplified56.2%
Taylor expanded in x around inf 45.3%
mul-1-neg45.3%
unsub-neg45.3%
Simplified45.3%
Taylor expanded in z around inf 31.1%
mul-1-neg31.1%
associate-/l*37.5%
distribute-neg-frac37.5%
Simplified37.5%
if -0.103999999999999995 < z < -1.45000000000000005e-155 or 1.35999999999999996e-117 < z < 2.49999999999999982e-40Initial program 60.9%
associate-*l/71.7%
Simplified71.7%
Taylor expanded in t around inf 43.7%
if -1.45000000000000005e-155 < z < 1.35999999999999996e-117 or 2.49999999999999982e-40 < z < 2.8499999999999999e57Initial program 74.9%
associate-*l/74.2%
Simplified74.2%
Taylor expanded in a around inf 48.3%
if 2.8499999999999999e57 < z < 8.0000000000000002e189Initial program 69.6%
associate-*l/87.8%
Simplified87.8%
Taylor expanded in z around -inf 58.1%
Taylor expanded in y around inf 42.3%
*-commutative42.3%
Simplified42.3%
Taylor expanded in a around inf 38.6%
associate-/l*46.1%
Simplified46.1%
Final simplification43.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (- 1.0 (/ z a)))))
(if (<= t -1e+160)
y
(if (<= t -1.4e+73)
(* (- z a) (/ x t))
(if (<= t -900.0)
t_1
(if (<= t -3.2e-149)
(* y (/ (- z t) a))
(if (<= t 8.5e+70) t_1 y)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (t <= -1e+160) {
tmp = y;
} else if (t <= -1.4e+73) {
tmp = (z - a) * (x / t);
} else if (t <= -900.0) {
tmp = t_1;
} else if (t <= -3.2e-149) {
tmp = y * ((z - t) / a);
} else if (t <= 8.5e+70) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = x * (1.0d0 - (z / a))
if (t <= (-1d+160)) then
tmp = y
else if (t <= (-1.4d+73)) then
tmp = (z - a) * (x / t)
else if (t <= (-900.0d0)) then
tmp = t_1
else if (t <= (-3.2d-149)) then
tmp = y * ((z - t) / a)
else if (t <= 8.5d+70) then
tmp = t_1
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (t <= -1e+160) {
tmp = y;
} else if (t <= -1.4e+73) {
tmp = (z - a) * (x / t);
} else if (t <= -900.0) {
tmp = t_1;
} else if (t <= -3.2e-149) {
tmp = y * ((z - t) / a);
} else if (t <= 8.5e+70) {
tmp = t_1;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (1.0 - (z / a)) tmp = 0 if t <= -1e+160: tmp = y elif t <= -1.4e+73: tmp = (z - a) * (x / t) elif t <= -900.0: tmp = t_1 elif t <= -3.2e-149: tmp = y * ((z - t) / a) elif t <= 8.5e+70: tmp = t_1 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(1.0 - Float64(z / a))) tmp = 0.0 if (t <= -1e+160) tmp = y; elseif (t <= -1.4e+73) tmp = Float64(Float64(z - a) * Float64(x / t)); elseif (t <= -900.0) tmp = t_1; elseif (t <= -3.2e-149) tmp = Float64(y * Float64(Float64(z - t) / a)); elseif (t <= 8.5e+70) tmp = t_1; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (1.0 - (z / a)); tmp = 0.0; if (t <= -1e+160) tmp = y; elseif (t <= -1.4e+73) tmp = (z - a) * (x / t); elseif (t <= -900.0) tmp = t_1; elseif (t <= -3.2e-149) tmp = y * ((z - t) / a); elseif (t <= 8.5e+70) tmp = t_1; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1e+160], y, If[LessEqual[t, -1.4e+73], N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -900.0], t$95$1, If[LessEqual[t, -3.2e-149], N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 8.5e+70], t$95$1, y]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{if}\;t \leq -1 \cdot 10^{+160}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -1.4 \cdot 10^{+73}:\\
\;\;\;\;\left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{elif}\;t \leq -900:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -3.2 \cdot 10^{-149}:\\
\;\;\;\;y \cdot \frac{z - t}{a}\\
\mathbf{elif}\;t \leq 8.5 \cdot 10^{+70}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -1.00000000000000001e160 or 8.4999999999999996e70 < t Initial program 43.7%
associate-*l/64.7%
Simplified64.7%
Taylor expanded in t around inf 49.2%
if -1.00000000000000001e160 < t < -1.40000000000000004e73Initial program 72.1%
associate-*l/71.7%
Simplified71.7%
Taylor expanded in t around inf 67.6%
associate--l+67.6%
associate-*r/67.6%
associate-*r/67.6%
div-sub67.6%
distribute-lft-out--67.6%
associate-*r/67.6%
mul-1-neg67.6%
unsub-neg67.6%
distribute-rgt-out--67.6%
associate-/l*72.1%
Simplified72.1%
Taylor expanded in y around 0 49.6%
associate-*l/53.7%
Simplified53.7%
if -1.40000000000000004e73 < t < -900 or -3.20000000000000002e-149 < t < 8.4999999999999996e70Initial program 88.3%
associate-*l/88.7%
Simplified88.7%
Taylor expanded in t around 0 71.8%
associate-/l*73.6%
Simplified73.6%
Taylor expanded in x around inf 61.6%
mul-1-neg61.6%
unsub-neg61.6%
Simplified61.6%
if -900 < t < -3.20000000000000002e-149Initial program 86.1%
+-commutative86.1%
associate-*l/89.4%
fma-def89.7%
Simplified89.7%
fma-udef89.4%
associate-/r/88.3%
div-inv86.0%
clear-num86.1%
Applied egg-rr86.1%
Taylor expanded in y around inf 64.2%
div-sub64.2%
Simplified64.2%
Taylor expanded in a around inf 48.9%
Final simplification56.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* z (/ y a))))
(if (<= t -1.55e+107)
y
(if (<= t -5.8e-22)
x
(if (<= t -9e-273)
t_1
(if (<= t 2.55e-183)
x
(if (<= t 2.1e-152) t_1 (if (<= t 2.7e+63) x y))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = z * (y / a);
double tmp;
if (t <= -1.55e+107) {
tmp = y;
} else if (t <= -5.8e-22) {
tmp = x;
} else if (t <= -9e-273) {
tmp = t_1;
} else if (t <= 2.55e-183) {
tmp = x;
} else if (t <= 2.1e-152) {
tmp = t_1;
} else if (t <= 2.7e+63) {
tmp = x;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = z * (y / a)
if (t <= (-1.55d+107)) then
tmp = y
else if (t <= (-5.8d-22)) then
tmp = x
else if (t <= (-9d-273)) then
tmp = t_1
else if (t <= 2.55d-183) then
tmp = x
else if (t <= 2.1d-152) then
tmp = t_1
else if (t <= 2.7d+63) then
tmp = x
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = z * (y / a);
double tmp;
if (t <= -1.55e+107) {
tmp = y;
} else if (t <= -5.8e-22) {
tmp = x;
} else if (t <= -9e-273) {
tmp = t_1;
} else if (t <= 2.55e-183) {
tmp = x;
} else if (t <= 2.1e-152) {
tmp = t_1;
} else if (t <= 2.7e+63) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = z * (y / a) tmp = 0 if t <= -1.55e+107: tmp = y elif t <= -5.8e-22: tmp = x elif t <= -9e-273: tmp = t_1 elif t <= 2.55e-183: tmp = x elif t <= 2.1e-152: tmp = t_1 elif t <= 2.7e+63: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(z * Float64(y / a)) tmp = 0.0 if (t <= -1.55e+107) tmp = y; elseif (t <= -5.8e-22) tmp = x; elseif (t <= -9e-273) tmp = t_1; elseif (t <= 2.55e-183) tmp = x; elseif (t <= 2.1e-152) tmp = t_1; elseif (t <= 2.7e+63) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = z * (y / a); tmp = 0.0; if (t <= -1.55e+107) tmp = y; elseif (t <= -5.8e-22) tmp = x; elseif (t <= -9e-273) tmp = t_1; elseif (t <= 2.55e-183) tmp = x; elseif (t <= 2.1e-152) tmp = t_1; elseif (t <= 2.7e+63) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.55e+107], y, If[LessEqual[t, -5.8e-22], x, If[LessEqual[t, -9e-273], t$95$1, If[LessEqual[t, 2.55e-183], x, If[LessEqual[t, 2.1e-152], t$95$1, If[LessEqual[t, 2.7e+63], x, y]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \frac{y}{a}\\
\mathbf{if}\;t \leq -1.55 \cdot 10^{+107}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -5.8 \cdot 10^{-22}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq -9 \cdot 10^{-273}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 2.55 \cdot 10^{-183}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 2.1 \cdot 10^{-152}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 2.7 \cdot 10^{+63}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -1.55000000000000013e107 or 2.70000000000000017e63 < t Initial program 47.0%
associate-*l/64.0%
Simplified64.0%
Taylor expanded in t around inf 45.2%
if -1.55000000000000013e107 < t < -5.8000000000000003e-22 or -8.99999999999999921e-273 < t < 2.55e-183 or 2.09999999999999999e-152 < t < 2.70000000000000017e63Initial program 87.6%
associate-*l/88.8%
Simplified88.8%
Taylor expanded in a around inf 40.9%
if -5.8000000000000003e-22 < t < -8.99999999999999921e-273 or 2.55e-183 < t < 2.09999999999999999e-152Initial program 86.9%
associate-*l/88.8%
Simplified88.8%
Taylor expanded in t around 0 70.6%
associate-/l*72.4%
Simplified72.4%
Taylor expanded in z around inf 59.1%
Taylor expanded in y around inf 46.0%
Final simplification43.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (/ (- x y) (/ t (- z a))))))
(if (<= t -1.35e+45)
t_1
(if (<= t 2.4e-69)
(- x (/ (- x y) (/ a (- z t))))
(if (<= t 8.8e-8)
(* (- z t) (/ y (- a t)))
(if (<= t 5e+49) (+ x (/ z (- (/ a x)))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + ((x - y) / (t / (z - a)));
double tmp;
if (t <= -1.35e+45) {
tmp = t_1;
} else if (t <= 2.4e-69) {
tmp = x - ((x - y) / (a / (z - t)));
} else if (t <= 8.8e-8) {
tmp = (z - t) * (y / (a - t));
} else if (t <= 5e+49) {
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / (z - a)))
if (t <= (-1.35d+45)) then
tmp = t_1
else if (t <= 2.4d-69) then
tmp = x - ((x - y) / (a / (z - t)))
else if (t <= 8.8d-8) then
tmp = (z - t) * (y / (a - t))
else if (t <= 5d+49) then
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / (z - a)));
double tmp;
if (t <= -1.35e+45) {
tmp = t_1;
} else if (t <= 2.4e-69) {
tmp = x - ((x - y) / (a / (z - t)));
} else if (t <= 8.8e-8) {
tmp = (z - t) * (y / (a - t));
} else if (t <= 5e+49) {
tmp = x + (z / -(a / x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + ((x - y) / (t / (z - a))) tmp = 0 if t <= -1.35e+45: tmp = t_1 elif t <= 2.4e-69: tmp = x - ((x - y) / (a / (z - t))) elif t <= 8.8e-8: tmp = (z - t) * (y / (a - t)) elif t <= 5e+49: tmp = x + (z / -(a / x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(Float64(x - y) / Float64(t / Float64(z - a)))) tmp = 0.0 if (t <= -1.35e+45) tmp = t_1; elseif (t <= 2.4e-69) tmp = Float64(x - Float64(Float64(x - y) / Float64(a / Float64(z - t)))); elseif (t <= 8.8e-8) tmp = Float64(Float64(z - t) * Float64(y / Float64(a - t))); elseif (t <= 5e+49) tmp = Float64(x + Float64(z / Float64(-Float64(a / x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + ((x - y) / (t / (z - a))); tmp = 0.0; if (t <= -1.35e+45) tmp = t_1; elseif (t <= 2.4e-69) tmp = x - ((x - y) / (a / (z - t))); elseif (t <= 8.8e-8) tmp = (z - t) * (y / (a - t)); elseif (t <= 5e+49) tmp = x + (z / -(a / x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(N[(x - y), $MachinePrecision] / N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.35e+45], t$95$1, If[LessEqual[t, 2.4e-69], N[(x - N[(N[(x - y), $MachinePrecision] / N[(a / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 8.8e-8], N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5e+49], N[(x + N[(z / (-N[(a / x), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \frac{x - y}{\frac{t}{z - a}}\\
\mathbf{if}\;t \leq -1.35 \cdot 10^{+45}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-69}:\\
\;\;\;\;x - \frac{x - y}{\frac{a}{z - t}}\\
\mathbf{elif}\;t \leq 8.8 \cdot 10^{-8}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+49}:\\
\;\;\;\;x + \frac{z}{-\frac{a}{x}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -1.34999999999999992e45 or 5.0000000000000004e49 < t Initial program 53.8%
associate-*l/69.5%
Simplified69.5%
Taylor expanded in t around inf 68.5%
associate--l+68.5%
associate-*r/68.5%
associate-*r/68.5%
div-sub68.5%
distribute-lft-out--68.5%
associate-*r/68.5%
mul-1-neg68.5%
unsub-neg68.5%
distribute-rgt-out--68.6%
associate-/l*81.9%
Simplified81.9%
if -1.34999999999999992e45 < t < 2.4000000000000001e-69Initial program 90.6%
associate-*l/89.6%
Simplified89.6%
Taylor expanded in a around inf 80.5%
associate-/l*84.0%
Simplified84.0%
if 2.4000000000000001e-69 < t < 8.7999999999999994e-8Initial program 73.2%
associate-*l/64.3%
Simplified64.3%
Taylor expanded in x around 0 73.2%
associate-/l*64.8%
associate-/r/73.5%
Simplified73.5%
if 8.7999999999999994e-8 < t < 5.0000000000000004e49Initial program 74.9%
associate-*l/93.3%
Simplified93.3%
Taylor expanded in t around 0 61.2%
associate-/l*73.8%
Simplified73.8%
Taylor expanded in y around 0 80.5%
associate-*r/80.5%
neg-mul-180.5%
Simplified80.5%
Final simplification82.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ y (/ (- x y) (/ t z)))))
(if (<= t -2.4e+45)
t_1
(if (<= t 3.6e-67)
(- x (/ (- x y) (/ a (- z t))))
(if (<= t 0.00022)
(* (- z t) (/ y (- a t)))
(if (<= t 6.5e+52) (+ x (/ z (- (/ a x)))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y + ((x - y) / (t / z));
double tmp;
if (t <= -2.4e+45) {
tmp = t_1;
} else if (t <= 3.6e-67) {
tmp = x - ((x - y) / (a / (z - t)));
} else if (t <= 0.00022) {
tmp = (z - t) * (y / (a - t));
} else if (t <= 6.5e+52) {
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / z))
if (t <= (-2.4d+45)) then
tmp = t_1
else if (t <= 3.6d-67) then
tmp = x - ((x - y) / (a / (z - t)))
else if (t <= 0.00022d0) then
tmp = (z - t) * (y / (a - t))
else if (t <= 6.5d+52) then
tmp = x + (z / -(a / 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 = y + ((x - y) / (t / z));
double tmp;
if (t <= -2.4e+45) {
tmp = t_1;
} else if (t <= 3.6e-67) {
tmp = x - ((x - y) / (a / (z - t)));
} else if (t <= 0.00022) {
tmp = (z - t) * (y / (a - t));
} else if (t <= 6.5e+52) {
tmp = x + (z / -(a / x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y + ((x - y) / (t / z)) tmp = 0 if t <= -2.4e+45: tmp = t_1 elif t <= 3.6e-67: tmp = x - ((x - y) / (a / (z - t))) elif t <= 0.00022: tmp = (z - t) * (y / (a - t)) elif t <= 6.5e+52: tmp = x + (z / -(a / x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y + Float64(Float64(x - y) / Float64(t / z))) tmp = 0.0 if (t <= -2.4e+45) tmp = t_1; elseif (t <= 3.6e-67) tmp = Float64(x - Float64(Float64(x - y) / Float64(a / Float64(z - t)))); elseif (t <= 0.00022) tmp = Float64(Float64(z - t) * Float64(y / Float64(a - t))); elseif (t <= 6.5e+52) tmp = Float64(x + Float64(z / Float64(-Float64(a / x)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y + ((x - y) / (t / z)); tmp = 0.0; if (t <= -2.4e+45) tmp = t_1; elseif (t <= 3.6e-67) tmp = x - ((x - y) / (a / (z - t))); elseif (t <= 0.00022) tmp = (z - t) * (y / (a - t)); elseif (t <= 6.5e+52) tmp = x + (z / -(a / x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y + N[(N[(x - y), $MachinePrecision] / N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.4e+45], t$95$1, If[LessEqual[t, 3.6e-67], N[(x - N[(N[(x - y), $MachinePrecision] / N[(a / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.00022], N[(N[(z - t), $MachinePrecision] * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.5e+52], N[(x + N[(z / (-N[(a / x), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \frac{x - y}{\frac{t}{z}}\\
\mathbf{if}\;t \leq -2.4 \cdot 10^{+45}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq 3.6 \cdot 10^{-67}:\\
\;\;\;\;x - \frac{x - y}{\frac{a}{z - t}}\\
\mathbf{elif}\;t \leq 0.00022:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a - t}\\
\mathbf{elif}\;t \leq 6.5 \cdot 10^{+52}:\\
\;\;\;\;x + \frac{z}{-\frac{a}{x}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -2.39999999999999989e45 or 6.49999999999999996e52 < t Initial program 53.8%
associate-*l/69.5%
Simplified69.5%
Taylor expanded in t around inf 68.5%
associate--l+68.5%
associate-*r/68.5%
associate-*r/68.5%
div-sub68.5%
distribute-lft-out--68.5%
associate-*r/68.5%
mul-1-neg68.5%
unsub-neg68.5%
distribute-rgt-out--68.6%
associate-/l*81.9%
Simplified81.9%
Taylor expanded in z around inf 73.8%
if -2.39999999999999989e45 < t < 3.59999999999999999e-67Initial program 90.6%
associate-*l/89.6%
Simplified89.6%
Taylor expanded in a around inf 80.5%
associate-/l*84.0%
Simplified84.0%
if 3.59999999999999999e-67 < t < 2.20000000000000008e-4Initial program 73.2%
associate-*l/64.3%
Simplified64.3%
Taylor expanded in x around 0 73.2%
associate-/l*64.8%
associate-/r/73.5%
Simplified73.5%
if 2.20000000000000008e-4 < t < 6.49999999999999996e52Initial program 74.9%
associate-*l/93.3%
Simplified93.3%
Taylor expanded in t around 0 61.2%
associate-/l*73.8%
Simplified73.8%
Taylor expanded in y around 0 80.5%
associate-*r/80.5%
neg-mul-180.5%
Simplified80.5%
Final simplification79.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- y) (+ (/ a t) -1.0))))
(if (<= t -3.2e+161)
t_1
(if (<= t -2e+64)
(* (- z a) (/ x t))
(if (<= t 6.1e-117)
(+ x (/ (* y z) a))
(if (<= t 7.2e+62) (* x (- 1.0 (/ z a))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = -y / ((a / t) + -1.0);
double tmp;
if (t <= -3.2e+161) {
tmp = t_1;
} else if (t <= -2e+64) {
tmp = (z - a) * (x / t);
} else if (t <= 6.1e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 7.2e+62) {
tmp = x * (1.0 - (z / 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 = -y / ((a / t) + (-1.0d0))
if (t <= (-3.2d+161)) then
tmp = t_1
else if (t <= (-2d+64)) then
tmp = (z - a) * (x / t)
else if (t <= 6.1d-117) then
tmp = x + ((y * z) / a)
else if (t <= 7.2d+62) then
tmp = x * (1.0d0 - (z / 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 = -y / ((a / t) + -1.0);
double tmp;
if (t <= -3.2e+161) {
tmp = t_1;
} else if (t <= -2e+64) {
tmp = (z - a) * (x / t);
} else if (t <= 6.1e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 7.2e+62) {
tmp = x * (1.0 - (z / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = -y / ((a / t) + -1.0) tmp = 0 if t <= -3.2e+161: tmp = t_1 elif t <= -2e+64: tmp = (z - a) * (x / t) elif t <= 6.1e-117: tmp = x + ((y * z) / a) elif t <= 7.2e+62: tmp = x * (1.0 - (z / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(-y) / Float64(Float64(a / t) + -1.0)) tmp = 0.0 if (t <= -3.2e+161) tmp = t_1; elseif (t <= -2e+64) tmp = Float64(Float64(z - a) * Float64(x / t)); elseif (t <= 6.1e-117) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= 7.2e+62) tmp = Float64(x * Float64(1.0 - Float64(z / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = -y / ((a / t) + -1.0); tmp = 0.0; if (t <= -3.2e+161) tmp = t_1; elseif (t <= -2e+64) tmp = (z - a) * (x / t); elseif (t <= 6.1e-117) tmp = x + ((y * z) / a); elseif (t <= 7.2e+62) tmp = x * (1.0 - (z / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[((-y) / N[(N[(a / t), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.2e+161], t$95$1, If[LessEqual[t, -2e+64], N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.1e-117], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 7.2e+62], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-y}{\frac{a}{t} + -1}\\
\mathbf{if}\;t \leq -3.2 \cdot 10^{+161}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -2 \cdot 10^{+64}:\\
\;\;\;\;\left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{elif}\;t \leq 6.1 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 7.2 \cdot 10^{+62}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -3.20000000000000002e161 or 7.2e62 < t Initial program 43.7%
+-commutative43.7%
associate-*l/64.7%
fma-def64.7%
Simplified64.7%
fma-udef64.7%
associate-/r/72.0%
div-inv72.0%
clear-num72.0%
Applied egg-rr72.0%
Taylor expanded in y around inf 68.4%
div-sub68.4%
Simplified68.4%
Taylor expanded in z around 0 37.4%
mul-1-neg37.4%
*-commutative37.4%
associate-/l*52.9%
div-sub52.9%
sub-neg52.9%
*-inverses52.9%
metadata-eval52.9%
Simplified52.9%
if -3.20000000000000002e161 < t < -2.00000000000000004e64Initial program 76.5%
associate-*l/76.1%
Simplified76.1%
Taylor expanded in t around inf 68.7%
associate--l+68.7%
associate-*r/68.7%
associate-*r/68.7%
div-sub68.7%
distribute-lft-out--68.7%
associate-*r/68.7%
mul-1-neg68.7%
unsub-neg68.7%
distribute-rgt-out--68.7%
associate-/l*72.5%
Simplified72.5%
Taylor expanded in y around 0 50.0%
associate-*l/53.4%
Simplified53.4%
if -2.00000000000000004e64 < t < 6.10000000000000002e-117Initial program 90.2%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in t around 0 76.8%
associate-/l*78.0%
Simplified78.0%
Taylor expanded in y around inf 67.1%
if 6.10000000000000002e-117 < t < 7.2e62Initial program 79.6%
associate-*l/84.1%
Simplified84.1%
Taylor expanded in t around 0 49.3%
associate-/l*54.0%
Simplified54.0%
Taylor expanded in x around inf 59.4%
mul-1-neg59.4%
unsub-neg59.4%
Simplified59.4%
Final simplification60.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (- y) (+ (/ a t) -1.0))))
(if (<= t -8.4e+189)
t_1
(if (<= t -1.2e+45)
(/ (- z) (/ t (- y x)))
(if (<= t 6e-117)
(+ x (/ (* y z) a))
(if (<= t 2.2e+71) (* x (- 1.0 (/ z a))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = -y / ((a / t) + -1.0);
double tmp;
if (t <= -8.4e+189) {
tmp = t_1;
} else if (t <= -1.2e+45) {
tmp = -z / (t / (y - x));
} else if (t <= 6e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 2.2e+71) {
tmp = x * (1.0 - (z / 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 = -y / ((a / t) + (-1.0d0))
if (t <= (-8.4d+189)) then
tmp = t_1
else if (t <= (-1.2d+45)) then
tmp = -z / (t / (y - x))
else if (t <= 6d-117) then
tmp = x + ((y * z) / a)
else if (t <= 2.2d+71) then
tmp = x * (1.0d0 - (z / 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 = -y / ((a / t) + -1.0);
double tmp;
if (t <= -8.4e+189) {
tmp = t_1;
} else if (t <= -1.2e+45) {
tmp = -z / (t / (y - x));
} else if (t <= 6e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 2.2e+71) {
tmp = x * (1.0 - (z / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = -y / ((a / t) + -1.0) tmp = 0 if t <= -8.4e+189: tmp = t_1 elif t <= -1.2e+45: tmp = -z / (t / (y - x)) elif t <= 6e-117: tmp = x + ((y * z) / a) elif t <= 2.2e+71: tmp = x * (1.0 - (z / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(-y) / Float64(Float64(a / t) + -1.0)) tmp = 0.0 if (t <= -8.4e+189) tmp = t_1; elseif (t <= -1.2e+45) tmp = Float64(Float64(-z) / Float64(t / Float64(y - x))); elseif (t <= 6e-117) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= 2.2e+71) tmp = Float64(x * Float64(1.0 - Float64(z / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = -y / ((a / t) + -1.0); tmp = 0.0; if (t <= -8.4e+189) tmp = t_1; elseif (t <= -1.2e+45) tmp = -z / (t / (y - x)); elseif (t <= 6e-117) tmp = x + ((y * z) / a); elseif (t <= 2.2e+71) tmp = x * (1.0 - (z / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[((-y) / N[(N[(a / t), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -8.4e+189], t$95$1, If[LessEqual[t, -1.2e+45], N[((-z) / N[(t / N[(y - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6e-117], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.2e+71], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-y}{\frac{a}{t} + -1}\\
\mathbf{if}\;t \leq -8.4 \cdot 10^{+189}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -1.2 \cdot 10^{+45}:\\
\;\;\;\;\frac{-z}{\frac{t}{y - x}}\\
\mathbf{elif}\;t \leq 6 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 2.2 \cdot 10^{+71}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -8.3999999999999997e189 or 2.19999999999999995e71 < t Initial program 43.9%
+-commutative43.9%
associate-*l/66.8%
fma-def66.8%
Simplified66.8%
fma-udef66.8%
associate-/r/73.8%
div-inv73.8%
clear-num73.8%
Applied egg-rr73.8%
Taylor expanded in y around inf 67.7%
div-sub67.7%
Simplified67.7%
Taylor expanded in z around 0 38.1%
mul-1-neg38.1%
*-commutative38.1%
associate-/l*55.1%
div-sub55.1%
sub-neg55.1%
*-inverses55.1%
metadata-eval55.1%
Simplified55.1%
if -8.3999999999999997e189 < t < -1.19999999999999995e45Initial program 69.7%
associate-*l/72.8%
Simplified72.8%
Taylor expanded in t around inf 67.5%
associate--l+67.5%
associate-*r/67.5%
associate-*r/67.5%
div-sub67.5%
distribute-lft-out--67.5%
associate-*r/67.5%
mul-1-neg67.5%
unsub-neg67.5%
distribute-rgt-out--67.5%
associate-/l*75.3%
Simplified75.3%
Taylor expanded in z around -inf 43.8%
mul-1-neg43.8%
associate-/l*49.0%
Simplified49.0%
if -1.19999999999999995e45 < t < 5.99999999999999982e-117Initial program 90.7%
associate-*l/89.7%
Simplified89.7%
Taylor expanded in t around 0 78.5%
associate-/l*78.9%
Simplified78.9%
Taylor expanded in y around inf 68.6%
if 5.99999999999999982e-117 < t < 2.19999999999999995e71Initial program 79.6%
associate-*l/84.1%
Simplified84.1%
Taylor expanded in t around 0 49.3%
associate-/l*54.0%
Simplified54.0%
Taylor expanded in x around inf 59.4%
mul-1-neg59.4%
unsub-neg59.4%
Simplified59.4%
Final simplification60.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (- 1.0 (/ z a)))))
(if (<= t -2.25e+138)
y
(if (<= t -240.0)
t_1
(if (<= t -5.5e-149) (* y (/ (- z t) a)) (if (<= t 1.6e+71) t_1 y))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (t <= -2.25e+138) {
tmp = y;
} else if (t <= -240.0) {
tmp = t_1;
} else if (t <= -5.5e-149) {
tmp = y * ((z - t) / a);
} else if (t <= 1.6e+71) {
tmp = t_1;
} else {
tmp = 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) :: t_1
real(8) :: tmp
t_1 = x * (1.0d0 - (z / a))
if (t <= (-2.25d+138)) then
tmp = y
else if (t <= (-240.0d0)) then
tmp = t_1
else if (t <= (-5.5d-149)) then
tmp = y * ((z - t) / a)
else if (t <= 1.6d+71) then
tmp = t_1
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (t <= -2.25e+138) {
tmp = y;
} else if (t <= -240.0) {
tmp = t_1;
} else if (t <= -5.5e-149) {
tmp = y * ((z - t) / a);
} else if (t <= 1.6e+71) {
tmp = t_1;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (1.0 - (z / a)) tmp = 0 if t <= -2.25e+138: tmp = y elif t <= -240.0: tmp = t_1 elif t <= -5.5e-149: tmp = y * ((z - t) / a) elif t <= 1.6e+71: tmp = t_1 else: tmp = y return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(1.0 - Float64(z / a))) tmp = 0.0 if (t <= -2.25e+138) tmp = y; elseif (t <= -240.0) tmp = t_1; elseif (t <= -5.5e-149) tmp = Float64(y * Float64(Float64(z - t) / a)); elseif (t <= 1.6e+71) tmp = t_1; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (1.0 - (z / a)); tmp = 0.0; if (t <= -2.25e+138) tmp = y; elseif (t <= -240.0) tmp = t_1; elseif (t <= -5.5e-149) tmp = y * ((z - t) / a); elseif (t <= 1.6e+71) tmp = t_1; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.25e+138], y, If[LessEqual[t, -240.0], t$95$1, If[LessEqual[t, -5.5e-149], N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.6e+71], t$95$1, y]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{if}\;t \leq -2.25 \cdot 10^{+138}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq -240:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -5.5 \cdot 10^{-149}:\\
\;\;\;\;y \cdot \frac{z - t}{a}\\
\mathbf{elif}\;t \leq 1.6 \cdot 10^{+71}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -2.24999999999999991e138 or 1.60000000000000012e71 < t Initial program 44.2%
associate-*l/63.7%
Simplified63.7%
Taylor expanded in t around inf 47.1%
if -2.24999999999999991e138 < t < -240 or -5.50000000000000043e-149 < t < 1.60000000000000012e71Initial program 87.5%
associate-*l/87.8%
Simplified87.8%
Taylor expanded in t around 0 68.7%
associate-/l*69.6%
Simplified69.6%
Taylor expanded in x around inf 58.8%
mul-1-neg58.8%
unsub-neg58.8%
Simplified58.8%
if -240 < t < -5.50000000000000043e-149Initial program 86.1%
+-commutative86.1%
associate-*l/89.4%
fma-def89.7%
Simplified89.7%
fma-udef89.4%
associate-/r/88.3%
div-inv86.0%
clear-num86.1%
Applied egg-rr86.1%
Taylor expanded in y around inf 64.2%
div-sub64.2%
Simplified64.2%
Taylor expanded in a around inf 48.9%
Final simplification54.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* y (/ (- t z) t))))
(if (<= t -2e+161)
t_1
(if (<= t -6.5e+63)
(* (- z a) (/ x t))
(if (<= t 4.5e-117)
(+ x (/ (* y z) a))
(if (<= t 1.55e+50) (* x (- 1.0 (/ z a))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * ((t - z) / t);
double tmp;
if (t <= -2e+161) {
tmp = t_1;
} else if (t <= -6.5e+63) {
tmp = (z - a) * (x / t);
} else if (t <= 4.5e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 1.55e+50) {
tmp = x * (1.0 - (z / 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 = y * ((t - z) / t)
if (t <= (-2d+161)) then
tmp = t_1
else if (t <= (-6.5d+63)) then
tmp = (z - a) * (x / t)
else if (t <= 4.5d-117) then
tmp = x + ((y * z) / a)
else if (t <= 1.55d+50) then
tmp = x * (1.0d0 - (z / 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 = y * ((t - z) / t);
double tmp;
if (t <= -2e+161) {
tmp = t_1;
} else if (t <= -6.5e+63) {
tmp = (z - a) * (x / t);
} else if (t <= 4.5e-117) {
tmp = x + ((y * z) / a);
} else if (t <= 1.55e+50) {
tmp = x * (1.0 - (z / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * ((t - z) / t) tmp = 0 if t <= -2e+161: tmp = t_1 elif t <= -6.5e+63: tmp = (z - a) * (x / t) elif t <= 4.5e-117: tmp = x + ((y * z) / a) elif t <= 1.55e+50: tmp = x * (1.0 - (z / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(Float64(t - z) / t)) tmp = 0.0 if (t <= -2e+161) tmp = t_1; elseif (t <= -6.5e+63) tmp = Float64(Float64(z - a) * Float64(x / t)); elseif (t <= 4.5e-117) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t <= 1.55e+50) tmp = Float64(x * Float64(1.0 - Float64(z / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * ((t - z) / t); tmp = 0.0; if (t <= -2e+161) tmp = t_1; elseif (t <= -6.5e+63) tmp = (z - a) * (x / t); elseif (t <= 4.5e-117) tmp = x + ((y * z) / a); elseif (t <= 1.55e+50) tmp = x * (1.0 - (z / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(t - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2e+161], t$95$1, If[LessEqual[t, -6.5e+63], N[(N[(z - a), $MachinePrecision] * N[(x / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.5e-117], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.55e+50], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{t - z}{t}\\
\mathbf{if}\;t \leq -2 \cdot 10^{+161}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t \leq -6.5 \cdot 10^{+63}:\\
\;\;\;\;\left(z - a\right) \cdot \frac{x}{t}\\
\mathbf{elif}\;t \leq 4.5 \cdot 10^{-117}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t \leq 1.55 \cdot 10^{+50}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if t < -2.0000000000000001e161 or 1.55000000000000001e50 < t Initial program 45.1%
+-commutative45.1%
associate-*l/65.6%
fma-def65.6%
Simplified65.6%
fma-udef65.6%
associate-/r/72.8%
div-inv72.7%
clear-num72.7%
Applied egg-rr72.7%
Taylor expanded in y around inf 69.2%
div-sub69.2%
Simplified69.2%
Taylor expanded in a around 0 64.2%
associate-*r/64.2%
neg-mul-164.2%
Simplified64.2%
if -2.0000000000000001e161 < t < -6.49999999999999992e63Initial program 76.5%
associate-*l/76.1%
Simplified76.1%
Taylor expanded in t around inf 68.7%
associate--l+68.7%
associate-*r/68.7%
associate-*r/68.7%
div-sub68.7%
distribute-lft-out--68.7%
associate-*r/68.7%
mul-1-neg68.7%
unsub-neg68.7%
distribute-rgt-out--68.7%
associate-/l*72.5%
Simplified72.5%
Taylor expanded in y around 0 50.0%
associate-*l/53.4%
Simplified53.4%
if -6.49999999999999992e63 < t < 4.49999999999999969e-117Initial program 90.2%
associate-*l/90.0%
Simplified90.0%
Taylor expanded in t around 0 76.8%
associate-/l*78.0%
Simplified78.0%
Taylor expanded in y around inf 67.1%
if 4.49999999999999969e-117 < t < 1.55000000000000001e50Initial program 78.4%
associate-*l/83.2%
Simplified83.2%
Taylor expanded in t around 0 49.3%
associate-/l*54.1%
Simplified54.1%
Taylor expanded in x around inf 59.9%
mul-1-neg59.9%
unsub-neg59.9%
Simplified59.9%
Final simplification63.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (- 1.0 (/ z a)))))
(if (<= x -1.55e-27)
t_1
(if (<= x 2.5e+38)
(* y (/ (- z t) (- a t)))
(if (<= x 1.25e+79)
(/ (- z) (/ t (- y x)))
(if (<= x 3.45e+115) (* y (/ (- t z) t)) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (x <= -1.55e-27) {
tmp = t_1;
} else if (x <= 2.5e+38) {
tmp = y * ((z - t) / (a - t));
} else if (x <= 1.25e+79) {
tmp = -z / (t / (y - x));
} else if (x <= 3.45e+115) {
tmp = y * ((t - z) / t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x * (1.0d0 - (z / a))
if (x <= (-1.55d-27)) then
tmp = t_1
else if (x <= 2.5d+38) then
tmp = y * ((z - t) / (a - t))
else if (x <= 1.25d+79) then
tmp = -z / (t / (y - x))
else if (x <= 3.45d+115) then
tmp = y * ((t - z) / t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (x <= -1.55e-27) {
tmp = t_1;
} else if (x <= 2.5e+38) {
tmp = y * ((z - t) / (a - t));
} else if (x <= 1.25e+79) {
tmp = -z / (t / (y - x));
} else if (x <= 3.45e+115) {
tmp = y * ((t - z) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (1.0 - (z / a)) tmp = 0 if x <= -1.55e-27: tmp = t_1 elif x <= 2.5e+38: tmp = y * ((z - t) / (a - t)) elif x <= 1.25e+79: tmp = -z / (t / (y - x)) elif x <= 3.45e+115: tmp = y * ((t - z) / t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(1.0 - Float64(z / a))) tmp = 0.0 if (x <= -1.55e-27) tmp = t_1; elseif (x <= 2.5e+38) tmp = Float64(y * Float64(Float64(z - t) / Float64(a - t))); elseif (x <= 1.25e+79) tmp = Float64(Float64(-z) / Float64(t / Float64(y - x))); elseif (x <= 3.45e+115) tmp = Float64(y * Float64(Float64(t - z) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (1.0 - (z / a)); tmp = 0.0; if (x <= -1.55e-27) tmp = t_1; elseif (x <= 2.5e+38) tmp = y * ((z - t) / (a - t)); elseif (x <= 1.25e+79) tmp = -z / (t / (y - x)); elseif (x <= 3.45e+115) tmp = y * ((t - z) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.55e-27], t$95$1, If[LessEqual[x, 2.5e+38], N[(y * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.25e+79], N[((-z) / N[(t / N[(y - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.45e+115], N[(y * N[(N[(t - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{if}\;x \leq -1.55 \cdot 10^{-27}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq 2.5 \cdot 10^{+38}:\\
\;\;\;\;y \cdot \frac{z - t}{a - t}\\
\mathbf{elif}\;x \leq 1.25 \cdot 10^{+79}:\\
\;\;\;\;\frac{-z}{\frac{t}{y - x}}\\
\mathbf{elif}\;x \leq 3.45 \cdot 10^{+115}:\\
\;\;\;\;y \cdot \frac{t - z}{t}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if x < -1.5499999999999999e-27 or 3.44999999999999983e115 < x Initial program 66.2%
associate-*l/77.6%
Simplified77.6%
Taylor expanded in t around 0 56.5%
associate-/l*61.8%
Simplified61.8%
Taylor expanded in x around inf 62.8%
mul-1-neg62.8%
unsub-neg62.8%
Simplified62.8%
if -1.5499999999999999e-27 < x < 2.49999999999999985e38Initial program 84.1%
+-commutative84.1%
associate-*l/86.4%
fma-def86.4%
Simplified86.4%
fma-udef86.4%
associate-/r/92.5%
div-inv91.9%
clear-num92.0%
Applied egg-rr92.0%
Taylor expanded in y around inf 74.9%
div-sub74.9%
Simplified74.9%
if 2.49999999999999985e38 < x < 1.25e79Initial program 42.7%
associate-*l/60.8%
Simplified60.8%
Taylor expanded in t around inf 62.0%
associate--l+62.0%
associate-*r/62.0%
associate-*r/62.0%
div-sub62.0%
distribute-lft-out--62.0%
associate-*r/62.0%
mul-1-neg62.0%
unsub-neg62.0%
distribute-rgt-out--62.0%
associate-/l*80.5%
Simplified80.5%
Taylor expanded in z around -inf 43.6%
mul-1-neg43.6%
associate-/l*62.1%
Simplified62.1%
if 1.25e79 < x < 3.44999999999999983e115Initial program 23.9%
+-commutative23.9%
associate-*l/22.2%
fma-def21.5%
Simplified21.5%
fma-udef22.2%
associate-/r/41.7%
div-inv41.7%
clear-num41.7%
Applied egg-rr41.7%
Taylor expanded in y around inf 81.6%
div-sub81.6%
Simplified81.6%
Taylor expanded in a around 0 81.6%
associate-*r/81.6%
neg-mul-181.6%
Simplified81.6%
Final simplification69.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* x (- 1.0 (/ z a)))))
(if (<= x -1.3e-27)
t_1
(if (<= x 1.1e+36)
(* y (/ (- z t) (- a t)))
(if (<= x 1.26e+81)
(* z (/ (- y x) (- a t)))
(if (<= x 2.1e+115) (* y (/ (- t z) t)) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (x <= -1.3e-27) {
tmp = t_1;
} else if (x <= 1.1e+36) {
tmp = y * ((z - t) / (a - t));
} else if (x <= 1.26e+81) {
tmp = z * ((y - x) / (a - t));
} else if (x <= 2.1e+115) {
tmp = y * ((t - z) / t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = x * (1.0d0 - (z / a))
if (x <= (-1.3d-27)) then
tmp = t_1
else if (x <= 1.1d+36) then
tmp = y * ((z - t) / (a - t))
else if (x <= 1.26d+81) then
tmp = z * ((y - x) / (a - t))
else if (x <= 2.1d+115) then
tmp = y * ((t - z) / t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x * (1.0 - (z / a));
double tmp;
if (x <= -1.3e-27) {
tmp = t_1;
} else if (x <= 1.1e+36) {
tmp = y * ((z - t) / (a - t));
} else if (x <= 1.26e+81) {
tmp = z * ((y - x) / (a - t));
} else if (x <= 2.1e+115) {
tmp = y * ((t - z) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x * (1.0 - (z / a)) tmp = 0 if x <= -1.3e-27: tmp = t_1 elif x <= 1.1e+36: tmp = y * ((z - t) / (a - t)) elif x <= 1.26e+81: tmp = z * ((y - x) / (a - t)) elif x <= 2.1e+115: tmp = y * ((t - z) / t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x * Float64(1.0 - Float64(z / a))) tmp = 0.0 if (x <= -1.3e-27) tmp = t_1; elseif (x <= 1.1e+36) tmp = Float64(y * Float64(Float64(z - t) / Float64(a - t))); elseif (x <= 1.26e+81) tmp = Float64(z * Float64(Float64(y - x) / Float64(a - t))); elseif (x <= 2.1e+115) tmp = Float64(y * Float64(Float64(t - z) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x * (1.0 - (z / a)); tmp = 0.0; if (x <= -1.3e-27) tmp = t_1; elseif (x <= 1.1e+36) tmp = y * ((z - t) / (a - t)); elseif (x <= 1.26e+81) tmp = z * ((y - x) / (a - t)); elseif (x <= 2.1e+115) tmp = y * ((t - z) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.3e-27], t$95$1, If[LessEqual[x, 1.1e+36], N[(y * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.26e+81], N[(z * N[(N[(y - x), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2.1e+115], N[(y * N[(N[(t - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{if}\;x \leq -1.3 \cdot 10^{-27}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq 1.1 \cdot 10^{+36}:\\
\;\;\;\;y \cdot \frac{z - t}{a - t}\\
\mathbf{elif}\;x \leq 1.26 \cdot 10^{+81}:\\
\;\;\;\;z \cdot \frac{y - x}{a - t}\\
\mathbf{elif}\;x \leq 2.1 \cdot 10^{+115}:\\
\;\;\;\;y \cdot \frac{t - z}{t}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if x < -1.30000000000000009e-27 or 2.10000000000000003e115 < x Initial program 66.2%
associate-*l/77.6%
Simplified77.6%
Taylor expanded in t around 0 56.5%
associate-/l*61.8%
Simplified61.8%
Taylor expanded in x around inf 62.8%
mul-1-neg62.8%
unsub-neg62.8%
Simplified62.8%
if -1.30000000000000009e-27 < x < 1.1e36Initial program 83.9%
+-commutative83.9%
associate-*l/86.2%
fma-def86.1%
Simplified86.1%
fma-udef86.2%
associate-/r/92.4%
div-inv91.8%
clear-num91.9%
Applied egg-rr91.9%
Taylor expanded in y around inf 74.5%
div-sub74.5%
Simplified74.5%
if 1.1e36 < x < 1.25999999999999996e81Initial program 52.3%
associate-*l/67.4%
Simplified67.4%
Taylor expanded in z around inf 75.5%
div-sub75.5%
Simplified75.5%
if 1.25999999999999996e81 < x < 2.10000000000000003e115Initial program 23.9%
+-commutative23.9%
associate-*l/22.2%
fma-def21.5%
Simplified21.5%
fma-udef22.2%
associate-/r/41.7%
div-inv41.7%
clear-num41.7%
Applied egg-rr41.7%
Taylor expanded in y around inf 81.6%
div-sub81.6%
Simplified81.6%
Taylor expanded in a around 0 81.6%
associate-*r/81.6%
neg-mul-181.6%
Simplified81.6%
Final simplification69.8%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -2.55e-196) (not (<= a 2.45e-127))) (- x (* (- z t) (/ (- x y) (- a t)))) (+ y (/ (- x y) (/ t (- z a))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -2.55e-196) || !(a <= 2.45e-127)) {
tmp = x - ((z - t) * ((x - y) / (a - t)));
} else {
tmp = y + ((x - y) / (t / (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 ((a <= (-2.55d-196)) .or. (.not. (a <= 2.45d-127))) then
tmp = x - ((z - t) * ((x - y) / (a - t)))
else
tmp = y + ((x - y) / (t / (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 ((a <= -2.55e-196) || !(a <= 2.45e-127)) {
tmp = x - ((z - t) * ((x - y) / (a - t)));
} else {
tmp = y + ((x - y) / (t / (z - a)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -2.55e-196) or not (a <= 2.45e-127): tmp = x - ((z - t) * ((x - y) / (a - t))) else: tmp = y + ((x - y) / (t / (z - a))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -2.55e-196) || !(a <= 2.45e-127)) tmp = Float64(x - Float64(Float64(z - t) * Float64(Float64(x - y) / Float64(a - t)))); else tmp = Float64(y + Float64(Float64(x - y) / Float64(t / Float64(z - a)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -2.55e-196) || ~((a <= 2.45e-127))) tmp = x - ((z - t) * ((x - y) / (a - t))); else tmp = y + ((x - y) / (t / (z - a))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -2.55e-196], N[Not[LessEqual[a, 2.45e-127]], $MachinePrecision]], N[(x - N[(N[(z - t), $MachinePrecision] * N[(N[(x - y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + N[(N[(x - y), $MachinePrecision] / N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.55 \cdot 10^{-196} \lor \neg \left(a \leq 2.45 \cdot 10^{-127}\right):\\
\;\;\;\;x - \left(z - t\right) \cdot \frac{x - y}{a - t}\\
\mathbf{else}:\\
\;\;\;\;y + \frac{x - y}{\frac{t}{z - a}}\\
\end{array}
\end{array}
if a < -2.55000000000000001e-196 or 2.45e-127 < a Initial program 77.1%
associate-*l/85.7%
Simplified85.7%
if -2.55000000000000001e-196 < a < 2.45e-127Initial program 62.6%
associate-*l/62.3%
Simplified62.3%
Taylor expanded in t around inf 88.4%
associate--l+88.4%
associate-*r/88.4%
associate-*r/88.4%
div-sub88.4%
distribute-lft-out--88.4%
associate-*r/88.4%
mul-1-neg88.4%
unsub-neg88.4%
distribute-rgt-out--88.4%
associate-/l*93.3%
Simplified93.3%
Final simplification87.4%
(FPCore (x y z t a) :precision binary64 (if (<= t -2.3e+138) y (if (<= t 2.6e+71) (* x (- 1.0 (/ z a))) y)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.3e+138) {
tmp = y;
} else if (t <= 2.6e+71) {
tmp = x * (1.0 - (z / a));
} else {
tmp = 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 <= (-2.3d+138)) then
tmp = y
else if (t <= 2.6d+71) then
tmp = x * (1.0d0 - (z / a))
else
tmp = 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 <= -2.3e+138) {
tmp = y;
} else if (t <= 2.6e+71) {
tmp = x * (1.0 - (z / a));
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.3e+138: tmp = y elif t <= 2.6e+71: tmp = x * (1.0 - (z / a)) else: tmp = y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.3e+138) tmp = y; elseif (t <= 2.6e+71) tmp = Float64(x * Float64(1.0 - Float64(z / a))); else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -2.3e+138) tmp = y; elseif (t <= 2.6e+71) tmp = x * (1.0 - (z / a)); else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.3e+138], y, If[LessEqual[t, 2.6e+71], N[(x * N[(1.0 - N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], y]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.3 \cdot 10^{+138}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq 2.6 \cdot 10^{+71}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{a}\right)\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -2.30000000000000008e138 or 2.59999999999999991e71 < t Initial program 44.2%
associate-*l/63.7%
Simplified63.7%
Taylor expanded in t around inf 47.1%
if -2.30000000000000008e138 < t < 2.59999999999999991e71Initial program 87.3%
associate-*l/88.1%
Simplified88.1%
Taylor expanded in t around 0 67.2%
associate-/l*68.5%
Simplified68.5%
Taylor expanded in x around inf 54.3%
mul-1-neg54.3%
unsub-neg54.3%
Simplified54.3%
Final simplification52.1%
(FPCore (x y z t a) :precision binary64 (if (<= t -2.8e+110) y (if (<= t 5.5e+66) x y)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.8e+110) {
tmp = y;
} else if (t <= 5.5e+66) {
tmp = x;
} else {
tmp = 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 <= (-2.8d+110)) then
tmp = y
else if (t <= 5.5d+66) then
tmp = x
else
tmp = 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 <= -2.8e+110) {
tmp = y;
} else if (t <= 5.5e+66) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.8e+110: tmp = y elif t <= 5.5e+66: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.8e+110) tmp = y; elseif (t <= 5.5e+66) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -2.8e+110) tmp = y; elseif (t <= 5.5e+66) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.8e+110], y, If[LessEqual[t, 5.5e+66], x, y]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.8 \cdot 10^{+110}:\\
\;\;\;\;y\\
\mathbf{elif}\;t \leq 5.5 \cdot 10^{+66}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if t < -2.79999999999999987e110 or 5.5e66 < t Initial program 47.0%
associate-*l/64.0%
Simplified64.0%
Taylor expanded in t around inf 45.2%
if -2.79999999999999987e110 < t < 5.5e66Initial program 87.4%
associate-*l/88.8%
Simplified88.8%
Taylor expanded in a around inf 34.9%
Final simplification38.4%
(FPCore (x y z t a) :precision binary64 x)
double code(double x, double y, double z, double t, double a) {
return x;
}
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
end function
public static double code(double x, double y, double z, double t, double a) {
return x;
}
def code(x, y, z, t, a): return x
function code(x, y, z, t, a) return x end
function tmp = code(x, y, z, t, a) tmp = x; end
code[x_, y_, z_, t_, a_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 73.8%
associate-*l/80.5%
Simplified80.5%
Taylor expanded in a around inf 26.1%
Final simplification26.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t))))))
(if (< a -1.6153062845442575e-142)
t_1
(if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((y - x) / 1.0) * ((z - t) / (a - t)));
double tmp;
if (a < -1.6153062845442575e-142) {
tmp = t_1;
} else if (a < 3.774403170083174e-182) {
tmp = y - ((z / t) * (y - 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 = x + (((y - x) / 1.0d0) * ((z - t) / (a - t)))
if (a < (-1.6153062845442575d-142)) then
tmp = t_1
else if (a < 3.774403170083174d-182) then
tmp = y - ((z / t) * (y - 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 - x) / 1.0) * ((z - t) / (a - t)));
double tmp;
if (a < -1.6153062845442575e-142) {
tmp = t_1;
} else if (a < 3.774403170083174e-182) {
tmp = y - ((z / t) * (y - x));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (((y - x) / 1.0) * ((z - t) / (a - t))) tmp = 0 if a < -1.6153062845442575e-142: tmp = t_1 elif a < 3.774403170083174e-182: tmp = y - ((z / t) * (y - x)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(Float64(y - x) / 1.0) * Float64(Float64(z - t) / Float64(a - t)))) tmp = 0.0 if (a < -1.6153062845442575e-142) tmp = t_1; elseif (a < 3.774403170083174e-182) tmp = Float64(y - Float64(Float64(z / t) * Float64(y - x))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (((y - x) / 1.0) * ((z - t) / (a - t))); tmp = 0.0; if (a < -1.6153062845442575e-142) tmp = t_1; elseif (a < 3.774403170083174e-182) tmp = y - ((z / t) * (y - x)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(N[(y - x), $MachinePrecision] / 1.0), $MachinePrecision] * N[(N[(z - t), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[a, -1.6153062845442575e-142], t$95$1, If[Less[a, 3.774403170083174e-182], N[(y - N[(N[(z / t), $MachinePrecision] * N[(y - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y - x}{1} \cdot \frac{z - t}{a - t}\\
\mathbf{if}\;a < -1.6153062845442575 \cdot 10^{-142}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;a < 3.774403170083174 \cdot 10^{-182}:\\
\;\;\;\;y - \frac{z}{t} \cdot \left(y - x\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
herbie shell --seed 2024023
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:linMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< a -1.6153062845442575e-142) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t)))) (if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t))))))
(+ x (/ (* (- y x) (- z t)) (- a t))))