
(FPCore (x.re x.im y.re y.im) :precision binary64 (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.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_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
code = ((x_46im * y_46re) - (x_46re * y_46im)) / ((y_46re * y_46re) + (y_46im * y_46im))
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): return ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))
function code(x_46_re, x_46_im, y_46_re, y_46_im) return Float64(Float64(Float64(x_46_im * y_46_re) - Float64(x_46_re * y_46_im)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(x$46$re * y$46$im), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x.im \cdot y.re - x.re \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 32 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x.re x.im y.re y.im) :precision binary64 (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.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_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
code = ((x_46im * y_46re) - (x_46re * y_46im)) / ((y_46re * y_46re) + (y_46im * y_46im))
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): return ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))
function code(x_46_re, x_46_im, y_46_re, y_46_im) return Float64(Float64(Float64(x_46_im * y_46_re) - Float64(x_46_re * y_46_im)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = ((x_46_im * y_46_re) - (x_46_re * y_46_im)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(x$46$re * y$46$im), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x.im \cdot y.re - x.re \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}
\end{array}
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
2e+247)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(* x.re (/ (/ (- y.im) (hypot y.im y.re)) (hypot y.im y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= 2e+247) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (x_46_re * ((-y_46_im / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= 2e+247) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(x_46_re * Float64(Float64(Float64(-y_46_im) / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re)))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+247], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(x$46$re * N[(N[((-y$46$im) / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq 2 \cdot 10^{+247}:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, x.re \cdot \frac{\frac{-y.im}{\mathsf{hypot}\left(y.im, y.re\right)}}{\mathsf{hypot}\left(y.im, y.re\right)}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < 1.9999999999999999e247Initial program 76.6%
*-un-lft-identity76.6%
add-sqr-sqrt76.6%
times-frac76.5%
hypot-define76.5%
fma-neg76.5%
distribute-rgt-neg-in76.5%
hypot-define96.9%
Applied egg-rr96.9%
associate-*r/97.1%
Applied egg-rr97.1%
if 1.9999999999999999e247 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 24.1%
div-sub18.1%
*-commutative18.1%
add-sqr-sqrt18.1%
times-frac22.3%
fma-neg22.3%
hypot-define22.3%
hypot-define56.3%
associate-/l*68.8%
add-sqr-sqrt68.8%
pow268.8%
hypot-define68.8%
Applied egg-rr68.8%
*-un-lft-identity68.8%
unpow268.8%
times-frac95.6%
hypot-undefine68.8%
+-commutative68.8%
hypot-define95.6%
hypot-undefine68.8%
+-commutative68.8%
hypot-define95.6%
Applied egg-rr95.6%
associate-*l/95.6%
*-lft-identity95.6%
Simplified95.6%
Final simplification96.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
2e+247)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
(* y.re (/ 1.0 (hypot y.im y.re)))
(/ x.im (hypot y.re y.im))
(* x.re (/ (/ (- y.im) (hypot y.im y.re)) (hypot y.im y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= 2e+247) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma((y_46_re * (1.0 / hypot(y_46_im, y_46_re))), (x_46_im / hypot(y_46_re, y_46_im)), (x_46_re * ((-y_46_im / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= 2e+247) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(Float64(y_46_re * Float64(1.0 / hypot(y_46_im, y_46_re))), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(x_46_re * Float64(Float64(Float64(-y_46_im) / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re)))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+247], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(N[(y$46$re * N[(1.0 / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(x$46$re * N[(N[((-y$46$im) / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq 2 \cdot 10^{+247}:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y.re \cdot \frac{1}{\mathsf{hypot}\left(y.im, y.re\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, x.re \cdot \frac{\frac{-y.im}{\mathsf{hypot}\left(y.im, y.re\right)}}{\mathsf{hypot}\left(y.im, y.re\right)}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < 1.9999999999999999e247Initial program 76.6%
*-un-lft-identity76.6%
add-sqr-sqrt76.6%
times-frac76.5%
hypot-define76.5%
fma-neg76.5%
distribute-rgt-neg-in76.5%
hypot-define96.9%
Applied egg-rr96.9%
associate-*r/97.1%
Applied egg-rr97.1%
if 1.9999999999999999e247 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 24.1%
div-sub18.1%
*-commutative18.1%
add-sqr-sqrt18.1%
times-frac22.3%
fma-neg22.3%
hypot-define22.3%
hypot-define56.3%
associate-/l*68.8%
add-sqr-sqrt68.8%
pow268.8%
hypot-define68.8%
Applied egg-rr68.8%
*-un-lft-identity68.8%
unpow268.8%
times-frac95.6%
hypot-undefine68.8%
+-commutative68.8%
hypot-define95.6%
hypot-undefine68.8%
+-commutative68.8%
hypot-define95.6%
Applied egg-rr95.6%
associate-*l/95.6%
*-lft-identity95.6%
Simplified95.6%
div-inv95.5%
hypot-undefine50.6%
+-commutative50.6%
hypot-define95.5%
Applied egg-rr95.5%
Final simplification96.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
2e+245)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(* x.re (/ (- y.im) (pow (hypot y.re y.im) 2.0))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= 2e+245) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (x_46_re * (-y_46_im / pow(hypot(y_46_re, y_46_im), 2.0))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= 2e+245) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(x_46_re * Float64(Float64(-y_46_im) / (hypot(y_46_re, y_46_im) ^ 2.0)))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+245], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(x$46$re * N[((-y$46$im) / N[Power[N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq 2 \cdot 10^{+245}:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, x.re \cdot \frac{-y.im}{{\left(\mathsf{hypot}\left(y.re, y.im\right)\right)}^{2}}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < 2.00000000000000009e245Initial program 76.5%
*-un-lft-identity76.5%
add-sqr-sqrt76.5%
times-frac76.4%
hypot-define76.4%
fma-neg76.4%
distribute-rgt-neg-in76.4%
hypot-define96.9%
Applied egg-rr96.9%
associate-*r/97.0%
Applied egg-rr97.0%
if 2.00000000000000009e245 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 25.2%
div-sub19.3%
*-commutative19.3%
add-sqr-sqrt19.3%
times-frac23.4%
fma-neg23.4%
hypot-define23.4%
hypot-define56.9%
associate-/l*69.2%
add-sqr-sqrt69.2%
pow269.2%
hypot-define69.2%
Applied egg-rr69.2%
Final simplification89.5%
(FPCore (x.re x.im y.re y.im) :precision binary64 (fma (/ y.re (hypot y.re y.im)) (/ x.im (hypot y.re y.im)) (/ -1.0 (/ (hypot y.im y.re) (* x.re (/ y.im (hypot y.im y.re)))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (-1.0 / (hypot(y_46_im, y_46_re) / (x_46_re * (y_46_im / hypot(y_46_im, y_46_re))))));
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) return fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(-1.0 / Float64(hypot(y_46_im, y_46_re) / Float64(x_46_re * Float64(y_46_im / hypot(y_46_im, y_46_re)))))) end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(-1.0 / N[(N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision] / N[(x$46$re * N[(y$46$im / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{-1}{\frac{\mathsf{hypot}\left(y.im, y.re\right)}{x.re \cdot \frac{y.im}{\mathsf{hypot}\left(y.im, y.re\right)}}}\right)
\end{array}
Initial program 62.7%
div-sub59.7%
*-commutative59.7%
add-sqr-sqrt59.7%
times-frac60.1%
fma-neg60.1%
hypot-define60.1%
hypot-define77.5%
associate-/l*80.4%
add-sqr-sqrt80.4%
pow280.4%
hypot-define80.4%
Applied egg-rr80.4%
*-un-lft-identity80.4%
unpow280.4%
times-frac94.3%
hypot-undefine80.4%
+-commutative80.4%
hypot-define94.3%
hypot-undefine80.4%
+-commutative80.4%
hypot-define94.3%
Applied egg-rr94.3%
associate-*l/94.3%
*-lft-identity94.3%
Simplified94.3%
associate-*r/96.6%
clear-num95.9%
Applied egg-rr95.9%
Final simplification95.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
INFINITY)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(/ 1.0 (/ (hypot y.im y.re) (* x.re (/ y.im y.re)))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= ((double) INFINITY)) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (1.0 / (hypot(y_46_im, y_46_re) / (x_46_re * (y_46_im / y_46_re)))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= Inf) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(1.0 / Float64(hypot(y_46_im, y_46_re) / Float64(x_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_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(1.0 / N[(N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision] / N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq \infty:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{1}{\frac{\mathsf{hypot}\left(y.im, y.re\right)}{x.re \cdot \frac{y.im}{y.re}}}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < +inf.0Initial program 75.7%
*-un-lft-identity75.7%
add-sqr-sqrt75.7%
times-frac75.6%
hypot-define75.6%
fma-neg75.6%
distribute-rgt-neg-in75.6%
hypot-define94.6%
Applied egg-rr94.6%
associate-*r/94.7%
Applied egg-rr94.7%
if +inf.0 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 0.0%
div-sub0.0%
*-commutative0.0%
add-sqr-sqrt0.0%
times-frac2.1%
fma-neg2.1%
hypot-define2.1%
hypot-define49.9%
associate-/l*58.3%
add-sqr-sqrt58.3%
pow258.3%
hypot-define58.3%
Applied egg-rr58.3%
*-un-lft-identity58.3%
unpow258.3%
times-frac95.5%
hypot-undefine58.3%
+-commutative58.3%
hypot-define95.5%
hypot-undefine58.3%
+-commutative58.3%
hypot-define95.5%
Applied egg-rr95.5%
associate-*l/95.5%
*-lft-identity95.5%
Simplified95.5%
associate-*r/99.8%
clear-num99.8%
Applied egg-rr99.8%
Taylor expanded in y.re around -inf 50.4%
mul-1-neg50.4%
associate-*r/59.6%
distribute-rgt-neg-in59.6%
distribute-neg-frac259.6%
Simplified59.6%
Final simplification88.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
5e+286)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
(/ y.re y.im)
(/ x.im (hypot y.re y.im))
(* x.re (/ (/ (- y.im) (hypot y.im y.re)) (hypot y.im y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= 5e+286) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma((y_46_re / y_46_im), (x_46_im / hypot(y_46_re, y_46_im)), (x_46_re * ((-y_46_im / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= 5e+286) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(Float64(y_46_re / y_46_im), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(x_46_re * Float64(Float64(Float64(-y_46_im) / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re)))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 5e+286], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(N[(y$46$re / y$46$im), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(x$46$re * N[(N[((-y$46$im) / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq 5 \cdot 10^{+286}:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{y.im}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, x.re \cdot \frac{\frac{-y.im}{\mathsf{hypot}\left(y.im, y.re\right)}}{\mathsf{hypot}\left(y.im, y.re\right)}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < 5.0000000000000004e286Initial program 77.2%
*-un-lft-identity77.2%
add-sqr-sqrt77.2%
times-frac77.1%
hypot-define77.1%
fma-neg77.1%
distribute-rgt-neg-in77.1%
hypot-define97.0%
Applied egg-rr97.0%
associate-*r/97.1%
Applied egg-rr97.1%
if 5.0000000000000004e286 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 18.1%
div-sub11.6%
*-commutative11.6%
add-sqr-sqrt11.6%
times-frac16.1%
fma-neg16.1%
hypot-define16.1%
hypot-define52.8%
associate-/l*66.3%
add-sqr-sqrt66.3%
pow266.3%
hypot-define66.3%
Applied egg-rr66.3%
*-un-lft-identity66.3%
unpow266.3%
times-frac95.2%
hypot-undefine66.3%
+-commutative66.3%
hypot-define95.2%
hypot-undefine66.3%
+-commutative66.3%
hypot-define95.2%
Applied egg-rr95.2%
associate-*l/95.3%
*-lft-identity95.3%
Simplified95.3%
Taylor expanded in y.re around 0 51.4%
Final simplification85.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<=
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))
INFINITY)
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))
(fma
1.0
(/ x.im (hypot y.re y.im))
(* x.re (/ (/ (- y.im) (hypot y.im y.re)) (hypot y.im y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))) <= ((double) INFINITY)) {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
} else {
tmp = fma(1.0, (x_46_im / hypot(y_46_re, y_46_im)), (x_46_re * ((-y_46_im / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re))));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) <= Inf) tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); else tmp = fma(1.0, Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(x_46_re * Float64(Float64(Float64(-y_46_im) / hypot(y_46_im, y_46_re)) / hypot(y_46_im, y_46_re)))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision], N[(1.0 * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[(x$46$re * N[(N[((-y$46$im) / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im} \leq \infty:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, x.re \cdot \frac{\frac{-y.im}{\mathsf{hypot}\left(y.im, y.re\right)}}{\mathsf{hypot}\left(y.im, y.re\right)}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < +inf.0Initial program 75.7%
*-un-lft-identity75.7%
add-sqr-sqrt75.7%
times-frac75.6%
hypot-define75.6%
fma-neg75.6%
distribute-rgt-neg-in75.6%
hypot-define94.6%
Applied egg-rr94.6%
associate-*r/94.7%
Applied egg-rr94.7%
if +inf.0 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 0.0%
div-sub0.0%
*-commutative0.0%
add-sqr-sqrt0.0%
times-frac2.1%
fma-neg2.1%
hypot-define2.1%
hypot-define49.9%
associate-/l*58.3%
add-sqr-sqrt58.3%
pow258.3%
hypot-define58.3%
Applied egg-rr58.3%
*-un-lft-identity58.3%
unpow258.3%
times-frac95.5%
hypot-undefine58.3%
+-commutative58.3%
hypot-define95.5%
hypot-undefine58.3%
+-commutative58.3%
hypot-define95.5%
Applied egg-rr95.5%
associate-*l/95.5%
*-lft-identity95.5%
Simplified95.5%
Taylor expanded in y.re around inf 49.2%
Final simplification86.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))))
(if (<= t_0 INFINITY)
t_0
(*
(/ y.im (hypot y.im y.re))
(/ (fma x.im (/ y.re y.im) x.re) (hypot y.im y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
double tmp;
if (t_0 <= ((double) INFINITY)) {
tmp = t_0;
} else {
tmp = (y_46_im / hypot(y_46_im, y_46_re)) * (fma(x_46_im, (y_46_re / y_46_im), x_46_re) / hypot(y_46_im, y_46_re));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) tmp = 0.0 if (t_0 <= Inf) tmp = t_0; else tmp = Float64(Float64(y_46_im / hypot(y_46_im, y_46_re)) * Float64(fma(x_46_im, Float64(y_46_re / y_46_im), x_46_re) / hypot(y_46_im, y_46_re))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, Infinity], t$95$0, N[(N[(y$46$im / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision] * N[(N[(x$46$im * N[(y$46$re / y$46$im), $MachinePrecision] + x$46$re), $MachinePrecision] / N[Sqrt[y$46$im ^ 2 + y$46$re ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;t\_0 \leq \infty:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im}{\mathsf{hypot}\left(y.im, y.re\right)} \cdot \frac{\mathsf{fma}\left(x.im, \frac{y.re}{y.im}, x.re\right)}{\mathsf{hypot}\left(y.im, y.re\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < +inf.0Initial program 75.7%
if +inf.0 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 0.0%
Taylor expanded in y.im around inf 0.0%
*-commutative0.0%
sub-neg0.0%
mul-1-neg0.0%
+-commutative0.0%
add-sqr-sqrt0.0%
hypot-undefine0.0%
hypot-undefine0.0%
times-frac34.7%
Applied egg-rr16.0%
Final simplification65.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))))
(if (<= t_0 1e+261)
t_0
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(/ (- x.re) y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
double tmp;
if (t_0 <= 1e+261) {
tmp = t_0;
} else {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (-x_46_re / y_46_im));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) tmp = 0.0 if (t_0 <= 1e+261) tmp = t_0; else tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(Float64(-x_46_re) / y_46_im)); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e+261], t$95$0, N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[((-x$46$re) / y$46$im), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;t\_0 \leq 10^{+261}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{-x.re}{y.im}\right)\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < 9.9999999999999993e260Initial program 76.8%
if 9.9999999999999993e260 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 21.8%
div-sub15.6%
*-commutative15.6%
add-sqr-sqrt15.6%
times-frac19.9%
fma-neg19.9%
hypot-define19.9%
hypot-define55.0%
associate-/l*67.8%
add-sqr-sqrt67.8%
pow267.8%
hypot-define67.8%
Applied egg-rr67.8%
Taylor expanded in y.im around inf 72.3%
Final simplification75.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im)))))
(if (<= t_0 INFINITY)
t_0
(pow (cbrt (+ (/ x.im y.re) (* x.re (* y.im (pow y.re -2.0))))) 3.0))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
double tmp;
if (t_0 <= ((double) INFINITY)) {
tmp = t_0;
} else {
tmp = pow(cbrt(((x_46_im / y_46_re) + (x_46_re * (y_46_im * pow(y_46_re, -2.0))))), 3.0);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
double tmp;
if (t_0 <= Double.POSITIVE_INFINITY) {
tmp = t_0;
} else {
tmp = Math.pow(Math.cbrt(((x_46_im / y_46_re) + (x_46_re * (y_46_im * Math.pow(y_46_re, -2.0))))), 3.0);
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) tmp = 0.0 if (t_0 <= Inf) tmp = t_0; else tmp = cbrt(Float64(Float64(x_46_im / y_46_re) + Float64(x_46_re * Float64(y_46_im * (y_46_re ^ -2.0))))) ^ 3.0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, Infinity], t$95$0, N[Power[N[Power[N[(N[(x$46$im / y$46$re), $MachinePrecision] + N[(x$46$re * N[(y$46$im * N[Power[y$46$re, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{if}\;t\_0 \leq \infty:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt[3]{\frac{x.im}{y.re} + x.re \cdot \left(y.im \cdot {y.re}^{-2}\right)}\right)}^{3}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) < +inf.0Initial program 75.7%
if +inf.0 < (/.f64 (-.f64 (*.f64 x.im y.re) (*.f64 x.re y.im)) (+.f64 (*.f64 y.re y.re) (*.f64 y.im y.im))) Initial program 0.0%
div-sub0.0%
*-commutative0.0%
add-sqr-sqrt0.0%
times-frac2.1%
fma-neg2.1%
hypot-define2.1%
hypot-define49.9%
associate-/l*58.3%
add-sqr-sqrt58.3%
pow258.3%
hypot-define58.3%
Applied egg-rr58.3%
Taylor expanded in y.im around 0 45.7%
+-commutative45.7%
mul-1-neg45.7%
unsub-neg45.7%
associate-/l*48.9%
Simplified48.9%
add-cube-cbrt47.7%
pow347.7%
cancel-sign-sub-inv47.7%
add-sqr-sqrt27.0%
sqrt-unprod40.5%
sqr-neg40.5%
sqrt-unprod20.6%
add-sqr-sqrt47.1%
pow247.1%
div-inv47.1%
pow247.1%
pow-flip47.1%
metadata-eval47.1%
Applied egg-rr47.1%
Final simplification70.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.im -3.2e+94)
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(/ (- x.re) y.im))
(*
(/ 1.0 (hypot y.re y.im))
(/ (fma x.im y.re (* y.im (- x.re))) (hypot y.re y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= -3.2e+94) {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (-x_46_re / y_46_im));
} else {
tmp = (1.0 / hypot(y_46_re, y_46_im)) * (fma(x_46_im, y_46_re, (y_46_im * -x_46_re)) / hypot(y_46_re, y_46_im));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= -3.2e+94) tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(Float64(-x_46_re) / y_46_im)); else tmp = Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * Float64(fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re))) / hypot(y_46_re, y_46_im))); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, -3.2e+94], N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[((-x$46$re) / y$46$im), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq -3.2 \cdot 10^{+94}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{-x.re}{y.im}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{\mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\end{array}
\end{array}
if y.im < -3.20000000000000014e94Initial program 39.1%
div-sub39.1%
*-commutative39.1%
add-sqr-sqrt39.1%
times-frac42.0%
fma-neg42.0%
hypot-define42.0%
hypot-define52.1%
associate-/l*64.0%
add-sqr-sqrt64.0%
pow264.0%
hypot-define64.0%
Applied egg-rr64.0%
Taylor expanded in y.im around inf 91.9%
if -3.20000000000000014e94 < y.im Initial program 66.4%
*-un-lft-identity66.4%
add-sqr-sqrt66.4%
times-frac66.3%
hypot-define66.3%
fma-neg66.3%
distribute-rgt-neg-in66.3%
hypot-define82.9%
Applied egg-rr82.9%
Final simplification84.1%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.im -1.65e+98)
(fma
(/ y.re (hypot y.re y.im))
(/ x.im (hypot y.re y.im))
(/ (- x.re) y.im))
(/
(* (/ 1.0 (hypot y.re y.im)) (fma x.im y.re (* y.im (- x.re))))
(hypot y.re y.im))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= -1.65e+98) {
tmp = fma((y_46_re / hypot(y_46_re, y_46_im)), (x_46_im / hypot(y_46_re, y_46_im)), (-x_46_re / y_46_im));
} else {
tmp = ((1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, (y_46_im * -x_46_re))) / hypot(y_46_re, y_46_im);
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= -1.65e+98) tmp = fma(Float64(y_46_re / hypot(y_46_re, y_46_im)), Float64(x_46_im / hypot(y_46_re, y_46_im)), Float64(Float64(-x_46_re) / y_46_im)); else tmp = Float64(Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * fma(x_46_im, y_46_re, Float64(y_46_im * Float64(-x_46_re)))) / hypot(y_46_re, y_46_im)); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, -1.65e+98], N[(N[(y$46$re / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] + N[((-x$46$re) / y$46$im), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x$46$im * y$46$re + N[(y$46$im * (-x$46$re)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq -1.65 \cdot 10^{+98}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y.re}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{x.im}{\mathsf{hypot}\left(y.re, y.im\right)}, \frac{-x.re}{y.im}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \mathsf{fma}\left(x.im, y.re, y.im \cdot \left(-x.re\right)\right)}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\end{array}
\end{array}
if y.im < -1.65000000000000014e98Initial program 39.1%
div-sub39.1%
*-commutative39.1%
add-sqr-sqrt39.1%
times-frac42.0%
fma-neg42.0%
hypot-define42.0%
hypot-define52.1%
associate-/l*64.0%
add-sqr-sqrt64.0%
pow264.0%
hypot-define64.0%
Applied egg-rr64.0%
Taylor expanded in y.im around inf 91.9%
if -1.65000000000000014e98 < y.im Initial program 66.4%
*-un-lft-identity66.4%
add-sqr-sqrt66.4%
times-frac66.3%
hypot-define66.3%
fma-neg66.3%
distribute-rgt-neg-in66.3%
hypot-define82.9%
Applied egg-rr82.9%
associate-*r/83.0%
Applied egg-rr83.0%
Final simplification84.2%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 1.35e+154) (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))) (* (/ 1.0 (hypot y.re y.im)) (/ (* y.re x.im) (hypot y.re y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.35e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = (1.0 / hypot(y_46_re, y_46_im)) * ((y_46_re * x_46_im) / hypot(y_46_re, y_46_im));
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.35e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = (1.0 / Math.hypot(y_46_re, y_46_im)) * ((y_46_re * x_46_im) / Math.hypot(y_46_re, y_46_im));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 1.35e+154: tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)) else: tmp = (1.0 / math.hypot(y_46_re, y_46_im)) * ((y_46_re * x_46_im) / math.hypot(y_46_re, y_46_im)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 1.35e+154) tmp = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); else tmp = Float64(Float64(1.0 / hypot(y_46_re, y_46_im)) * Float64(Float64(y_46_re * x_46_im) / hypot(y_46_re, y_46_im))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 1.35e+154) tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); else tmp = (1.0 / hypot(y_46_re, y_46_im)) * ((y_46_re * x_46_im) / hypot(y_46_re, y_46_im)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 1.35e+154], N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(N[(y$46$re * x$46$im), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{y.re \cdot x.im}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\end{array}
\end{array}
if y.re < 1.35000000000000003e154Initial program 68.5%
if 1.35000000000000003e154 < y.re Initial program 35.4%
*-un-lft-identity35.4%
add-sqr-sqrt35.4%
times-frac35.4%
hypot-define35.4%
fma-neg35.4%
distribute-rgt-neg-in35.4%
hypot-define67.4%
Applied egg-rr67.4%
Taylor expanded in x.im around inf 61.5%
*-commutative61.5%
Simplified61.5%
Final simplification67.2%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 1.85e+132) (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))) (* (/ (* y.im x.re) (hypot y.re y.im)) (/ -1.0 (hypot y.re y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.85e+132) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = ((y_46_im * x_46_re) / hypot(y_46_re, y_46_im)) * (-1.0 / hypot(y_46_re, y_46_im));
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.85e+132) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = ((y_46_im * x_46_re) / Math.hypot(y_46_re, y_46_im)) * (-1.0 / Math.hypot(y_46_re, y_46_im));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 1.85e+132: tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)) else: tmp = ((y_46_im * x_46_re) / math.hypot(y_46_re, y_46_im)) * (-1.0 / math.hypot(y_46_re, y_46_im)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 1.85e+132) tmp = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); else tmp = Float64(Float64(Float64(y_46_im * x_46_re) / hypot(y_46_re, y_46_im)) * Float64(-1.0 / hypot(y_46_re, y_46_im))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 1.85e+132) tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); else tmp = ((y_46_im * x_46_re) / hypot(y_46_re, y_46_im)) * (-1.0 / hypot(y_46_re, y_46_im)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 1.85e+132], N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(y$46$im * x$46$re), $MachinePrecision] / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision] * N[(-1.0 / N[Sqrt[y$46$re ^ 2 + y$46$im ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 1.85 \cdot 10^{+132}:\\
\;\;\;\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot x.re}{\mathsf{hypot}\left(y.re, y.im\right)} \cdot \frac{-1}{\mathsf{hypot}\left(y.re, y.im\right)}\\
\end{array}
\end{array}
if y.im < 1.85000000000000005e132Initial program 66.8%
if 1.85000000000000005e132 < y.im Initial program 34.5%
*-un-lft-identity34.5%
add-sqr-sqrt34.5%
times-frac34.5%
hypot-define34.5%
fma-neg34.5%
distribute-rgt-neg-in34.5%
hypot-define70.4%
Applied egg-rr70.4%
Taylor expanded in x.im around 0 59.2%
mul-1-neg59.2%
distribute-rgt-neg-out59.2%
Simplified59.2%
Final simplification65.8%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 1.66e+115) (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))) (pow (cbrt (/ (- x.im (* x.re (/ y.im y.re))) y.re)) 3.0)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.66e+115) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = pow(cbrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)), 3.0);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.66e+115) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = Math.pow(Math.cbrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)), 3.0);
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 1.66e+115) tmp = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); else tmp = cbrt(Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re)) ^ 3.0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 1.66e+115], N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[Power[N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 1.66 \cdot 10^{+115}:\\
\;\;\;\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt[3]{\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}}\right)}^{3}\\
\end{array}
\end{array}
if y.re < 1.66e115Initial program 68.5%
if 1.66e115 < y.re Initial program 40.8%
Taylor expanded in y.re around inf 83.8%
mul-1-neg83.8%
unsub-neg83.8%
associate-/l*89.4%
Simplified89.4%
add-cube-cbrt87.8%
pow387.8%
Applied egg-rr87.8%
Final simplification72.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 1e+154) (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))) (pow (sqrt (/ (- x.im (* x.re (/ y.im y.re))) y.re)) 2.0)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = pow(sqrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)), 2.0);
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 1d+154) then
tmp = ((y_46re * x_46im) - (y_46im * x_46re)) / ((y_46re * y_46re) + (y_46im * y_46im))
else
tmp = sqrt(((x_46im - (x_46re * (y_46im / y_46re))) / y_46re)) ** 2.0d0
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = Math.pow(Math.sqrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)), 2.0);
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 1e+154: tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)) else: tmp = math.pow(math.sqrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)), 2.0) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 1e+154) tmp = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); else tmp = sqrt(Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re)) ^ 2.0; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 1e+154) tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); else tmp = sqrt(((x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re)) ^ 2.0; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 1e+154], N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[Sqrt[N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 10^{+154}:\\
\;\;\;\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt{\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}}\right)}^{2}\\
\end{array}
\end{array}
if y.re < 1.00000000000000004e154Initial program 68.5%
if 1.00000000000000004e154 < y.re Initial program 35.4%
Taylor expanded in y.re around inf 84.8%
mul-1-neg84.8%
unsub-neg84.8%
associate-/l*91.5%
Simplified91.5%
add-sqr-sqrt57.8%
pow257.8%
Applied egg-rr57.8%
Final simplification66.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 2.7e+154) (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))) (/ (log1p (expm1 (- x.im (* x.re (/ y.im y.re))))) y.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 2.7e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = log1p(expm1((x_46_im - (x_46_re * (y_46_im / y_46_re))))) / y_46_re;
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 2.7e+154) {
tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
} else {
tmp = Math.log1p(Math.expm1((x_46_im - (x_46_re * (y_46_im / y_46_re))))) / y_46_re;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 2.7e+154: tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)) else: tmp = math.log1p(math.expm1((x_46_im - (x_46_re * (y_46_im / y_46_re))))) / y_46_re return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 2.7e+154) tmp = Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); else tmp = Float64(log1p(expm1(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))))) / y_46_re); end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 2.7e+154], N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Log[1 + N[(Exp[N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision] / y$46$re), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 2.7 \cdot 10^{+154}:\\
\;\;\;\;\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{log1p}\left(\mathsf{expm1}\left(x.im - x.re \cdot \frac{y.im}{y.re}\right)\right)}{y.re}\\
\end{array}
\end{array}
if y.re < 2.70000000000000006e154Initial program 68.2%
if 2.70000000000000006e154 < y.re Initial program 36.1%
Taylor expanded in y.re around inf 84.5%
mul-1-neg84.5%
unsub-neg84.5%
associate-/l*91.3%
Simplified91.3%
log1p-expm1-u53.5%
Applied egg-rr53.5%
Final simplification65.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 4.6e-55) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (* y.im x.re) (- (* y.re (- y.re)) (* y.im y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 4.6e-55) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im));
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46im <= 4.6d-55) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = (y_46im * x_46re) / ((y_46re * -y_46re) - (y_46im * y_46im))
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 4.6e-55) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 4.6e-55: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 4.6e-55) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(y_46_im * x_46_re) / Float64(Float64(y_46_re * Float64(-y_46_re)) - Float64(y_46_im * y_46_im))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 4.6e-55) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 4.6e-55], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(y$46$im * x$46$re), $MachinePrecision] / N[(N[(y$46$re * (-y$46$re)), $MachinePrecision] - N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 4.6 \cdot 10^{-55}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot x.re}{y.re \cdot \left(-y.re\right) - y.im \cdot y.im}\\
\end{array}
\end{array}
if y.im < 4.60000000000000023e-55Initial program 64.6%
Taylor expanded in y.re around inf 62.8%
mul-1-neg62.8%
unsub-neg62.8%
associate-/l*63.3%
Simplified63.3%
if 4.60000000000000023e-55 < y.im Initial program 58.0%
Taylor expanded in x.im around 0 50.5%
mul-1-neg50.5%
distribute-rgt-neg-out50.5%
Simplified50.5%
Final simplification59.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 2.8e-55) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (* y.im x.re) (- (* y.re (- y.re)) (* y.im y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 2.8e-55) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im));
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46im <= 2.8d-55) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = (y_46im * x_46re) / ((y_46re * -y_46re) - (y_46im * y_46im))
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 2.8e-55) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 2.8e-55: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 2.8e-55) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(y_46_im * x_46_re) / Float64(Float64(y_46_re * Float64(-y_46_re)) - Float64(y_46_im * y_46_im))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 2.8e-55) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = (y_46_im * x_46_re) / ((y_46_re * -y_46_re) - (y_46_im * y_46_im)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 2.8e-55], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(y$46$im * x$46$re), $MachinePrecision] / N[(N[(y$46$re * (-y$46$re)), $MachinePrecision] - N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 2.8 \cdot 10^{-55}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot x.re}{y.re \cdot \left(-y.re\right) - y.im \cdot y.im}\\
\end{array}
\end{array}
if y.im < 2.79999999999999984e-55Initial program 64.6%
Taylor expanded in y.re around inf 62.8%
mul-1-neg62.8%
unsub-neg62.8%
associate-/l*63.3%
Simplified63.3%
if 2.79999999999999984e-55 < y.im Initial program 58.0%
Taylor expanded in y.im around inf 58.0%
Taylor expanded in y.im around inf 50.5%
neg-mul-150.5%
*-commutative50.5%
distribute-rgt-neg-in50.5%
Simplified50.5%
Final simplification59.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 29000000000000.0) (/ (- (/ x.im (/ y.im y.re)) x.re) y.im) (- (/ x.im y.re) (* x.re (/ (/ y.im y.re) y.re)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 29000000000000.0) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (x_46_im / y_46_re) - (x_46_re * ((y_46_im / y_46_re) / y_46_re));
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 29000000000000.0d0) then
tmp = ((x_46im / (y_46im / y_46re)) - x_46re) / y_46im
else
tmp = (x_46im / y_46re) - (x_46re * ((y_46im / y_46re) / y_46re))
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 29000000000000.0) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (x_46_im / y_46_re) - (x_46_re * ((y_46_im / y_46_re) / y_46_re));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 29000000000000.0: tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im else: tmp = (x_46_im / y_46_re) - (x_46_re * ((y_46_im / y_46_re) / y_46_re)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 29000000000000.0) tmp = Float64(Float64(Float64(x_46_im / Float64(y_46_im / y_46_re)) - x_46_re) / y_46_im); else tmp = Float64(Float64(x_46_im / y_46_re) - Float64(x_46_re * Float64(Float64(y_46_im / y_46_re) / y_46_re))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 29000000000000.0) tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im; else tmp = (x_46_im / y_46_re) - (x_46_re * ((y_46_im / y_46_re) / y_46_re)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 29000000000000.0], N[(N[(N[(x$46$im / N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], N[(N[(x$46$im / y$46$re), $MachinePrecision] - N[(x$46$re * N[(N[(y$46$im / y$46$re), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 29000000000000:\\
\;\;\;\;\frac{\frac{x.im}{\frac{y.im}{y.re}} - x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.re} - x.re \cdot \frac{\frac{y.im}{y.re}}{y.re}\\
\end{array}
\end{array}
if y.re < 2.9e13Initial program 65.6%
Taylor expanded in y.re around 0 56.5%
+-commutative56.5%
mul-1-neg56.5%
unsub-neg56.5%
unpow256.5%
associate-/r*60.6%
div-sub62.9%
associate-/l*65.6%
Simplified65.6%
clear-num65.6%
un-div-inv65.6%
Applied egg-rr65.6%
if 2.9e13 < y.re Initial program 55.4%
div-sub55.3%
*-commutative55.3%
add-sqr-sqrt55.3%
times-frac57.2%
fma-neg57.2%
hypot-define57.2%
hypot-define85.2%
associate-/l*86.4%
add-sqr-sqrt86.4%
pow286.4%
hypot-define86.4%
Applied egg-rr86.4%
Taylor expanded in y.im around 0 80.3%
+-commutative80.3%
mul-1-neg80.3%
unsub-neg80.3%
associate-/l*78.9%
Simplified78.9%
*-un-lft-identity78.9%
pow278.9%
times-frac82.6%
Applied egg-rr82.6%
associate-*l/82.7%
*-lft-identity82.7%
Simplified82.7%
Final simplification70.5%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 4.4e+21) (/ (- (/ x.im (/ y.im y.re)) x.re) y.im) (/ (* y.im (- (/ x.im y.im) (/ x.re y.re))) y.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 4.4e+21) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (y_46_im * ((x_46_im / y_46_im) - (x_46_re / y_46_re))) / y_46_re;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 4.4d+21) then
tmp = ((x_46im / (y_46im / y_46re)) - x_46re) / y_46im
else
tmp = (y_46im * ((x_46im / y_46im) - (x_46re / y_46re))) / y_46re
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 4.4e+21) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (y_46_im * ((x_46_im / y_46_im) - (x_46_re / y_46_re))) / y_46_re;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 4.4e+21: tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im else: tmp = (y_46_im * ((x_46_im / y_46_im) - (x_46_re / y_46_re))) / y_46_re return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 4.4e+21) tmp = Float64(Float64(Float64(x_46_im / Float64(y_46_im / y_46_re)) - x_46_re) / y_46_im); else tmp = Float64(Float64(y_46_im * Float64(Float64(x_46_im / y_46_im) - Float64(x_46_re / y_46_re))) / y_46_re); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 4.4e+21) tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im; else tmp = (y_46_im * ((x_46_im / y_46_im) - (x_46_re / y_46_re))) / y_46_re; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 4.4e+21], N[(N[(N[(x$46$im / N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], N[(N[(y$46$im * N[(N[(x$46$im / y$46$im), $MachinePrecision] - N[(x$46$re / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 4.4 \cdot 10^{+21}:\\
\;\;\;\;\frac{\frac{x.im}{\frac{y.im}{y.re}} - x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot \left(\frac{x.im}{y.im} - \frac{x.re}{y.re}\right)}{y.re}\\
\end{array}
\end{array}
if y.re < 4.4e21Initial program 65.8%
Taylor expanded in y.re around 0 56.3%
+-commutative56.3%
mul-1-neg56.3%
unsub-neg56.3%
unpow256.3%
associate-/r*60.3%
div-sub62.6%
associate-/l*65.3%
Simplified65.3%
clear-num65.3%
un-div-inv65.3%
Applied egg-rr65.3%
if 4.4e21 < y.re Initial program 54.8%
Taylor expanded in y.re around inf 84.0%
mul-1-neg84.0%
unsub-neg84.0%
associate-/l*86.1%
Simplified86.1%
Taylor expanded in y.im around inf 77.6%
Final simplification68.8%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 4.2e+172) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (* y.re x.im) (+ (* y.re y.re) (* y.im y.im)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 4.2e+172) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_re * x_46_im) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46im <= 4.2d+172) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = (y_46re * x_46im) / ((y_46re * y_46re) + (y_46im * y_46im))
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 4.2e+172) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (y_46_re * x_46_im) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 4.2e+172: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = (y_46_re * x_46_im) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 4.2e+172) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(y_46_re * x_46_im) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 4.2e+172) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = (y_46_re * x_46_im) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 4.2e+172], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(y$46$re * x$46$im), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 4.2 \cdot 10^{+172}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{y.re \cdot x.im}{y.re \cdot y.re + y.im \cdot y.im}\\
\end{array}
\end{array}
if y.im < 4.2000000000000003e172Initial program 64.8%
Taylor expanded in y.re around inf 59.5%
mul-1-neg59.5%
unsub-neg59.5%
associate-/l*60.8%
Simplified60.8%
if 4.2000000000000003e172 < y.im Initial program 42.0%
Taylor expanded in x.im around inf 42.6%
*-commutative42.6%
Simplified42.6%
Final simplification59.1%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 1.05e+90) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (* (- x.re (* x.im (/ y.re y.im))) (/ -1.0 y.im))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.05e+90) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (x_46_re - (x_46_im * (y_46_re / y_46_im))) * (-1.0 / y_46_im);
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46im <= 1.05d+90) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = (x_46re - (x_46im * (y_46re / y_46im))) * ((-1.0d0) / y_46im)
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.05e+90) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (x_46_re - (x_46_im * (y_46_re / y_46_im))) * (-1.0 / y_46_im);
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 1.05e+90: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = (x_46_re - (x_46_im * (y_46_re / y_46_im))) * (-1.0 / y_46_im) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 1.05e+90) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(x_46_re - Float64(x_46_im * Float64(y_46_re / y_46_im))) * Float64(-1.0 / y_46_im)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 1.05e+90) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = (x_46_re - (x_46_im * (y_46_re / y_46_im))) * (-1.0 / y_46_im); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 1.05e+90], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(x$46$re - N[(x$46$im * N[(y$46$re / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-1.0 / y$46$im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 1.05 \cdot 10^{+90}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\left(x.re - x.im \cdot \frac{y.re}{y.im}\right) \cdot \frac{-1}{y.im}\\
\end{array}
\end{array}
if y.im < 1.0499999999999999e90Initial program 66.2%
Taylor expanded in y.re around inf 61.7%
mul-1-neg61.7%
unsub-neg61.7%
associate-/l*62.6%
Simplified62.6%
if 1.0499999999999999e90 < y.im Initial program 43.6%
Taylor expanded in y.re around 0 67.1%
+-commutative67.1%
mul-1-neg67.1%
unsub-neg67.1%
unpow267.1%
associate-/r*72.1%
div-sub72.1%
associate-/l*79.4%
Simplified79.4%
frac-2neg79.4%
sub-neg79.4%
associate-*r/72.1%
mul-1-neg72.1%
+-commutative72.1%
div-inv72.1%
Applied egg-rr79.4%
Final simplification65.2%
(FPCore (x.re x.im y.re y.im) :precision binary64 (/ (- (* y.re x.im) (* y.im x.re)) (+ (* y.re y.re) (* y.im y.im))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
code = ((y_46re * x_46im) - (y_46im * x_46re)) / ((y_46re * y_46re) + (y_46im * y_46im))
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im));
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): return ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im))
function code(x_46_re, x_46_im, y_46_re, y_46_im) return Float64(Float64(Float64(y_46_re * x_46_im) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_re * y_46_re) + Float64(y_46_im * y_46_im))) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = ((y_46_re * x_46_im) - (y_46_im * x_46_re)) / ((y_46_re * y_46_re) + (y_46_im * y_46_im)); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$re * y$46$re), $MachinePrecision] + N[(y$46$im * y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{y.re \cdot x.im - y.im \cdot x.re}{y.re \cdot y.re + y.im \cdot y.im}
\end{array}
Initial program 62.7%
Final simplification62.7%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 1.4e-47) (/ (- x.re) y.im) (/ (- x.im (* x.re (/ y.im y.re))) y.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.4e-47) {
tmp = -x_46_re / y_46_im;
} else {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 1.4d-47) then
tmp = -x_46re / y_46im
else
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 1.4e-47) {
tmp = -x_46_re / y_46_im;
} else {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 1.4e-47: tmp = -x_46_re / y_46_im else: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 1.4e-47) tmp = Float64(Float64(-x_46_re) / y_46_im); else tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 1.4e-47) tmp = -x_46_re / y_46_im; else tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 1.4e-47], N[((-x$46$re) / y$46$im), $MachinePrecision], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 1.4 \cdot 10^{-47}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\end{array}
\end{array}
if y.re < 1.39999999999999996e-47Initial program 64.5%
Taylor expanded in y.re around 0 51.6%
associate-*r/51.6%
neg-mul-151.6%
Simplified51.6%
if 1.39999999999999996e-47 < y.re Initial program 58.7%
Taylor expanded in y.re around inf 81.8%
mul-1-neg81.8%
unsub-neg81.8%
associate-/l*83.8%
Simplified83.8%
Final simplification61.7%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 8e-41) (/ (- x.re) y.im) (/ (- x.im (/ (* y.im x.re) y.re)) y.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 8e-41) {
tmp = -x_46_re / y_46_im;
} else {
tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 8d-41) then
tmp = -x_46re / y_46im
else
tmp = (x_46im - ((y_46im * x_46re) / y_46re)) / y_46re
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 8e-41) {
tmp = -x_46_re / y_46_im;
} else {
tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 8e-41: tmp = -x_46_re / y_46_im else: tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 8e-41) tmp = Float64(Float64(-x_46_re) / y_46_im); else tmp = Float64(Float64(x_46_im - Float64(Float64(y_46_im * x_46_re) / y_46_re)) / y_46_re); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 8e-41) tmp = -x_46_re / y_46_im; else tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 8e-41], N[((-x$46$re) / y$46$im), $MachinePrecision], N[(N[(x$46$im - N[(N[(y$46$im * x$46$re), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 8 \cdot 10^{-41}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im - \frac{y.im \cdot x.re}{y.re}}{y.re}\\
\end{array}
\end{array}
if y.re < 8.00000000000000005e-41Initial program 64.5%
Taylor expanded in y.re around 0 51.6%
associate-*r/51.6%
neg-mul-151.6%
Simplified51.6%
if 8.00000000000000005e-41 < y.re Initial program 58.7%
div-sub58.6%
*-commutative58.6%
add-sqr-sqrt58.6%
times-frac60.4%
fma-neg60.4%
hypot-define60.4%
hypot-define86.3%
associate-/l*87.4%
add-sqr-sqrt87.4%
pow287.4%
hypot-define87.4%
Applied egg-rr87.4%
Taylor expanded in y.re around inf 81.8%
mul-1-neg81.8%
unsub-neg81.8%
*-commutative81.8%
Simplified81.8%
Final simplification61.1%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re -4.7e-19) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (- (* x.im (/ y.re y.im)) x.re) y.im)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= -4.7e-19) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = ((x_46_im * (y_46_re / y_46_im)) - x_46_re) / y_46_im;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= (-4.7d-19)) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = ((x_46im * (y_46re / y_46im)) - x_46re) / y_46im
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= -4.7e-19) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = ((x_46_im * (y_46_re / y_46_im)) - x_46_re) / y_46_im;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= -4.7e-19: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = ((x_46_im * (y_46_re / y_46_im)) - x_46_re) / y_46_im return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= -4.7e-19) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(Float64(x_46_im * Float64(y_46_re / y_46_im)) - x_46_re) / y_46_im); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= -4.7e-19) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = ((x_46_im * (y_46_re / y_46_im)) - x_46_re) / y_46_im; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, -4.7e-19], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(N[(x$46$im * N[(y$46$re / y$46$im), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -4.7 \cdot 10^{-19}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im \cdot \frac{y.re}{y.im} - x.re}{y.im}\\
\end{array}
\end{array}
if y.re < -4.7e-19Initial program 51.9%
Taylor expanded in y.re around inf 78.2%
mul-1-neg78.2%
unsub-neg78.2%
associate-/l*78.3%
Simplified78.3%
if -4.7e-19 < y.re Initial program 66.0%
Taylor expanded in y.re around 0 53.6%
+-commutative53.6%
mul-1-neg53.6%
unsub-neg53.6%
unpow253.6%
associate-/r*57.2%
div-sub59.3%
associate-/l*60.3%
Simplified60.3%
Final simplification64.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re -3.3e-18) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (- (/ x.im (/ y.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) {
double tmp;
if (y_46_re <= -3.3e-18) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= (-3.3d-18)) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = ((x_46im / (y_46im / y_46re)) - x_46re) / y_46im
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= -3.3e-18) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= -3.3e-18: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = ((x_46_im / (y_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) tmp = 0.0 if (y_46_re <= -3.3e-18) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(Float64(x_46_im / Float64(y_46_im / y_46_re)) - x_46_re) / y_46_im); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= -3.3e-18) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, -3.3e-18], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(N[(x$46$im / N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -3.3 \cdot 10^{-18}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x.im}{\frac{y.im}{y.re}} - x.re}{y.im}\\
\end{array}
\end{array}
if y.re < -3.3000000000000002e-18Initial program 51.9%
Taylor expanded in y.re around inf 78.2%
mul-1-neg78.2%
unsub-neg78.2%
associate-/l*78.3%
Simplified78.3%
if -3.3000000000000002e-18 < y.re Initial program 66.0%
Taylor expanded in y.re around 0 53.6%
+-commutative53.6%
mul-1-neg53.6%
unsub-neg53.6%
unpow253.6%
associate-/r*57.2%
div-sub59.3%
associate-/l*60.3%
Simplified60.3%
clear-num60.3%
un-div-inv60.3%
Applied egg-rr60.3%
Final simplification64.6%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re 200000000000.0) (/ (- (/ x.im (/ y.im y.re)) x.re) y.im) (/ (- x.im (/ (* y.im x.re) y.re)) y.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 200000000000.0) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= 200000000000.0d0) then
tmp = ((x_46im / (y_46im / y_46re)) - x_46re) / y_46im
else
tmp = (x_46im - ((y_46im * x_46re) / y_46re)) / y_46re
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= 200000000000.0) {
tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im;
} else {
tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= 200000000000.0: tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im else: tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= 200000000000.0) tmp = Float64(Float64(Float64(x_46_im / Float64(y_46_im / y_46_re)) - x_46_re) / y_46_im); else tmp = Float64(Float64(x_46_im - Float64(Float64(y_46_im * x_46_re) / y_46_re)) / y_46_re); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= 200000000000.0) tmp = ((x_46_im / (y_46_im / y_46_re)) - x_46_re) / y_46_im; else tmp = (x_46_im - ((y_46_im * x_46_re) / y_46_re)) / y_46_re; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, 200000000000.0], N[(N[(N[(x$46$im / N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], N[(N[(x$46$im - N[(N[(y$46$im * x$46$re), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq 200000000000:\\
\;\;\;\;\frac{\frac{x.im}{\frac{y.im}{y.re}} - x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im - \frac{y.im \cdot x.re}{y.re}}{y.re}\\
\end{array}
\end{array}
if y.re < 2e11Initial program 65.6%
Taylor expanded in y.re around 0 56.5%
+-commutative56.5%
mul-1-neg56.5%
unsub-neg56.5%
unpow256.5%
associate-/r*60.6%
div-sub62.9%
associate-/l*65.6%
Simplified65.6%
clear-num65.6%
un-div-inv65.6%
Applied egg-rr65.6%
if 2e11 < y.re Initial program 55.4%
Taylor expanded in y.re around inf 84.2%
mul-1-neg84.2%
unsub-neg84.2%
unsub-neg84.2%
distribute-frac-neg84.2%
distribute-rgt-neg-out84.2%
Simplified84.2%
Final simplification71.0%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im 1.1e+92) (/ (- x.im (* x.re (/ y.im y.re))) y.re) (/ (- (/ (* y.re x.im) y.im) x.re) y.im)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.1e+92) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (((y_46_re * x_46_im) / y_46_im) - x_46_re) / y_46_im;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46im <= 1.1d+92) then
tmp = (x_46im - (x_46re * (y_46im / y_46re))) / y_46re
else
tmp = (((y_46re * x_46im) / y_46im) - x_46re) / y_46im
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_im <= 1.1e+92) {
tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re;
} else {
tmp = (((y_46_re * x_46_im) / y_46_im) - x_46_re) / y_46_im;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= 1.1e+92: tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re else: tmp = (((y_46_re * x_46_im) / y_46_im) - x_46_re) / y_46_im return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= 1.1e+92) tmp = Float64(Float64(x_46_im - Float64(x_46_re * Float64(y_46_im / y_46_re))) / y_46_re); else tmp = Float64(Float64(Float64(Float64(y_46_re * x_46_im) / y_46_im) - x_46_re) / y_46_im); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_im <= 1.1e+92) tmp = (x_46_im - (x_46_re * (y_46_im / y_46_re))) / y_46_re; else tmp = (((y_46_re * x_46_im) / y_46_im) - x_46_re) / y_46_im; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, 1.1e+92], N[(N[(x$46$im - N[(x$46$re * N[(y$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], N[(N[(N[(N[(y$46$re * x$46$im), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq 1.1 \cdot 10^{+92}:\\
\;\;\;\;\frac{x.im - x.re \cdot \frac{y.im}{y.re}}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{y.re \cdot x.im}{y.im} - x.re}{y.im}\\
\end{array}
\end{array}
if y.im < 1.09999999999999996e92Initial program 66.2%
Taylor expanded in y.re around inf 61.7%
mul-1-neg61.7%
unsub-neg61.7%
associate-/l*62.6%
Simplified62.6%
if 1.09999999999999996e92 < y.im Initial program 43.6%
Taylor expanded in y.im around inf 72.1%
Final simplification64.1%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re -4.6e-18) (/ 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) {
double tmp;
if (y_46_re <= -4.6e-18) {
tmp = x_46_im / y_46_re;
} else {
tmp = -x_46_re / y_46_im;
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: tmp
if (y_46re <= (-4.6d-18)) then
tmp = x_46im / y_46re
else
tmp = -x_46re / y_46im
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (y_46_re <= -4.6e-18) {
tmp = x_46_im / y_46_re;
} else {
tmp = -x_46_re / y_46_im;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_re <= -4.6e-18: tmp = x_46_im / y_46_re else: tmp = -x_46_re / y_46_im return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_re <= -4.6e-18) tmp = Float64(x_46_im / y_46_re); else tmp = Float64(Float64(-x_46_re) / y_46_im); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0; if (y_46_re <= -4.6e-18) tmp = x_46_im / y_46_re; else tmp = -x_46_re / y_46_im; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$re, -4.6e-18], N[(x$46$im / y$46$re), $MachinePrecision], N[((-x$46$re) / y$46$im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -4.6 \cdot 10^{-18}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\end{array}
\end{array}
if y.re < -4.6000000000000002e-18Initial program 51.9%
Taylor expanded in y.re around inf 65.2%
if -4.6000000000000002e-18 < y.re Initial program 66.0%
Taylor expanded in y.re around 0 46.5%
associate-*r/46.5%
neg-mul-146.5%
Simplified46.5%
Final simplification50.9%
(FPCore (x.re x.im y.re y.im) :precision binary64 (/ x.im y.re))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return x_46_im / y_46_re;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
code = x_46im / y_46re
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return x_46_im / y_46_re;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): return x_46_im / y_46_re
function code(x_46_re, x_46_im, y_46_re, y_46_im) return Float64(x_46_im / y_46_re) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = x_46_im / y_46_re; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(x$46$im / y$46$re), $MachinePrecision]
\begin{array}{l}
\\
\frac{x.im}{y.re}
\end{array}
Initial program 62.7%
Taylor expanded in y.re around inf 42.7%
Final simplification42.7%
herbie shell --seed 2024066
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, imaginary part"
:precision binary64
(/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))))