| Alternative 1 | |
|---|---|
| Accuracy | 78.3% |
| Cost | 20560 |
(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))))
(if (<= y.im -7.4e+118)
(+ (/ x.im y.im) (/ x.re (/ (* y.im y.im) y.re)))
(if (<= y.im -6.2e-155)
(* (/ 1.0 (hypot y.re y.im)) t_0)
(if (<= y.im 8.5e-155)
(+ (/ x.re y.re) (/ (/ (* y.im x.im) y.re) y.re))
(if (<= y.im 2.5e+43)
(/ 1.0 (/ (hypot y.re y.im) t_0))
(+ (/ x.im y.im) (* y.re (/ (/ x.re y.im) 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);
double tmp;
if (y_46_im <= -7.4e+118) {
tmp = (x_46_im / y_46_im) + (x_46_re / ((y_46_im * y_46_im) / y_46_re));
} else if (y_46_im <= -6.2e-155) {
tmp = (1.0 / hypot(y_46_re, y_46_im)) * t_0;
} else if (y_46_im <= 8.5e-155) {
tmp = (x_46_re / y_46_re) + (((y_46_im * x_46_im) / y_46_re) / y_46_re);
} else if (y_46_im <= 2.5e+43) {
tmp = 1.0 / (hypot(y_46_re, y_46_im) / t_0);
} else {
tmp = (x_46_im / y_46_im) + (y_46_re * ((x_46_re / y_46_im) / 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(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_im <= -7.4e+118) tmp = Float64(Float64(x_46_im / y_46_im) + Float64(x_46_re / Float64(Float64(y_46_im * y_46_im) / y_46_re))); elseif (y_46_im <= -6.2e-155) tmp = Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * t_0); elseif (y_46_im <= 8.5e-155) tmp = Float64(Float64(x_46_re / y_46_re) + Float64(Float64(Float64(y_46_im * x_46_im) / y_46_re) / y_46_re)); elseif (y_46_im <= 2.5e+43) tmp = Float64(1.0 / Float64(hypot(y_46_re, y_46_im) / t_0)); else tmp = Float64(Float64(x_46_im / y_46_im) + Float64(y_46_re * Float64(Float64(x_46_re / y_46_im) / 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[(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]}, If[LessEqual[y$46$im, -7.4e+118], N[(N[(x$46$im / y$46$im), $MachinePrecision] + N[(x$46$re / N[(N[(y$46$im * y$46$im), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, -6.2e-155], N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[y$46$im, 8.5e-155], N[(N[(x$46$re / y$46$re), $MachinePrecision] + N[(N[(N[(y$46$im * x$46$im), $MachinePrecision] / y$46$re), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, 2.5e+43], N[(1.0 / N[(N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$im / y$46$im), $MachinePrecision] + N[(y$46$re * N[(N[(x$46$re / y$46$im), $MachinePrecision] / y$46$im), $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{\mathsf{fma}\left(x.re, y.re, y.im \cdot x.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{if}\;y.im \leq -7.4 \cdot 10^{+118}:\\
\;\;\;\;\frac{x.im}{y.im} + \frac{x.re}{\frac{y.im \cdot y.im}{y.re}}\\
\mathbf{elif}\;y.im \leq -6.2 \cdot 10^{-155}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot t_0\\
\mathbf{elif}\;y.im \leq 8.5 \cdot 10^{-155}:\\
\;\;\;\;\frac{x.re}{y.re} + \frac{\frac{y.im \cdot x.im}{y.re}}{y.re}\\
\mathbf{elif}\;y.im \leq 2.5 \cdot 10^{+43}:\\
\;\;\;\;\frac{1}{\frac{\mathsf{hypot}\left(y.re, y.im\right)}{t_0}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.im} + y.re \cdot \frac{\frac{x.re}{y.im}}{y.im}\\
\end{array}
if y.im < -7.39999999999999973e118Initial program 35.5%
Applied egg-rr18.1%
[Start]35.5 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-commutative [=>]35.5 | \[ \frac{\color{blue}{y.re \cdot x.re} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
add-sqr-sqrt [=>]18.1 | \[ \frac{y.re \cdot \color{blue}{\left(\sqrt{x.re} \cdot \sqrt{x.re}\right)} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
associate-*r* [=>]18.1 | \[ \frac{\color{blue}{\left(y.re \cdot \sqrt{x.re}\right) \cdot \sqrt{x.re}} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
fma-def [=>]18.1 | \[ \frac{\color{blue}{\mathsf{fma}\left(y.re \cdot \sqrt{x.re}, \sqrt{x.re}, x.im \cdot y.im\right)}}{y.re \cdot y.re + y.im \cdot y.im}
\] |
Taylor expanded in y.re around 0 76.3%
Simplified77.8%
[Start]76.3 | \[ \frac{x.re \cdot y.re}{{y.im}^{2}} + \frac{x.im}{y.im}
\] |
|---|---|
+-commutative [=>]76.3 | \[ \color{blue}{\frac{x.im}{y.im} + \frac{x.re \cdot y.re}{{y.im}^{2}}}
\] |
associate-/l* [=>]77.8 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{x.re}{\frac{{y.im}^{2}}{y.re}}}
\] |
unpow2 [=>]77.8 | \[ \frac{x.im}{y.im} + \frac{x.re}{\frac{\color{blue}{y.im \cdot y.im}}{y.re}}
\] |
if -7.39999999999999973e118 < y.im < -6.2e-155Initial program 71.6%
Applied egg-rr78.5%
[Start]71.6 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-un-lft-identity [=>]71.6 | \[ \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 [=>]71.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 [=>]71.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 [=>]71.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 [=>]71.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 [=>]78.5 | \[ \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 -6.2e-155 < y.im < 8.4999999999999996e-155Initial program 53.6%
Taylor expanded in y.re around inf 73.4%
Simplified72.1%
[Start]73.4 | \[ \frac{x.re}{y.re} + \frac{y.im \cdot x.im}{{y.re}^{2}}
\] |
|---|---|
associate-/l* [=>]69.9 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{\frac{{y.re}^{2}}{x.im}}}
\] |
associate-/r/ [=>]72.1 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im}{{y.re}^{2}} \cdot x.im}
\] |
unpow2 [=>]72.1 | \[ \frac{x.re}{y.re} + \frac{y.im}{\color{blue}{y.re \cdot y.re}} \cdot x.im
\] |
Applied egg-rr78.1%
[Start]72.1 | \[ \frac{x.re}{y.re} + \frac{y.im}{y.re \cdot y.re} \cdot x.im
\] |
|---|---|
associate-*l/ [=>]73.4 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{y.im \cdot x.im}{y.re \cdot y.re}}
\] |
associate-/r* [=>]78.1 | \[ \frac{x.re}{y.re} + \color{blue}{\frac{\frac{y.im \cdot x.im}{y.re}}{y.re}}
\] |
if 8.4999999999999996e-155 < y.im < 2.5000000000000002e43Initial program 72.8%
Applied egg-rr80.4%
[Start]72.8 | \[ \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.8 | \[ \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.8 | \[ \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.7 | \[ \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.7 | \[ \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.7 | \[ \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 [=>]80.4 | \[ \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-rr80.0%
[Start]80.4 | \[ \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/ [=>]80.6 | \[ \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)}}
\] |
associate-/l* [=>]80.0 | \[ \color{blue}{\frac{1}{\frac{\mathsf{hypot}\left(y.re, y.im\right)}{\frac{\mathsf{fma}\left(x.re, y.re, x.im \cdot y.im\right)}{\mathsf{hypot}\left(y.re, y.im\right)}}}}
\] |
if 2.5000000000000002e43 < y.im Initial program 44.4%
Applied egg-rr23.2%
[Start]44.4 | \[ \frac{x.re \cdot y.re + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
|---|---|
*-commutative [=>]44.4 | \[ \frac{\color{blue}{y.re \cdot x.re} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
add-sqr-sqrt [=>]23.2 | \[ \frac{y.re \cdot \color{blue}{\left(\sqrt{x.re} \cdot \sqrt{x.re}\right)} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
associate-*r* [=>]23.2 | \[ \frac{\color{blue}{\left(y.re \cdot \sqrt{x.re}\right) \cdot \sqrt{x.re}} + x.im \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\] |
fma-def [=>]23.2 | \[ \frac{\color{blue}{\mathsf{fma}\left(y.re \cdot \sqrt{x.re}, \sqrt{x.re}, x.im \cdot y.im\right)}}{y.re \cdot y.re + y.im \cdot y.im}
\] |
Taylor expanded in y.re around 0 70.9%
Simplified77.5%
[Start]70.9 | \[ \frac{x.re \cdot y.re}{{y.im}^{2}} + \frac{x.im}{y.im}
\] |
|---|---|
+-commutative [=>]70.9 | \[ \color{blue}{\frac{x.im}{y.im} + \frac{x.re \cdot y.re}{{y.im}^{2}}}
\] |
associate-/l* [=>]71.4 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{x.re}{\frac{{y.im}^{2}}{y.re}}}
\] |
associate-/r/ [=>]73.0 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{x.re}{{y.im}^{2}} \cdot y.re}
\] |
unpow2 [=>]73.0 | \[ \frac{x.im}{y.im} + \frac{x.re}{\color{blue}{y.im \cdot y.im}} \cdot y.re
\] |
associate-/r* [=>]77.5 | \[ \frac{x.im}{y.im} + \color{blue}{\frac{\frac{x.re}{y.im}}{y.im}} \cdot y.re
\] |
Final simplification78.3%
| Alternative 1 | |
|---|---|
| Accuracy | 78.3% |
| Cost | 20560 |
| Alternative 2 | |
|---|---|
| Accuracy | 67.1% |
| Cost | 1496 |
| Alternative 3 | |
|---|---|
| Accuracy | 75.6% |
| Cost | 1488 |
| Alternative 4 | |
|---|---|
| Accuracy | 68.2% |
| Cost | 1232 |
| Alternative 5 | |
|---|---|
| Accuracy | 69.0% |
| Cost | 1232 |
| Alternative 6 | |
|---|---|
| Accuracy | 71.7% |
| Cost | 1228 |
| Alternative 7 | |
|---|---|
| Accuracy | 71.3% |
| Cost | 1100 |
| Alternative 8 | |
|---|---|
| Accuracy | 67.2% |
| Cost | 969 |
| Alternative 9 | |
|---|---|
| Accuracy | 66.1% |
| Cost | 969 |
| Alternative 10 | |
|---|---|
| Accuracy | 60.4% |
| Cost | 720 |
| Alternative 11 | |
|---|---|
| Accuracy | 61.4% |
| Cost | 456 |
| Alternative 12 | |
|---|---|
| Accuracy | 38.8% |
| Cost | 192 |
herbie shell --seed 2023153
(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))))