| Alternative 1 | |
|---|---|
| Accuracy | 78.7% |
| Cost | 7236 |
(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 (* x.im y.im)) (hypot y.re y.im))
(hypot y.re y.im))))
(if (<= y.re -6.2e+147)
(+ (/ x.re y.re) (/ (/ x.im y.re) (/ y.re y.im)))
(if (<= y.re -2e-175)
t_0
(if (<= y.re 8e-127)
(+ (/ x.im y.im) (* x.re (/ (/ y.re y.im) y.im)))
(if (<= y.re 2.6e+121)
t_0
(+ (/ x.re y.re) (* (/ x.im y.re) (/ y.im y.re)))))))))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, (x_46_im * y_46_im)) / hypot(y_46_re, y_46_im)) / hypot(y_46_re, y_46_im);
double tmp;
if (y_46_re <= -6.2e+147) {
tmp = (x_46_re / y_46_re) + ((x_46_im / y_46_re) / (y_46_re / y_46_im));
} else if (y_46_re <= -2e-175) {
tmp = t_0;
} else if (y_46_re <= 8e-127) {
tmp = (x_46_im / y_46_im) + (x_46_re * ((y_46_re / y_46_im) / y_46_im));
} else if (y_46_re <= 2.6e+121) {
tmp = t_0;
} else {
tmp = (x_46_re / y_46_re) + ((x_46_im / y_46_re) * (y_46_im / y_46_re));
}
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(x_46_im * y_46_im)) / hypot(y_46_re, y_46_im)) / hypot(y_46_re, y_46_im)) tmp = 0.0 if (y_46_re <= -6.2e+147) tmp = Float64(Float64(x_46_re / y_46_re) + Float64(Float64(x_46_im / y_46_re) / Float64(y_46_re / y_46_im))); elseif (y_46_re <= -2e-175) tmp = t_0; elseif (y_46_re <= 8e-127) tmp = Float64(Float64(x_46_im / y_46_im) + Float64(x_46_re * Float64(Float64(y_46_re / y_46_im) / y_46_im))); elseif (y_46_re <= 2.6e+121) tmp = t_0; else tmp = Float64(Float64(x_46_re / y_46_re) + Float64(Float64(x_46_im / y_46_re) * Float64(y_46_im / y_46_re))); 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[(x$46$im * y$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$re, -6.2e+147], N[(N[(x$46$re / y$46$re), $MachinePrecision] + N[(N[(x$46$im / y$46$re), $MachinePrecision] / N[(y$46$re / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -2e-175], t$95$0, If[LessEqual[y$46$re, 8e-127], N[(N[(x$46$im / y$46$im), $MachinePrecision] + N[(x$46$re * N[(N[(y$46$re / y$46$im), $MachinePrecision] / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 2.6e+121], t$95$0, N[(N[(x$46$re / y$46$re), $MachinePrecision] + N[(N[(x$46$im / y$46$re), $MachinePrecision] * N[(y$46$im / y$46$re), $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, x.im \cdot y.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{if}\;y.re \leq -6.2 \cdot 10^{+147}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{\frac{x.im}{y.re}}{\frac{y.re}{y.im}}\\
\mathbf{elif}\;y.re \leq -2 \cdot 10^{-175}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y.re \leq 8 \cdot 10^{-127}:\\
\;\;\;\;\frac{x.im}{y.im} + x.re \cdot \frac{\frac{y.re}{y.im}}{y.im}\\
\mathbf{elif}\;y.re \leq 2.6 \cdot 10^{+121}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{x.im}{y.re} \cdot \frac{y.im}{y.re}\\
\end{array}
if y.re < -6.2000000000000001e147Initial program 29.0%
Taylor expanded in y.re around inf 75.6%
Simplified88.5%
[Start]75.6 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{{y.re}^{2}}
\] |
|---|---|
unpow2 [=>]75.6 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{\color{blue}{y.re \cdot y.re}}
\] |
times-frac [=>]88.5 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{y.re} \cdot \frac{x.im}{y.re}}
\] |
Applied egg-rr88.3%
[Start]88.5 | \[ \frac{x.re}{y.re} + \frac{y.im}{y.re} \cdot \frac{x.im}{y.re}
\] |
|---|---|
*-commutative [=>]88.5 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{x.im}{y.re} \cdot \frac{y.im}{y.re}}
\] |
clear-num [=>]88.3 | \[ \frac{x.re}{y.re} + \frac{x.im}{y.re} \cdot \color{blue}{\frac{1}{\frac{y.re}{y.im}}}
\] |
un-div-inv [=>]88.3 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{\frac{x.im}{y.re}}{\frac{y.re}{y.im}}}
\] |
if -6.2000000000000001e147 < y.re < -2e-175 or 8.0000000000000002e-127 < y.re < 2.5999999999999999e121Initial program 75.0%
Applied egg-rr82.7%
[Start]75.0 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-un-lft-identity [=>]75.0 | \[ \frac{\color{blue}{1 \cdot \left(x.re \cdot y.re + x.im \cdot y.im\right)}}{y.re \cdot y.re + y.im \cdot y.im}
\] |
add-sqr-sqrt [=>]75.0 | \[ \frac{1 \cdot \left(x.re \cdot y.re + x.im \cdot y.im\right)}{\color{blue}{\sqrt{y.re \cdot y.re + y.im \cdot y.im} \cdot \sqrt{y.re \cdot y.re + y.im \cdot y.im}}}
\] |
times-frac [=>]75.0 | \[ \color{blue}{\frac{1}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}} \cdot \frac{x.re \cdot y.re + x.im \cdot y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}
\] |
hypot-def [=>]75.0 | \[ \frac{1}{\color{blue}{\mathsf{hypot}\left(y.re, y.im\right)}} \cdot \frac{x.re \cdot y.re + x.im \cdot y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}
\] |
fma-def [=>]75.0 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{\color{blue}{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}
\] |
hypot-def [=>]82.7 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}{\color{blue}{\mathsf{hypot}\left(y.re, y.im\right)}}
\] |
Applied egg-rr82.9%
[Start]82.7 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}
\] |
|---|---|
associate-*l/ [=>]82.9 | \[ \color{blue}{\frac{1 \cdot \frac{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}{\mathsf{hypot}\left(y.re, y.im\right)}}
\] |
*-un-lft-identity [<=]82.9 | \[ \frac{\color{blue}{\frac{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}}{\mathsf{hypot}\left(y.re, y.im\right)}
\] |
if -2e-175 < y.re < 8.0000000000000002e-127Initial program 62.6%
Taylor expanded in y.re around 0 84.8%
Simplified83.0%
[Start]84.8 | \[ \frac{x.re \cdot y.re}{{y.im}^{2}} + \frac{x.im}{y.im}
\] |
|---|---|
+-commutative [=>]84.8 | \[ \color{blue}{\frac{x.im}{y.im} + \frac{x.re \cdot y.re}{{y.im}^{2}}}
\] |
unpow2 [=>]84.8 | \[ \frac{x.im}{y.im} + \frac{x.re \cdot y.re}{\color{blue}{y.im \cdot y.im}}
\] |
associate-/l* [=>]83.0 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{x.re}{\frac{y.im \cdot y.im}{y.re}}}
\] |
Applied egg-rr89.7%
[Start]83.0 | \[ \frac{x.im}{y.im} + \frac{x.re}{\frac{y.im \cdot y.im}{y.re}}
\] |
|---|---|
clear-num [=>]83.0 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{1}{\frac{\frac{y.im \cdot y.im}{y.re}}{x.re}}}
\] |
associate-/r/ [=>]83.0 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{1}{\frac{y.im \cdot y.im}{y.re}} \cdot x.re}
\] |
clear-num [<=]83.2 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{y.re}{y.im \cdot y.im}} \cdot x.re
\] |
associate-/r* [=>]89.7 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{\frac{y.re}{y.im}}{y.im}} \cdot x.re
\] |
if 2.5999999999999999e121 < y.re Initial program 35.9%
Taylor expanded in y.re around inf 73.5%
Simplified82.9%
[Start]73.5 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{{y.re}^{2}}
\] |
|---|---|
unpow2 [=>]73.5 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{\color{blue}{y.re \cdot y.re}}
\] |
times-frac [=>]82.9 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{y.re} \cdot \frac{x.im}{y.re}}
\] |
Final simplification85.2%
| Alternative 1 | |
|---|---|
| Accuracy | 78.7% |
| Cost | 7236 |
| Alternative 2 | |
|---|---|
| Accuracy | 72.5% |
| Cost | 1497 |
| Alternative 3 | |
|---|---|
| Accuracy | 72.3% |
| Cost | 1496 |
| Alternative 4 | |
|---|---|
| Accuracy | 72.3% |
| Cost | 1496 |
| Alternative 5 | |
|---|---|
| Accuracy | 72.2% |
| Cost | 1496 |
| Alternative 6 | |
|---|---|
| Accuracy | 78.1% |
| Cost | 1488 |
| Alternative 7 | |
|---|---|
| Accuracy | 69.3% |
| Cost | 969 |
| Alternative 8 | |
|---|---|
| Accuracy | 63.2% |
| Cost | 456 |
| Alternative 9 | |
|---|---|
| Accuracy | 41.8% |
| Cost | 324 |
| Alternative 10 | |
|---|---|
| Accuracy | 41.4% |
| Cost | 192 |
herbie shell --seed 2023147
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, real part"
:precision binary64
(/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))))