
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (log (sqrt (+ (* x.re x.re) (* x.im x.im))))))
(*
(exp (- (* t_0 y.re) (* (atan2 x.im x.re) y.im)))
(sin (+ (* t_0 y.im) (* (atan2 x.im x.re) y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
return exp(((t_0 * y_46_re) - (atan2(x_46_im, x_46_re) * y_46_im))) * sin(((t_0 * y_46_im) + (atan2(x_46_im, x_46_re) * 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
real(8) :: t_0
t_0 = log(sqrt(((x_46re * x_46re) + (x_46im * x_46im))))
code = exp(((t_0 * y_46re) - (atan2(x_46im, x_46re) * y_46im))) * sin(((t_0 * y_46im) + (atan2(x_46im, x_46re) * y_46re)))
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 = Math.log(Math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
return Math.exp(((t_0 * y_46_re) - (Math.atan2(x_46_im, x_46_re) * y_46_im))) * Math.sin(((t_0 * y_46_im) + (Math.atan2(x_46_im, x_46_re) * y_46_re)));
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.log(math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))) return math.exp(((t_0 * y_46_re) - (math.atan2(x_46_im, x_46_re) * y_46_im))) * math.sin(((t_0 * y_46_im) + (math.atan2(x_46_im, x_46_re) * y_46_re)))
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = log(sqrt(Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)))) return Float64(exp(Float64(Float64(t_0 * y_46_re) - Float64(atan(x_46_im, x_46_re) * y_46_im))) * sin(Float64(Float64(t_0 * y_46_im) + Float64(atan(x_46_im, x_46_re) * y_46_re)))) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))); tmp = exp(((t_0 * y_46_re) - (atan2(x_46_im, x_46_re) * y_46_im))) * sin(((t_0 * y_46_im) + (atan2(x_46_im, x_46_re) * y_46_re))); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Log[N[Sqrt[N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]}, N[(N[Exp[N[(N[(t$95$0 * y$46$re), $MachinePrecision] - N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Sin[N[(N[(t$95$0 * y$46$im), $MachinePrecision] + N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$re), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right)\\
e^{t\_0 \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \sin \left(t\_0 \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (log (sqrt (+ (* x.re x.re) (* x.im x.im))))))
(*
(exp (- (* t_0 y.re) (* (atan2 x.im x.re) y.im)))
(sin (+ (* t_0 y.im) (* (atan2 x.im x.re) y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
return exp(((t_0 * y_46_re) - (atan2(x_46_im, x_46_re) * y_46_im))) * sin(((t_0 * y_46_im) + (atan2(x_46_im, x_46_re) * 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
real(8) :: t_0
t_0 = log(sqrt(((x_46re * x_46re) + (x_46im * x_46im))))
code = exp(((t_0 * y_46re) - (atan2(x_46im, x_46re) * y_46im))) * sin(((t_0 * y_46im) + (atan2(x_46im, x_46re) * y_46re)))
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 = Math.log(Math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
return Math.exp(((t_0 * y_46_re) - (Math.atan2(x_46_im, x_46_re) * y_46_im))) * Math.sin(((t_0 * y_46_im) + (Math.atan2(x_46_im, x_46_re) * y_46_re)));
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.log(math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))) return math.exp(((t_0 * y_46_re) - (math.atan2(x_46_im, x_46_re) * y_46_im))) * math.sin(((t_0 * y_46_im) + (math.atan2(x_46_im, x_46_re) * y_46_re)))
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = log(sqrt(Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)))) return Float64(exp(Float64(Float64(t_0 * y_46_re) - Float64(atan(x_46_im, x_46_re) * y_46_im))) * sin(Float64(Float64(t_0 * y_46_im) + Float64(atan(x_46_im, x_46_re) * y_46_re)))) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))); tmp = exp(((t_0 * y_46_re) - (atan2(x_46_im, x_46_re) * y_46_im))) * sin(((t_0 * y_46_im) + (atan2(x_46_im, x_46_re) * y_46_re))); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Log[N[Sqrt[N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]}, N[(N[Exp[N[(N[(t$95$0 * y$46$re), $MachinePrecision] - N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Sin[N[(N[(t$95$0 * y$46$im), $MachinePrecision] + N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$re), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right)\\
e^{t\_0 \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \sin \left(t\_0 \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)
\end{array}
\end{array}
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* (atan2 x.im x.re) y.im))
(t_1 (log (hypot x.im x.re)))
(t_2 (* y.re (atan2 x.im x.re)))
(t_3
(*
(exp (- (* y.re (log (sqrt (+ (* x.re x.re) (* x.im x.im))))) t_0))
(+ (sin t_2) (* y.im (* t_1 (cos t_2)))))))
(if (<= y.re -3.6e+85)
t_3
(if (<= y.re 7.6)
(/
(sin (* y.im (+ t_1 (* y.re (/ (atan2 x.im x.re) y.im)))))
(/ (exp t_0) (pow (hypot x.re x.im) y.re)))
t_3))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = log(hypot(x_46_im, x_46_re));
double t_2 = y_46_re * atan2(x_46_im, x_46_re);
double t_3 = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * (sin(t_2) + (y_46_im * (t_1 * cos(t_2))));
double tmp;
if (y_46_re <= -3.6e+85) {
tmp = t_3;
} else if (y_46_re <= 7.6) {
tmp = sin((y_46_im * (t_1 + (y_46_re * (atan2(x_46_im, x_46_re) / y_46_im))))) / (exp(t_0) / pow(hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_3;
}
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 = Math.atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = Math.log(Math.hypot(x_46_im, x_46_re));
double t_2 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_3 = Math.exp(((y_46_re * Math.log(Math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * (Math.sin(t_2) + (y_46_im * (t_1 * Math.cos(t_2))));
double tmp;
if (y_46_re <= -3.6e+85) {
tmp = t_3;
} else if (y_46_re <= 7.6) {
tmp = Math.sin((y_46_im * (t_1 + (y_46_re * (Math.atan2(x_46_im, x_46_re) / y_46_im))))) / (Math.exp(t_0) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_3;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.atan2(x_46_im, x_46_re) * y_46_im t_1 = math.log(math.hypot(x_46_im, x_46_re)) t_2 = y_46_re * math.atan2(x_46_im, x_46_re) t_3 = math.exp(((y_46_re * math.log(math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * (math.sin(t_2) + (y_46_im * (t_1 * math.cos(t_2)))) tmp = 0 if y_46_re <= -3.6e+85: tmp = t_3 elif y_46_re <= 7.6: tmp = math.sin((y_46_im * (t_1 + (y_46_re * (math.atan2(x_46_im, x_46_re) / y_46_im))))) / (math.exp(t_0) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) else: tmp = t_3 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(atan(x_46_im, x_46_re) * y_46_im) t_1 = log(hypot(x_46_im, x_46_re)) t_2 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_3 = Float64(exp(Float64(Float64(y_46_re * log(sqrt(Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im))))) - t_0)) * Float64(sin(t_2) + Float64(y_46_im * Float64(t_1 * cos(t_2))))) tmp = 0.0 if (y_46_re <= -3.6e+85) tmp = t_3; elseif (y_46_re <= 7.6) tmp = Float64(sin(Float64(y_46_im * Float64(t_1 + Float64(y_46_re * Float64(atan(x_46_im, x_46_re) / y_46_im))))) / Float64(exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re))); else tmp = t_3; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = atan2(x_46_im, x_46_re) * y_46_im; t_1 = log(hypot(x_46_im, x_46_re)); t_2 = y_46_re * atan2(x_46_im, x_46_re); t_3 = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * (sin(t_2) + (y_46_im * (t_1 * cos(t_2)))); tmp = 0.0; if (y_46_re <= -3.6e+85) tmp = t_3; elseif (y_46_re <= 7.6) tmp = sin((y_46_im * (t_1 + (y_46_re * (atan2(x_46_im, x_46_re) / y_46_im))))) / (exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re)); else tmp = t_3; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]}, Block[{t$95$1 = N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Exp[N[(N[(y$46$re * N[Log[N[Sqrt[N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision] * N[(N[Sin[t$95$2], $MachinePrecision] + N[(y$46$im * N[(t$95$1 * N[Cos[t$95$2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -3.6e+85], t$95$3, If[LessEqual[y$46$re, 7.6], N[(N[Sin[N[(y$46$im * N[(t$95$1 + N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[Exp[t$95$0], $MachinePrecision] / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$3]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im\\
t_1 := \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
t_2 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_3 := e^{y.re \cdot \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) - t\_0} \cdot \left(\sin t\_2 + y.im \cdot \left(t\_1 \cdot \cos t\_2\right)\right)\\
\mathbf{if}\;y.re \leq -3.6 \cdot 10^{+85}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y.re \leq 7.6:\\
\;\;\;\;\frac{\sin \left(y.im \cdot \left(t\_1 + y.re \cdot \frac{\tan^{-1}_* \frac{x.im}{x.re}}{y.im}\right)\right)}{\frac{e^{t\_0}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if y.re < -3.5999999999999998e85 or 7.5999999999999996 < y.re Initial program 43.6%
Taylor expanded in y.im around 0
+-lowering-+.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6474.6%
Simplified74.6%
if -3.5999999999999998e85 < y.re < 7.5999999999999996Initial program 34.3%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.0%
Taylor expanded in y.im around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6479.6%
Simplified79.6%
Final simplification77.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (pow (hypot x.im x.re) y.re)) (t_1 (log (hypot x.im x.re))))
(if (<= y.re -3e+86)
(* t_0 (sin (/ 1.0 (/ (/ 1.0 (atan2 x.im x.re)) y.re))))
(if (<= y.re 2300000000.0)
(/
(sin (* y.im (+ t_1 (* y.re (/ (atan2 x.im x.re) y.im)))))
(/ (exp (* (atan2 x.im x.re) y.im)) (pow (hypot x.re x.im) y.re)))
(* t_0 (* y.im t_1))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(hypot(x_46_im, x_46_re), y_46_re);
double t_1 = log(hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -3e+86) {
tmp = t_0 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 2300000000.0) {
tmp = sin((y_46_im * (t_1 + (y_46_re * (atan2(x_46_im, x_46_re) / y_46_im))))) / (exp((atan2(x_46_im, x_46_re) * y_46_im)) / pow(hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_0 * (y_46_im * t_1);
}
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 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double t_1 = Math.log(Math.hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -3e+86) {
tmp = t_0 * Math.sin((1.0 / ((1.0 / Math.atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 2300000000.0) {
tmp = Math.sin((y_46_im * (t_1 + (y_46_re * (Math.atan2(x_46_im, x_46_re) / y_46_im))))) / (Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im)) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_0 * (y_46_im * t_1);
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) t_1 = math.log(math.hypot(x_46_im, x_46_re)) tmp = 0 if y_46_re <= -3e+86: tmp = t_0 * math.sin((1.0 / ((1.0 / math.atan2(x_46_im, x_46_re)) / y_46_re))) elif y_46_re <= 2300000000.0: tmp = math.sin((y_46_im * (t_1 + (y_46_re * (math.atan2(x_46_im, x_46_re) / y_46_im))))) / (math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) else: tmp = t_0 * (y_46_im * t_1) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re t_1 = log(hypot(x_46_im, x_46_re)) tmp = 0.0 if (y_46_re <= -3e+86) tmp = Float64(t_0 * sin(Float64(1.0 / Float64(Float64(1.0 / atan(x_46_im, x_46_re)) / y_46_re)))); elseif (y_46_re <= 2300000000.0) tmp = Float64(sin(Float64(y_46_im * Float64(t_1 + Float64(y_46_re * Float64(atan(x_46_im, x_46_re) / y_46_im))))) / Float64(exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re))); else tmp = Float64(t_0 * Float64(y_46_im * t_1)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re; t_1 = log(hypot(x_46_im, x_46_re)); tmp = 0.0; if (y_46_re <= -3e+86) tmp = t_0 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re))); elseif (y_46_re <= 2300000000.0) tmp = sin((y_46_im * (t_1 + (y_46_re * (atan2(x_46_im, x_46_re) / y_46_im))))) / (exp((atan2(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re)); else tmp = t_0 * (y_46_im * t_1); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$1 = N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -3e+86], N[(t$95$0 * N[Sin[N[(1.0 / N[(N[(1.0 / N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 2300000000.0], N[(N[Sin[N[(y$46$im * N[(t$95$1 + N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] / y$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision] / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[(y$46$im * t$95$1), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
t_1 := \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{if}\;y.re \leq -3 \cdot 10^{+86}:\\
\;\;\;\;t\_0 \cdot \sin \left(\frac{1}{\frac{\frac{1}{\tan^{-1}_* \frac{x.im}{x.re}}}{y.re}}\right)\\
\mathbf{elif}\;y.re \leq 2300000000:\\
\;\;\;\;\frac{\sin \left(y.im \cdot \left(t\_1 + y.re \cdot \frac{\tan^{-1}_* \frac{x.im}{x.re}}{y.im}\right)\right)}{\frac{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \left(y.im \cdot t\_1\right)\\
\end{array}
\end{array}
if y.re < -2.99999999999999977e86Initial program 48.2%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.7%
Simplified78.7%
remove-double-divN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6480.5%
Applied egg-rr80.5%
if -2.99999999999999977e86 < y.re < 2.3e9Initial program 34.1%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.1%
Taylor expanded in y.im around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6479.8%
Simplified79.8%
if 2.3e9 < y.re Initial program 39.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.7%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.8%
Simplified52.8%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.8%
Applied egg-rr52.8%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.3%
Simplified64.3%
Final simplification76.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (pow (hypot x.im x.re) y.re)))
(if (<= y.re -27.0)
(* t_0 (sin (/ 1.0 (/ (/ 1.0 (atan2 x.im x.re)) y.re))))
(if (<= y.re 4500000000.0)
(/
(sin (+ (* y.re (atan2 x.im x.re)) (* y.im (log (hypot x.re x.im)))))
(/ (exp (* (atan2 x.im x.re) y.im)) (pow (hypot x.re x.im) y.re)))
(* t_0 (* y.im (log (hypot x.im x.re))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(hypot(x_46_im, x_46_re), y_46_re);
double tmp;
if (y_46_re <= -27.0) {
tmp = t_0 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 4500000000.0) {
tmp = sin(((y_46_re * atan2(x_46_im, x_46_re)) + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp((atan2(x_46_im, x_46_re) * y_46_im)) / pow(hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_0 * (y_46_im * log(hypot(x_46_im, x_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 t_0 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double tmp;
if (y_46_re <= -27.0) {
tmp = t_0 * Math.sin((1.0 / ((1.0 / Math.atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 4500000000.0) {
tmp = Math.sin(((y_46_re * Math.atan2(x_46_im, x_46_re)) + (y_46_im * Math.log(Math.hypot(x_46_re, x_46_im))))) / (Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im)) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_0 * (y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) tmp = 0 if y_46_re <= -27.0: tmp = t_0 * math.sin((1.0 / ((1.0 / math.atan2(x_46_im, x_46_re)) / y_46_re))) elif y_46_re <= 4500000000.0: tmp = math.sin(((y_46_re * math.atan2(x_46_im, x_46_re)) + (y_46_im * math.log(math.hypot(x_46_re, x_46_im))))) / (math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) else: tmp = t_0 * (y_46_im * math.log(math.hypot(x_46_im, x_46_re))) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re tmp = 0.0 if (y_46_re <= -27.0) tmp = Float64(t_0 * sin(Float64(1.0 / Float64(Float64(1.0 / atan(x_46_im, x_46_re)) / y_46_re)))); elseif (y_46_re <= 4500000000.0) tmp = Float64(sin(Float64(Float64(y_46_re * atan(x_46_im, x_46_re)) + Float64(y_46_im * log(hypot(x_46_re, x_46_im))))) / Float64(exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re))); else tmp = Float64(t_0 * Float64(y_46_im * log(hypot(x_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) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re; tmp = 0.0; if (y_46_re <= -27.0) tmp = t_0 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re))); elseif (y_46_re <= 4500000000.0) tmp = sin(((y_46_re * atan2(x_46_im, x_46_re)) + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp((atan2(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re)); else tmp = t_0 * (y_46_im * log(hypot(x_46_im, x_46_re))); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, If[LessEqual[y$46$re, -27.0], N[(t$95$0 * N[Sin[N[(1.0 / N[(N[(1.0 / N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 4500000000.0], N[(N[Sin[N[(N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] + N[(y$46$im * N[Log[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision] / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
\mathbf{if}\;y.re \leq -27:\\
\;\;\;\;t\_0 \cdot \sin \left(\frac{1}{\frac{\frac{1}{\tan^{-1}_* \frac{x.im}{x.re}}}{y.re}}\right)\\
\mathbf{elif}\;y.re \leq 4500000000:\\
\;\;\;\;\frac{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + y.im \cdot \log \left(\mathsf{hypot}\left(x.re, x.im\right)\right)\right)}{\frac{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
\end{array}
\end{array}
if y.re < -27Initial program 42.0%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.4%
Simplified78.4%
remove-double-divN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6479.8%
Applied egg-rr79.8%
if -27 < y.re < 4.5e9Initial program 35.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.3%
if 4.5e9 < y.re Initial program 39.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.7%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.8%
Simplified52.8%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.8%
Applied egg-rr52.8%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.3%
Simplified64.3%
Final simplification76.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (sin (* y.im (log (hypot x.im x.re)))))
(t_1
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im))))))
(t_2 (pow (hypot x.re x.im) y.re))
(t_3 (/ (exp (* (atan2 x.im x.re) y.im)) t_2)))
(if (<= y.im -1.2e+206)
t_1
(if (<= y.im -3.1e-172)
(/ t_0 t_3)
(if (<= y.im 3.3e-122)
(/ (sin (* y.re (atan2 x.im x.re))) t_3)
(if (<= y.im 1.05e+100)
(/ t_0 (/ (exp (/ (atan2 x.im x.re) (/ 1.0 y.im))) t_2))
t_1))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = sin((y_46_im * log(hypot(x_46_im, x_46_re))));
double t_1 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double t_2 = pow(hypot(x_46_re, x_46_im), y_46_re);
double t_3 = exp((atan2(x_46_im, x_46_re) * y_46_im)) / t_2;
double tmp;
if (y_46_im <= -1.2e+206) {
tmp = t_1;
} else if (y_46_im <= -3.1e-172) {
tmp = t_0 / t_3;
} else if (y_46_im <= 3.3e-122) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / t_3;
} else if (y_46_im <= 1.05e+100) {
tmp = t_0 / (exp((atan2(x_46_im, x_46_re) / (1.0 / y_46_im))) / t_2);
} else {
tmp = t_1;
}
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 = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))));
double t_1 = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double t_2 = Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re);
double t_3 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im)) / t_2;
double tmp;
if (y_46_im <= -1.2e+206) {
tmp = t_1;
} else if (y_46_im <= -3.1e-172) {
tmp = t_0 / t_3;
} else if (y_46_im <= 3.3e-122) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) / t_3;
} else if (y_46_im <= 1.05e+100) {
tmp = t_0 / (Math.exp((Math.atan2(x_46_im, x_46_re) / (1.0 / y_46_im))) / t_2);
} else {
tmp = t_1;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) t_1 = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) t_2 = math.pow(math.hypot(x_46_re, x_46_im), y_46_re) t_3 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) / t_2 tmp = 0 if y_46_im <= -1.2e+206: tmp = t_1 elif y_46_im <= -3.1e-172: tmp = t_0 / t_3 elif y_46_im <= 3.3e-122: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) / t_3 elif y_46_im <= 1.05e+100: tmp = t_0 / (math.exp((math.atan2(x_46_im, x_46_re) / (1.0 / y_46_im))) / t_2) else: tmp = t_1 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) t_1 = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))) t_2 = hypot(x_46_re, x_46_im) ^ y_46_re t_3 = Float64(exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) / t_2) tmp = 0.0 if (y_46_im <= -1.2e+206) tmp = t_1; elseif (y_46_im <= -3.1e-172) tmp = Float64(t_0 / t_3); elseif (y_46_im <= 3.3e-122) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) / t_3); elseif (y_46_im <= 1.05e+100) tmp = Float64(t_0 / Float64(exp(Float64(atan(x_46_im, x_46_re) / Float64(1.0 / y_46_im))) / t_2)); else tmp = t_1; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = sin((y_46_im * log(hypot(x_46_im, x_46_re)))); t_1 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); t_2 = hypot(x_46_re, x_46_im) ^ y_46_re; t_3 = exp((atan2(x_46_im, x_46_re) * y_46_im)) / t_2; tmp = 0.0; if (y_46_im <= -1.2e+206) tmp = t_1; elseif (y_46_im <= -3.1e-172) tmp = t_0 / t_3; elseif (y_46_im <= 3.3e-122) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / t_3; elseif (y_46_im <= 1.05e+100) tmp = t_0 / (exp((atan2(x_46_im, x_46_re) / (1.0 / y_46_im))) / t_2); else tmp = t_1; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Sin[N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$3 = N[(N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision] / t$95$2), $MachinePrecision]}, If[LessEqual[y$46$im, -1.2e+206], t$95$1, If[LessEqual[y$46$im, -3.1e-172], N[(t$95$0 / t$95$3), $MachinePrecision], If[LessEqual[y$46$im, 3.3e-122], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t$95$3), $MachinePrecision], If[LessEqual[y$46$im, 1.05e+100], N[(t$95$0 / N[(N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] / N[(1.0 / y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
t_1 := y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
t_2 := {\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}\\
t_3 := \frac{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}{t\_2}\\
\mathbf{if}\;y.im \leq -1.2 \cdot 10^{+206}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq -3.1 \cdot 10^{-172}:\\
\;\;\;\;\frac{t\_0}{t\_3}\\
\mathbf{elif}\;y.im \leq 3.3 \cdot 10^{-122}:\\
\;\;\;\;\frac{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)}{t\_3}\\
\mathbf{elif}\;y.im \leq 1.05 \cdot 10^{+100}:\\
\;\;\;\;\frac{t\_0}{\frac{e^{\frac{\tan^{-1}_* \frac{x.im}{x.re}}{\frac{1}{y.im}}}}{t\_2}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y.im < -1.2e206 or 1.0499999999999999e100 < y.im Initial program 27.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6452.6%
Simplified52.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6469.3%
Simplified69.3%
if -1.2e206 < y.im < -3.1000000000000003e-172Initial program 39.4%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified71.2%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6467.0%
Simplified67.0%
if -3.1000000000000003e-172 < y.im < 3.29999999999999999e-122Initial program 50.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified93.1%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
if 3.29999999999999999e-122 < y.im < 1.0499999999999999e100Initial program 37.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified69.7%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6477.9%
Simplified77.9%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6477.9%
Applied egg-rr77.9%
Final simplification73.3%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im))))))
(t_1 (/ (exp (* (atan2 x.im x.re) y.im)) (pow (hypot x.re x.im) y.re)))
(t_2 (/ (sin (* y.im (log (hypot x.im x.re)))) t_1)))
(if (<= y.im -3.8e+201)
t_0
(if (<= y.im -7.5e-175)
t_2
(if (<= y.im 3.2e-122)
(/ (sin (* y.re (atan2 x.im x.re))) t_1)
(if (<= y.im 2.05e+100) t_2 t_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 * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double t_1 = exp((atan2(x_46_im, x_46_re) * y_46_im)) / pow(hypot(x_46_re, x_46_im), y_46_re);
double t_2 = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1;
double tmp;
if (y_46_im <= -3.8e+201) {
tmp = t_0;
} else if (y_46_im <= -7.5e-175) {
tmp = t_2;
} else if (y_46_im <= 3.2e-122) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / t_1;
} else if (y_46_im <= 2.05e+100) {
tmp = t_2;
} else {
tmp = t_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 * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double t_1 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im)) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re);
double t_2 = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)))) / t_1;
double tmp;
if (y_46_im <= -3.8e+201) {
tmp = t_0;
} else if (y_46_im <= -7.5e-175) {
tmp = t_2;
} else if (y_46_im <= 3.2e-122) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) / t_1;
} else if (y_46_im <= 2.05e+100) {
tmp = t_2;
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) t_1 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re) t_2 = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) / t_1 tmp = 0 if y_46_im <= -3.8e+201: tmp = t_0 elif y_46_im <= -7.5e-175: tmp = t_2 elif y_46_im <= 3.2e-122: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) / t_1 elif y_46_im <= 2.05e+100: tmp = t_2 else: tmp = t_0 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))) t_1 = Float64(exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re)) t_2 = Float64(sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1) tmp = 0.0 if (y_46_im <= -3.8e+201) tmp = t_0; elseif (y_46_im <= -7.5e-175) tmp = t_2; elseif (y_46_im <= 3.2e-122) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) / t_1); elseif (y_46_im <= 2.05e+100) tmp = t_2; 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 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); t_1 = exp((atan2(x_46_im, x_46_re) * y_46_im)) / (hypot(x_46_re, x_46_im) ^ y_46_re); t_2 = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1; tmp = 0.0; if (y_46_im <= -3.8e+201) tmp = t_0; elseif (y_46_im <= -7.5e-175) tmp = t_2; elseif (y_46_im <= 3.2e-122) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / t_1; elseif (y_46_im <= 2.05e+100) tmp = t_2; 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[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision] / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Sin[N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t$95$1), $MachinePrecision]}, If[LessEqual[y$46$im, -3.8e+201], t$95$0, If[LessEqual[y$46$im, -7.5e-175], t$95$2, If[LessEqual[y$46$im, 3.2e-122], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t$95$1), $MachinePrecision], If[LessEqual[y$46$im, 2.05e+100], t$95$2, t$95$0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
t_1 := \frac{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}\\
t_2 := \frac{\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)}{t\_1}\\
\mathbf{if}\;y.im \leq -3.8 \cdot 10^{+201}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.im \leq -7.5 \cdot 10^{-175}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{-122}:\\
\;\;\;\;\frac{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)}{t\_1}\\
\mathbf{elif}\;y.im \leq 2.05 \cdot 10^{+100}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.im < -3.79999999999999995e201 or 2.0500000000000001e100 < y.im Initial program 27.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6452.6%
Simplified52.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6469.3%
Simplified69.3%
if -3.79999999999999995e201 < y.im < -7.50000000000000053e-175 or 3.2000000000000002e-122 < y.im < 2.0500000000000001e100Initial program 38.7%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified70.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6471.3%
Simplified71.3%
if -7.50000000000000053e-175 < y.im < 3.2000000000000002e-122Initial program 50.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified93.1%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
Final simplification73.3%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.im (log (hypot x.im x.re))))
(t_1 (* (pow (hypot x.im x.re) y.re) t_0))
(t_2 (exp (* (atan2 x.im x.re) y.im))))
(if (<= y.im -7.6e-6)
(/ (sin t_0) t_2)
(if (<= y.im -1.8e-172)
t_1
(if (<= y.im 3.2e-122)
(/
(sin (* y.re (atan2 x.im x.re)))
(/ t_2 (pow (hypot x.re x.im) y.re)))
(if (<= y.im 1.6e+15)
t_1
(*
y.re
(*
(atan2 x.im x.re)
(exp (* (atan2 x.im x.re) (- 0.0 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_im * log(hypot(x_46_im, x_46_re));
double t_1 = pow(hypot(x_46_im, x_46_re), y_46_re) * t_0;
double t_2 = exp((atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_im <= -7.6e-6) {
tmp = sin(t_0) / t_2;
} else if (y_46_im <= -1.8e-172) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / (t_2 / pow(hypot(x_46_re, x_46_im), y_46_re));
} else if (y_46_im <= 1.6e+15) {
tmp = t_1;
} else {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - 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 t_0 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double t_1 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re) * t_0;
double t_2 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_im <= -7.6e-6) {
tmp = Math.sin(t_0) / t_2;
} else if (y_46_im <= -1.8e-172) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) / (t_2 / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else if (y_46_im <= 1.6e+15) {
tmp = t_1;
} else {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) t_1 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) * t_0 t_2 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) tmp = 0 if y_46_im <= -7.6e-6: tmp = math.sin(t_0) / t_2 elif y_46_im <= -1.8e-172: tmp = t_1 elif y_46_im <= 3.2e-122: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) / (t_2 / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) elif y_46_im <= 1.6e+15: tmp = t_1 else: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) t_1 = Float64((hypot(x_46_im, x_46_re) ^ y_46_re) * t_0) t_2 = exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) tmp = 0.0 if (y_46_im <= -7.6e-6) tmp = Float64(sin(t_0) / t_2); elseif (y_46_im <= -1.8e-172) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) / Float64(t_2 / (hypot(x_46_re, x_46_im) ^ y_46_re))); elseif (y_46_im <= 1.6e+15) tmp = t_1; else tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_im * log(hypot(x_46_im, x_46_re)); t_1 = (hypot(x_46_im, x_46_re) ^ y_46_re) * t_0; t_2 = exp((atan2(x_46_im, x_46_re) * y_46_im)); tmp = 0.0; if (y_46_im <= -7.6e-6) tmp = sin(t_0) / t_2; elseif (y_46_im <= -1.8e-172) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) / (t_2 / (hypot(x_46_re, x_46_im) ^ y_46_re)); elseif (y_46_im <= 1.6e+15) tmp = t_1; else tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision] * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$im, -7.6e-6], N[(N[Sin[t$95$0], $MachinePrecision] / t$95$2), $MachinePrecision], If[LessEqual[y$46$im, -1.8e-172], t$95$1, If[LessEqual[y$46$im, 3.2e-122], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(t$95$2 / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, 1.6e+15], t$95$1, N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
t_1 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re} \cdot t\_0\\
t_2 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}\\
\mathbf{if}\;y.im \leq -7.6 \cdot 10^{-6}:\\
\;\;\;\;\frac{\sin t\_0}{t\_2}\\
\mathbf{elif}\;y.im \leq -1.8 \cdot 10^{-172}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{-122}:\\
\;\;\;\;\frac{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)}{\frac{t\_2}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{elif}\;y.im \leq 1.6 \cdot 10^{+15}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\end{array}
\end{array}
if y.im < -7.6000000000000001e-6Initial program 28.5%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified55.9%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6462.5%
Simplified62.5%
if -7.6000000000000001e-6 < y.im < -1.80000000000000007e-172 or 3.2000000000000002e-122 < y.im < 1.6e15Initial program 50.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified82.2%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6482.0%
Simplified82.0%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6482.0%
Applied egg-rr82.0%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6480.2%
Simplified80.2%
if -1.80000000000000007e-172 < y.im < 3.2000000000000002e-122Initial program 50.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified93.1%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
if 1.6e15 < y.im Initial program 22.7%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6450.7%
Simplified50.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6461.6%
Simplified61.6%
Final simplification71.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (pow (hypot x.im x.re) y.re))
(t_1 (* y.im (log (hypot x.im x.re))))
(t_2 (* t_0 t_1)))
(if (<= y.im -1.85e-5)
(/ (sin t_1) (exp (* (atan2 x.im x.re) y.im)))
(if (<= y.im -2.7e-170)
t_2
(if (<= y.im 3.2e-122)
(* (sin (* y.re (atan2 x.im x.re))) t_0)
(if (<= y.im 1.45e+15)
t_2
(*
y.re
(*
(atan2 x.im x.re)
(exp (* (atan2 x.im x.re) (- 0.0 y.im)))))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(hypot(x_46_im, x_46_re), y_46_re);
double t_1 = y_46_im * log(hypot(x_46_im, x_46_re));
double t_2 = t_0 * t_1;
double tmp;
if (y_46_im <= -1.85e-5) {
tmp = sin(t_1) / exp((atan2(x_46_im, x_46_re) * y_46_im));
} else if (y_46_im <= -2.7e-170) {
tmp = t_2;
} else if (y_46_im <= 3.2e-122) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * t_0;
} else if (y_46_im <= 1.45e+15) {
tmp = t_2;
} else {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - 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 t_0 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double t_1 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double t_2 = t_0 * t_1;
double tmp;
if (y_46_im <= -1.85e-5) {
tmp = Math.sin(t_1) / Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
} else if (y_46_im <= -2.7e-170) {
tmp = t_2;
} else if (y_46_im <= 3.2e-122) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) * t_0;
} else if (y_46_im <= 1.45e+15) {
tmp = t_2;
} else {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) t_1 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) t_2 = t_0 * t_1 tmp = 0 if y_46_im <= -1.85e-5: tmp = math.sin(t_1) / math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) elif y_46_im <= -2.7e-170: tmp = t_2 elif y_46_im <= 3.2e-122: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) * t_0 elif y_46_im <= 1.45e+15: tmp = t_2 else: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re t_1 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) t_2 = Float64(t_0 * t_1) tmp = 0.0 if (y_46_im <= -1.85e-5) tmp = Float64(sin(t_1) / exp(Float64(atan(x_46_im, x_46_re) * y_46_im))); elseif (y_46_im <= -2.7e-170) tmp = t_2; elseif (y_46_im <= 3.2e-122) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) * t_0); elseif (y_46_im <= 1.45e+15) tmp = t_2; else tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re; t_1 = y_46_im * log(hypot(x_46_im, x_46_re)); t_2 = t_0 * t_1; tmp = 0.0; if (y_46_im <= -1.85e-5) tmp = sin(t_1) / exp((atan2(x_46_im, x_46_re) * y_46_im)); elseif (y_46_im <= -2.7e-170) tmp = t_2; elseif (y_46_im <= 3.2e-122) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * t_0; elseif (y_46_im <= 1.45e+15) tmp = t_2; else tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$1 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 * t$95$1), $MachinePrecision]}, If[LessEqual[y$46$im, -1.85e-5], N[(N[Sin[t$95$1], $MachinePrecision] / N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$im, -2.7e-170], t$95$2, If[LessEqual[y$46$im, 3.2e-122], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[y$46$im, 1.45e+15], t$95$2, N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
t_1 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
t_2 := t\_0 \cdot t\_1\\
\mathbf{if}\;y.im \leq -1.85 \cdot 10^{-5}:\\
\;\;\;\;\frac{\sin t\_1}{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}\\
\mathbf{elif}\;y.im \leq -2.7 \cdot 10^{-170}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{-122}:\\
\;\;\;\;\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot t\_0\\
\mathbf{elif}\;y.im \leq 1.45 \cdot 10^{+15}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\end{array}
\end{array}
if y.im < -1.84999999999999991e-5Initial program 28.5%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified55.9%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6462.5%
Simplified62.5%
if -1.84999999999999991e-5 < y.im < -2.6999999999999999e-170 or 3.2000000000000002e-122 < y.im < 1.45e15Initial program 50.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified82.2%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6482.0%
Simplified82.0%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6482.0%
Applied egg-rr82.0%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6480.2%
Simplified80.2%
if -2.6999999999999999e-170 < y.im < 3.2000000000000002e-122Initial program 50.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
if 1.45e15 < y.im Initial program 22.7%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6450.7%
Simplified50.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6461.6%
Simplified61.6%
Final simplification71.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (pow (hypot x.im x.re) y.re))
(t_1 (* t_0 (* y.im (log (hypot x.im x.re)))))
(t_2
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))))
(if (<= y.im -5.5e+60)
t_2
(if (<= y.im -1.9e-174)
t_1
(if (<= y.im 3.2e-122)
(* (sin (* y.re (atan2 x.im x.re))) t_0)
(if (<= y.im 2.3e+15) t_1 t_2))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(hypot(x_46_im, x_46_re), y_46_re);
double t_1 = t_0 * (y_46_im * log(hypot(x_46_im, x_46_re)));
double t_2 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double tmp;
if (y_46_im <= -5.5e+60) {
tmp = t_2;
} else if (y_46_im <= -1.9e-174) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * t_0;
} else if (y_46_im <= 2.3e+15) {
tmp = t_1;
} else {
tmp = t_2;
}
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 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double t_1 = t_0 * (y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)));
double t_2 = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double tmp;
if (y_46_im <= -5.5e+60) {
tmp = t_2;
} else if (y_46_im <= -1.9e-174) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) * t_0;
} else if (y_46_im <= 2.3e+15) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) t_1 = t_0 * (y_46_im * math.log(math.hypot(x_46_im, x_46_re))) t_2 = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) tmp = 0 if y_46_im <= -5.5e+60: tmp = t_2 elif y_46_im <= -1.9e-174: tmp = t_1 elif y_46_im <= 3.2e-122: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) * t_0 elif y_46_im <= 2.3e+15: tmp = t_1 else: tmp = t_2 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re t_1 = Float64(t_0 * Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) t_2 = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))) tmp = 0.0 if (y_46_im <= -5.5e+60) tmp = t_2; elseif (y_46_im <= -1.9e-174) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) * t_0); elseif (y_46_im <= 2.3e+15) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re; t_1 = t_0 * (y_46_im * log(hypot(x_46_im, x_46_re))); t_2 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); tmp = 0.0; if (y_46_im <= -5.5e+60) tmp = t_2; elseif (y_46_im <= -1.9e-174) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * t_0; elseif (y_46_im <= 2.3e+15) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$im, -5.5e+60], t$95$2, If[LessEqual[y$46$im, -1.9e-174], t$95$1, If[LessEqual[y$46$im, 3.2e-122], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[y$46$im, 2.3e+15], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
t_1 := t\_0 \cdot \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
t_2 := y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\mathbf{if}\;y.im \leq -5.5 \cdot 10^{+60}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.im \leq -1.9 \cdot 10^{-174}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{-122}:\\
\;\;\;\;\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot t\_0\\
\mathbf{elif}\;y.im \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.im < -5.5000000000000001e60 or 2.3e15 < y.im Initial program 24.0%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6450.7%
Simplified50.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6461.7%
Simplified61.7%
if -5.5000000000000001e60 < y.im < -1.9000000000000001e-174 or 3.2000000000000002e-122 < y.im < 2.3e15Initial program 47.4%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified76.8%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6477.8%
Simplified77.8%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6477.8%
Applied egg-rr77.8%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6473.9%
Simplified73.9%
if -1.9000000000000001e-174 < y.im < 3.2000000000000002e-122Initial program 50.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
Final simplification70.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (pow (hypot x.im x.re) y.re))
(t_1 (* t_0 (* y.im (log (hypot x.im x.re)))))
(t_2
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))))
(if (<= y.im -1.2e+60)
t_2
(if (<= y.im -3.9e-174)
t_1
(if (<= y.im 3.2e-122)
(* (* y.re (atan2 x.im x.re)) t_0)
(if (<= y.im 1.75e+15) t_1 t_2))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(hypot(x_46_im, x_46_re), y_46_re);
double t_1 = t_0 * (y_46_im * log(hypot(x_46_im, x_46_re)));
double t_2 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double tmp;
if (y_46_im <= -1.2e+60) {
tmp = t_2;
} else if (y_46_im <= -3.9e-174) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = (y_46_re * atan2(x_46_im, x_46_re)) * t_0;
} else if (y_46_im <= 1.75e+15) {
tmp = t_1;
} else {
tmp = t_2;
}
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 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double t_1 = t_0 * (y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)));
double t_2 = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
double tmp;
if (y_46_im <= -1.2e+60) {
tmp = t_2;
} else if (y_46_im <= -3.9e-174) {
tmp = t_1;
} else if (y_46_im <= 3.2e-122) {
tmp = (y_46_re * Math.atan2(x_46_im, x_46_re)) * t_0;
} else if (y_46_im <= 1.75e+15) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) t_1 = t_0 * (y_46_im * math.log(math.hypot(x_46_im, x_46_re))) t_2 = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) tmp = 0 if y_46_im <= -1.2e+60: tmp = t_2 elif y_46_im <= -3.9e-174: tmp = t_1 elif y_46_im <= 3.2e-122: tmp = (y_46_re * math.atan2(x_46_im, x_46_re)) * t_0 elif y_46_im <= 1.75e+15: tmp = t_1 else: tmp = t_2 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re t_1 = Float64(t_0 * Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) t_2 = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))) tmp = 0.0 if (y_46_im <= -1.2e+60) tmp = t_2; elseif (y_46_im <= -3.9e-174) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = Float64(Float64(y_46_re * atan(x_46_im, x_46_re)) * t_0); elseif (y_46_im <= 1.75e+15) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = hypot(x_46_im, x_46_re) ^ y_46_re; t_1 = t_0 * (y_46_im * log(hypot(x_46_im, x_46_re))); t_2 = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); tmp = 0.0; if (y_46_im <= -1.2e+60) tmp = t_2; elseif (y_46_im <= -3.9e-174) tmp = t_1; elseif (y_46_im <= 3.2e-122) tmp = (y_46_re * atan2(x_46_im, x_46_re)) * t_0; elseif (y_46_im <= 1.75e+15) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$im, -1.2e+60], t$95$2, If[LessEqual[y$46$im, -3.9e-174], t$95$1, If[LessEqual[y$46$im, 3.2e-122], N[(N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[y$46$im, 1.75e+15], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
t_1 := t\_0 \cdot \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
t_2 := y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\mathbf{if}\;y.im \leq -1.2 \cdot 10^{+60}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.im \leq -3.9 \cdot 10^{-174}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 3.2 \cdot 10^{-122}:\\
\;\;\;\;\left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot t\_0\\
\mathbf{elif}\;y.im \leq 1.75 \cdot 10^{+15}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.im < -1.2e60 or 1.75e15 < y.im Initial program 24.0%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6450.7%
Simplified50.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6461.7%
Simplified61.7%
if -1.2e60 < y.im < -3.8999999999999999e-174 or 3.2000000000000002e-122 < y.im < 1.75e15Initial program 47.4%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified76.8%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6477.8%
Simplified77.8%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6477.8%
Applied egg-rr77.8%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6473.9%
Simplified73.9%
if -3.8999999999999999e-174 < y.im < 3.2000000000000002e-122Initial program 50.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.3%
Simplified82.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6477.3%
Simplified77.3%
Final simplification69.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.im (log (hypot x.im x.re))))
(t_1 (pow (hypot x.im x.re) y.re))
(t_2 (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(if (<= y.re -1.66e-10)
(* t_1 (sin (/ 1.0 (/ (/ 1.0 (atan2 x.im x.re)) y.re))))
(if (<= y.re 9.5e-247)
(* t_0 t_2)
(if (<= y.re 5.3e+20)
(* y.re (* (atan2 x.im x.re) t_2))
(* 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 = y_46_im * log(hypot(x_46_im, x_46_re));
double t_1 = pow(hypot(x_46_im, x_46_re), y_46_re);
double t_2 = exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)));
double tmp;
if (y_46_re <= -1.66e-10) {
tmp = t_1 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 9.5e-247) {
tmp = t_0 * t_2;
} else if (y_46_re <= 5.3e+20) {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * t_2);
} else {
tmp = t_1 * t_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_im * Math.log(Math.hypot(x_46_im, x_46_re));
double t_1 = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
double t_2 = Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)));
double tmp;
if (y_46_re <= -1.66e-10) {
tmp = t_1 * Math.sin((1.0 / ((1.0 / Math.atan2(x_46_im, x_46_re)) / y_46_re)));
} else if (y_46_re <= 9.5e-247) {
tmp = t_0 * t_2;
} else if (y_46_re <= 5.3e+20) {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * t_2);
} else {
tmp = t_1 * t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) t_1 = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) t_2 = math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))) tmp = 0 if y_46_re <= -1.66e-10: tmp = t_1 * math.sin((1.0 / ((1.0 / math.atan2(x_46_im, x_46_re)) / y_46_re))) elif y_46_re <= 9.5e-247: tmp = t_0 * t_2 elif y_46_re <= 5.3e+20: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * t_2) else: tmp = t_1 * t_0 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) t_1 = hypot(x_46_im, x_46_re) ^ y_46_re t_2 = exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))) tmp = 0.0 if (y_46_re <= -1.66e-10) tmp = Float64(t_1 * sin(Float64(1.0 / Float64(Float64(1.0 / atan(x_46_im, x_46_re)) / y_46_re)))); elseif (y_46_re <= 9.5e-247) tmp = Float64(t_0 * t_2); elseif (y_46_re <= 5.3e+20) tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * t_2)); else tmp = Float64(t_1 * t_0); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_im * log(hypot(x_46_im, x_46_re)); t_1 = hypot(x_46_im, x_46_re) ^ y_46_re; t_2 = exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))); tmp = 0.0; if (y_46_re <= -1.66e-10) tmp = t_1 * sin((1.0 / ((1.0 / atan2(x_46_im, x_46_re)) / y_46_re))); elseif (y_46_re <= 9.5e-247) tmp = t_0 * t_2; elseif (y_46_re <= 5.3e+20) tmp = y_46_re * (atan2(x_46_im, x_46_re) * t_2); else tmp = t_1 * 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[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]}, Block[{t$95$2 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -1.66e-10], N[(t$95$1 * N[Sin[N[(1.0 / N[(N[(1.0 / N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] / y$46$re), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[(t$95$0 * t$95$2), $MachinePrecision], If[LessEqual[y$46$re, 5.3e+20], N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * t$95$0), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
t_1 := {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
t_2 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\\
\mathbf{if}\;y.re \leq -1.66 \cdot 10^{-10}:\\
\;\;\;\;t\_1 \cdot \sin \left(\frac{1}{\frac{\frac{1}{\tan^{-1}_* \frac{x.im}{x.re}}}{y.re}}\right)\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;t\_0 \cdot t\_2\\
\mathbf{elif}\;y.re \leq 5.3 \cdot 10^{+20}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot t\_2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot t\_0\\
\end{array}
\end{array}
if y.re < -1.66e-10Initial program 41.4%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.7%
Simplified78.7%
remove-double-divN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6480.1%
Applied egg-rr80.1%
if -1.66e-10 < y.re < 9.49999999999999939e-247Initial program 35.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6467.3%
Simplified67.3%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6467.1%
Simplified67.1%
if 9.49999999999999939e-247 < y.re < 5.3e20Initial program 37.5%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.9%
Simplified36.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6466.0%
Simplified66.0%
if 5.3e20 < y.re Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.9%
Simplified52.9%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.9%
Applied egg-rr52.9%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.9%
Simplified64.9%
Final simplification70.0%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.im (log (hypot x.im x.re))))
(t_1 (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(if (<= y.re -0.98)
(*
(sin (* y.re (atan2 x.im x.re)))
(pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))
(if (<= y.re 9.5e-247)
(* t_0 t_1)
(if (<= y.re 4.4e+19)
(* y.re (* (atan2 x.im x.re) t_1))
(* (pow (hypot x.im x.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 = y_46_im * log(hypot(x_46_im, x_46_re));
double t_1 = exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)));
double tmp;
if (y_46_re <= -0.98) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = t_0 * t_1;
} else if (y_46_re <= 4.4e+19) {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * t_1);
} else {
tmp = pow(hypot(x_46_im, x_46_re), y_46_re) * t_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_im * Math.log(Math.hypot(x_46_im, x_46_re));
double t_1 = Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)));
double tmp;
if (y_46_re <= -0.98) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = t_0 * t_1;
} else if (y_46_re <= 4.4e+19) {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * t_1);
} else {
tmp = Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re) * t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) t_1 = math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))) tmp = 0 if y_46_re <= -0.98: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) elif y_46_re <= 9.5e-247: tmp = t_0 * t_1 elif y_46_re <= 4.4e+19: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * t_1) else: tmp = math.pow(math.hypot(x_46_im, x_46_re), y_46_re) * t_0 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) t_1 = exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))) tmp = 0.0 if (y_46_re <= -0.98) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); elseif (y_46_re <= 9.5e-247) tmp = Float64(t_0 * t_1); elseif (y_46_re <= 4.4e+19) tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * t_1)); else tmp = Float64((hypot(x_46_im, x_46_re) ^ y_46_re) * t_0); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_im * log(hypot(x_46_im, x_46_re)); t_1 = exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))); tmp = 0.0; if (y_46_re <= -0.98) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re); elseif (y_46_re <= 9.5e-247) tmp = t_0 * t_1; elseif (y_46_re <= 4.4e+19) tmp = y_46_re * (atan2(x_46_im, x_46_re) * t_1); else tmp = (hypot(x_46_im, x_46_re) ^ y_46_re) * 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[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -0.98], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[N[(x$46$im + N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] * 0.5), $MachinePrecision] / x$46$im), $MachinePrecision]), $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[(t$95$0 * t$95$1), $MachinePrecision], If[LessEqual[y$46$re, 4.4e+19], N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], N[(N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision] * t$95$0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
t_1 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\\
\mathbf{if}\;y.re \leq -0.98:\\
\;\;\;\;\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;t\_0 \cdot t\_1\\
\mathbf{elif}\;y.re \leq 4.4 \cdot 10^{+19}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re} \cdot t\_0\\
\end{array}
\end{array}
if y.re < -0.97999999999999998Initial program 41.4%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.7%
Simplified78.7%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.8%
Simplified78.8%
if -0.97999999999999998 < y.re < 9.49999999999999939e-247Initial program 35.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6467.3%
Simplified67.3%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6467.1%
Simplified67.1%
if 9.49999999999999939e-247 < y.re < 4.4e19Initial program 37.5%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.9%
Simplified36.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6466.0%
Simplified66.0%
if 4.4e19 < y.re Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.9%
Simplified52.9%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.9%
Applied egg-rr52.9%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.9%
Simplified64.9%
Final simplification69.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im))))
(if (<= y.re -6.7e-12)
(*
(sin (* y.re (atan2 x.im x.re)))
(pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))
(if (<= y.re 2e+20)
(* y.re (* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(/ (* y.im (pow t_0 (/ y.re 2.0))) (/ 2.0 (log 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 * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -6.7e-12) {
tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
} else if (y_46_re <= 2e+20) {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
} else {
tmp = (y_46_im * pow(t_0, (y_46_re / 2.0))) / (2.0 / log(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 * x_46re) + (x_46im * x_46im)
if (y_46re <= (-6.7d-12)) then
tmp = sin((y_46re * atan2(x_46im, x_46re))) * ((x_46im + (((x_46re * x_46re) * 0.5d0) / x_46im)) ** y_46re)
else if (y_46re <= 2d+20) then
tmp = y_46re * (atan2(x_46im, x_46re) * exp((atan2(x_46im, x_46re) * (0.0d0 - y_46im))))
else
tmp = (y_46im * (t_0 ** (y_46re / 2.0d0))) / (2.0d0 / log(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 * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -6.7e-12) {
tmp = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re))) * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
} else if (y_46_re <= 2e+20) {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
} else {
tmp = (y_46_im * Math.pow(t_0, (y_46_re / 2.0))) / (2.0 / Math.log(t_0));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im) tmp = 0 if y_46_re <= -6.7e-12: tmp = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) elif y_46_re <= 2e+20: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) else: tmp = (y_46_im * math.pow(t_0, (y_46_re / 2.0))) / (2.0 / math.log(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 * x_46_re) + Float64(x_46_im * x_46_im)) tmp = 0.0 if (y_46_re <= -6.7e-12) tmp = Float64(sin(Float64(y_46_re * atan(x_46_im, x_46_re))) * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); elseif (y_46_re <= 2e+20) tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))); else tmp = Float64(Float64(y_46_im * (t_0 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(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 * x_46_re) + (x_46_im * x_46_im); tmp = 0.0; if (y_46_re <= -6.7e-12) tmp = sin((y_46_re * atan2(x_46_im, x_46_re))) * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re); elseif (y_46_re <= 2e+20) tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); else tmp = (y_46_im * (t_0 ^ (y_46_re / 2.0))) / (2.0 / log(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$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -6.7e-12], N[(N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[N[(x$46$im + N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] * 0.5), $MachinePrecision] / x$46$im), $MachinePrecision]), $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 2e+20], N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y$46$im * N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
\mathbf{if}\;y.re \leq -6.7 \cdot 10^{-12}:\\
\;\;\;\;\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\mathbf{elif}\;y.re \leq 2 \cdot 10^{+20}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_0}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_0}}\\
\end{array}
\end{array}
if y.re < -6.7000000000000001e-12Initial program 40.3%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6476.6%
Simplified76.6%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.7%
Simplified76.7%
if -6.7000000000000001e-12 < y.re < 2e20Initial program 36.9%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6434.1%
Simplified34.1%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6455.5%
Simplified55.5%
if 2e20 < y.re Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.9%
Simplified52.9%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.9%
Applied egg-rr52.9%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.9%
Simplified64.9%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr60.9%
Final simplification62.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im))))
(if (<= y.re -6.7e-12)
(* (* y.re (atan2 x.im x.re)) (pow (hypot x.im x.re) y.re))
(if (<= y.re 2.85e+19)
(* y.re (* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(/ (* y.im (pow t_0 (/ y.re 2.0))) (/ 2.0 (log 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 * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -6.7e-12) {
tmp = (y_46_re * atan2(x_46_im, x_46_re)) * pow(hypot(x_46_im, x_46_re), y_46_re);
} else if (y_46_re <= 2.85e+19) {
tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
} else {
tmp = (y_46_im * pow(t_0, (y_46_re / 2.0))) / (2.0 / log(t_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 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -6.7e-12) {
tmp = (y_46_re * Math.atan2(x_46_im, x_46_re)) * Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
} else if (y_46_re <= 2.85e+19) {
tmp = y_46_re * (Math.atan2(x_46_im, x_46_re) * Math.exp((Math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im))));
} else {
tmp = (y_46_im * Math.pow(t_0, (y_46_re / 2.0))) / (2.0 / Math.log(t_0));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im) tmp = 0 if y_46_re <= -6.7e-12: tmp = (y_46_re * math.atan2(x_46_im, x_46_re)) * math.pow(math.hypot(x_46_im, x_46_re), y_46_re) elif y_46_re <= 2.85e+19: tmp = y_46_re * (math.atan2(x_46_im, x_46_re) * math.exp((math.atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))) else: tmp = (y_46_im * math.pow(t_0, (y_46_re / 2.0))) / (2.0 / math.log(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 * x_46_re) + Float64(x_46_im * x_46_im)) tmp = 0.0 if (y_46_re <= -6.7e-12) tmp = Float64(Float64(y_46_re * atan(x_46_im, x_46_re)) * (hypot(x_46_im, x_46_re) ^ y_46_re)); elseif (y_46_re <= 2.85e+19) tmp = Float64(y_46_re * Float64(atan(x_46_im, x_46_re) * exp(Float64(atan(x_46_im, x_46_re) * Float64(0.0 - y_46_im))))); else tmp = Float64(Float64(y_46_im * (t_0 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(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 * x_46_re) + (x_46_im * x_46_im); tmp = 0.0; if (y_46_re <= -6.7e-12) tmp = (y_46_re * atan2(x_46_im, x_46_re)) * (hypot(x_46_im, x_46_re) ^ y_46_re); elseif (y_46_re <= 2.85e+19) tmp = y_46_re * (atan2(x_46_im, x_46_re) * exp((atan2(x_46_im, x_46_re) * (0.0 - y_46_im)))); else tmp = (y_46_im * (t_0 ^ (y_46_re / 2.0))) / (2.0 / log(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$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -6.7e-12], N[(N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 2.85e+19], N[(y$46$re * N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * N[(0.0 - y$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y$46$im * N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
\mathbf{if}\;y.re \leq -6.7 \cdot 10^{-12}:\\
\;\;\;\;\left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
\mathbf{elif}\;y.re \leq 2.85 \cdot 10^{+19}:\\
\;\;\;\;y.re \cdot \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot \left(0 - y.im\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_0}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_0}}\\
\end{array}
\end{array}
if y.re < -6.7000000000000001e-12Initial program 40.3%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6476.6%
Simplified76.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6474.1%
Simplified74.1%
if -6.7000000000000001e-12 < y.re < 2.85e19Initial program 36.9%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6434.1%
Simplified34.1%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6455.5%
Simplified55.5%
if 2.85e19 < y.re Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified54.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.9%
Simplified52.9%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6452.9%
Applied egg-rr52.9%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6464.9%
Simplified64.9%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr60.9%
Final simplification61.8%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im)))
(t_1 (* y.re (atan2 x.im x.re)))
(t_2 (* y.im (log (hypot x.im x.re)))))
(if (<= y.re -2.35e+122)
(* t_2 (pow x.im y.re))
(if (<= y.re -6.2e-10)
(* (sin t_1) (pow x.re y.re))
(if (<= y.re 9.5e-247)
(sin t_2)
(if (<= y.re 1.35e-34)
t_1
(/ (* y.im (pow t_0 (/ y.re 2.0))) (/ 2.0 (log 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 * x_46_re) + (x_46_im * x_46_im);
double t_1 = y_46_re * atan2(x_46_im, x_46_re);
double t_2 = y_46_im * log(hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -2.35e+122) {
tmp = t_2 * pow(x_46_im, y_46_re);
} else if (y_46_re <= -6.2e-10) {
tmp = sin(t_1) * pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = sin(t_2);
} else if (y_46_re <= 1.35e-34) {
tmp = t_1;
} else {
tmp = (y_46_im * pow(t_0, (y_46_re / 2.0))) / (2.0 / log(t_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 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double t_1 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_2 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -2.35e+122) {
tmp = t_2 * Math.pow(x_46_im, y_46_re);
} else if (y_46_re <= -6.2e-10) {
tmp = Math.sin(t_1) * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = Math.sin(t_2);
} else if (y_46_re <= 1.35e-34) {
tmp = t_1;
} else {
tmp = (y_46_im * Math.pow(t_0, (y_46_re / 2.0))) / (2.0 / Math.log(t_0));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im) t_1 = y_46_re * math.atan2(x_46_im, x_46_re) t_2 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) tmp = 0 if y_46_re <= -2.35e+122: tmp = t_2 * math.pow(x_46_im, y_46_re) elif y_46_re <= -6.2e-10: tmp = math.sin(t_1) * math.pow(x_46_re, y_46_re) elif y_46_re <= 9.5e-247: tmp = math.sin(t_2) elif y_46_re <= 1.35e-34: tmp = t_1 else: tmp = (y_46_im * math.pow(t_0, (y_46_re / 2.0))) / (2.0 / math.log(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 * x_46_re) + Float64(x_46_im * x_46_im)) t_1 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_2 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) tmp = 0.0 if (y_46_re <= -2.35e+122) tmp = Float64(t_2 * (x_46_im ^ y_46_re)); elseif (y_46_re <= -6.2e-10) tmp = Float64(sin(t_1) * (x_46_re ^ y_46_re)); elseif (y_46_re <= 9.5e-247) tmp = sin(t_2); elseif (y_46_re <= 1.35e-34) tmp = t_1; else tmp = Float64(Float64(y_46_im * (t_0 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(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 * x_46_re) + (x_46_im * x_46_im); t_1 = y_46_re * atan2(x_46_im, x_46_re); t_2 = y_46_im * log(hypot(x_46_im, x_46_re)); tmp = 0.0; if (y_46_re <= -2.35e+122) tmp = t_2 * (x_46_im ^ y_46_re); elseif (y_46_re <= -6.2e-10) tmp = sin(t_1) * (x_46_re ^ y_46_re); elseif (y_46_re <= 9.5e-247) tmp = sin(t_2); elseif (y_46_re <= 1.35e-34) tmp = t_1; else tmp = (y_46_im * (t_0 ^ (y_46_re / 2.0))) / (2.0 / log(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$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -2.35e+122], N[(t$95$2 * N[Power[x$46$im, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -6.2e-10], N[(N[Sin[t$95$1], $MachinePrecision] * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[Sin[t$95$2], $MachinePrecision], If[LessEqual[y$46$re, 1.35e-34], t$95$1, N[(N[(y$46$im * N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
t_1 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_2 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{if}\;y.re \leq -2.35 \cdot 10^{+122}:\\
\;\;\;\;t\_2 \cdot {x.im}^{y.re}\\
\mathbf{elif}\;y.re \leq -6.2 \cdot 10^{-10}:\\
\;\;\;\;\sin t\_1 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;\sin t\_2\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-34}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_0}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_0}}\\
\end{array}
\end{array}
if y.re < -2.35000000000000012e122Initial program 47.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified63.0%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6471.7%
Simplified71.7%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6471.7%
Applied egg-rr71.7%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6480.5%
Simplified80.5%
Taylor expanded in x.im around inf
Simplified65.5%
if -2.35000000000000012e122 < y.re < -6.2000000000000003e-10Initial program 29.2%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6483.5%
Simplified83.5%
Taylor expanded in x.im around 0
pow-lowering-pow.f6467.2%
Simplified67.2%
if -6.2000000000000003e-10 < y.re < 9.49999999999999939e-247Initial program 35.7%
Taylor expanded in y.im around 0
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6422.9%
Simplified22.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6432.7%
Simplified32.7%
if 9.49999999999999939e-247 < y.re < 1.35000000000000008e-34Initial program 38.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
if 1.35000000000000008e-34 < y.re Initial program 37.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified58.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.2%
Simplified55.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6455.2%
Applied egg-rr55.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.5%
Simplified60.5%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr58.6%
Final simplification48.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im)))
(t_1 (* y.re (atan2 x.im x.re)))
(t_2 (* y.im (log (hypot x.im x.re)))))
(if (<= y.re -5e+169)
(* t_2 (pow x.im y.re))
(if (<= y.re -6.8e-8)
(* t_1 (pow x.re y.re))
(if (<= y.re 9.5e-247)
(sin t_2)
(if (<= y.re 1.2e-33)
t_1
(/ (* y.im (pow t_0 (/ y.re 2.0))) (/ 2.0 (log 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 * x_46_re) + (x_46_im * x_46_im);
double t_1 = y_46_re * atan2(x_46_im, x_46_re);
double t_2 = y_46_im * log(hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -5e+169) {
tmp = t_2 * pow(x_46_im, y_46_re);
} else if (y_46_re <= -6.8e-8) {
tmp = t_1 * pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = sin(t_2);
} else if (y_46_re <= 1.2e-33) {
tmp = t_1;
} else {
tmp = (y_46_im * pow(t_0, (y_46_re / 2.0))) / (2.0 / log(t_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 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double t_1 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_2 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double tmp;
if (y_46_re <= -5e+169) {
tmp = t_2 * Math.pow(x_46_im, y_46_re);
} else if (y_46_re <= -6.8e-8) {
tmp = t_1 * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = Math.sin(t_2);
} else if (y_46_re <= 1.2e-33) {
tmp = t_1;
} else {
tmp = (y_46_im * Math.pow(t_0, (y_46_re / 2.0))) / (2.0 / Math.log(t_0));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im) t_1 = y_46_re * math.atan2(x_46_im, x_46_re) t_2 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) tmp = 0 if y_46_re <= -5e+169: tmp = t_2 * math.pow(x_46_im, y_46_re) elif y_46_re <= -6.8e-8: tmp = t_1 * math.pow(x_46_re, y_46_re) elif y_46_re <= 9.5e-247: tmp = math.sin(t_2) elif y_46_re <= 1.2e-33: tmp = t_1 else: tmp = (y_46_im * math.pow(t_0, (y_46_re / 2.0))) / (2.0 / math.log(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 * x_46_re) + Float64(x_46_im * x_46_im)) t_1 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_2 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) tmp = 0.0 if (y_46_re <= -5e+169) tmp = Float64(t_2 * (x_46_im ^ y_46_re)); elseif (y_46_re <= -6.8e-8) tmp = Float64(t_1 * (x_46_re ^ y_46_re)); elseif (y_46_re <= 9.5e-247) tmp = sin(t_2); elseif (y_46_re <= 1.2e-33) tmp = t_1; else tmp = Float64(Float64(y_46_im * (t_0 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(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 * x_46_re) + (x_46_im * x_46_im); t_1 = y_46_re * atan2(x_46_im, x_46_re); t_2 = y_46_im * log(hypot(x_46_im, x_46_re)); tmp = 0.0; if (y_46_re <= -5e+169) tmp = t_2 * (x_46_im ^ y_46_re); elseif (y_46_re <= -6.8e-8) tmp = t_1 * (x_46_re ^ y_46_re); elseif (y_46_re <= 9.5e-247) tmp = sin(t_2); elseif (y_46_re <= 1.2e-33) tmp = t_1; else tmp = (y_46_im * (t_0 ^ (y_46_re / 2.0))) / (2.0 / log(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$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -5e+169], N[(t$95$2 * N[Power[x$46$im, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -6.8e-8], N[(t$95$1 * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[Sin[t$95$2], $MachinePrecision], If[LessEqual[y$46$re, 1.2e-33], t$95$1, N[(N[(y$46$im * N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
t_1 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_2 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{if}\;y.re \leq -5 \cdot 10^{+169}:\\
\;\;\;\;t\_2 \cdot {x.im}^{y.re}\\
\mathbf{elif}\;y.re \leq -6.8 \cdot 10^{-8}:\\
\;\;\;\;t\_1 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;\sin t\_2\\
\mathbf{elif}\;y.re \leq 1.2 \cdot 10^{-33}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_0}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_0}}\\
\end{array}
\end{array}
if y.re < -5.00000000000000017e169Initial program 52.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified71.1%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6473.7%
Simplified73.7%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6473.7%
Applied egg-rr73.7%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6484.3%
Simplified84.3%
Taylor expanded in x.im around inf
Simplified71.3%
if -5.00000000000000017e169 < y.re < -6.8e-8Initial program 28.1%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6468.9%
Simplified68.9%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6468.9%
Applied egg-rr68.9%
Taylor expanded in x.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6456.7%
Simplified56.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6459.9%
Simplified59.9%
if -6.8e-8 < y.re < 9.49999999999999939e-247Initial program 35.7%
Taylor expanded in y.im around 0
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6422.9%
Simplified22.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6432.7%
Simplified32.7%
if 9.49999999999999939e-247 < y.re < 1.2e-33Initial program 38.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
if 1.2e-33 < y.re Initial program 37.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified58.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.2%
Simplified55.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6455.2%
Applied egg-rr55.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.5%
Simplified60.5%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr58.6%
Final simplification48.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))) (t_1 (+ (* x.re x.re) (* x.im x.im))))
(if (<= y.re -1e-105)
(* t_0 (pow (hypot x.im x.re) y.re))
(if (<= y.re 9.5e-247)
(sin (* y.im (log (hypot x.im x.re))))
(if (<= y.re 9.5e-34)
t_0
(/ (* y.im (pow t_1 (/ y.re 2.0))) (/ 2.0 (log t_1))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -1e-105) {
tmp = t_0 * pow(hypot(x_46_im, x_46_re), y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 9.5e-34) {
tmp = t_0;
} else {
tmp = (y_46_im * pow(t_1, (y_46_re / 2.0))) / (2.0 / log(t_1));
}
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 * Math.atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -1e-105) {
tmp = t_0 * Math.pow(Math.hypot(x_46_im, x_46_re), y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 9.5e-34) {
tmp = t_0;
} else {
tmp = (y_46_im * Math.pow(t_1, (y_46_re / 2.0))) / (2.0 / Math.log(t_1));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im) tmp = 0 if y_46_re <= -1e-105: tmp = t_0 * math.pow(math.hypot(x_46_im, x_46_re), y_46_re) elif y_46_re <= 9.5e-247: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) elif y_46_re <= 9.5e-34: tmp = t_0 else: tmp = (y_46_im * math.pow(t_1, (y_46_re / 2.0))) / (2.0 / math.log(t_1)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)) tmp = 0.0 if (y_46_re <= -1e-105) tmp = Float64(t_0 * (hypot(x_46_im, x_46_re) ^ y_46_re)); elseif (y_46_re <= 9.5e-247) tmp = sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 9.5e-34) tmp = t_0; else tmp = Float64(Float64(y_46_im * (t_1 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(t_1))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im); tmp = 0.0; if (y_46_re <= -1e-105) tmp = t_0 * (hypot(x_46_im, x_46_re) ^ y_46_re); elseif (y_46_re <= 9.5e-247) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 9.5e-34) tmp = t_0; else tmp = (y_46_im * (t_1 ^ (y_46_re / 2.0))) / (2.0 / log(t_1)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1e-105], N[(t$95$0 * N[Power[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[Sin[N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[y$46$re, 9.5e-34], t$95$0, N[(N[(y$46$im * N[Power[t$95$1, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := x.re \cdot x.re + x.im \cdot x.im\\
\mathbf{if}\;y.re \leq -1 \cdot 10^{-105}:\\
\;\;\;\;t\_0 \cdot {\left(\mathsf{hypot}\left(x.im, x.re\right)\right)}^{y.re}\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-34}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_1}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_1}}\\
\end{array}
\end{array}
if y.re < -9.99999999999999965e-106Initial program 42.1%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6468.2%
Simplified68.2%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6466.2%
Simplified66.2%
if -9.99999999999999965e-106 < y.re < 9.49999999999999939e-247Initial program 33.3%
Taylor expanded in y.im around 0
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6422.2%
Simplified22.2%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6436.9%
Simplified36.9%
if 9.49999999999999939e-247 < y.re < 9.49999999999999985e-34Initial program 38.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
if 9.49999999999999985e-34 < y.re Initial program 37.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified58.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.2%
Simplified55.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6455.2%
Applied egg-rr55.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.5%
Simplified60.5%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr58.6%
Final simplification51.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))) (t_1 (+ (* x.re x.re) (* x.im x.im))))
(if (<= y.re -0.0031)
(* t_0 (pow x.re y.re))
(if (<= y.re 9.5e-247)
(sin (* y.im (log (hypot x.im x.re))))
(if (<= y.re 1.25e-33)
t_0
(/ (* y.im (pow t_1 (/ y.re 2.0))) (/ 2.0 (log t_1))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -0.0031) {
tmp = t_0 * pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 1.25e-33) {
tmp = t_0;
} else {
tmp = (y_46_im * pow(t_1, (y_46_re / 2.0))) / (2.0 / log(t_1));
}
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 * Math.atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -0.0031) {
tmp = t_0 * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 9.5e-247) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 1.25e-33) {
tmp = t_0;
} else {
tmp = (y_46_im * Math.pow(t_1, (y_46_re / 2.0))) / (2.0 / Math.log(t_1));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im) tmp = 0 if y_46_re <= -0.0031: tmp = t_0 * math.pow(x_46_re, y_46_re) elif y_46_re <= 9.5e-247: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) elif y_46_re <= 1.25e-33: tmp = t_0 else: tmp = (y_46_im * math.pow(t_1, (y_46_re / 2.0))) / (2.0 / math.log(t_1)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)) tmp = 0.0 if (y_46_re <= -0.0031) tmp = Float64(t_0 * (x_46_re ^ y_46_re)); elseif (y_46_re <= 9.5e-247) tmp = sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 1.25e-33) tmp = t_0; else tmp = Float64(Float64(y_46_im * (t_1 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(t_1))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im); tmp = 0.0; if (y_46_re <= -0.0031) tmp = t_0 * (x_46_re ^ y_46_re); elseif (y_46_re <= 9.5e-247) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 1.25e-33) tmp = t_0; else tmp = (y_46_im * (t_1 ^ (y_46_re / 2.0))) / (2.0 / log(t_1)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -0.0031], N[(t$95$0 * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9.5e-247], N[Sin[N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[y$46$re, 1.25e-33], t$95$0, N[(N[(y$46$im * N[Power[t$95$1, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := x.re \cdot x.re + x.im \cdot x.im\\
\mathbf{if}\;y.re \leq -0.0031:\\
\;\;\;\;t\_0 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
\mathbf{elif}\;y.re \leq 1.25 \cdot 10^{-33}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_1}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_1}}\\
\end{array}
\end{array}
if y.re < -0.00309999999999999989Initial program 41.4%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.7%
Simplified78.7%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6478.7%
Applied egg-rr78.7%
Taylor expanded in x.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6454.8%
Simplified54.8%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6453.4%
Simplified53.4%
if -0.00309999999999999989 < y.re < 9.49999999999999939e-247Initial program 35.7%
Taylor expanded in y.im around 0
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6422.9%
Simplified22.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6432.7%
Simplified32.7%
if 9.49999999999999939e-247 < y.re < 1.25000000000000007e-33Initial program 38.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
if 1.25000000000000007e-33 < y.re Initial program 37.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified58.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.2%
Simplified55.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6455.2%
Applied egg-rr55.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.5%
Simplified60.5%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr58.6%
Final simplification44.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))) (t_1 (+ (* x.re x.re) (* x.im x.im))))
(if (<= y.re -1.36e-102)
(* t_0 (pow x.re y.re))
(if (<= y.re 9e-247)
(* y.im (log (hypot x.im x.re)))
(if (<= y.re 1.25e-33)
t_0
(/ (* y.im (pow t_1 (/ y.re 2.0))) (/ 2.0 (log t_1))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -1.36e-102) {
tmp = t_0 * pow(x_46_re, y_46_re);
} else if (y_46_re <= 9e-247) {
tmp = y_46_im * log(hypot(x_46_im, x_46_re));
} else if (y_46_re <= 1.25e-33) {
tmp = t_0;
} else {
tmp = (y_46_im * pow(t_1, (y_46_re / 2.0))) / (2.0 / log(t_1));
}
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 * Math.atan2(x_46_im, x_46_re);
double t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double tmp;
if (y_46_re <= -1.36e-102) {
tmp = t_0 * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 9e-247) {
tmp = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
} else if (y_46_re <= 1.25e-33) {
tmp = t_0;
} else {
tmp = (y_46_im * Math.pow(t_1, (y_46_re / 2.0))) / (2.0 / Math.log(t_1));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im) tmp = 0 if y_46_re <= -1.36e-102: tmp = t_0 * math.pow(x_46_re, y_46_re) elif y_46_re <= 9e-247: tmp = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) elif y_46_re <= 1.25e-33: tmp = t_0 else: tmp = (y_46_im * math.pow(t_1, (y_46_re / 2.0))) / (2.0 / math.log(t_1)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)) tmp = 0.0 if (y_46_re <= -1.36e-102) tmp = Float64(t_0 * (x_46_re ^ y_46_re)); elseif (y_46_re <= 9e-247) tmp = Float64(y_46_im * log(hypot(x_46_im, x_46_re))); elseif (y_46_re <= 1.25e-33) tmp = t_0; else tmp = Float64(Float64(y_46_im * (t_1 ^ Float64(y_46_re / 2.0))) / Float64(2.0 / log(t_1))); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); t_1 = (x_46_re * x_46_re) + (x_46_im * x_46_im); tmp = 0.0; if (y_46_re <= -1.36e-102) tmp = t_0 * (x_46_re ^ y_46_re); elseif (y_46_re <= 9e-247) tmp = y_46_im * log(hypot(x_46_im, x_46_re)); elseif (y_46_re <= 1.25e-33) tmp = t_0; else tmp = (y_46_im * (t_1 ^ (y_46_re / 2.0))) / (2.0 / log(t_1)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.36e-102], N[(t$95$0 * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 9e-247], N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.25e-33], t$95$0, N[(N[(y$46$im * N[Power[t$95$1, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[Log[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := x.re \cdot x.re + x.im \cdot x.im\\
\mathbf{if}\;y.re \leq -1.36 \cdot 10^{-102}:\\
\;\;\;\;t\_0 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 9 \cdot 10^{-247}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{elif}\;y.re \leq 1.25 \cdot 10^{-33}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{y.im \cdot {t\_1}^{\left(\frac{y.re}{2}\right)}}{\frac{2}{\log t\_1}}\\
\end{array}
\end{array}
if y.re < -1.36000000000000001e-102Initial program 42.1%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6468.2%
Simplified68.2%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6468.2%
Applied egg-rr68.2%
Taylor expanded in x.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6446.8%
Simplified46.8%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6445.6%
Simplified45.6%
if -1.36000000000000001e-102 < y.re < 9.0000000000000005e-247Initial program 33.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6471.1%
Simplified71.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6436.4%
Simplified36.4%
if 9.0000000000000005e-247 < y.re < 1.25000000000000007e-33Initial program 38.8%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6436.6%
Simplified36.6%
if 1.25000000000000007e-33 < y.re Initial program 37.9%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified58.6%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.2%
Simplified55.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6455.2%
Applied egg-rr55.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.5%
Simplified60.5%
associate-*l*N/A
pow1/2N/A
log-powN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
associate-/r/N/A
un-div-invN/A
Applied egg-rr58.6%
Final simplification44.6%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* (* y.re (atan2 x.im x.re)) (pow x.re y.re))))
(if (<= y.re -1.35e-102)
t_0
(if (<= y.re 9.5e-247)
(* y.im (log (hypot x.im x.re)))
(if (<= y.re 0.011)
(sin (* y.re (/ 1.0 (/ 1.0 (atan2 x.im x.re)))))
t_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 * atan2(x_46_im, x_46_re)) * pow(x_46_re, y_46_re);
double tmp;
if (y_46_re <= -1.35e-102) {
tmp = t_0;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * log(hypot(x_46_im, x_46_re));
} else if (y_46_re <= 0.011) {
tmp = sin((y_46_re * (1.0 / (1.0 / atan2(x_46_im, x_46_re)))));
} else {
tmp = t_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 * Math.atan2(x_46_im, x_46_re)) * Math.pow(x_46_re, y_46_re);
double tmp;
if (y_46_re <= -1.35e-102) {
tmp = t_0;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
} else if (y_46_re <= 0.011) {
tmp = Math.sin((y_46_re * (1.0 / (1.0 / Math.atan2(x_46_im, x_46_re)))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = (y_46_re * math.atan2(x_46_im, x_46_re)) * math.pow(x_46_re, y_46_re) tmp = 0 if y_46_re <= -1.35e-102: tmp = t_0 elif y_46_re <= 9.5e-247: tmp = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) elif y_46_re <= 0.011: tmp = math.sin((y_46_re * (1.0 / (1.0 / math.atan2(x_46_im, x_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(y_46_re * atan(x_46_im, x_46_re)) * (x_46_re ^ y_46_re)) tmp = 0.0 if (y_46_re <= -1.35e-102) tmp = t_0; elseif (y_46_re <= 9.5e-247) tmp = Float64(y_46_im * log(hypot(x_46_im, x_46_re))); elseif (y_46_re <= 0.011) tmp = sin(Float64(y_46_re * Float64(1.0 / Float64(1.0 / atan(x_46_im, x_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 = (y_46_re * atan2(x_46_im, x_46_re)) * (x_46_re ^ y_46_re); tmp = 0.0; if (y_46_re <= -1.35e-102) tmp = t_0; elseif (y_46_re <= 9.5e-247) tmp = y_46_im * log(hypot(x_46_im, x_46_re)); elseif (y_46_re <= 0.011) tmp = sin((y_46_re * (1.0 / (1.0 / atan2(x_46_im, x_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[(N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision] * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.35e-102], t$95$0, If[LessEqual[y$46$re, 9.5e-247], N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 0.011], N[Sin[N[(y$46$re * N[(1.0 / N[(1.0 / N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot {x.re}^{y.re}\\
\mathbf{if}\;y.re \leq -1.35 \cdot 10^{-102}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{elif}\;y.re \leq 0.011:\\
\;\;\;\;\sin \left(y.re \cdot \frac{1}{\frac{1}{\tan^{-1}_* \frac{x.im}{x.re}}}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -1.35e-102 or 0.010999999999999999 < y.re Initial program 40.9%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6463.4%
Simplified63.4%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6463.4%
Applied egg-rr63.4%
Taylor expanded in x.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6447.5%
Simplified47.5%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6445.4%
Simplified45.4%
if -1.35e-102 < y.re < 9.49999999999999939e-247Initial program 33.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6471.1%
Simplified71.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6436.4%
Simplified36.4%
if 9.49999999999999939e-247 < y.re < 0.010999999999999999Initial program 37.7%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6433.9%
Simplified33.9%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6433.9%
Applied egg-rr33.9%
Taylor expanded in y.re around 0
Simplified33.9%
Final simplification40.8%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))) (t_1 (* t_0 (pow x.re y.re))))
(if (<= y.re -1.36e-102)
t_1
(if (<= y.re 9.5e-247)
(* y.im (log (hypot x.im x.re)))
(if (<= y.re 3.1e-21) 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 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = t_0 * pow(x_46_re, y_46_re);
double tmp;
if (y_46_re <= -1.36e-102) {
tmp = t_1;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * log(hypot(x_46_im, x_46_re));
} else if (y_46_re <= 3.1e-21) {
tmp = t_0;
} else {
tmp = t_1;
}
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 * Math.atan2(x_46_im, x_46_re);
double t_1 = t_0 * Math.pow(x_46_re, y_46_re);
double tmp;
if (y_46_re <= -1.36e-102) {
tmp = t_1;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
} else if (y_46_re <= 3.1e-21) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = t_0 * math.pow(x_46_re, y_46_re) tmp = 0 if y_46_re <= -1.36e-102: tmp = t_1 elif y_46_re <= 9.5e-247: tmp = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) elif y_46_re <= 3.1e-21: tmp = t_0 else: tmp = t_1 return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(t_0 * (x_46_re ^ y_46_re)) tmp = 0.0 if (y_46_re <= -1.36e-102) tmp = t_1; elseif (y_46_re <= 9.5e-247) tmp = Float64(y_46_im * log(hypot(x_46_im, x_46_re))); elseif (y_46_re <= 3.1e-21) tmp = t_0; else tmp = t_1; end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); t_1 = t_0 * (x_46_re ^ y_46_re); tmp = 0.0; if (y_46_re <= -1.36e-102) tmp = t_1; elseif (y_46_re <= 9.5e-247) tmp = y_46_im * log(hypot(x_46_im, x_46_re)); elseif (y_46_re <= 3.1e-21) tmp = t_0; else tmp = t_1; end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := Block[{t$95$0 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.36e-102], t$95$1, If[LessEqual[y$46$re, 9.5e-247], N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 3.1e-21], t$95$0, t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := t\_0 \cdot {x.re}^{y.re}\\
\mathbf{if}\;y.re \leq -1.36 \cdot 10^{-102}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{elif}\;y.re \leq 3.1 \cdot 10^{-21}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y.re < -1.36000000000000001e-102 or 3.0999999999999998e-21 < y.re Initial program 40.6%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6463.0%
Simplified63.0%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f6463.0%
Applied egg-rr63.0%
Taylor expanded in x.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6447.2%
Simplified47.2%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6445.1%
Simplified45.1%
if -1.36000000000000001e-102 < y.re < 9.49999999999999939e-247Initial program 33.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6471.1%
Simplified71.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6436.4%
Simplified36.4%
if 9.49999999999999939e-247 < y.re < 3.0999999999999998e-21Initial program 38.4%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6434.5%
Simplified34.5%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6434.5%
Simplified34.5%
Final simplification40.8%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* (pow x.im y.re) (* y.im (log x.im)))))
(if (<= y.re -2.85e+87)
t_0
(if (<= y.re 9.5e-247)
(* y.im (log (hypot x.im x.re)))
(if (<= y.re 1.6e+35) (* y.re (atan2 x.im x.re)) t_0)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = pow(x_46_im, y_46_re) * (y_46_im * log(x_46_im));
double tmp;
if (y_46_re <= -2.85e+87) {
tmp = t_0;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * log(hypot(x_46_im, x_46_re));
} else if (y_46_re <= 1.6e+35) {
tmp = y_46_re * atan2(x_46_im, x_46_re);
} else {
tmp = t_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 = Math.pow(x_46_im, y_46_re) * (y_46_im * Math.log(x_46_im));
double tmp;
if (y_46_re <= -2.85e+87) {
tmp = t_0;
} else if (y_46_re <= 9.5e-247) {
tmp = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
} else if (y_46_re <= 1.6e+35) {
tmp = y_46_re * Math.atan2(x_46_im, x_46_re);
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = math.pow(x_46_im, y_46_re) * (y_46_im * math.log(x_46_im)) tmp = 0 if y_46_re <= -2.85e+87: tmp = t_0 elif y_46_re <= 9.5e-247: tmp = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) elif y_46_re <= 1.6e+35: tmp = y_46_re * math.atan2(x_46_im, x_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((x_46_im ^ y_46_re) * Float64(y_46_im * log(x_46_im))) tmp = 0.0 if (y_46_re <= -2.85e+87) tmp = t_0; elseif (y_46_re <= 9.5e-247) tmp = Float64(y_46_im * log(hypot(x_46_im, x_46_re))); elseif (y_46_re <= 1.6e+35) tmp = Float64(y_46_re * atan(x_46_im, x_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_im ^ y_46_re) * (y_46_im * log(x_46_im)); tmp = 0.0; if (y_46_re <= -2.85e+87) tmp = t_0; elseif (y_46_re <= 9.5e-247) tmp = y_46_im * log(hypot(x_46_im, x_46_re)); elseif (y_46_re <= 1.6e+35) tmp = y_46_re * atan2(x_46_im, x_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[(N[Power[x$46$im, y$46$re], $MachinePrecision] * N[(y$46$im * N[Log[x$46$im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -2.85e+87], t$95$0, If[LessEqual[y$46$re, 9.5e-247], N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.6e+35], N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x.im}^{y.re} \cdot \left(y.im \cdot \log x.im\right)\\
\mathbf{if}\;y.re \leq -2.85 \cdot 10^{+87}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 9.5 \cdot 10^{-247}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{elif}\;y.re \leq 1.6 \cdot 10^{+35}:\\
\;\;\;\;y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -2.85000000000000019e87 or 1.59999999999999991e35 < y.re Initial program 44.7%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified60.2%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.2%
Simplified60.2%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6460.2%
Applied egg-rr60.2%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6472.0%
Simplified72.0%
Taylor expanded in x.re around 0
associate-*r*N/A
remove-double-negN/A
log-recN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
pow-lowering-pow.f6434.2%
Simplified34.2%
if -2.85000000000000019e87 < y.re < 9.49999999999999939e-247Initial program 33.2%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6461.5%
Simplified61.5%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6427.2%
Simplified27.2%
if 9.49999999999999939e-247 < y.re < 1.59999999999999991e35Initial program 35.5%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6435.6%
Simplified35.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6430.5%
Simplified30.5%
Final simplification30.8%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re)))
(t_1 (* y.im (log (hypot x.im x.re)))))
(if (<= y.im -1.15e-185)
t_1
(if (<= y.im 3.75e-115)
(sin t_0)
(if (<= y.im 1950000000000.0) 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 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = y_46_im * log(hypot(x_46_im, x_46_re));
double tmp;
if (y_46_im <= -1.15e-185) {
tmp = t_1;
} else if (y_46_im <= 3.75e-115) {
tmp = sin(t_0);
} else if (y_46_im <= 1950000000000.0) {
tmp = t_1;
} else {
tmp = t_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 * Math.atan2(x_46_im, x_46_re);
double t_1 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double tmp;
if (y_46_im <= -1.15e-185) {
tmp = t_1;
} else if (y_46_im <= 3.75e-115) {
tmp = Math.sin(t_0);
} else if (y_46_im <= 1950000000000.0) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) tmp = 0 if y_46_im <= -1.15e-185: tmp = t_1 elif y_46_im <= 3.75e-115: tmp = math.sin(t_0) elif y_46_im <= 1950000000000.0: 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(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) tmp = 0.0 if (y_46_im <= -1.15e-185) tmp = t_1; elseif (y_46_im <= 3.75e-115) tmp = sin(t_0); elseif (y_46_im <= 1950000000000.0) tmp = t_1; 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 = y_46_re * atan2(x_46_im, x_46_re); t_1 = y_46_im * log(hypot(x_46_im, x_46_re)); tmp = 0.0; if (y_46_im <= -1.15e-185) tmp = t_1; elseif (y_46_im <= 3.75e-115) tmp = sin(t_0); elseif (y_46_im <= 1950000000000.0) tmp = t_1; 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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$im, -1.15e-185], t$95$1, If[LessEqual[y$46$im, 3.75e-115], N[Sin[t$95$0], $MachinePrecision], If[LessEqual[y$46$im, 1950000000000.0], t$95$1, t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{if}\;y.im \leq -1.15 \cdot 10^{-185}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 3.75 \cdot 10^{-115}:\\
\;\;\;\;\sin t\_0\\
\mathbf{elif}\;y.im \leq 1950000000000:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.im < -1.15e-185 or 3.75000000000000019e-115 < y.im < 1.95e12Initial program 39.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6452.6%
Simplified52.6%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6421.9%
Simplified21.9%
if -1.15e-185 < y.im < 3.75000000000000019e-115Initial program 52.6%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6481.7%
Simplified81.7%
Taylor expanded in y.re around 0
Simplified39.3%
if 1.95e12 < y.im Initial program 23.2%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6429.3%
Simplified29.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6410.7%
Simplified10.7%
Final simplification22.8%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re)))
(t_1 (* y.im (log (hypot x.im x.re)))))
(if (<= y.im -1e-182)
t_1
(if (<= y.im 4.2e-115) t_0 (if (<= y.im 2200000000000.0) 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 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = y_46_im * log(hypot(x_46_im, x_46_re));
double tmp;
if (y_46_im <= -1e-182) {
tmp = t_1;
} else if (y_46_im <= 4.2e-115) {
tmp = t_0;
} else if (y_46_im <= 2200000000000.0) {
tmp = t_1;
} else {
tmp = t_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 * Math.atan2(x_46_im, x_46_re);
double t_1 = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
double tmp;
if (y_46_im <= -1e-182) {
tmp = t_1;
} else if (y_46_im <= 4.2e-115) {
tmp = t_0;
} else if (y_46_im <= 2200000000000.0) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): t_0 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) tmp = 0 if y_46_im <= -1e-182: tmp = t_1 elif y_46_im <= 4.2e-115: tmp = t_0 elif y_46_im <= 2200000000000.0: 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(y_46_re * atan(x_46_im, x_46_re)) t_1 = Float64(y_46_im * log(hypot(x_46_im, x_46_re))) tmp = 0.0 if (y_46_im <= -1e-182) tmp = t_1; elseif (y_46_im <= 4.2e-115) tmp = t_0; elseif (y_46_im <= 2200000000000.0) tmp = t_1; 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 = y_46_re * atan2(x_46_im, x_46_re); t_1 = y_46_im * log(hypot(x_46_im, x_46_re)); tmp = 0.0; if (y_46_im <= -1e-182) tmp = t_1; elseif (y_46_im <= 4.2e-115) tmp = t_0; elseif (y_46_im <= 2200000000000.0) tmp = t_1; 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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$im, -1e-182], t$95$1, If[LessEqual[y$46$im, 4.2e-115], t$95$0, If[LessEqual[y$46$im, 2200000000000.0], t$95$1, t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{if}\;y.im \leq -1 \cdot 10^{-182}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.im \leq 4.2 \cdot 10^{-115}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.im \leq 2200000000000:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.im < -1e-182 or 4.20000000000000003e-115 < y.im < 2.2e12Initial program 40.2%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6453.0%
Simplified53.0%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6422.1%
Simplified22.1%
if -1e-182 < y.im < 4.20000000000000003e-115 or 2.2e12 < y.im Initial program 36.4%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6453.0%
Simplified53.0%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6423.7%
Simplified23.7%
(FPCore (x.re x.im y.re y.im) :precision binary64 (if (<= x.re 11.5) (* (pow x.im y.re) (* y.im (log x.im))) (* (pow x.re y.re) (* y.im (log x.re)))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if (x_46_re <= 11.5) {
tmp = pow(x_46_im, y_46_re) * (y_46_im * log(x_46_im));
} else {
tmp = pow(x_46_re, y_46_re) * (y_46_im * log(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 (x_46re <= 11.5d0) then
tmp = (x_46im ** y_46re) * (y_46im * log(x_46im))
else
tmp = (x_46re ** y_46re) * (y_46im * log(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 (x_46_re <= 11.5) {
tmp = Math.pow(x_46_im, y_46_re) * (y_46_im * Math.log(x_46_im));
} else {
tmp = Math.pow(x_46_re, y_46_re) * (y_46_im * Math.log(x_46_re));
}
return tmp;
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): tmp = 0 if x_46_re <= 11.5: tmp = math.pow(x_46_im, y_46_re) * (y_46_im * math.log(x_46_im)) else: tmp = math.pow(x_46_re, y_46_re) * (y_46_im * math.log(x_46_re)) return tmp
function code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = 0.0 if (x_46_re <= 11.5) tmp = Float64((x_46_im ^ y_46_re) * Float64(y_46_im * log(x_46_im))); else tmp = Float64((x_46_re ^ y_46_re) * Float64(y_46_im * log(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 (x_46_re <= 11.5) tmp = (x_46_im ^ y_46_re) * (y_46_im * log(x_46_im)); else tmp = (x_46_re ^ y_46_re) * (y_46_im * log(x_46_re)); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := If[LessEqual[x$46$re, 11.5], N[(N[Power[x$46$im, y$46$re], $MachinePrecision] * N[(y$46$im * N[Log[x$46$im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[x$46$re, y$46$re], $MachinePrecision] * N[(y$46$im * N[Log[x$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x.re \leq 11.5:\\
\;\;\;\;{x.im}^{y.re} \cdot \left(y.im \cdot \log x.im\right)\\
\mathbf{else}:\\
\;\;\;\;{x.re}^{y.re} \cdot \left(y.im \cdot \log x.re\right)\\
\end{array}
\end{array}
if x.re < 11.5Initial program 42.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified71.9%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6460.8%
Simplified60.8%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6460.8%
Applied egg-rr60.8%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6445.2%
Simplified45.2%
Taylor expanded in x.re around 0
associate-*r*N/A
remove-double-negN/A
log-recN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
pow-lowering-pow.f6421.7%
Simplified21.7%
if 11.5 < x.re Initial program 27.7%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified67.8%
Taylor expanded in y.re around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6462.0%
Simplified62.0%
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
atan2-lowering-atan2.f64N/A
/-lowering-/.f6462.0%
Applied egg-rr62.0%
Taylor expanded in y.im around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f6455.5%
Simplified55.5%
Taylor expanded in x.im around 0
associate-*r*N/A
remove-double-negN/A
log-recN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
pow-lowering-pow.f6454.2%
Simplified54.2%
Final simplification30.3%
(FPCore (x.re x.im y.re y.im) :precision binary64 (* y.re (atan2 x.im x.re)))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return y_46_re * atan2(x_46_im, x_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 = y_46re * atan2(x_46im, x_46re)
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 * Math.atan2(x_46_im, x_46_re);
}
def code(x_46_re, x_46_im, y_46_re, y_46_im): return y_46_re * math.atan2(x_46_im, x_46_re)
function code(x_46_re, x_46_im, y_46_re, y_46_im) return Float64(y_46_re * atan(x_46_im, x_46_re)) end
function tmp = code(x_46_re, x_46_im, y_46_re, y_46_im) tmp = y_46_re * atan2(x_46_im, x_46_re); end
code[x$46$re_, x$46$im_, y$46$re_, y$46$im_] := N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}
\end{array}
Initial program 38.3%
Taylor expanded in y.im around 0
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
unpow2N/A
unpow2N/A
hypot-defineN/A
hypot-lowering-hypot.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6444.6%
Simplified44.6%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6414.7%
Simplified14.7%
herbie shell --seed 2024164
(FPCore (x.re x.im y.re y.im)
:name "powComplex, imaginary part"
:precision binary64
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))