
(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 17 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 (* y.re (atan2 x.im x.re)))
(t_2 (sin t_1))
(t_3 (log (sqrt (+ (* x.re x.re) (* x.im x.im)))))
(t_4 (exp (- (* y.re t_3) t_0))))
(if (<= y.re -1.08e-47)
(* t_4 t_2)
(if (<= y.re 1e+110)
(/
(sin (+ t_1 (* y.im (log (hypot x.re x.im)))))
(/ (exp t_0) (pow (hypot x.re x.im) y.re)))
(* t_4 (+ t_2 (* t_3 (* y.im (cos t_1)))))))))
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 = y_46_re * atan2(x_46_im, x_46_re);
double t_2 = sin(t_1);
double t_3 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
double t_4 = exp(((y_46_re * t_3) - t_0));
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = t_4 * t_2;
} else if (y_46_re <= 1e+110) {
tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp(t_0) / pow(hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_4 * (t_2 + (t_3 * (y_46_im * cos(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.atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_2 = Math.sin(t_1);
double t_3 = Math.log(Math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))));
double t_4 = Math.exp(((y_46_re * t_3) - t_0));
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = t_4 * t_2;
} else if (y_46_re <= 1e+110) {
tmp = Math.sin((t_1 + (y_46_im * Math.log(Math.hypot(x_46_re, x_46_im))))) / (Math.exp(t_0) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_4 * (t_2 + (t_3 * (y_46_im * Math.cos(t_1))));
}
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 = y_46_re * math.atan2(x_46_im, x_46_re) t_2 = math.sin(t_1) t_3 = math.log(math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))) t_4 = math.exp(((y_46_re * t_3) - t_0)) tmp = 0 if y_46_re <= -1.08e-47: tmp = t_4 * t_2 elif y_46_re <= 1e+110: tmp = math.sin((t_1 + (y_46_im * math.log(math.hypot(x_46_re, x_46_im))))) / (math.exp(t_0) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) else: tmp = t_4 * (t_2 + (t_3 * (y_46_im * math.cos(t_1)))) 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 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_2 = sin(t_1) t_3 = log(sqrt(Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im)))) t_4 = exp(Float64(Float64(y_46_re * t_3) - t_0)) tmp = 0.0 if (y_46_re <= -1.08e-47) tmp = Float64(t_4 * t_2); elseif (y_46_re <= 1e+110) tmp = Float64(sin(Float64(t_1 + Float64(y_46_im * log(hypot(x_46_re, x_46_im))))) / Float64(exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re))); else tmp = Float64(t_4 * Float64(t_2 + Float64(t_3 * Float64(y_46_im * cos(t_1))))); 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 = y_46_re * atan2(x_46_im, x_46_re); t_2 = sin(t_1); t_3 = log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im)))); t_4 = exp(((y_46_re * t_3) - t_0)); tmp = 0.0; if (y_46_re <= -1.08e-47) tmp = t_4 * t_2; elseif (y_46_re <= 1e+110) tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re)); else tmp = t_4 * (t_2 + (t_3 * (y_46_im * cos(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[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Sin[t$95$1], $MachinePrecision]}, Block[{t$95$3 = N[Log[N[Sqrt[N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$4 = N[Exp[N[(N[(y$46$re * t$95$3), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -1.08e-47], N[(t$95$4 * t$95$2), $MachinePrecision], If[LessEqual[y$46$re, 1e+110], N[(N[Sin[N[(t$95$1 + N[(y$46$im * N[Log[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $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], N[(t$95$4 * N[(t$95$2 + N[(t$95$3 * N[(y$46$im * N[Cos[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im\\
t_1 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_2 := \sin t\_1\\
t_3 := \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right)\\
t_4 := e^{y.re \cdot t\_3 - t\_0}\\
\mathbf{if}\;y.re \leq -1.08 \cdot 10^{-47}:\\
\;\;\;\;t\_4 \cdot t\_2\\
\mathbf{elif}\;y.re \leq 10^{+110}:\\
\;\;\;\;\frac{\sin \left(t\_1 + y.im \cdot \log \left(\mathsf{hypot}\left(x.re, x.im\right)\right)\right)}{\frac{e^{t\_0}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{else}:\\
\;\;\;\;t\_4 \cdot \left(t\_2 + t\_3 \cdot \left(y.im \cdot \cos t\_1\right)\right)\\
\end{array}
\end{array}
if y.re < -1.08000000000000005e-47Initial program 45.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6489.4%
Simplified89.4%
if -1.08000000000000005e-47 < y.re < 1e110Initial program 42.5%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified82.0%
if 1e110 < y.re Initial program 30.4%
Taylor expanded in y.im around 0
+-lowering-+.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
log-lowering-log.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
Final simplification80.0%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* (atan2 x.im x.re) y.im)) (t_1 (* y.re (atan2 x.im x.re))))
(if (<= y.re -1.08e-47)
(*
(exp (- (* y.re (log (sqrt (+ (* x.re x.re) (* x.im x.im))))) t_0))
(sin t_1))
(if (<= y.re 5.5e+245)
(/
(sin (+ t_1 (* y.im (log (hypot x.re x.im)))))
(/ (exp t_0) (pow (hypot x.re x.im) y.re)))
(* t_1 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = y_46_re * atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * sin(t_1);
} else if (y_46_re <= 5.5e+245) {
tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp(t_0) / pow(hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_1 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = Math.atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = y_46_re * Math.atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = 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_1);
} else if (y_46_re <= 5.5e+245) {
tmp = Math.sin((t_1 + (y_46_im * Math.log(Math.hypot(x_46_re, x_46_im))))) / (Math.exp(t_0) / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else {
tmp = t_1 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
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 = y_46_re * math.atan2(x_46_im, x_46_re) tmp = 0 if y_46_re <= -1.08e-47: tmp = 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_1) elif y_46_re <= 5.5e+245: tmp = math.sin((t_1 + (y_46_im * math.log(math.hypot(x_46_re, x_46_im))))) / (math.exp(t_0) / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) else: tmp = t_1 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) 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 = Float64(y_46_re * atan(x_46_im, x_46_re)) tmp = 0.0 if (y_46_re <= -1.08e-47) tmp = 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)) * sin(t_1)); elseif (y_46_re <= 5.5e+245) tmp = Float64(sin(Float64(t_1 + Float64(y_46_im * log(hypot(x_46_re, x_46_im))))) / Float64(exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re))); else tmp = Float64(t_1 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = atan2(x_46_im, x_46_re) * y_46_im; t_1 = y_46_re * atan2(x_46_im, x_46_re); tmp = 0.0; if (y_46_re <= -1.08e-47) tmp = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * sin(t_1); elseif (y_46_re <= 5.5e+245) tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / (exp(t_0) / (hypot(x_46_re, x_46_im) ^ y_46_re)); else tmp = t_1 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_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[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.08e-47], 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[Sin[t$95$1], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 5.5e+245], N[(N[Sin[N[(t$95$1 + N[(y$46$im * N[Log[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $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], N[(t$95$1 * 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]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im\\
t_1 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -1.08 \cdot 10^{-47}:\\
\;\;\;\;e^{y.re \cdot \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) - t\_0} \cdot \sin t\_1\\
\mathbf{elif}\;y.re \leq 5.5 \cdot 10^{+245}:\\
\;\;\;\;\frac{\sin \left(t\_1 + y.im \cdot \log \left(\mathsf{hypot}\left(x.re, x.im\right)\right)\right)}{\frac{e^{t\_0}}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -1.08000000000000005e-47Initial program 45.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6489.4%
Simplified89.4%
if -1.08000000000000005e-47 < y.re < 5.4999999999999997e245Initial program 41.0%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified76.5%
if 5.4999999999999997e245 < y.re Initial program 21.4%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.9%
Simplified42.9%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.9%
Simplified42.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6457.1%
Simplified57.1%
Final simplification79.2%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* (atan2 x.im x.re) y.im))
(t_1 (* y.re (atan2 x.im x.re)))
(t_2
(*
(exp (- (* y.re (log (sqrt (+ (* x.re x.re) (* x.im x.im))))) t_0))
(sin t_1))))
(if (<= y.re -1.08e-47)
t_2
(if (<= y.re 2e+55)
(/ (sin (+ t_1 (* y.im (log (hypot x.re x.im))))) (exp t_0))
t_2))))
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 = y_46_re * atan2(x_46_im, x_46_re);
double t_2 = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * sin(t_1);
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = t_2;
} else if (y_46_re <= 2e+55) {
tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / exp(t_0);
} 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.atan2(x_46_im, x_46_re) * y_46_im;
double t_1 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_2 = 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_1);
double tmp;
if (y_46_re <= -1.08e-47) {
tmp = t_2;
} else if (y_46_re <= 2e+55) {
tmp = Math.sin((t_1 + (y_46_im * Math.log(Math.hypot(x_46_re, x_46_im))))) / Math.exp(t_0);
} else {
tmp = t_2;
}
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 = y_46_re * math.atan2(x_46_im, x_46_re) t_2 = 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_1) tmp = 0 if y_46_re <= -1.08e-47: tmp = t_2 elif y_46_re <= 2e+55: tmp = math.sin((t_1 + (y_46_im * math.log(math.hypot(x_46_re, x_46_im))))) / math.exp(t_0) else: tmp = t_2 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 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_2 = 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)) * sin(t_1)) tmp = 0.0 if (y_46_re <= -1.08e-47) tmp = t_2; elseif (y_46_re <= 2e+55) tmp = Float64(sin(Float64(t_1 + Float64(y_46_im * log(hypot(x_46_re, x_46_im))))) / exp(t_0)); 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 = atan2(x_46_im, x_46_re) * y_46_im; t_1 = y_46_re * atan2(x_46_im, x_46_re); t_2 = exp(((y_46_re * log(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))))) - t_0)) * sin(t_1); tmp = 0.0; if (y_46_re <= -1.08e-47) tmp = t_2; elseif (y_46_re <= 2e+55) tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / exp(t_0); 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[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]}, Block[{t$95$1 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = 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[Sin[t$95$1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.08e-47], t$95$2, If[LessEqual[y$46$re, 2e+55], N[(N[Sin[N[(t$95$1 + N[(y$46$im * N[Log[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Exp[t$95$0], $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im\\
t_1 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_2 := e^{y.re \cdot \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) - t\_0} \cdot \sin t\_1\\
\mathbf{if}\;y.re \leq -1.08 \cdot 10^{-47}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.re \leq 2 \cdot 10^{+55}:\\
\;\;\;\;\frac{\sin \left(t\_1 + y.im \cdot \log \left(\mathsf{hypot}\left(x.re, x.im\right)\right)\right)}{e^{t\_0}}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.re < -1.08000000000000005e-47 or 2.00000000000000002e55 < y.re Initial program 40.6%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6474.5%
Simplified74.5%
if -1.08000000000000005e-47 < y.re < 2.00000000000000002e55Initial program 41.8%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified82.7%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6482.7%
Simplified82.7%
Final simplification78.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))))
(if (<= y.re -8e-17)
(* (sin t_1) (pow (* t_0 t_0) (* y.re 0.25)))
(if (<= y.re 1.85e+59)
(/
(sin (+ t_1 (* y.im (log (hypot x.re x.im)))))
(exp (* (atan2 x.im x.re) y.im)))
(* t_1 (pow (sqrt t_0) y.re))))))
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 tmp;
if (y_46_re <= -8e-17) {
tmp = sin(t_1) * pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= 1.85e+59) {
tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / exp((atan2(x_46_im, x_46_re) * y_46_im));
} else {
tmp = t_1 * pow(sqrt(t_0), y_46_re);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double 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 tmp;
if (y_46_re <= -8e-17) {
tmp = Math.sin(t_1) * Math.pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= 1.85e+59) {
tmp = Math.sin((t_1 + (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));
} else {
tmp = t_1 * Math.pow(Math.sqrt(t_0), y_46_re);
}
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) tmp = 0 if y_46_re <= -8e-17: tmp = math.sin(t_1) * math.pow((t_0 * t_0), (y_46_re * 0.25)) elif y_46_re <= 1.85e+59: tmp = math.sin((t_1 + (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)) else: tmp = t_1 * math.pow(math.sqrt(t_0), y_46_re) 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)) tmp = 0.0 if (y_46_re <= -8e-17) tmp = Float64(sin(t_1) * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))); elseif (y_46_re <= 1.85e+59) tmp = Float64(sin(Float64(t_1 + Float64(y_46_im * log(hypot(x_46_re, x_46_im))))) / exp(Float64(atan(x_46_im, x_46_re) * y_46_im))); else tmp = Float64(t_1 * (sqrt(t_0) ^ y_46_re)); 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); tmp = 0.0; if (y_46_re <= -8e-17) tmp = sin(t_1) * ((t_0 * t_0) ^ (y_46_re * 0.25)); elseif (y_46_re <= 1.85e+59) tmp = sin((t_1 + (y_46_im * log(hypot(x_46_re, x_46_im))))) / exp((atan2(x_46_im, x_46_re) * y_46_im)); else tmp = t_1 * (sqrt(t_0) ^ y_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[(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]}, If[LessEqual[y$46$re, -8e-17], N[(N[Sin[t$95$1], $MachinePrecision] * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 1.85e+59], N[(N[Sin[N[(t$95$1 + N[(y$46$im * N[Log[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[Power[N[Sqrt[t$95$0], $MachinePrecision], y$46$re], $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}\\
\mathbf{if}\;y.re \leq -8 \cdot 10^{-17}:\\
\;\;\;\;\sin t\_1 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
\mathbf{elif}\;y.re \leq 1.85 \cdot 10^{+59}:\\
\;\;\;\;\frac{\sin \left(t\_1 + y.im \cdot \log \left(\mathsf{hypot}\left(x.re, x.im\right)\right)\right)}{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot {\left(\sqrt{t\_0}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -8.00000000000000057e-17Initial program 45.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.9%
Simplified81.9%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval81.9%
Applied egg-rr81.9%
if -8.00000000000000057e-17 < y.re < 1.84999999999999999e59Initial program 41.7%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified82.6%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6481.7%
Simplified81.7%
if 1.84999999999999999e59 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6454.5%
Simplified54.5%
Final simplification75.7%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im)))
(t_1 (sin (* y.re (atan2 x.im x.re))))
(t_2 (* t_1 (pow (* t_0 t_0) (* y.re 0.25))))
(t_3 (exp (* (atan2 x.im x.re) y.im))))
(if (<= y.re -38000.0)
t_2
(if (<= y.re -6.5e-135)
(/ t_1 (/ t_3 (pow (hypot x.re x.im) y.re)))
(if (<= y.re 3.55e-65)
(/ (sin (* y.im (log (hypot x.im x.re)))) t_3)
(if (<= y.re 1.18e+55)
(/ t_1 (/ t_3 (pow t_0 (/ y.re 2.0))))
t_2))))))
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 = sin((y_46_re * atan2(x_46_im, x_46_re)));
double t_2 = t_1 * pow((t_0 * t_0), (y_46_re * 0.25));
double t_3 = exp((atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_re <= -38000.0) {
tmp = t_2;
} else if (y_46_re <= -6.5e-135) {
tmp = t_1 / (t_3 / pow(hypot(x_46_re, x_46_im), y_46_re));
} else if (y_46_re <= 3.55e-65) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_3;
} else if (y_46_re <= 1.18e+55) {
tmp = t_1 / (t_3 / pow(t_0, (y_46_re / 2.0)));
} 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 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double t_1 = Math.sin((y_46_re * Math.atan2(x_46_im, x_46_re)));
double t_2 = t_1 * Math.pow((t_0 * t_0), (y_46_re * 0.25));
double t_3 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_re <= -38000.0) {
tmp = t_2;
} else if (y_46_re <= -6.5e-135) {
tmp = t_1 / (t_3 / Math.pow(Math.hypot(x_46_re, x_46_im), y_46_re));
} else if (y_46_re <= 3.55e-65) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)))) / t_3;
} else if (y_46_re <= 1.18e+55) {
tmp = t_1 / (t_3 / Math.pow(t_0, (y_46_re / 2.0)));
} else {
tmp = t_2;
}
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 = math.sin((y_46_re * math.atan2(x_46_im, x_46_re))) t_2 = t_1 * math.pow((t_0 * t_0), (y_46_re * 0.25)) t_3 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) tmp = 0 if y_46_re <= -38000.0: tmp = t_2 elif y_46_re <= -6.5e-135: tmp = t_1 / (t_3 / math.pow(math.hypot(x_46_re, x_46_im), y_46_re)) elif y_46_re <= 3.55e-65: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) / t_3 elif y_46_re <= 1.18e+55: tmp = t_1 / (t_3 / math.pow(t_0, (y_46_re / 2.0))) else: tmp = t_2 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 = sin(Float64(y_46_re * atan(x_46_im, x_46_re))) t_2 = Float64(t_1 * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))) t_3 = exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) tmp = 0.0 if (y_46_re <= -38000.0) tmp = t_2; elseif (y_46_re <= -6.5e-135) tmp = Float64(t_1 / Float64(t_3 / (hypot(x_46_re, x_46_im) ^ y_46_re))); elseif (y_46_re <= 3.55e-65) tmp = Float64(sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) / t_3); elseif (y_46_re <= 1.18e+55) tmp = Float64(t_1 / Float64(t_3 / (t_0 ^ Float64(y_46_re / 2.0)))); 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 = (x_46_re * x_46_re) + (x_46_im * x_46_im); t_1 = sin((y_46_re * atan2(x_46_im, x_46_re))); t_2 = t_1 * ((t_0 * t_0) ^ (y_46_re * 0.25)); t_3 = exp((atan2(x_46_im, x_46_re) * y_46_im)); tmp = 0.0; if (y_46_re <= -38000.0) tmp = t_2; elseif (y_46_re <= -6.5e-135) tmp = t_1 / (t_3 / (hypot(x_46_re, x_46_im) ^ y_46_re)); elseif (y_46_re <= 3.55e-65) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_3; elseif (y_46_re <= 1.18e+55) tmp = t_1 / (t_3 / (t_0 ^ (y_46_re / 2.0))); 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[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -38000.0], t$95$2, If[LessEqual[y$46$re, -6.5e-135], N[(t$95$1 / N[(t$95$3 / N[Power[N[Sqrt[x$46$re ^ 2 + x$46$im ^ 2], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 3.55e-65], 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$3), $MachinePrecision], If[LessEqual[y$46$re, 1.18e+55], N[(t$95$1 / N[(t$95$3 / N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
t_1 := \sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)\\
t_2 := t\_1 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
t_3 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}\\
\mathbf{if}\;y.re \leq -38000:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.re \leq -6.5 \cdot 10^{-135}:\\
\;\;\;\;\frac{t\_1}{\frac{t\_3}{{\left(\mathsf{hypot}\left(x.re, x.im\right)\right)}^{y.re}}}\\
\mathbf{elif}\;y.re \leq 3.55 \cdot 10^{-65}:\\
\;\;\;\;\frac{\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)}{t\_3}\\
\mathbf{elif}\;y.re \leq 1.18 \cdot 10^{+55}:\\
\;\;\;\;\frac{t\_1}{\frac{t\_3}{{t\_0}^{\left(\frac{y.re}{2}\right)}}}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.re < -38000 or 1.1799999999999999e55 < y.re Initial program 40.5%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.2%
Simplified69.2%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval69.2%
Applied egg-rr69.2%
if -38000 < y.re < -6.50000000000000056e-135Initial program 45.0%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified81.1%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6480.2%
Simplified80.2%
if -6.50000000000000056e-135 < y.re < 3.55000000000000014e-65Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.4%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.4%
Simplified85.4%
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.f6470.0%
Simplified70.0%
if 3.55000000000000014e-65 < y.re < 1.1799999999999999e55Initial program 47.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified70.9%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6448.9%
Simplified48.9%
frac-2negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-lowering-neg.f64N/A
neg-sub0N/A
metadata-evalN/A
metadata-evalN/A
--lowering--.f64N/A
metadata-evalN/A
metadata-evalN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sqrt-pow2N/A
Applied egg-rr63.4%
Final simplification68.0%
(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 (sin t_1))
(t_3 (* t_2 (pow (* t_0 t_0) (* y.re 0.25))))
(t_4 (exp (* (atan2 x.im x.re) y.im))))
(if (<= y.re -1.35)
t_3
(if (<= y.re -2.8e-137)
(/ t_1 t_4)
(if (<= y.re 5.2e-65)
(/ (sin (* y.im (log (hypot x.im x.re)))) t_4)
(if (<= y.re 1.18e+55)
(/ t_2 (/ t_4 (pow t_0 (/ y.re 2.0))))
t_3))))))
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 = sin(t_1);
double t_3 = t_2 * pow((t_0 * t_0), (y_46_re * 0.25));
double t_4 = exp((atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_re <= -1.35) {
tmp = t_3;
} else if (y_46_re <= -2.8e-137) {
tmp = t_1 / t_4;
} else if (y_46_re <= 5.2e-65) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_4;
} else if (y_46_re <= 1.18e+55) {
tmp = t_2 / (t_4 / pow(t_0, (y_46_re / 2.0)));
} 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 = (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 = Math.sin(t_1);
double t_3 = t_2 * Math.pow((t_0 * t_0), (y_46_re * 0.25));
double t_4 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
double tmp;
if (y_46_re <= -1.35) {
tmp = t_3;
} else if (y_46_re <= -2.8e-137) {
tmp = t_1 / t_4;
} else if (y_46_re <= 5.2e-65) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)))) / t_4;
} else if (y_46_re <= 1.18e+55) {
tmp = t_2 / (t_4 / Math.pow(t_0, (y_46_re / 2.0)));
} else {
tmp = t_3;
}
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 = math.sin(t_1) t_3 = t_2 * math.pow((t_0 * t_0), (y_46_re * 0.25)) t_4 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) tmp = 0 if y_46_re <= -1.35: tmp = t_3 elif y_46_re <= -2.8e-137: tmp = t_1 / t_4 elif y_46_re <= 5.2e-65: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) / t_4 elif y_46_re <= 1.18e+55: tmp = t_2 / (t_4 / math.pow(t_0, (y_46_re / 2.0))) else: tmp = t_3 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 = sin(t_1) t_3 = Float64(t_2 * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))) t_4 = exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) tmp = 0.0 if (y_46_re <= -1.35) tmp = t_3; elseif (y_46_re <= -2.8e-137) tmp = Float64(t_1 / t_4); elseif (y_46_re <= 5.2e-65) tmp = Float64(sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) / t_4); elseif (y_46_re <= 1.18e+55) tmp = Float64(t_2 / Float64(t_4 / (t_0 ^ Float64(y_46_re / 2.0)))); 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 = (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 = sin(t_1); t_3 = t_2 * ((t_0 * t_0) ^ (y_46_re * 0.25)); t_4 = exp((atan2(x_46_im, x_46_re) * y_46_im)); tmp = 0.0; if (y_46_re <= -1.35) tmp = t_3; elseif (y_46_re <= -2.8e-137) tmp = t_1 / t_4; elseif (y_46_re <= 5.2e-65) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_4; elseif (y_46_re <= 1.18e+55) tmp = t_2 / (t_4 / (t_0 ^ (y_46_re / 2.0))); 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[(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[Sin[t$95$1], $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y$46$re, -1.35], t$95$3, If[LessEqual[y$46$re, -2.8e-137], N[(t$95$1 / t$95$4), $MachinePrecision], If[LessEqual[y$46$re, 5.2e-65], 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$4), $MachinePrecision], If[LessEqual[y$46$re, 1.18e+55], N[(t$95$2 / N[(t$95$4 / N[Power[t$95$0, N[(y$46$re / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$3]]]]]]]]]
\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 := \sin t\_1\\
t_3 := t\_2 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
t_4 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}\\
\mathbf{if}\;y.re \leq -1.35:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y.re \leq -2.8 \cdot 10^{-137}:\\
\;\;\;\;\frac{t\_1}{t\_4}\\
\mathbf{elif}\;y.re \leq 5.2 \cdot 10^{-65}:\\
\;\;\;\;\frac{\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)}{t\_4}\\
\mathbf{elif}\;y.re \leq 1.18 \cdot 10^{+55}:\\
\;\;\;\;\frac{t\_2}{\frac{t\_4}{{t\_0}^{\left(\frac{y.re}{2}\right)}}}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if y.re < -1.3500000000000001 or 1.1799999999999999e55 < y.re Initial program 40.6%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.7%
Simplified69.7%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval69.7%
Applied egg-rr69.7%
if -1.3500000000000001 < y.re < -2.7999999999999999e-137Initial program 44.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.8%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.8%
Simplified78.8%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.5%
Simplified78.5%
if -2.7999999999999999e-137 < y.re < 5.20000000000000019e-65Initial program 39.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.4%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.4%
Simplified85.4%
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.f6470.0%
Simplified70.0%
if 5.20000000000000019e-65 < y.re < 1.1799999999999999e55Initial program 47.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified70.9%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6448.9%
Simplified48.9%
frac-2negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-lowering-neg.f64N/A
neg-sub0N/A
metadata-evalN/A
metadata-evalN/A
--lowering--.f64N/A
metadata-evalN/A
metadata-evalN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
sqrt-pow2N/A
Applied egg-rr63.4%
Final simplification67.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im)))
(t_1 (exp (* (atan2 x.im x.re) y.im)))
(t_2 (* y.re (atan2 x.im x.re))))
(if (<= y.re -1.65)
(* (sin t_2) (pow (* t_0 t_0) (* y.re 0.25)))
(if (<= y.re -1.45e-136)
(/ t_2 t_1)
(if (<= y.re 4.2e-87)
(/ (sin (* y.im (log (hypot x.im x.re)))) t_1)
(if (<= y.re 4.8e+58)
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(* t_2 (pow (sqrt t_0) y.re))))))))
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 = exp((atan2(x_46_im, x_46_re) * y_46_im));
double t_2 = y_46_re * atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -1.65) {
tmp = sin(t_2) * pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= -1.45e-136) {
tmp = t_2 / t_1;
} else if (y_46_re <= 4.2e-87) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1;
} else if (y_46_re <= 4.8e+58) {
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 = t_2 * pow(sqrt(t_0), y_46_re);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double t_1 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
double t_2 = y_46_re * Math.atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -1.65) {
tmp = Math.sin(t_2) * Math.pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= -1.45e-136) {
tmp = t_2 / t_1;
} else if (y_46_re <= 4.2e-87) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re)))) / t_1;
} else if (y_46_re <= 4.8e+58) {
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 = t_2 * Math.pow(Math.sqrt(t_0), y_46_re);
}
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 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) t_2 = y_46_re * math.atan2(x_46_im, x_46_re) tmp = 0 if y_46_re <= -1.65: tmp = math.sin(t_2) * math.pow((t_0 * t_0), (y_46_re * 0.25)) elif y_46_re <= -1.45e-136: tmp = t_2 / t_1 elif y_46_re <= 4.2e-87: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) / t_1 elif y_46_re <= 4.8e+58: 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 = t_2 * math.pow(math.sqrt(t_0), y_46_re) 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 = exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) t_2 = Float64(y_46_re * atan(x_46_im, x_46_re)) tmp = 0.0 if (y_46_re <= -1.65) tmp = Float64(sin(t_2) * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))); elseif (y_46_re <= -1.45e-136) tmp = Float64(t_2 / t_1); elseif (y_46_re <= 4.2e-87) tmp = Float64(sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1); elseif (y_46_re <= 4.8e+58) 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(t_2 * (sqrt(t_0) ^ y_46_re)); 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 = exp((atan2(x_46_im, x_46_re) * y_46_im)); t_2 = y_46_re * atan2(x_46_im, x_46_re); tmp = 0.0; if (y_46_re <= -1.65) tmp = sin(t_2) * ((t_0 * t_0) ^ (y_46_re * 0.25)); elseif (y_46_re <= -1.45e-136) tmp = t_2 / t_1; elseif (y_46_re <= 4.2e-87) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))) / t_1; elseif (y_46_re <= 4.8e+58) 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 = t_2 * (sqrt(t_0) ^ y_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[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.65], N[(N[Sin[t$95$2], $MachinePrecision] * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -1.45e-136], N[(t$95$2 / t$95$1), $MachinePrecision], If[LessEqual[y$46$re, 4.2e-87], 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$re, 4.8e+58], 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[(t$95$2 * N[Power[N[Sqrt[t$95$0], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
t_1 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}\\
t_2 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -1.65:\\
\;\;\;\;\sin t\_2 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
\mathbf{elif}\;y.re \leq -1.45 \cdot 10^{-136}:\\
\;\;\;\;\frac{t\_2}{t\_1}\\
\mathbf{elif}\;y.re \leq 4.2 \cdot 10^{-87}:\\
\;\;\;\;\frac{\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)}{t\_1}\\
\mathbf{elif}\;y.re \leq 4.8 \cdot 10^{+58}:\\
\;\;\;\;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}:\\
\;\;\;\;t\_2 \cdot {\left(\sqrt{t\_0}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -1.6499999999999999Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval83.0%
Applied egg-rr83.0%
if -1.6499999999999999 < y.re < -1.44999999999999997e-136Initial program 44.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.8%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.8%
Simplified78.8%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.5%
Simplified78.5%
if -1.44999999999999997e-136 < y.re < 4.20000000000000014e-87Initial program 40.1%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.5%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.5%
Simplified85.5%
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.f6471.6%
Simplified71.6%
if 4.20000000000000014e-87 < y.re < 4.8e58Initial program 41.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6438.9%
Simplified38.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6456.8%
Simplified56.8%
if 4.8e58 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6454.5%
Simplified54.5%
Final simplification70.1%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (+ (* x.re x.re) (* x.im x.im)))
(t_1 (exp (* (atan2 x.im x.re) y.im)))
(t_2 (* y.re (atan2 x.im x.re))))
(if (<= y.re -4.4)
(* (sin t_2) (pow (* t_0 t_0) (* y.re 0.25)))
(if (<= y.re -2.6e-137)
(/ t_2 t_1)
(if (<= y.re 1.35e-104)
(/ (* y.im (log (hypot x.im x.re))) t_1)
(if (<= y.re 7e+58)
(*
y.re
(* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(* t_2 (pow (sqrt t_0) y.re))))))))
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 = exp((atan2(x_46_im, x_46_re) * y_46_im));
double t_2 = y_46_re * atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -4.4) {
tmp = sin(t_2) * pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= -2.6e-137) {
tmp = t_2 / t_1;
} else if (y_46_re <= 1.35e-104) {
tmp = (y_46_im * log(hypot(x_46_im, x_46_re))) / t_1;
} else if (y_46_re <= 7e+58) {
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 = t_2 * pow(sqrt(t_0), y_46_re);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (x_46_re * x_46_re) + (x_46_im * x_46_im);
double t_1 = Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
double t_2 = y_46_re * Math.atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -4.4) {
tmp = Math.sin(t_2) * Math.pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= -2.6e-137) {
tmp = t_2 / t_1;
} else if (y_46_re <= 1.35e-104) {
tmp = (y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))) / t_1;
} else if (y_46_re <= 7e+58) {
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 = t_2 * Math.pow(Math.sqrt(t_0), y_46_re);
}
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 = math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) t_2 = y_46_re * math.atan2(x_46_im, x_46_re) tmp = 0 if y_46_re <= -4.4: tmp = math.sin(t_2) * math.pow((t_0 * t_0), (y_46_re * 0.25)) elif y_46_re <= -2.6e-137: tmp = t_2 / t_1 elif y_46_re <= 1.35e-104: tmp = (y_46_im * math.log(math.hypot(x_46_im, x_46_re))) / t_1 elif y_46_re <= 7e+58: 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 = t_2 * math.pow(math.sqrt(t_0), y_46_re) 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 = exp(Float64(atan(x_46_im, x_46_re) * y_46_im)) t_2 = Float64(y_46_re * atan(x_46_im, x_46_re)) tmp = 0.0 if (y_46_re <= -4.4) tmp = Float64(sin(t_2) * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))); elseif (y_46_re <= -2.6e-137) tmp = Float64(t_2 / t_1); elseif (y_46_re <= 1.35e-104) tmp = Float64(Float64(y_46_im * log(hypot(x_46_im, x_46_re))) / t_1); elseif (y_46_re <= 7e+58) 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(t_2 * (sqrt(t_0) ^ y_46_re)); 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 = exp((atan2(x_46_im, x_46_re) * y_46_im)); t_2 = y_46_re * atan2(x_46_im, x_46_re); tmp = 0.0; if (y_46_re <= -4.4) tmp = sin(t_2) * ((t_0 * t_0) ^ (y_46_re * 0.25)); elseif (y_46_re <= -2.6e-137) tmp = t_2 / t_1; elseif (y_46_re <= 1.35e-104) tmp = (y_46_im * log(hypot(x_46_im, x_46_re))) / t_1; elseif (y_46_re <= 7e+58) 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 = t_2 * (sqrt(t_0) ^ y_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[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -4.4], N[(N[Sin[t$95$2], $MachinePrecision] * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -2.6e-137], N[(t$95$2 / t$95$1), $MachinePrecision], If[LessEqual[y$46$re, 1.35e-104], N[(N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], If[LessEqual[y$46$re, 7e+58], 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[(t$95$2 * N[Power[N[Sqrt[t$95$0], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x.re \cdot x.re + x.im \cdot x.im\\
t_1 := e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}\\
t_2 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -4.4:\\
\;\;\;\;\sin t\_2 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
\mathbf{elif}\;y.re \leq -2.6 \cdot 10^{-137}:\\
\;\;\;\;\frac{t\_2}{t\_1}\\
\mathbf{elif}\;y.re \leq 1.35 \cdot 10^{-104}:\\
\;\;\;\;\frac{y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)}{t\_1}\\
\mathbf{elif}\;y.re \leq 7 \cdot 10^{+58}:\\
\;\;\;\;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}:\\
\;\;\;\;t\_2 \cdot {\left(\sqrt{t\_0}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -4.4000000000000004Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval83.0%
Applied egg-rr83.0%
if -4.4000000000000004 < y.re < -2.6e-137Initial program 44.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified79.8%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.8%
Simplified78.8%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.5%
Simplified78.5%
if -2.6e-137 < y.re < 1.3499999999999999e-104Initial program 40.7%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified86.4%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6486.4%
Simplified86.4%
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.f6472.4%
Simplified72.4%
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.f6470.9%
Simplified70.9%
if 1.3499999999999999e-104 < y.re < 6.9999999999999995e58Initial program 39.8%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6437.5%
Simplified37.5%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6454.8%
Simplified54.8%
if 6.9999999999999995e58 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6454.5%
Simplified54.5%
Final simplification69.7%
(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))))
(if (<= y.re -2.6)
(* (sin t_1) (pow (* t_0 t_0) (* y.re 0.25)))
(if (<= y.re 4.1e+58)
(* y.re (* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(* t_1 (pow (sqrt t_0) y.re))))))
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 tmp;
if (y_46_re <= -2.6) {
tmp = sin(t_1) * pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= 4.1e+58) {
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 = t_1 * pow(sqrt(t_0), y_46_re);
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x_46re * x_46re) + (x_46im * x_46im)
t_1 = y_46re * atan2(x_46im, x_46re)
if (y_46re <= (-2.6d0)) then
tmp = sin(t_1) * ((t_0 * t_0) ** (y_46re * 0.25d0))
else if (y_46re <= 4.1d+58) then
tmp = y_46re * (atan2(x_46im, x_46re) * exp((atan2(x_46im, x_46re) * (0.0d0 - y_46im))))
else
tmp = t_1 * (sqrt(t_0) ** y_46re)
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double 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 tmp;
if (y_46_re <= -2.6) {
tmp = Math.sin(t_1) * Math.pow((t_0 * t_0), (y_46_re * 0.25));
} else if (y_46_re <= 4.1e+58) {
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 = t_1 * Math.pow(Math.sqrt(t_0), y_46_re);
}
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) tmp = 0 if y_46_re <= -2.6: tmp = math.sin(t_1) * math.pow((t_0 * t_0), (y_46_re * 0.25)) elif y_46_re <= 4.1e+58: 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 = t_1 * math.pow(math.sqrt(t_0), y_46_re) 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)) tmp = 0.0 if (y_46_re <= -2.6) tmp = Float64(sin(t_1) * (Float64(t_0 * t_0) ^ Float64(y_46_re * 0.25))); elseif (y_46_re <= 4.1e+58) 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(t_1 * (sqrt(t_0) ^ y_46_re)); 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); tmp = 0.0; if (y_46_re <= -2.6) tmp = sin(t_1) * ((t_0 * t_0) ^ (y_46_re * 0.25)); elseif (y_46_re <= 4.1e+58) 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 = t_1 * (sqrt(t_0) ^ y_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[(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]}, If[LessEqual[y$46$re, -2.6], N[(N[Sin[t$95$1], $MachinePrecision] * N[Power[N[(t$95$0 * t$95$0), $MachinePrecision], N[(y$46$re * 0.25), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 4.1e+58], 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[(t$95$1 * N[Power[N[Sqrt[t$95$0], $MachinePrecision], y$46$re], $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}\\
\mathbf{if}\;y.re \leq -2.6:\\
\;\;\;\;\sin t\_1 \cdot {\left(t\_0 \cdot t\_0\right)}^{\left(y.re \cdot 0.25\right)}\\
\mathbf{elif}\;y.re \leq 4.1 \cdot 10^{+58}:\\
\;\;\;\;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}:\\
\;\;\;\;t\_1 \cdot {\left(\sqrt{t\_0}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -2.60000000000000009Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
pow1/2N/A
sqr-powN/A
pow-prod-downN/A
pow-unpowN/A
associate-/r/N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval83.0%
Applied egg-rr83.0%
if -2.60000000000000009 < y.re < 4.1e58Initial program 41.4%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6438.1%
Simplified38.1%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6458.7%
Simplified58.7%
if 4.1e58 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6454.5%
Simplified54.5%
Final simplification64.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0
(*
(* y.re (atan2 x.im x.re))
(pow (sqrt (+ (* x.re x.re) (* x.im x.im))) y.re))))
(if (<= y.re -1.8)
t_0
(if (<= y.re 4.6e+60)
(* y.re (* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
t_0))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (y_46_re * atan2(x_46_im, x_46_re)) * pow(sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))), y_46_re);
double tmp;
if (y_46_re <= -1.8) {
tmp = t_0;
} else if (y_46_re <= 4.6e+60) {
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 = 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 = (y_46re * atan2(x_46im, x_46re)) * (sqrt(((x_46re * x_46re) + (x_46im * x_46im))) ** y_46re)
if (y_46re <= (-1.8d0)) then
tmp = t_0
else if (y_46re <= 4.6d+60) then
tmp = y_46re * (atan2(x_46im, x_46re) * exp((atan2(x_46im, x_46re) * (0.0d0 - y_46im))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = (y_46_re * Math.atan2(x_46_im, x_46_re)) * Math.pow(Math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))), y_46_re);
double tmp;
if (y_46_re <= -1.8) {
tmp = t_0;
} else if (y_46_re <= 4.6e+60) {
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 = 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(math.sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))), y_46_re) tmp = 0 if y_46_re <= -1.8: tmp = t_0 elif y_46_re <= 4.6e+60: 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 = 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)) * (sqrt(Float64(Float64(x_46_re * x_46_re) + Float64(x_46_im * x_46_im))) ^ y_46_re)) tmp = 0.0 if (y_46_re <= -1.8) tmp = t_0; elseif (y_46_re <= 4.6e+60) 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 = 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)) * (sqrt(((x_46_re * x_46_re) + (x_46_im * x_46_im))) ^ y_46_re); tmp = 0.0; if (y_46_re <= -1.8) tmp = t_0; elseif (y_46_re <= 4.6e+60) 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 = 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[N[Sqrt[N[(N[(x$46$re * x$46$re), $MachinePrecision] + N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], y$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.8], t$95$0, If[LessEqual[y$46$re, 4.6e+60], 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], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right) \cdot {\left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right)}^{y.re}\\
\mathbf{if}\;y.re \leq -1.8:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 4.6 \cdot 10^{+60}:\\
\;\;\;\;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}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -1.80000000000000004 or 4.60000000000000034e60 < y.re Initial program 40.9%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.4%
Simplified69.4%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6467.8%
Simplified67.8%
if -1.80000000000000004 < y.re < 4.60000000000000034e60Initial program 41.4%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6438.1%
Simplified38.1%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6458.7%
Simplified58.7%
Final simplification63.2%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))))
(if (<= y.re -7.8)
(* (sin t_0) (pow x.re y.re))
(if (<= y.re 5.2e+59)
(* y.re (* (atan2 x.im x.re) (exp (* (atan2 x.im x.re) (- 0.0 y.im)))))
(* t_0 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -7.8) {
tmp = sin(t_0) * pow(x_46_re, y_46_re);
} else if (y_46_re <= 5.2e+59) {
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 = t_0 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: t_0
real(8) :: tmp
t_0 = y_46re * atan2(x_46im, x_46re)
if (y_46re <= (-7.8d0)) then
tmp = sin(t_0) * (x_46re ** y_46re)
else if (y_46re <= 5.2d+59) then
tmp = y_46re * (atan2(x_46im, x_46re) * exp((atan2(x_46im, x_46re) * (0.0d0 - y_46im))))
else
tmp = t_0 * ((x_46im + (((x_46re * x_46re) * 0.5d0) / x_46im)) ** y_46re)
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * Math.atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -7.8) {
tmp = Math.sin(t_0) * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 5.2e+59) {
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 = t_0 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
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) tmp = 0 if y_46_re <= -7.8: tmp = math.sin(t_0) * math.pow(x_46_re, y_46_re) elif y_46_re <= 5.2e+59: 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 = t_0 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) 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)) tmp = 0.0 if (y_46_re <= -7.8) tmp = Float64(sin(t_0) * (x_46_re ^ y_46_re)); elseif (y_46_re <= 5.2e+59) 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(t_0 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); tmp = 0.0; if (y_46_re <= -7.8) tmp = sin(t_0) * (x_46_re ^ y_46_re); elseif (y_46_re <= 5.2e+59) 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 = t_0 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -7.8], N[(N[Sin[t$95$0], $MachinePrecision] * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 5.2e+59], 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[(t$95$0 * 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]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -7.8:\\
\;\;\;\;\sin t\_0 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 5.2 \cdot 10^{+59}:\\
\;\;\;\;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}:\\
\;\;\;\;t\_0 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -7.79999999999999982Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
Taylor expanded in x.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f6473.1%
Simplified73.1%
if -7.79999999999999982 < y.re < 5.19999999999999999e59Initial program 41.4%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6438.1%
Simplified38.1%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
atan2-lowering-atan2.f6458.7%
Simplified58.7%
if 5.19999999999999999e59 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.3%
Simplified49.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6451.0%
Simplified51.0%
Final simplification60.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))) (t_1 (sin t_0)))
(if (<= y.re -0.034)
(* t_1 (pow x.re y.re))
(if (<= y.re -9e-139)
t_1
(if (<= y.re 1.95e-87)
(sin (* y.im (log (hypot x.im x.re))))
(if (<= y.re 2.4e+14)
t_1
(* t_0 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = sin(t_0);
double tmp;
if (y_46_re <= -0.034) {
tmp = t_1 * pow(x_46_re, y_46_re);
} else if (y_46_re <= -9e-139) {
tmp = t_1;
} else if (y_46_re <= 1.95e-87) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 2.4e+14) {
tmp = t_1;
} else {
tmp = t_0 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
return tmp;
}
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_1 = Math.sin(t_0);
double tmp;
if (y_46_re <= -0.034) {
tmp = t_1 * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= -9e-139) {
tmp = t_1;
} else if (y_46_re <= 1.95e-87) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 2.4e+14) {
tmp = t_1;
} else {
tmp = t_0 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
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 = math.sin(t_0) tmp = 0 if y_46_re <= -0.034: tmp = t_1 * math.pow(x_46_re, y_46_re) elif y_46_re <= -9e-139: tmp = t_1 elif y_46_re <= 1.95e-87: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) elif y_46_re <= 2.4e+14: tmp = t_1 else: tmp = t_0 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) 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 = sin(t_0) tmp = 0.0 if (y_46_re <= -0.034) tmp = Float64(t_1 * (x_46_re ^ y_46_re)); elseif (y_46_re <= -9e-139) tmp = t_1; elseif (y_46_re <= 1.95e-87) tmp = sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 2.4e+14) tmp = t_1; else tmp = Float64(t_0 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); t_1 = sin(t_0); tmp = 0.0; if (y_46_re <= -0.034) tmp = t_1 * (x_46_re ^ y_46_re); elseif (y_46_re <= -9e-139) tmp = t_1; elseif (y_46_re <= 1.95e-87) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 2.4e+14) tmp = t_1; else tmp = t_0 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, If[LessEqual[y$46$re, -0.034], N[(t$95$1 * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, -9e-139], t$95$1, If[LessEqual[y$46$re, 1.95e-87], 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, 2.4e+14], t$95$1, N[(t$95$0 * 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]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := \sin t\_0\\
\mathbf{if}\;y.re \leq -0.034:\\
\;\;\;\;t\_1 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq -9 \cdot 10^{-139}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 1.95 \cdot 10^{-87}:\\
\;\;\;\;\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
\mathbf{elif}\;y.re \leq 2.4 \cdot 10^{+14}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -0.034000000000000002Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
Taylor expanded in x.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f6473.1%
Simplified73.1%
if -0.034000000000000002 < y.re < -9.00000000000000046e-139 or 1.9499999999999999e-87 < y.re < 2.4e14Initial program 47.2%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified80.9%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6478.8%
Simplified78.8%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6440.5%
Simplified40.5%
if -9.00000000000000046e-139 < y.re < 1.9499999999999999e-87Initial program 40.1%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.5%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.5%
Simplified85.5%
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.f6471.6%
Simplified71.6%
Taylor expanded in y.im around 0
Simplified37.1%
if 2.4e14 < y.re Initial program 32.3%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.1%
Simplified50.1%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.9%
Simplified46.9%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6448.5%
Simplified48.5%
Final simplification50.4%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re)))
(t_1 (sin t_0))
(t_2 (* t_0 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))
(if (<= y.re -15.0)
t_2
(if (<= y.re -2.8e-138)
t_1
(if (<= y.re 1.45e-84)
(sin (* y.im (log (hypot x.im x.re))))
(if (<= y.re 2.4e+14) 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 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = sin(t_0);
double t_2 = t_0 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
double tmp;
if (y_46_re <= -15.0) {
tmp = t_2;
} else if (y_46_re <= -2.8e-138) {
tmp = t_1;
} else if (y_46_re <= 1.45e-84) {
tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 2.4e+14) {
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 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_1 = Math.sin(t_0);
double t_2 = t_0 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
double tmp;
if (y_46_re <= -15.0) {
tmp = t_2;
} else if (y_46_re <= -2.8e-138) {
tmp = t_1;
} else if (y_46_re <= 1.45e-84) {
tmp = Math.sin((y_46_im * Math.log(Math.hypot(x_46_im, x_46_re))));
} else if (y_46_re <= 2.4e+14) {
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 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = math.sin(t_0) t_2 = t_0 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) tmp = 0 if y_46_re <= -15.0: tmp = t_2 elif y_46_re <= -2.8e-138: tmp = t_1 elif y_46_re <= 1.45e-84: tmp = math.sin((y_46_im * math.log(math.hypot(x_46_im, x_46_re)))) elif y_46_re <= 2.4e+14: 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 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = sin(t_0) t_2 = Float64(t_0 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)) tmp = 0.0 if (y_46_re <= -15.0) tmp = t_2; elseif (y_46_re <= -2.8e-138) tmp = t_1; elseif (y_46_re <= 1.45e-84) tmp = sin(Float64(y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 2.4e+14) 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 = y_46_re * atan2(x_46_im, x_46_re); t_1 = sin(t_0); t_2 = t_0 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re); tmp = 0.0; if (y_46_re <= -15.0) tmp = t_2; elseif (y_46_re <= -2.8e-138) tmp = t_1; elseif (y_46_re <= 1.45e-84) tmp = sin((y_46_im * log(hypot(x_46_im, x_46_re)))); elseif (y_46_re <= 2.4e+14) 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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 * 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, -15.0], t$95$2, If[LessEqual[y$46$re, -2.8e-138], t$95$1, If[LessEqual[y$46$re, 1.45e-84], 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, 2.4e+14], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := \sin t\_0\\
t_2 := t\_0 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\mathbf{if}\;y.re \leq -15:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.re \leq -2.8 \cdot 10^{-138}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 1.45 \cdot 10^{-84}:\\
\;\;\;\;\sin \left(y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\right)\\
\mathbf{elif}\;y.re \leq 2.4 \cdot 10^{+14}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.re < -15 or 2.4e14 < y.re Initial program 38.9%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.3%
Simplified67.3%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.8%
Simplified59.8%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6458.2%
Simplified58.2%
if -15 < y.re < -2.80000000000000001e-138 or 1.4500000000000001e-84 < y.re < 2.4e14Initial program 48.3%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified81.3%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6477.3%
Simplified77.3%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6439.8%
Simplified39.8%
if -2.80000000000000001e-138 < y.re < 1.4500000000000001e-84Initial program 40.1%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.5%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.5%
Simplified85.5%
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.f6471.6%
Simplified71.6%
Taylor expanded in y.im around 0
Simplified37.1%
Final simplification48.5%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))))
(if (<= y.re -14.5)
(* (sin t_0) (pow x.re y.re))
(if (<= y.re 7.5e+59)
(/ t_0 (exp (* (atan2 x.im x.re) y.im)))
(* t_0 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))))
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -14.5) {
tmp = sin(t_0) * pow(x_46_re, y_46_re);
} else if (y_46_re <= 7.5e+59) {
tmp = t_0 / exp((atan2(x_46_im, x_46_re) * y_46_im));
} else {
tmp = t_0 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
return tmp;
}
real(8) function code(x_46re, x_46im, y_46re, y_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8), intent (in) :: y_46re
real(8), intent (in) :: y_46im
real(8) :: t_0
real(8) :: tmp
t_0 = y_46re * atan2(x_46im, x_46re)
if (y_46re <= (-14.5d0)) then
tmp = sin(t_0) * (x_46re ** y_46re)
else if (y_46re <= 7.5d+59) then
tmp = t_0 / exp((atan2(x_46im, x_46re) * y_46im))
else
tmp = t_0 * ((x_46im + (((x_46re * x_46re) * 0.5d0) / x_46im)) ** y_46re)
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double t_0 = y_46_re * Math.atan2(x_46_im, x_46_re);
double tmp;
if (y_46_re <= -14.5) {
tmp = Math.sin(t_0) * Math.pow(x_46_re, y_46_re);
} else if (y_46_re <= 7.5e+59) {
tmp = t_0 / Math.exp((Math.atan2(x_46_im, x_46_re) * y_46_im));
} else {
tmp = t_0 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
}
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) tmp = 0 if y_46_re <= -14.5: tmp = math.sin(t_0) * math.pow(x_46_re, y_46_re) elif y_46_re <= 7.5e+59: tmp = t_0 / math.exp((math.atan2(x_46_im, x_46_re) * y_46_im)) else: tmp = t_0 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) 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)) tmp = 0.0 if (y_46_re <= -14.5) tmp = Float64(sin(t_0) * (x_46_re ^ y_46_re)); elseif (y_46_re <= 7.5e+59) tmp = Float64(t_0 / exp(Float64(atan(x_46_im, x_46_re) * y_46_im))); else tmp = Float64(t_0 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)); end return tmp end
function tmp_2 = code(x_46_re, x_46_im, y_46_re, y_46_im) t_0 = y_46_re * atan2(x_46_im, x_46_re); tmp = 0.0; if (y_46_re <= -14.5) tmp = sin(t_0) * (x_46_re ^ y_46_re); elseif (y_46_re <= 7.5e+59) tmp = t_0 / exp((atan2(x_46_im, x_46_re) * y_46_im)); else tmp = t_0 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -14.5], N[(N[Sin[t$95$0], $MachinePrecision] * N[Power[x$46$re, y$46$re], $MachinePrecision]), $MachinePrecision], If[LessEqual[y$46$re, 7.5e+59], N[(t$95$0 / N[Exp[N[(N[ArcTan[x$46$im / x$46$re], $MachinePrecision] * y$46$im), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$0 * 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]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -14.5:\\
\;\;\;\;\sin t\_0 \cdot {x.re}^{y.re}\\
\mathbf{elif}\;y.re \leq 7.5 \cdot 10^{+59}:\\
\;\;\;\;\frac{t\_0}{e^{\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im}}\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\end{array}
\end{array}
if y.re < -14.5Initial program 45.7%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.0%
Simplified83.0%
Taylor expanded in x.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f6473.1%
Simplified73.1%
if -14.5 < y.re < 7.4999999999999996e59Initial program 41.4%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified81.9%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6458.7%
Simplified58.7%
Taylor expanded in y.re around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6458.7%
Simplified58.7%
if 7.4999999999999996e59 < y.re Initial program 35.1%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.3%
Simplified49.3%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6451.0%
Simplified51.0%
Final simplification60.9%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re)))
(t_1 (sin t_0))
(t_2 (* t_0 (pow (+ x.im (/ (* (* x.re x.re) 0.5) x.im)) y.re))))
(if (<= y.re -15.0)
t_2
(if (<= y.re -1.9e-136)
t_1
(if (<= y.re 5e-133)
(* y.im (log (hypot x.im x.re)))
(if (<= y.re 2.4e+14) 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 = y_46_re * atan2(x_46_im, x_46_re);
double t_1 = sin(t_0);
double t_2 = t_0 * pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
double tmp;
if (y_46_re <= -15.0) {
tmp = t_2;
} else if (y_46_re <= -1.9e-136) {
tmp = t_1;
} else if (y_46_re <= 5e-133) {
tmp = y_46_im * log(hypot(x_46_im, x_46_re));
} else if (y_46_re <= 2.4e+14) {
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 = y_46_re * Math.atan2(x_46_im, x_46_re);
double t_1 = Math.sin(t_0);
double t_2 = t_0 * Math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re);
double tmp;
if (y_46_re <= -15.0) {
tmp = t_2;
} else if (y_46_re <= -1.9e-136) {
tmp = t_1;
} else if (y_46_re <= 5e-133) {
tmp = y_46_im * Math.log(Math.hypot(x_46_im, x_46_re));
} else if (y_46_re <= 2.4e+14) {
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 = y_46_re * math.atan2(x_46_im, x_46_re) t_1 = math.sin(t_0) t_2 = t_0 * math.pow((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)), y_46_re) tmp = 0 if y_46_re <= -15.0: tmp = t_2 elif y_46_re <= -1.9e-136: tmp = t_1 elif y_46_re <= 5e-133: tmp = y_46_im * math.log(math.hypot(x_46_im, x_46_re)) elif y_46_re <= 2.4e+14: 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 = Float64(y_46_re * atan(x_46_im, x_46_re)) t_1 = sin(t_0) t_2 = Float64(t_0 * (Float64(x_46_im + Float64(Float64(Float64(x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re)) tmp = 0.0 if (y_46_re <= -15.0) tmp = t_2; elseif (y_46_re <= -1.9e-136) tmp = t_1; elseif (y_46_re <= 5e-133) tmp = Float64(y_46_im * log(hypot(x_46_im, x_46_re))); elseif (y_46_re <= 2.4e+14) 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 = y_46_re * atan2(x_46_im, x_46_re); t_1 = sin(t_0); t_2 = t_0 * ((x_46_im + (((x_46_re * x_46_re) * 0.5) / x_46_im)) ^ y_46_re); tmp = 0.0; if (y_46_re <= -15.0) tmp = t_2; elseif (y_46_re <= -1.9e-136) tmp = t_1; elseif (y_46_re <= 5e-133) tmp = y_46_im * log(hypot(x_46_im, x_46_re)); elseif (y_46_re <= 2.4e+14) 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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 * 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, -15.0], t$95$2, If[LessEqual[y$46$re, -1.9e-136], t$95$1, If[LessEqual[y$46$re, 5e-133], 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.4e+14], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
t_1 := \sin t\_0\\
t_2 := t\_0 \cdot {\left(x.im + \frac{\left(x.re \cdot x.re\right) \cdot 0.5}{x.im}\right)}^{y.re}\\
\mathbf{if}\;y.re \leq -15:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y.re \leq -1.9 \cdot 10^{-136}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y.re \leq 5 \cdot 10^{-133}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{elif}\;y.re \leq 2.4 \cdot 10^{+14}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y.re < -15 or 2.4e14 < y.re Initial program 38.9%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.3%
Simplified67.3%
Taylor expanded in x.re around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.8%
Simplified59.8%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6458.2%
Simplified58.2%
if -15 < y.re < -1.9000000000000001e-136 or 4.9999999999999999e-133 < y.re < 2.4e14Initial program 48.4%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified81.1%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6477.5%
Simplified77.5%
Taylor expanded in y.im around 0
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6438.9%
Simplified38.9%
if -1.9000000000000001e-136 < y.re < 4.9999999999999999e-133Initial program 39.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.8%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.8%
Simplified85.8%
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.f6472.6%
Simplified72.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.f6436.0%
Simplified36.0%
Final simplification48.0%
(FPCore (x.re x.im y.re y.im)
:precision binary64
(let* ((t_0 (* y.re (atan2 x.im x.re))))
(if (<= y.re -1.5e-135)
t_0
(if (<= y.re 5e-133) (* y.im (log (hypot 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);
double tmp;
if (y_46_re <= -1.5e-135) {
tmp = t_0;
} else if (y_46_re <= 5e-133) {
tmp = y_46_im * log(hypot(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);
double tmp;
if (y_46_re <= -1.5e-135) {
tmp = t_0;
} else if (y_46_re <= 5e-133) {
tmp = y_46_im * Math.log(Math.hypot(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) tmp = 0 if y_46_re <= -1.5e-135: tmp = t_0 elif y_46_re <= 5e-133: tmp = y_46_im * math.log(math.hypot(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(y_46_re * atan(x_46_im, x_46_re)) tmp = 0.0 if (y_46_re <= -1.5e-135) tmp = t_0; elseif (y_46_re <= 5e-133) tmp = Float64(y_46_im * log(hypot(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); tmp = 0.0; if (y_46_re <= -1.5e-135) tmp = t_0; elseif (y_46_re <= 5e-133) tmp = y_46_im * log(hypot(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[(y$46$re * N[ArcTan[x$46$im / x$46$re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y$46$re, -1.5e-135], t$95$0, If[LessEqual[y$46$re, 5e-133], N[(y$46$im * N[Log[N[Sqrt[x$46$im ^ 2 + x$46$re ^ 2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\\
\mathbf{if}\;y.re \leq -1.5 \cdot 10^{-135}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y.re \leq 5 \cdot 10^{-133}:\\
\;\;\;\;y.im \cdot \log \left(\mathsf{hypot}\left(x.im, x.re\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y.re < -1.50000000000000006e-135 or 4.9999999999999999e-133 < y.re Initial program 41.8%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.2%
Simplified55.2%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6417.6%
Simplified17.6%
if -1.50000000000000006e-135 < y.re < 4.9999999999999999e-133Initial program 39.6%
exp-diffN/A
associate-*l/N/A
associate-/l*N/A
*-commutativeN/A
associate-/r/N/A
exp-diffN/A
Simplified85.8%
Taylor expanded in y.re around 0
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f6485.8%
Simplified85.8%
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.f6472.6%
Simplified72.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.f6436.0%
Simplified36.0%
(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 41.2%
Taylor expanded in y.im around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
atan2-lowering-atan2.f64N/A
pow-lowering-pow.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6443.5%
Simplified43.5%
Taylor expanded in y.re around 0
*-lowering-*.f64N/A
atan2-lowering-atan2.f6416.7%
Simplified16.7%
herbie shell --seed 2024144
(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)))))