| Alternative 1 | |
|---|---|
| Accuracy | 80.5% |
| Cost | 7828 |
(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 (/ 1.0 (hypot y.re y.im)))
(t_1 (* t_0 (/ (fma x.re y.re (* y.im x.im)) (hypot y.re y.im)))))
(if (<= y.re -4.6e+73)
(+ (/ x.re y.re) (/ (* y.im (/ x.im y.re)) y.re))
(if (<= y.re -8.6e-160)
t_1
(if (<= y.re 4.1e-262)
(+ (/ x.im y.im) (/ x.re (* y.im (/ y.im y.re))))
(if (<= y.re 5.8e+97)
t_1
(* t_0 (+ x.re (* x.im (/ 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 = 1.0 / hypot(y_46_re, y_46_im);
double t_1 = t_0 * (fma(x_46_re, y_46_re, (y_46_im * x_46_im)) / hypot(y_46_re, y_46_im));
double tmp;
if (y_46_re <= -4.6e+73) {
tmp = (x_46_re / y_46_re) + ((y_46_im * (x_46_im / y_46_re)) / y_46_re);
} else if (y_46_re <= -8.6e-160) {
tmp = t_1;
} else if (y_46_re <= 4.1e-262) {
tmp = (x_46_im / y_46_im) + (x_46_re / (y_46_im * (y_46_im / y_46_re)));
} else if (y_46_re <= 5.8e+97) {
tmp = t_1;
} else {
tmp = t_0 * (x_46_re + (x_46_im * (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(1.0 / hypot(y_46_re, y_46_im)) t_1 = Float64(t_0 * Float64(fma(x_46_re, y_46_re, Float64(y_46_im * x_46_im)) / hypot(y_46_re, y_46_im))) tmp = 0.0 if (y_46_re <= -4.6e+73) tmp = Float64(Float64(x_46_re / y_46_re) + Float64(Float64(y_46_im * Float64(x_46_im / y_46_re)) / y_46_re)); elseif (y_46_re <= -8.6e-160) tmp = t_1; elseif (y_46_re <= 4.1e-262) tmp = Float64(Float64(x_46_im / y_46_im) + Float64(x_46_re / Float64(y_46_im * Float64(y_46_im / y_46_re)))); elseif (y_46_re <= 5.8e+97) tmp = t_1; else tmp = Float64(t_0 * Float64(x_46_re + Float64(x_46_im * 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[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * 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]), $MachinePrecision]}, If[LessEqual[y$46$re, -4.6e+73], N[(N[(x$46$re / y$46$re), $MachinePrecision] + N[(N[(y$46$im * N[(x$46$im / y$46$re), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -8.6e-160], t$95$1, If[LessEqual[y$46$re, 4.1e-262], N[(N[(x$46$im / y$46$im), $MachinePrecision] + N[(x$46$re / N[(y$46$im * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 5.8e+97], t$95$1, N[(t$95$0 * N[(x$46$re + N[(x$46$im * N[(y$46$im / y$46$re), $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{1}{\mathsf{hypot}\left(y.re, y.im\right)}\\
t_1 := t_0 \cdot \frac{\mathsf{fma}\left(x.re, y.re, y.im \cdot x.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{if}\;y.re \leq -4.6 \cdot 10^{+73}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{y.im \cdot \frac{x.im}{y.re}}{y.re}\\
\mathbf{elif}\;y.re \leq -8.6 \cdot 10^{-160}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y.re \leq 4.1 \cdot 10^{-262}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{y.im \cdot \frac{y.im}{y.re}}\\
\mathbf{elif}\;y.re \leq 5.8 \cdot 10^{+97}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot \left(x.re + x.im \cdot \frac{y.im}{y.re}\right)\\
\end{array}
if y.re < -4.6e73Initial program 40.3%
Taylor expanded in y.re around inf 72.9%
Simplified73.0%
[Start]72.9 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{{y.re}^{2}}
\] |
|---|---|
associate-/l* [=>]74.9 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{\frac{{y.re}^{2}}{x.im}}}
\] |
associate-/r/ [=>]73.0 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{{y.re}^{2}} \cdot x.im}
\] |
unpow2 [=>]73.0 | \[ \frac{x.re}{y.re} + \frac{y.im}{\color{blue}{y.re \cdot y.re}} \cdot x.im
\] |
Applied egg-rr81.6%
[Start]73.0 | \[ \frac{x.re}{y.re} + \frac{y.im}{y.re \cdot y.re} \cdot x.im
\] |
|---|---|
associate-/r* [=>]75.8 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{\frac{y.im}{y.re}}{y.re}} \cdot x.im
\] |
associate-*l/ [=>]81.6 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{\frac{y.im}{y.re} \cdot x.im}{y.re}}
\] |
Taylor expanded in y.im around 0 76.8%
Simplified82.0%
[Start]76.8 | \[ \frac{x.re}{y.re} + \frac{\frac{y.im \cdot x.im}{y.re}}{y.re}
\] |
|---|---|
associate-*r/ [<=]82.0 | \[ \frac{x.re}{y.re} + \frac{\color{blue}{y.im \cdot \frac{x.im}{y.re}}}{y.re}
\] |
if -4.6e73 < y.re < -8.60000000000000028e-160 or 4.10000000000000026e-262 < y.re < 5.79999999999999974e97Initial program 72.3%
Applied egg-rr81.2%
[Start]72.3 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-un-lft-identity [=>]72.3 | \[ \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 [=>]72.3 | \[ \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 [=>]72.3 | \[ \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 [=>]72.3 | \[ \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 [=>]72.3 | \[ \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 [=>]81.2 | \[ \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)}}
\] |
if -8.60000000000000028e-160 < y.re < 4.10000000000000026e-262Initial program 62.2%
Taylor expanded in y.re around 0 85.5%
Simplified89.8%
[Start]85.5 | \[ \frac{x.re \cdot y.re}{{y.im}^{2}} + \frac{x.im}{y.im}
\] |
|---|---|
+-commutative [=>]85.5 | \[ \color{blue}{\frac{x.im}{y.im} + \frac{x.re \cdot y.re}{{y.im}^{2}}}
\] |
*-commutative [=>]85.5 | \[ \frac{x.im}{y.im} + \frac{\color{blue}{y.re \cdot x.re}}{{y.im}^{2}}
\] |
unpow2 [=>]85.5 | \[ \frac{x.im}{y.im} + \frac{y.re \cdot x.re}{\color{blue}{y.im \cdot y.im}}
\] |
times-frac [=>]89.8 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{y.re}{y.im} \cdot \frac{x.re}{y.im}}
\] |
Applied egg-rr92.4%
[Start]89.8 | \[ \frac{x.im}{y.im} + \frac{y.re}{y.im} \cdot \frac{x.re}{y.im}
\] |
|---|---|
clear-num [=>]89.6 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{1}{\frac{y.im}{y.re}}} \cdot \frac{x.re}{y.im}
\] |
frac-times [=>]92.4 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{1 \cdot x.re}{\frac{y.im}{y.re} \cdot y.im}}
\] |
*-un-lft-identity [<=]92.4 | \[ \frac{x.im}{y.im} + \frac{\color{blue}{x.re}}{\frac{y.im}{y.re} \cdot y.im}
\] |
if 5.79999999999999974e97 < y.re Initial program 39.5%
Applied egg-rr59.2%
[Start]39.5 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-un-lft-identity [=>]39.5 | \[ \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 [=>]39.5 | \[ \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 [=>]39.5 | \[ \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 [=>]39.5 | \[ \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 [=>]39.5 | \[ \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 [=>]59.2 | \[ \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)}}
\] |
Taylor expanded in y.re around inf 77.9%
Simplified83.6%
[Start]77.9 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.re + \frac{y.im \cdot x.im}{y.re}\right)
\] |
|---|---|
associate-/l* [=>]84.3 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.re + \color{blue}{\frac{y.im}{\frac{y.re}{x.im}}}\right)
\] |
associate-/r/ [=>]83.6 | \[ \frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \left(x.re + \color{blue}{\frac{y.im}{y.re} \cdot x.im}\right)
\] |
Final simplification83.4%
| Alternative 1 | |
|---|---|
| Accuracy | 80.5% |
| Cost | 7828 |
| Alternative 2 | |
|---|---|
| Accuracy | 80.5% |
| Cost | 7636 |
| Alternative 3 | |
|---|---|
| Accuracy | 80.6% |
| Cost | 1488 |
| Alternative 4 | |
|---|---|
| Accuracy | 69.0% |
| Cost | 1233 |
| Alternative 5 | |
|---|---|
| Accuracy | 75.9% |
| Cost | 1232 |
| Alternative 6 | |
|---|---|
| Accuracy | 73.0% |
| Cost | 969 |
| Alternative 7 | |
|---|---|
| Accuracy | 74.8% |
| Cost | 969 |
| Alternative 8 | |
|---|---|
| Accuracy | 76.1% |
| Cost | 969 |
| Alternative 9 | |
|---|---|
| Accuracy | 63.2% |
| Cost | 456 |
| Alternative 10 | |
|---|---|
| Accuracy | 41.8% |
| Cost | 192 |
herbie shell --seed 2023151
(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))))