\[x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\]
↓
\[\begin{array}{l}
t_1 := \mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
t_2 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{if}\;t_2 \leq -1 \cdot 10^{-296}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-262}:\\
\;\;\;\;t + \frac{a - y}{\frac{-z}{x}}\\
\mathbf{elif}\;t_2 \leq 4 \cdot 10^{-147}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-144}:\\
\;\;\;\;t\\
\mathbf{elif}\;t_2 \leq 2 \cdot 10^{+199}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\]
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
↓
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (- t x) (/ (- y z) (- a z)) x))
(t_2 (+ x (* (- y z) (/ (- t x) (- a z))))))
(if (<= t_2 -1e-296)
t_1
(if (<= t_2 5e-262)
(+ t (/ (- a y) (/ (- z) x)))
(if (<= t_2 4e-147)
t_1
(if (<= t_2 5e-144) t (if (<= t_2 2e+199) t_2 t_1)))))))double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
↓
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((t - x), ((y - z) / (a - z)), x);
double t_2 = x + ((y - z) * ((t - x) / (a - z)));
double tmp;
if (t_2 <= -1e-296) {
tmp = t_1;
} else if (t_2 <= 5e-262) {
tmp = t + ((a - y) / (-z / x));
} else if (t_2 <= 4e-147) {
tmp = t_1;
} else if (t_2 <= 5e-144) {
tmp = t;
} else if (t_2 <= 2e+199) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a)
return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z))))
end
↓
function code(x, y, z, t, a)
t_1 = fma(Float64(t - x), Float64(Float64(y - z) / Float64(a - z)), x)
t_2 = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z))))
tmp = 0.0
if (t_2 <= -1e-296)
tmp = t_1;
elseif (t_2 <= 5e-262)
tmp = Float64(t + Float64(Float64(a - y) / Float64(Float64(-z) / x)));
elseif (t_2 <= 4e-147)
tmp = t_1;
elseif (t_2 <= 5e-144)
tmp = t;
elseif (t_2 <= 2e+199)
tmp = t_2;
else
tmp = t_1;
end
return tmp
end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-296], t$95$1, If[LessEqual[t$95$2, 5e-262], N[(t + N[(N[(a - y), $MachinePrecision] / N[((-z) / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 4e-147], t$95$1, If[LessEqual[t$95$2, 5e-144], t, If[LessEqual[t$95$2, 2e+199], t$95$2, t$95$1]]]]]]]
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
↓
\begin{array}{l}
t_1 := \mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
t_2 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{if}\;t_2 \leq -1 \cdot 10^{-296}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-262}:\\
\;\;\;\;t + \frac{a - y}{\frac{-z}{x}}\\
\mathbf{elif}\;t_2 \leq 4 \cdot 10^{-147}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-144}:\\
\;\;\;\;t\\
\mathbf{elif}\;t_2 \leq 2 \cdot 10^{+199}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
