\[\frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\]
↓
\[\begin{array}{l}
t_0 := \frac{\frac{\mathsf{fma}\left(x.re, y.re, y.im \cdot x.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{if}\;y.im \leq -3.4 \cdot 10^{+103}:\\
\;\;\;\;\left(x.im + \frac{x.re}{\frac{y.im}{y.re}}\right) \cdot \frac{-1}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{elif}\;y.im \leq -1.02 \cdot 10^{-196}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 2.15 \cdot 10^{-181}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{elif}\;y.im \leq 1.1 \cdot 10^{+96}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.im + y.re \cdot \frac{x.re}{y.im}\right)\\
\end{array}
\]
(FPCore (x.re x.im y.re y.im)
:precision binary64
(/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))))
↓
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(/
(/ (fma x.re y.re (* y.im x.im)) (hypot y.re y.im))
(hypot y.re y.im))))
(if (<= y.im -3.4e+103)
(* (+ x.im (/ x.re (/ y.im y.re))) (/ -1.0 (hypot y.re y.im)))
(if (<= y.im -1.02e-196)
t_0
(if (<= y.im 2.15e-181)
(+ (/ x.re y.re) (* x.im (* (/ y.im y.re) (/ 1.0 y.re))))
(if (<= y.im 1.1e+96)
t_0
(* (/ 1.0 (hypot y.re y.im)) (+ x.im (* y.re (/ x.re y.im))))))))))double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((x_46_re * y_46_re) + (x_46_im * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
↓
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (fma(x_46_re, y_46_re, (y_46_im * x_46_im)) / hypot(y_46_re, y_46_im)) / hypot(y_46_re, y_46_im);
double tmp;
if (y_46_im <= -3.4e+103) {
tmp = (x_46_im + (x_46_re / (y_46_im / y_46_re))) * (-1.0 / hypot(y_46_re, y_46_im));
} else if (y_46_im <= -1.02e-196) {
tmp = t_0;
} else if (y_46_im <= 2.15e-181) {
tmp = (x_46_re / y_46_re) + (x_46_im * ((y_46_im / y_46_re) * (1.0 / y_46_re)));
} else if (y_46_im <= 1.1e+96) {
tmp = t_0;
} else {
tmp = (1.0 / hypot(y_46_re, y_46_im)) * (x_46_im + (y_46_re * (x_46_re / y_46_im)));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im)
return Float64(Float64(Float64(x_46_re * y_46_re) + Float64(x_46_im * y_46_im)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im)))
end
↓
function code(x_46_re, x_46_im, y_46_re, y_46_im)
t_0 = Float64(Float64(fma(x_46_re, y_46_re, Float64(y_46_im * x_46_im)) / hypot(y_46_re, y_46_im)) / hypot(y_46_re, y_46_im))
tmp = 0.0
if (y_46_im <= -3.4e+103)
tmp = Float64(Float64(x_46_im + Float64(x_46_re / Float64(y_46_im / y_46_re))) * Float64(-1.0 / hypot(y_46_re, y_46_im)));
elseif (y_46_im <= -1.02e-196)
tmp = t_0;
elseif (y_46_im <= 2.15e-181)
tmp = Float64(Float64(x_46_re / y_46_re) + Float64(x_46_im * Float64(Float64(y_46_im / y_46_re) * Float64(1.0 / y_46_re))));
elseif (y_46_im <= 1.1e+96)
tmp = t_0;
else
tmp = Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * Float64(x_46_im + Float64(y_46_re * Float64(x_46_re / y_46_im))));
end
return tmp
end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(N[(N[(x$46$re * y$46$re), $MachinePrecision] + N[(x$46$im * y$46$im), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(N[(x$46$re * y$46$re + N[(y$46$im * x$46$im), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$im, -3.4e+103], N[(N[(x$46$im + N[(x$46$re / N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, -1.02e-196], t$95$0, If[LessEqual[y$46$im, 2.15e-181], N[(N[(x$46$re / y$46$re), $MachinePrecision] + N[(x$46$im * N[(N[(y$46$im / y$46$re), $MachinePrecision] * N[(1.0 / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, 1.1e+96], t$95$0, N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im + N[(y$46$re * N[(x$46$re / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
↓
\begin{array}{l}
t_0 := \frac{\frac{\mathsf{fma}\left(x.re, y.re, y.im \cdot x.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{if}\;y.im \leq -3.4 \cdot 10^{+103}:\\
\;\;\;\;\left(x.im + \frac{x.re}{\frac{y.im}{y.re}}\right) \cdot \frac{-1}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{elif}\;y.im \leq -1.02 \cdot 10^{-196}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 2.15 \cdot 10^{-181}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{elif}\;y.im \leq 1.1 \cdot 10^{+96}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.im + y.re \cdot \frac{x.re}{y.im}\right)\\
\end{array}
Alternatives
| Alternative 1 |
|---|
| Error | 13.8 |
|---|
| Cost | 7828 |
|---|
\[\begin{array}{l}
t_0 := \frac{y.im \cdot x.im + y.re \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;y.im \leq -2.4 \cdot 10^{+126}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -1.95 \cdot 10^{-123}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 4.2 \cdot 10^{-181}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{elif}\;y.im \leq 1.06 \cdot 10^{-97}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 3 \cdot 10^{+67}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{\frac{x.im}{y.re}}{\frac{y.re}{y.im}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.im + y.re \cdot \frac{x.re}{y.im}\right)\\
\end{array}
\]
| Alternative 2 |
|---|
| Error | 13.6 |
|---|
| Cost | 7828 |
|---|
\[\begin{array}{l}
t_0 := \frac{y.im \cdot x.im + y.re \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;y.im \leq -2.9 \cdot 10^{+94}:\\
\;\;\;\;\left(x.im + \frac{x.re}{\frac{y.im}{y.re}}\right) \cdot \frac{-1}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{elif}\;y.im \leq -2.05 \cdot 10^{-122}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 4.2 \cdot 10^{-181}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{elif}\;y.im \leq 6 \cdot 10^{-98}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{+67}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{\frac{x.im}{y.re}}{\frac{y.re}{y.im}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.im + y.re \cdot \frac{x.re}{y.im}\right)\\
\end{array}
\]
| Alternative 3 |
|---|
| Error | 13.9 |
|---|
| Cost | 1488 |
|---|
\[\begin{array}{l}
t_0 := \frac{y.im \cdot x.im + y.re \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;y.im \leq -3.5 \cdot 10^{+115}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -2.5 \cdot 10^{-122}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 4.2 \cdot 10^{-181}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{elif}\;y.im \leq 7.5 \cdot 10^{-98}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{+67}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{\frac{x.im}{y.re}}{\frac{y.re}{y.im}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\end{array}
\]
| Alternative 4 |
|---|
| Error | 18.0 |
|---|
| Cost | 1360 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.im \leq -2.8 \cdot 10^{+162}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -1.1 \cdot 10^{+72}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{y.im}{y.re} \cdot \frac{x.im}{y.re}\\
\mathbf{elif}\;y.im \leq -1.9 \cdot 10^{+16}:\\
\;\;\;\;\frac{y.im \cdot x.im}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{elif}\;y.im \leq 1.75 \cdot 10^{+61}:\\
\;\;\;\;\frac{x.re}{y.re} + x.im \cdot \left(\frac{y.im}{y.re} \cdot \frac{1}{y.re}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\end{array}
\]
| Alternative 5 |
|---|
| Error | 22.6 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
t_0 := y.re \cdot y.re + y.im \cdot y.im\\
\mathbf{if}\;y.re \leq -3.8 \cdot 10^{+92}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\mathbf{elif}\;y.re \leq -1.5 \cdot 10^{-95}:\\
\;\;\;\;x.im \cdot \frac{y.im}{t_0}\\
\mathbf{elif}\;y.re \leq 8 \cdot 10^{-47}:\\
\;\;\;\;\frac{x.im}{y.im}\\
\mathbf{elif}\;y.re \leq 1.3 \cdot 10^{+58}:\\
\;\;\;\;y.re \cdot \frac{x.re}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\end{array}
\]
| Alternative 6 |
|---|
| Error | 19.6 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
t_0 := \frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\mathbf{if}\;y.re \leq -5.4 \cdot 10^{+118}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\mathbf{elif}\;y.re \leq 6.6 \cdot 10^{-65}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.re \leq 1400000000000:\\
\;\;\;\;y.re \cdot \frac{x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{elif}\;y.re \leq 1.3 \cdot 10^{+23}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\end{array}
\]
| Alternative 7 |
|---|
| Error | 19.0 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.re \leq -4.6 \cdot 10^{+117}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\mathbf{elif}\;y.re \leq 3.9 \cdot 10^{-64}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re \cdot \frac{y.re}{y.im}}{y.im}\\
\mathbf{elif}\;y.re \leq 1500000000000:\\
\;\;\;\;y.re \cdot \frac{x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{elif}\;y.re \leq 2.15 \cdot 10^{+23}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\end{array}
\]
| Alternative 8 |
|---|
| Error | 18.2 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
t_0 := \frac{x.re}{y.re} + \frac{y.im}{y.re} \cdot \frac{x.im}{y.re}\\
\mathbf{if}\;y.im \leq -2.8 \cdot 10^{+162}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -4.4 \cdot 10^{+71}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.im \leq -4 \cdot 10^{+28}:\\
\;\;\;\;\frac{x.im}{y.im}\\
\mathbf{elif}\;y.im \leq 3.4 \cdot 10^{+67}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\end{array}
\]
| Alternative 9 |
|---|
| Error | 18.3 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.im \leq -2.8 \cdot 10^{+162}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -2.9 \cdot 10^{+72}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{y.im}{y.re} \cdot \frac{x.im}{y.re}\\
\mathbf{elif}\;y.im \leq -1.65 \cdot 10^{+16}:\\
\;\;\;\;\frac{x.im}{y.im}\\
\mathbf{elif}\;y.im \leq 1.05 \cdot 10^{+61}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{x.im}{y.re \cdot \frac{y.re}{y.im}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\end{array}
\]
| Alternative 10 |
|---|
| Error | 18.1 |
|---|
| Cost | 1232 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.im \leq -2.8 \cdot 10^{+162}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{\frac{x.re}{y.im}}{\frac{y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -4 \cdot 10^{+72}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{y.im}{y.re} \cdot \frac{x.im}{y.re}\\
\mathbf{elif}\;y.im \leq -1.05 \cdot 10^{+20}:\\
\;\;\;\;\frac{y.im \cdot x.im}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{elif}\;y.im \leq 7.5 \cdot 10^{+62}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{x.im}{y.re \cdot \frac{y.re}{y.im}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im} \cdot \frac{y.re}{y.im}\\
\end{array}
\]
| Alternative 11 |
|---|
| Error | 22.6 |
|---|
| Cost | 1100 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.re \leq -5 \cdot 10^{+92}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\mathbf{elif}\;y.re \leq 3.2 \cdot 10^{-43}:\\
\;\;\;\;\frac{x.im}{y.im}\\
\mathbf{elif}\;y.re \leq 1.65 \cdot 10^{+58}:\\
\;\;\;\;y.re \cdot \frac{x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\end{array}
\]
| Alternative 12 |
|---|
| Error | 25.1 |
|---|
| Cost | 722 |
|---|
\[\begin{array}{l}
\mathbf{if}\;y.im \leq -2.8 \cdot 10^{+162} \lor \neg \left(y.im \leq -4.2 \cdot 10^{+72}\right) \land \left(y.im \leq -5.1 \cdot 10^{+19} \lor \neg \left(y.im \leq 2.5 \cdot 10^{+77}\right)\right):\\
\;\;\;\;\frac{x.im}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re}\\
\end{array}
\]
| Alternative 13 |
|---|
| Error | 36.9 |
|---|
| Cost | 192 |
|---|
\[\frac{x.im}{y.im}
\]