
(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 16 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
(let* ((t_0 (fma y.im y.im (* y.re y.re)))
(t_1 (fma (/ y.re t_0) x.im (* (/ x.re t_0) (- y.im)))))
(if (<= y.re -3.1e+125)
(/ (- x.im (* (/ x.re y.re) y.im)) y.re)
(if (<= y.re -1.1e+108)
(/
(-
(fma (/ y.re y.im) x.im (* (/ x.re y.im) (/ (* y.re y.re) y.im)))
x.re)
y.im)
(if (<= y.re -1.05e-104)
t_1
(if (<= y.re 1.85e-148)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 4.4e+91)
t_1
(fma (/ (- x.re) y.re) (/ y.im y.re) (/ x.im y.re)))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = fma(y_46_im, y_46_im, (y_46_re * y_46_re));
double t_1 = fma((y_46_re / t_0), x_46_im, ((x_46_re / t_0) * -y_46_im));
double tmp;
if (y_46_re <= -3.1e+125) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else if (y_46_re <= -1.1e+108) {
tmp = (fma((y_46_re / y_46_im), x_46_im, ((x_46_re / y_46_im) * ((y_46_re * y_46_re) / y_46_im))) - x_46_re) / y_46_im;
} else if (y_46_re <= -1.05e-104) {
tmp = t_1;
} else if (y_46_re <= 1.85e-148) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 4.4e+91) {
tmp = t_1;
} else {
tmp = fma((-x_46_re / y_46_re), (y_46_im / y_46_re), (x_46_im / y_46_re));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re)) t_1 = fma(Float64(y_46_re / t_0), x_46_im, Float64(Float64(x_46_re / t_0) * Float64(-y_46_im))) tmp = 0.0 if (y_46_re <= -3.1e+125) tmp = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re); elseif (y_46_re <= -1.1e+108) tmp = Float64(Float64(fma(Float64(y_46_re / y_46_im), x_46_im, Float64(Float64(x_46_re / y_46_im) * Float64(Float64(y_46_re * y_46_re) / y_46_im))) - x_46_re) / y_46_im); elseif (y_46_re <= -1.05e-104) tmp = t_1; elseif (y_46_re <= 1.85e-148) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 4.4e+91) tmp = t_1; else tmp = fma(Float64(Float64(-x_46_re) / y_46_re), Float64(y_46_im / y_46_re), Float64(x_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[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(y$46$re / t$95$0), $MachinePrecision] * x$46$im + N[(N[(x$46$re / t$95$0), $MachinePrecision] * (-y$46$im)), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -3.1e+125], N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(N[(y$46$re / y$46$im), $MachinePrecision] * x$46$im + N[(N[(x$46$re / y$46$im), $MachinePrecision] * N[(N[(y$46$re * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -1.05e-104], t$95$1, If[LessEqual[y$46$re, 1.85e-148], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 4.4e+91], t$95$1, N[(N[((-x$46$re) / y$46$re), $MachinePrecision] * N[(y$46$im / y$46$re), $MachinePrecision] + N[(x$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)\\
t_1 := \mathsf{fma}\left(\frac{y.re}{t\_0}, x.im, \frac{x.re}{t\_0} \cdot \left(-y.im\right)\right)\\
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;\frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{y.re}{y.im}, x.im, \frac{x.re}{y.im} \cdot \frac{y.re \cdot y.re}{y.im}\right) - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq -1.05 \cdot 10^{-104}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 1.85 \cdot 10^{-148}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 4.4 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{-x.re}{y.re}, \frac{y.im}{y.re}, \frac{x.im}{y.re}\right)\\
\end{array}
\end{array}
if y.re < -3.1e125Initial program 47.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6447.0
Applied rewrites47.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
frac-subN/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites0.3%
Taylor expanded in y.im around -inf
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -1.04999999999999999e-104 or 1.85000000000000017e-148 < y.re < 4.39999999999999999e91Initial program 77.0%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
sub-negN/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied rewrites83.7%
if -1.04999999999999999e-104 < y.re < 1.85000000000000017e-148Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6494.3
Applied rewrites94.3%
if 4.39999999999999999e91 < y.re Initial program 37.1%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6437.2
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6437.2
Applied rewrites37.2%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6491.0
Applied rewrites91.0%
Final simplification89.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (fma y.im y.im (* y.re y.re)))
(t_1 (fma (/ y.re t_0) x.im (* (/ x.re t_0) (- y.im)))))
(if (<= y.re -3.1e+125)
(/ (- x.im (* (/ x.re y.re) y.im)) y.re)
(if (<= y.re -1.1e+108)
(/ (fma (/ x.im y.im) y.re (- x.re)) y.im)
(if (<= y.re -1.05e-104)
t_1
(if (<= y.re 1.85e-148)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 4.4e+91)
t_1
(fma (/ (- x.re) y.re) (/ y.im y.re) (/ x.im y.re)))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = fma(y_46_im, y_46_im, (y_46_re * y_46_re));
double t_1 = fma((y_46_re / t_0), x_46_im, ((x_46_re / t_0) * -y_46_im));
double tmp;
if (y_46_re <= -3.1e+125) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else if (y_46_re <= -1.1e+108) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else if (y_46_re <= -1.05e-104) {
tmp = t_1;
} else if (y_46_re <= 1.85e-148) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 4.4e+91) {
tmp = t_1;
} else {
tmp = fma((-x_46_re / y_46_re), (y_46_im / y_46_re), (x_46_im / y_46_re));
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re)) t_1 = fma(Float64(y_46_re / t_0), x_46_im, Float64(Float64(x_46_re / t_0) * Float64(-y_46_im))) tmp = 0.0 if (y_46_re <= -3.1e+125) tmp = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re); elseif (y_46_re <= -1.1e+108) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); elseif (y_46_re <= -1.05e-104) tmp = t_1; elseif (y_46_re <= 1.85e-148) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 4.4e+91) tmp = t_1; else tmp = fma(Float64(Float64(-x_46_re) / y_46_re), Float64(y_46_im / y_46_re), Float64(x_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[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(y$46$re / t$95$0), $MachinePrecision] * x$46$im + N[(N[(x$46$re / t$95$0), $MachinePrecision] * (-y$46$im)), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -3.1e+125], N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -1.05e-104], t$95$1, If[LessEqual[y$46$re, 1.85e-148], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 4.4e+91], t$95$1, N[(N[((-x$46$re) / y$46$re), $MachinePrecision] * N[(y$46$im / y$46$re), $MachinePrecision] + N[(x$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)\\
t_1 := \mathsf{fma}\left(\frac{y.re}{t\_0}, x.im, \frac{x.re}{t\_0} \cdot \left(-y.im\right)\right)\\
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;\frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{elif}\;y.re \leq -1.05 \cdot 10^{-104}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 1.85 \cdot 10^{-148}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 4.4 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{-x.re}{y.re}, \frac{y.im}{y.re}, \frac{x.im}{y.re}\right)\\
\end{array}
\end{array}
if y.re < -3.1e125Initial program 47.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6447.0
Applied rewrites47.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6413.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6413.0
Applied rewrites13.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6417.0
Applied rewrites17.0%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6487.4
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -1.04999999999999999e-104 or 1.85000000000000017e-148 < y.re < 4.39999999999999999e91Initial program 77.0%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
sub-negN/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied rewrites83.7%
if -1.04999999999999999e-104 < y.re < 1.85000000000000017e-148Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6494.3
Applied rewrites94.3%
if 4.39999999999999999e91 < y.re Initial program 37.1%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6437.2
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6437.2
Applied rewrites37.2%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6491.0
Applied rewrites91.0%
Final simplification89.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.re -3.1e+125)
(/ (- x.im (* (/ x.re y.re) y.im)) y.re)
(if (<= y.re -1.1e+108)
(/ (fma (/ x.im y.im) y.re (- x.re)) y.im)
(if (<= y.re -6.5e-80)
(*
(/ -1.0 (fma y.im y.im (* y.re y.re)))
(fma (- x.im) y.re (* y.im x.re)))
(if (<= y.re 1.35e-71)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 2.5e+85)
(/ (- (* x.im y.re) (* y.im x.re)) (+ (* y.im y.im) (* y.re y.re)))
(fma (/ (- x.re) y.re) (/ y.im y.re) (/ x.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 <= -3.1e+125) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else if (y_46_re <= -1.1e+108) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else if (y_46_re <= -6.5e-80) {
tmp = (-1.0 / fma(y_46_im, y_46_im, (y_46_re * y_46_re))) * fma(-x_46_im, y_46_re, (y_46_im * x_46_re));
} else if (y_46_re <= 1.35e-71) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 2.5e+85) {
tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / ((y_46_im * y_46_im) + (y_46_re * y_46_re));
} else {
tmp = fma((-x_46_re / y_46_re), (y_46_im / y_46_re), (x_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 (y_46_re <= -3.1e+125) tmp = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re); elseif (y_46_re <= -1.1e+108) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); elseif (y_46_re <= -6.5e-80) tmp = Float64(Float64(-1.0 / fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re))) * fma(Float64(-x_46_im), y_46_re, Float64(y_46_im * x_46_re))); elseif (y_46_re <= 1.35e-71) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 2.5e+85) tmp = Float64(Float64(Float64(x_46_im * y_46_re) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_im * y_46_im) + Float64(y_46_re * y_46_re))); else tmp = fma(Float64(Float64(-x_46_re) / y_46_re), Float64(y_46_im / y_46_re), Float64(x_46_im / 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, -3.1e+125], N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -6.5e-80], N[(N[(-1.0 / N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[((-x$46$im) * y$46$re + N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.35e-71], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 2.5e+85], N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$im * y$46$im), $MachinePrecision] + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[((-x$46$re) / y$46$re), $MachinePrecision] * N[(y$46$im / y$46$re), $MachinePrecision] + N[(x$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;\frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{elif}\;y.re \leq -6.5 \cdot 10^{-80}:\\
\;\;\;\;\frac{-1}{\mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)} \cdot \mathsf{fma}\left(-x.im, y.re, y.im \cdot x.re\right)\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-71}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 2.5 \cdot 10^{+85}:\\
\;\;\;\;\frac{x.im \cdot y.re - y.im \cdot x.re}{y.im \cdot y.im + y.re \cdot y.re}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{-x.re}{y.re}, \frac{y.im}{y.re}, \frac{x.im}{y.re}\right)\\
\end{array}
\end{array}
if y.re < -3.1e125Initial program 47.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6447.0
Applied rewrites47.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6413.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6413.0
Applied rewrites13.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6417.0
Applied rewrites17.0%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6487.4
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -6.49999999999999984e-80Initial program 75.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6475.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6475.0
Applied rewrites75.0%
if -6.49999999999999984e-80 < y.re < 1.3500000000000001e-71Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
if 1.3500000000000001e-71 < y.re < 2.5e85Initial program 84.8%
if 2.5e85 < y.re Initial program 36.4%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6436.4
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6436.4
Applied rewrites36.4%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6489.0
Applied rewrites89.0%
Final simplification87.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.re -3.1e+125)
(/ (- x.im (* (/ x.re y.re) y.im)) y.re)
(if (<= y.re -1.1e+108)
(/ (fma (/ x.im y.im) y.re (- x.re)) y.im)
(if (<= y.re -6.5e-80)
(*
(/ -1.0 (fma y.im y.im (* y.re y.re)))
(fma (- x.im) y.re (* y.im x.re)))
(if (<= y.re 1.35e-71)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 2.5e+85)
(/ (- (* x.im y.re) (* y.im x.re)) (+ (* y.im y.im) (* y.re y.re)))
(fma (- x.re) (/ (/ y.im y.re) y.re) (/ x.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 <= -3.1e+125) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else if (y_46_re <= -1.1e+108) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else if (y_46_re <= -6.5e-80) {
tmp = (-1.0 / fma(y_46_im, y_46_im, (y_46_re * y_46_re))) * fma(-x_46_im, y_46_re, (y_46_im * x_46_re));
} else if (y_46_re <= 1.35e-71) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 2.5e+85) {
tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / ((y_46_im * y_46_im) + (y_46_re * y_46_re));
} else {
tmp = fma(-x_46_re, ((y_46_im / y_46_re) / y_46_re), (x_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 (y_46_re <= -3.1e+125) tmp = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re); elseif (y_46_re <= -1.1e+108) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); elseif (y_46_re <= -6.5e-80) tmp = Float64(Float64(-1.0 / fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re))) * fma(Float64(-x_46_im), y_46_re, Float64(y_46_im * x_46_re))); elseif (y_46_re <= 1.35e-71) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 2.5e+85) tmp = Float64(Float64(Float64(x_46_im * y_46_re) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_im * y_46_im) + Float64(y_46_re * y_46_re))); else tmp = fma(Float64(-x_46_re), Float64(Float64(y_46_im / y_46_re) / y_46_re), Float64(x_46_im / 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, -3.1e+125], N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -6.5e-80], N[(N[(-1.0 / N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[((-x$46$im) * y$46$re + N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.35e-71], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 2.5e+85], N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$im * y$46$im), $MachinePrecision] + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[((-x$46$re) * N[(N[(y$46$im / y$46$re), $MachinePrecision] / y$46$re), $MachinePrecision] + N[(x$46$im / y$46$re), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;\frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{elif}\;y.re \leq -6.5 \cdot 10^{-80}:\\
\;\;\;\;\frac{-1}{\mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)} \cdot \mathsf{fma}\left(-x.im, y.re, y.im \cdot x.re\right)\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-71}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 2.5 \cdot 10^{+85}:\\
\;\;\;\;\frac{x.im \cdot y.re - y.im \cdot x.re}{y.im \cdot y.im + y.re \cdot y.re}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-x.re, \frac{\frac{y.im}{y.re}}{y.re}, \frac{x.im}{y.re}\right)\\
\end{array}
\end{array}
if y.re < -3.1e125Initial program 47.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6447.0
Applied rewrites47.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6495.1
Applied rewrites95.1%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6413.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6413.0
Applied rewrites13.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6417.0
Applied rewrites17.0%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6487.4
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -6.49999999999999984e-80Initial program 75.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6475.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6475.0
Applied rewrites75.0%
if -6.49999999999999984e-80 < y.re < 1.3500000000000001e-71Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
if 1.3500000000000001e-71 < y.re < 2.5e85Initial program 84.8%
if 2.5e85 < y.re Initial program 36.4%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6436.4
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6436.4
Applied rewrites36.4%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6489.0
Applied rewrites89.0%
Applied rewrites89.0%
Final simplification87.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.im (* (/ x.re y.re) y.im)) y.re)))
(if (<= y.re -3.1e+125)
t_0
(if (<= y.re -1.1e+108)
(/ (fma (/ x.im y.im) y.re (- x.re)) y.im)
(if (<= y.re -6.5e-80)
(*
(/ -1.0 (fma y.im y.im (* y.re y.re)))
(fma (- x.im) y.re (* y.im x.re)))
(if (<= y.re 1.35e-71)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 2.5e+85)
(/
(- (* x.im y.re) (* y.im x.re))
(+ (* y.im y.im) (* y.re y.re)))
t_0)))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
double tmp;
if (y_46_re <= -3.1e+125) {
tmp = t_0;
} else if (y_46_re <= -1.1e+108) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else if (y_46_re <= -6.5e-80) {
tmp = (-1.0 / fma(y_46_im, y_46_im, (y_46_re * y_46_re))) * fma(-x_46_im, y_46_re, (y_46_im * x_46_re));
} else if (y_46_re <= 1.35e-71) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 2.5e+85) {
tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / ((y_46_im * y_46_im) + (y_46_re * y_46_re));
} else {
tmp = t_0;
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re) tmp = 0.0 if (y_46_re <= -3.1e+125) tmp = t_0; elseif (y_46_re <= -1.1e+108) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); elseif (y_46_re <= -6.5e-80) tmp = Float64(Float64(-1.0 / fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re))) * fma(Float64(-x_46_im), y_46_re, Float64(y_46_im * x_46_re))); elseif (y_46_re <= 1.35e-71) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 2.5e+85) tmp = Float64(Float64(Float64(x_46_im * y_46_re) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_im * y_46_im) + Float64(y_46_re * y_46_re))); else tmp = t_0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]}, If[LessEqual[y$46$re, -3.1e+125], t$95$0, If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -6.5e-80], N[(N[(-1.0 / N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[((-x$46$im) * y$46$re + N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.35e-71], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 2.5e+85], N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$im * y$46$im), $MachinePrecision] + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{elif}\;y.re \leq -6.5 \cdot 10^{-80}:\\
\;\;\;\;\frac{-1}{\mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)} \cdot \mathsf{fma}\left(-x.im, y.re, y.im \cdot x.re\right)\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-71}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 2.5 \cdot 10^{+85}:\\
\;\;\;\;\frac{x.im \cdot y.re - y.im \cdot x.re}{y.im \cdot y.im + y.re \cdot y.re}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -3.1e125 or 2.5e85 < y.re Initial program 41.4%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6441.5
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6441.5
Applied rewrites41.5%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6491.9
Applied rewrites91.9%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6491.9
Applied rewrites91.9%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6413.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6413.0
Applied rewrites13.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6417.0
Applied rewrites17.0%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6487.4
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -6.49999999999999984e-80Initial program 75.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6475.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6475.0
Applied rewrites75.0%
if -6.49999999999999984e-80 < y.re < 1.3500000000000001e-71Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
if 1.3500000000000001e-71 < y.re < 2.5e85Initial program 84.8%
Final simplification87.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.im (* (/ x.re y.re) y.im)) y.re))
(t_1
(/ (- (* x.im y.re) (* y.im x.re)) (+ (* y.im y.im) (* y.re y.re)))))
(if (<= y.re -3.1e+125)
t_0
(if (<= y.re -1.1e+108)
(/ (fma (/ x.im y.im) y.re (- x.re)) y.im)
(if (<= y.re -6.5e-80)
t_1
(if (<= y.re 1.35e-71)
(/ (- (/ (* x.im y.re) y.im) x.re) y.im)
(if (<= y.re 2.5e+85) t_1 t_0)))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
double t_1 = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / ((y_46_im * y_46_im) + (y_46_re * y_46_re));
double tmp;
if (y_46_re <= -3.1e+125) {
tmp = t_0;
} else if (y_46_re <= -1.1e+108) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else if (y_46_re <= -6.5e-80) {
tmp = t_1;
} else if (y_46_re <= 1.35e-71) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else if (y_46_re <= 2.5e+85) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re) t_1 = Float64(Float64(Float64(x_46_im * y_46_re) - Float64(y_46_im * x_46_re)) / Float64(Float64(y_46_im * y_46_im) + Float64(y_46_re * y_46_re))) tmp = 0.0 if (y_46_re <= -3.1e+125) tmp = t_0; elseif (y_46_re <= -1.1e+108) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); elseif (y_46_re <= -6.5e-80) tmp = t_1; elseif (y_46_re <= 1.35e-71) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); elseif (y_46_re <= 2.5e+85) tmp = t_1; else tmp = t_0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(N[(y$46$im * y$46$im), $MachinePrecision] + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -3.1e+125], t$95$0, If[LessEqual[y$46$re, -1.1e+108], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, -6.5e-80], t$95$1, If[LessEqual[y$46$re, 1.35e-71], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 2.5e+85], t$95$1, t$95$0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
t_1 := \frac{x.im \cdot y.re - y.im \cdot x.re}{y.im \cdot y.im + y.re \cdot y.re}\\
\mathbf{if}\;y.re \leq -3.1 \cdot 10^{+125}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq -1.1 \cdot 10^{+108}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{elif}\;y.re \leq -6.5 \cdot 10^{-80}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-71}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 2.5 \cdot 10^{+85}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -3.1e125 or 2.5e85 < y.re Initial program 41.4%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6441.5
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6441.5
Applied rewrites41.5%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6491.9
Applied rewrites91.9%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6491.9
Applied rewrites91.9%
if -3.1e125 < y.re < -1.1000000000000001e108Initial program 12.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6413.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6413.0
Applied rewrites13.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6417.0
Applied rewrites17.0%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6487.4
Applied rewrites87.4%
if -1.1000000000000001e108 < y.re < -6.49999999999999984e-80 or 1.3500000000000001e-71 < y.re < 2.5e85Initial program 79.8%
if -6.49999999999999984e-80 < y.re < 1.3500000000000001e-71Initial program 66.8%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
Final simplification87.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (fma y.im y.im (* y.re y.re))) (t_1 (/ (- x.re) y.im)))
(if (<= y.im -8.5e+170)
t_1
(if (<= y.im -1.2e-80)
(* (/ y.im t_0) (- x.re))
(if (<= y.im 1.52e-64)
(/ x.im y.re)
(if (<= y.im 3.6e+73) (/ (* (- x.re) y.im) t_0) t_1))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = fma(y_46_im, y_46_im, (y_46_re * y_46_re));
double t_1 = -x_46_re / y_46_im;
double tmp;
if (y_46_im <= -8.5e+170) {
tmp = t_1;
} else if (y_46_im <= -1.2e-80) {
tmp = (y_46_im / t_0) * -x_46_re;
} else if (y_46_im <= 1.52e-64) {
tmp = x_46_im / y_46_re;
} else if (y_46_im <= 3.6e+73) {
tmp = (-x_46_re * y_46_im) / t_0;
} else {
tmp = t_1;
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re)) t_1 = Float64(Float64(-x_46_re) / y_46_im) tmp = 0.0 if (y_46_im <= -8.5e+170) tmp = t_1; elseif (y_46_im <= -1.2e-80) tmp = Float64(Float64(y_46_im / t_0) * Float64(-x_46_re)); elseif (y_46_im <= 1.52e-64) tmp = Float64(x_46_im / y_46_re); elseif (y_46_im <= 3.6e+73) tmp = Float64(Float64(Float64(-x_46_re) * y_46_im) / t_0); else tmp = t_1; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[((-x$46$re) / y$46$im), $MachinePrecision]}, If[LessEqual[y$46$im, -8.5e+170], t$95$1, If[LessEqual[y$46$im, -1.2e-80], N[(N[(y$46$im / t$95$0), $MachinePrecision] * (-x$46$re)), $MachinePrecision], If[LessEqual[y$46$im, 1.52e-64], N[(x$46$im / y$46$re), $MachinePrecision], If[LessEqual[y$46$im, 3.6e+73], N[(N[((-x$46$re) * y$46$im), $MachinePrecision] / t$95$0), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)\\
t_1 := \frac{-x.re}{y.im}\\
\mathbf{if}\;y.im \leq -8.5 \cdot 10^{+170}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq -1.2 \cdot 10^{-80}:\\
\;\;\;\;\frac{y.im}{t\_0} \cdot \left(-x.re\right)\\
\mathbf{elif}\;y.im \leq 1.52 \cdot 10^{-64}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{elif}\;y.im \leq 3.6 \cdot 10^{+73}:\\
\;\;\;\;\frac{\left(-x.re\right) \cdot y.im}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y.im < -8.5000000000000004e170 or 3.5999999999999999e73 < y.im Initial program 33.9%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6468.8
Applied rewrites68.8%
if -8.5000000000000004e170 < y.im < -1.2e-80Initial program 65.9%
Taylor expanded in x.re around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6459.9
Applied rewrites59.9%
if -1.2e-80 < y.im < 1.5200000000000001e-64Initial program 75.7%
Taylor expanded in y.re around inf
lower-/.f6478.6
Applied rewrites78.6%
if 1.5200000000000001e-64 < y.im < 3.5999999999999999e73Initial program 85.0%
Taylor expanded in x.re around inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6463.9
Applied rewrites63.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f6463.9
Applied rewrites63.9%
Final simplification70.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.re) y.im))
(t_1 (* (/ y.im (fma y.im y.im (* y.re y.re))) (- x.re))))
(if (<= y.im -8.5e+170)
t_0
(if (<= y.im -1.2e-80)
t_1
(if (<= y.im 1.52e-64) (/ x.im y.re) (if (<= y.im 3.8e+77) t_1 t_0))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = -x_46_re / y_46_im;
double t_1 = (y_46_im / fma(y_46_im, y_46_im, (y_46_re * y_46_re))) * -x_46_re;
double tmp;
if (y_46_im <= -8.5e+170) {
tmp = t_0;
} else if (y_46_im <= -1.2e-80) {
tmp = t_1;
} else if (y_46_im <= 1.52e-64) {
tmp = x_46_im / y_46_re;
} else if (y_46_im <= 3.8e+77) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(-x_46_re) / y_46_im) t_1 = Float64(Float64(y_46_im / fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re))) * Float64(-x_46_re)) tmp = 0.0 if (y_46_im <= -8.5e+170) tmp = t_0; elseif (y_46_im <= -1.2e-80) tmp = t_1; elseif (y_46_im <= 1.52e-64) tmp = Float64(x_46_im / y_46_re); elseif (y_46_im <= 3.8e+77) tmp = t_1; else tmp = t_0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[((-x$46$re) / y$46$im), $MachinePrecision]}, Block[{t$95$1 = N[(N[(y$46$im / N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * (-x$46$re)), $MachinePrecision]}, If[LessEqual[y$46$im, -8.5e+170], t$95$0, If[LessEqual[y$46$im, -1.2e-80], t$95$1, If[LessEqual[y$46$im, 1.52e-64], N[(x$46$im / y$46$re), $MachinePrecision], If[LessEqual[y$46$im, 3.8e+77], t$95$1, t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-x.re}{y.im}\\
t_1 := \frac{y.im}{\mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)} \cdot \left(-x.re\right)\\
\mathbf{if}\;y.im \leq -8.5 \cdot 10^{+170}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.im \leq -1.2 \cdot 10^{-80}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 1.52 \cdot 10^{-64}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{elif}\;y.im \leq 3.8 \cdot 10^{+77}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.im < -8.5000000000000004e170 or 3.8000000000000001e77 < y.im Initial program 34.7%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6469.2
Applied rewrites69.2%
if -8.5000000000000004e170 < y.im < -1.2e-80 or 1.5200000000000001e-64 < y.im < 3.8000000000000001e77Initial program 71.1%
Taylor expanded in x.re around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.0
Applied rewrites61.0%
if -1.2e-80 < y.im < 1.5200000000000001e-64Initial program 75.7%
Taylor expanded in y.re around inf
lower-/.f6478.6
Applied rewrites78.6%
Final simplification70.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.re -1.1e+46)
(/ x.im y.re)
(if (<= y.re 1.75e-70)
(/ (- x.re) y.im)
(if (<= y.re 4e+83)
(/ (- (* x.im y.re) (* y.im x.re)) (* y.re y.re))
(/ x.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 <= -1.1e+46) {
tmp = x_46_im / y_46_re;
} else if (y_46_re <= 1.75e-70) {
tmp = -x_46_re / y_46_im;
} else if (y_46_re <= 4e+83) {
tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / (y_46_re * y_46_re);
} else {
tmp = x_46_im / 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.1d+46)) then
tmp = x_46im / y_46re
else if (y_46re <= 1.75d-70) then
tmp = -x_46re / y_46im
else if (y_46re <= 4d+83) then
tmp = ((x_46im * y_46re) - (y_46im * x_46re)) / (y_46re * y_46re)
else
tmp = x_46im / 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.1e+46) {
tmp = x_46_im / y_46_re;
} else if (y_46_re <= 1.75e-70) {
tmp = -x_46_re / y_46_im;
} else if (y_46_re <= 4e+83) {
tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / (y_46_re * y_46_re);
} else {
tmp = x_46_im / 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.1e+46: tmp = x_46_im / y_46_re elif y_46_re <= 1.75e-70: tmp = -x_46_re / y_46_im elif y_46_re <= 4e+83: tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / (y_46_re * y_46_re) else: tmp = x_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 (y_46_re <= -1.1e+46) tmp = Float64(x_46_im / y_46_re); elseif (y_46_re <= 1.75e-70) tmp = Float64(Float64(-x_46_re) / y_46_im); elseif (y_46_re <= 4e+83) tmp = Float64(Float64(Float64(x_46_im * y_46_re) - Float64(y_46_im * x_46_re)) / Float64(y_46_re * y_46_re)); else tmp = Float64(x_46_im / 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.1e+46) tmp = x_46_im / y_46_re; elseif (y_46_re <= 1.75e-70) tmp = -x_46_re / y_46_im; elseif (y_46_re <= 4e+83) tmp = ((x_46_im * y_46_re) - (y_46_im * x_46_re)) / (y_46_re * y_46_re); else tmp = x_46_im / 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.1e+46], N[(x$46$im / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, 1.75e-70], N[((-x$46$re) / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 4e+83], N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] - N[(y$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] / N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision], N[(x$46$im / y$46$re), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -1.1 \cdot 10^{+46}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{elif}\;y.re \leq 1.75 \cdot 10^{-70}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 4 \cdot 10^{+83}:\\
\;\;\;\;\frac{x.im \cdot y.re - y.im \cdot x.re}{y.re \cdot y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\end{array}
\end{array}
if y.re < -1.1e46 or 4.00000000000000012e83 < y.re Initial program 45.1%
Taylor expanded in y.re around inf
lower-/.f6475.1
Applied rewrites75.1%
if -1.1e46 < y.re < 1.74999999999999987e-70Initial program 69.0%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6466.9
Applied rewrites66.9%
if 1.74999999999999987e-70 < y.re < 4.00000000000000012e83Initial program 83.6%
Taylor expanded in y.re around inf
unpow2N/A
lower-*.f6468.6
Applied rewrites68.6%
Final simplification70.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= y.re -1.1e+46)
(/ x.im y.re)
(if (<= y.re 1.12e-68)
(/ (- x.re) y.im)
(if (<= y.re 1.5e+86)
(* (/ y.re (fma y.im y.im (* y.re y.re))) x.im)
(/ x.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 <= -1.1e+46) {
tmp = x_46_im / y_46_re;
} else if (y_46_re <= 1.12e-68) {
tmp = -x_46_re / y_46_im;
} else if (y_46_re <= 1.5e+86) {
tmp = (y_46_re / fma(y_46_im, y_46_im, (y_46_re * y_46_re))) * x_46_im;
} else {
tmp = x_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 (y_46_re <= -1.1e+46) tmp = Float64(x_46_im / y_46_re); elseif (y_46_re <= 1.12e-68) tmp = Float64(Float64(-x_46_re) / y_46_im); elseif (y_46_re <= 1.5e+86) tmp = Float64(Float64(y_46_re / fma(y_46_im, y_46_im, Float64(y_46_re * y_46_re))) * x_46_im); else tmp = Float64(x_46_im / 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, -1.1e+46], N[(x$46$im / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, 1.12e-68], N[((-x$46$re) / y$46$im), $MachinePrecision], If[LessEqual[y$46$re, 1.5e+86], N[(N[(y$46$re / N[(y$46$im * y$46$im + N[(y$46$re * y$46$re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * x$46$im), $MachinePrecision], N[(x$46$im / y$46$re), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -1.1 \cdot 10^{+46}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{elif}\;y.re \leq 1.12 \cdot 10^{-68}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{elif}\;y.re \leq 1.5 \cdot 10^{+86}:\\
\;\;\;\;\frac{y.re}{\mathsf{fma}\left(y.im, y.im, y.re \cdot y.re\right)} \cdot x.im\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\end{array}
\end{array}
if y.re < -1.1e46 or 1.49999999999999988e86 < y.re Initial program 43.5%
Taylor expanded in y.re around inf
lower-/.f6475.2
Applied rewrites75.2%
if -1.1e46 < y.re < 1.11999999999999992e-68Initial program 69.3%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6467.2
Applied rewrites67.2%
if 1.11999999999999992e-68 < y.re < 1.49999999999999988e86Initial program 84.4%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
sub-negN/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied rewrites86.9%
Taylor expanded in x.re around 0
*-commutativeN/A
associate-*l/N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Final simplification68.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.im (* (/ x.re y.re) y.im)) y.re)))
(if (<= y.re -2.2e+40)
t_0
(if (<= y.re 7.8e-53) (/ (- (/ (* x.im y.re) y.im) x.re) y.im) t_0))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
double tmp;
if (y_46_re <= -2.2e+40) {
tmp = t_0;
} else if (y_46_re <= 7.8e-53) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else {
tmp = t_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) :: t_0
real(8) :: tmp
t_0 = (x_46im - ((x_46re / y_46re) * y_46im)) / y_46re
if (y_46re <= (-2.2d+40)) then
tmp = t_0
else if (y_46re <= 7.8d-53) then
tmp = (((x_46im * y_46re) / y_46im) - x_46re) / y_46im
else
tmp = t_0
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 t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
double tmp;
if (y_46_re <= -2.2e+40) {
tmp = t_0;
} else if (y_46_re <= 7.8e-53) {
tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im;
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re tmp = 0 if y_46_re <= -2.2e+40: tmp = t_0 elif y_46_re <= 7.8e-53: tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im else: tmp = t_0 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re) tmp = 0.0 if (y_46_re <= -2.2e+40) tmp = t_0; elseif (y_46_re <= 7.8e-53) tmp = Float64(Float64(Float64(Float64(x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im); else tmp = t_0; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re; tmp = 0.0; if (y_46_re <= -2.2e+40) tmp = t_0; elseif (y_46_re <= 7.8e-53) tmp = (((x_46_im * y_46_re) / y_46_im) - x_46_re) / y_46_im; else tmp = t_0; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]}, If[LessEqual[y$46$re, -2.2e+40], t$95$0, If[LessEqual[y$46$re, 7.8e-53], N[(N[(N[(N[(x$46$im * y$46$re), $MachinePrecision] / y$46$im), $MachinePrecision] - x$46$re), $MachinePrecision] / y$46$im), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{if}\;y.re \leq -2.2 \cdot 10^{+40}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 7.8 \cdot 10^{-53}:\\
\;\;\;\;\frac{\frac{x.im \cdot y.re}{y.im} - x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -2.1999999999999999e40 or 7.8000000000000004e-53 < y.re Initial program 53.5%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6453.5
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6453.5
Applied rewrites53.5%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6481.0
Applied rewrites81.0%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6480.9
Applied rewrites80.9%
if -2.1999999999999999e40 < y.re < 7.8000000000000004e-53Initial program 70.1%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-*.f6484.3
Applied rewrites84.3%
Final simplification82.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.im (* (/ x.re y.re) y.im)) y.re)))
(if (<= y.re -2.2e+40)
t_0
(if (<= y.re 7.8e-53) (/ (fma (/ x.im y.im) y.re (- x.re)) y.im) t_0))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
double tmp;
if (y_46_re <= -2.2e+40) {
tmp = t_0;
} else if (y_46_re <= 7.8e-53) {
tmp = fma((x_46_im / y_46_im), y_46_re, -x_46_re) / y_46_im;
} else {
tmp = t_0;
}
return tmp;
}
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re) tmp = 0.0 if (y_46_re <= -2.2e+40) tmp = t_0; elseif (y_46_re <= 7.8e-53) tmp = Float64(fma(Float64(x_46_im / y_46_im), y_46_re, Float64(-x_46_re)) / y_46_im); else tmp = t_0; end return tmp end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]}, If[LessEqual[y$46$re, -2.2e+40], t$95$0, If[LessEqual[y$46$re, 7.8e-53], N[(N[(N[(x$46$im / y$46$im), $MachinePrecision] * y$46$re + (-x$46$re)), $MachinePrecision] / y$46$im), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{if}\;y.re \leq -2.2 \cdot 10^{+40}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 7.8 \cdot 10^{-53}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{x.im}{y.im}, y.re, -x.re\right)}{y.im}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -2.1999999999999999e40 or 7.8000000000000004e-53 < y.re Initial program 53.5%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6453.5
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6453.5
Applied rewrites53.5%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6481.0
Applied rewrites81.0%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6480.9
Applied rewrites80.9%
if -2.1999999999999999e40 < y.re < 7.8000000000000004e-53Initial program 70.1%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6470.0
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6470.0
Applied rewrites70.0%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6424.5
Applied rewrites24.5%
Taylor expanded in y.re around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
unpow2N/A
associate-/r*N/A
div-subN/A
lower-/.f64N/A
sub-negN/A
associate-*l/N/A
mul-1-negN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6482.1
Applied rewrites82.1%
Final simplification81.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (/ (- x.re) y.im)))
(if (<= y.im -4.3e+52)
t_0
(if (<= y.im 1.85e+84) (/ (- x.im (* (/ x.re y.re) y.im)) y.re) t_0))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = -x_46_re / y_46_im;
double tmp;
if (y_46_im <= -4.3e+52) {
tmp = t_0;
} else if (y_46_im <= 1.85e+84) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else {
tmp = t_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) :: t_0
real(8) :: tmp
t_0 = -x_46re / y_46im
if (y_46im <= (-4.3d+52)) then
tmp = t_0
else if (y_46im <= 1.85d+84) then
tmp = (x_46im - ((x_46re / y_46re) * y_46im)) / y_46re
else
tmp = t_0
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 t_0 = -x_46_re / y_46_im;
double tmp;
if (y_46_im <= -4.3e+52) {
tmp = t_0;
} else if (y_46_im <= 1.85e+84) {
tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re;
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = -x_46_re / y_46_im tmp = 0 if y_46_im <= -4.3e+52: tmp = t_0 elif y_46_im <= 1.85e+84: tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re else: tmp = t_0 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(Float64(-x_46_re) / y_46_im) tmp = 0.0 if (y_46_im <= -4.3e+52) tmp = t_0; elseif (y_46_im <= 1.85e+84) tmp = Float64(Float64(x_46_im - Float64(Float64(x_46_re / y_46_re) * y_46_im)) / y_46_re); else tmp = t_0; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = -x_46_re / y_46_im; tmp = 0.0; if (y_46_im <= -4.3e+52) tmp = t_0; elseif (y_46_im <= 1.85e+84) tmp = (x_46_im - ((x_46_re / y_46_re) * y_46_im)) / y_46_re; else tmp = t_0; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[((-x$46$re) / y$46$im), $MachinePrecision]}, If[LessEqual[y$46$im, -4.3e+52], t$95$0, If[LessEqual[y$46$im, 1.85e+84], N[(N[(x$46$im - N[(N[(x$46$re / y$46$re), $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-x.re}{y.im}\\
\mathbf{if}\;y.im \leq -4.3 \cdot 10^{+52}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.im \leq 1.85 \cdot 10^{+84}:\\
\;\;\;\;\frac{x.im - \frac{x.re}{y.re} \cdot y.im}{y.re}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.im < -4.3e52 or 1.85e84 < y.im Initial program 38.7%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6470.1
Applied rewrites70.1%
if -4.3e52 < y.im < 1.85e84Initial program 74.9%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
lift--.f64N/A
sub-negN/A
distribute-neg-inN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
remove-double-negN/A
lower-fma.f64N/A
lower-neg.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6474.9
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6474.9
Applied rewrites74.9%
Taylor expanded in y.im around 0
associate-*r/N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-/.f6475.8
Applied rewrites75.8%
Taylor expanded in y.re around inf
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6477.1
Applied rewrites77.1%
Final simplification74.4%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.im -1.6e+51) (/ (- x.re) y.im) (if (<= y.im 3.6e-63) (/ x.im y.re) (/ -1.0 (/ y.im x.re)))))
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.6e+51) {
tmp = -x_46_re / y_46_im;
} else if (y_46_im <= 3.6e-63) {
tmp = x_46_im / y_46_re;
} else {
tmp = -1.0 / (y_46_im / x_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_46im <= (-1.6d+51)) then
tmp = -x_46re / y_46im
else if (y_46im <= 3.6d-63) then
tmp = x_46im / y_46re
else
tmp = (-1.0d0) / (y_46im / x_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_im <= -1.6e+51) {
tmp = -x_46_re / y_46_im;
} else if (y_46_im <= 3.6e-63) {
tmp = x_46_im / y_46_re;
} else {
tmp = -1.0 / (y_46_im / x_46_re);
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if y_46_im <= -1.6e+51: tmp = -x_46_re / y_46_im elif y_46_im <= 3.6e-63: tmp = x_46_im / y_46_re else: tmp = -1.0 / (y_46_im / x_46_re) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (y_46_im <= -1.6e+51) tmp = Float64(Float64(-x_46_re) / y_46_im); elseif (y_46_im <= 3.6e-63) tmp = Float64(x_46_im / y_46_re); else tmp = Float64(-1.0 / Float64(y_46_im / x_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_im <= -1.6e+51) tmp = -x_46_re / y_46_im; elseif (y_46_im <= 3.6e-63) tmp = x_46_im / y_46_re; else tmp = -1.0 / (y_46_im / x_46_re); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[y$46$im, -1.6e+51], N[((-x$46$re) / y$46$im), $MachinePrecision], If[LessEqual[y$46$im, 3.6e-63], N[(x$46$im / y$46$re), $MachinePrecision], N[(-1.0 / N[(y$46$im / x$46$re), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.im \leq -1.6 \cdot 10^{+51}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{elif}\;y.im \leq 3.6 \cdot 10^{-63}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{\frac{y.im}{x.re}}\\
\end{array}
\end{array}
if y.im < -1.6000000000000001e51Initial program 31.1%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6470.7
Applied rewrites70.7%
if -1.6000000000000001e51 < y.im < 3.60000000000000008e-63Initial program 74.1%
Taylor expanded in y.re around inf
lower-/.f6470.5
Applied rewrites70.5%
if 3.60000000000000008e-63 < y.im Initial program 57.6%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6458.8
Applied rewrites58.8%
Applied rewrites59.5%
Final simplification66.9%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= y.re -1.1e+46) (/ x.im y.re) (if (<= y.re 1.95e-53) (/ (- x.re) y.im) (/ x.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 <= -1.1e+46) {
tmp = x_46_im / y_46_re;
} else if (y_46_re <= 1.95e-53) {
tmp = -x_46_re / y_46_im;
} else {
tmp = x_46_im / 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.1d+46)) then
tmp = x_46im / y_46re
else if (y_46re <= 1.95d-53) then
tmp = -x_46re / y_46im
else
tmp = x_46im / 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.1e+46) {
tmp = x_46_im / y_46_re;
} else if (y_46_re <= 1.95e-53) {
tmp = -x_46_re / y_46_im;
} else {
tmp = x_46_im / 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.1e+46: tmp = x_46_im / y_46_re elif y_46_re <= 1.95e-53: tmp = -x_46_re / y_46_im else: tmp = x_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 (y_46_re <= -1.1e+46) tmp = Float64(x_46_im / y_46_re); elseif (y_46_re <= 1.95e-53) tmp = Float64(Float64(-x_46_re) / y_46_im); else tmp = Float64(x_46_im / 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.1e+46) tmp = x_46_im / y_46_re; elseif (y_46_re <= 1.95e-53) tmp = -x_46_re / y_46_im; else tmp = x_46_im / 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.1e+46], N[(x$46$im / y$46$re), $MachinePrecision], If[LessEqual[y$46$re, 1.95e-53], N[((-x$46$re) / y$46$im), $MachinePrecision], N[(x$46$im / y$46$re), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y.re \leq -1.1 \cdot 10^{+46}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\mathbf{elif}\;y.re \leq 1.95 \cdot 10^{-53}:\\
\;\;\;\;\frac{-x.re}{y.im}\\
\mathbf{else}:\\
\;\;\;\;\frac{x.im}{y.re}\\
\end{array}
\end{array}
if y.re < -1.1e46 or 1.9500000000000001e-53 < y.re Initial program 53.5%
Taylor expanded in y.re around inf
lower-/.f6467.1
Applied rewrites67.1%
if -1.1e46 < y.re < 1.9500000000000001e-53Initial program 70.1%
Taylor expanded in y.re around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6466.4
Applied rewrites66.4%
Final simplification66.7%
(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 61.1%
Taylor expanded in y.re around inf
lower-/.f6445.7
Applied rewrites45.7%
herbie shell --seed 2024304
(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))))