\[\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}
\]
↓
\[\begin{array}{l}
t_0 := \frac{b}{d} + c \cdot \frac{\frac{a}{d}}{d}\\
t_1 := \frac{\frac{\mathsf{fma}\left(a, c, d \cdot b\right)}{\mathsf{hypot}\left(c, d\right)}}{\mathsf{hypot}\left(c, d\right)}\\
\mathbf{if}\;c \leq -4.490116382232771 \cdot 10^{+114}:\\
\;\;\;\;\frac{a}{c} + \frac{d}{c} \cdot \frac{b}{c}\\
\mathbf{elif}\;c \leq -5.612379604288589 \cdot 10^{+70}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;c \leq -4.567322487019083 \cdot 10^{+29}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;c \leq -261.0441113546419:\\
\;\;\;\;t_0\\
\mathbf{elif}\;c \leq -1 \cdot 10^{-120}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;c \leq 10^{-152}:\\
\;\;\;\;\frac{b}{d} + \frac{\frac{c \cdot \left(a - c \cdot \frac{b}{d}\right)}{d}}{d}\\
\mathbf{elif}\;c \leq 8.79419289432676 \cdot 10^{+153}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{c} + \frac{\frac{d}{c}}{\frac{c}{b - \frac{a}{\frac{c}{d}}}}\\
\end{array}
\]
double code(double a, double b, double c, double d) {
return ((a * c) + (b * d)) / ((c * c) + (d * d));
}
↓
double code(double a, double b, double c, double d) {
double t_0 = (b / d) + (c * ((a / d) / d));
double t_1 = (fma(a, c, (d * b)) / hypot(c, d)) / hypot(c, d);
double tmp;
if (c <= -4.490116382232771e+114) {
tmp = (a / c) + ((d / c) * (b / c));
} else if (c <= -5.612379604288589e+70) {
tmp = t_0;
} else if (c <= -4.567322487019083e+29) {
tmp = t_1;
} else if (c <= -261.0441113546419) {
tmp = t_0;
} else if (c <= -1e-120) {
tmp = t_1;
} else if (c <= 1e-152) {
tmp = (b / d) + (((c * (a - (c * (b / d)))) / d) / d);
} else if (c <= 8.79419289432676e+153) {
tmp = t_1;
} else {
tmp = (a / c) + ((d / c) / (c / (b - (a / (c / d)))));
}
return tmp;
}
function code(a, b, c, d)
return Float64(Float64(Float64(a * c) + Float64(b * d)) / Float64(Float64(c * c) + Float64(d * d)))
end
↓
function code(a, b, c, d)
t_0 = Float64(Float64(b / d) + Float64(c * Float64(Float64(a / d) / d)))
t_1 = Float64(Float64(fma(a, c, Float64(d * b)) / hypot(c, d)) / hypot(c, d))
tmp = 0.0
if (c <= -4.490116382232771e+114)
tmp = Float64(Float64(a / c) + Float64(Float64(d / c) * Float64(b / c)));
elseif (c <= -5.612379604288589e+70)
tmp = t_0;
elseif (c <= -4.567322487019083e+29)
tmp = t_1;
elseif (c <= -261.0441113546419)
tmp = t_0;
elseif (c <= -1e-120)
tmp = t_1;
elseif (c <= 1e-152)
tmp = Float64(Float64(b / d) + Float64(Float64(Float64(c * Float64(a - Float64(c * Float64(b / d)))) / d) / d));
elseif (c <= 8.79419289432676e+153)
tmp = t_1;
else
tmp = Float64(Float64(a / c) + Float64(Float64(d / c) / Float64(c / Float64(b - Float64(a / Float64(c / d))))));
end
return tmp
end
code[a_, b_, c_, d_] := N[(N[(N[(a * c), $MachinePrecision] + N[(b * d), $MachinePrecision]), $MachinePrecision] / N[(N[(c * c), $MachinePrecision] + N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[a_, b_, c_, d_] := Block[{t$95$0 = N[(N[(b / d), $MachinePrecision] + N[(c * N[(N[(a / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(a * c + N[(d * b), $MachinePrecision]), $MachinePrecision] / N[Sqrt[c ^ 2 + d ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[c ^ 2 + d ^ 2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -4.490116382232771e+114], N[(N[(a / c), $MachinePrecision] + N[(N[(d / c), $MachinePrecision] * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, -5.612379604288589e+70], t$95$0, If[LessEqual[c, -4.567322487019083e+29], t$95$1, If[LessEqual[c, -261.0441113546419], t$95$0, If[LessEqual[c, -1e-120], t$95$1, If[LessEqual[c, 1e-152], N[(N[(b / d), $MachinePrecision] + N[(N[(N[(c * N[(a - N[(c * N[(b / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 8.79419289432676e+153], t$95$1, N[(N[(a / c), $MachinePrecision] + N[(N[(d / c), $MachinePrecision] / N[(c / N[(b - N[(a / N[(c / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]
\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}
↓
\begin{array}{l}
t_0 := \frac{b}{d} + c \cdot \frac{\frac{a}{d}}{d}\\
t_1 := \frac{\frac{\mathsf{fma}\left(a, c, d \cdot b\right)}{\mathsf{hypot}\left(c, d\right)}}{\mathsf{hypot}\left(c, d\right)}\\
\mathbf{if}\;c \leq -4.490116382232771 \cdot 10^{+114}:\\
\;\;\;\;\frac{a}{c} + \frac{d}{c} \cdot \frac{b}{c}\\
\mathbf{elif}\;c \leq -5.612379604288589 \cdot 10^{+70}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;c \leq -4.567322487019083 \cdot 10^{+29}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;c \leq -261.0441113546419:\\
\;\;\;\;t_0\\
\mathbf{elif}\;c \leq -1 \cdot 10^{-120}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;c \leq 10^{-152}:\\
\;\;\;\;\frac{b}{d} + \frac{\frac{c \cdot \left(a - c \cdot \frac{b}{d}\right)}{d}}{d}\\
\mathbf{elif}\;c \leq 8.79419289432676 \cdot 10^{+153}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{c} + \frac{\frac{d}{c}}{\frac{c}{b - \frac{a}{\frac{c}{d}}}}\\
\end{array}
