
(FPCore (x.re x.im) :precision binary64 (- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))
double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
real(8) function code(x_46re, x_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
code = (((x_46re * x_46re) - (x_46im * x_46im)) * x_46re) - (((x_46re * x_46im) + (x_46im * x_46re)) * x_46im)
end function
public static double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
def code(x_46_re, x_46_im): return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im)
function code(x_46_re, x_46_im) return Float64(Float64(Float64(Float64(x_46_re * x_46_re) - Float64(x_46_im * x_46_im)) * x_46_re) - Float64(Float64(Float64(x_46_re * x_46_im) + Float64(x_46_im * x_46_re)) * x_46_im)) end
function tmp = code(x_46_re, x_46_im) tmp = (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im); end
code[x$46$re_, x$46$im_] := N[(N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] - N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision] * x$46$re), $MachinePrecision] - N[(N[(N[(x$46$re * x$46$im), $MachinePrecision] + N[(x$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] * x$46$im), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x.re x.im) :precision binary64 (- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))
double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
real(8) function code(x_46re, x_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
code = (((x_46re * x_46re) - (x_46im * x_46im)) * x_46re) - (((x_46re * x_46im) + (x_46im * x_46re)) * x_46im)
end function
public static double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
def code(x_46_re, x_46_im): return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im)
function code(x_46_re, x_46_im) return Float64(Float64(Float64(Float64(x_46_re * x_46_re) - Float64(x_46_im * x_46_im)) * x_46_re) - Float64(Float64(Float64(x_46_re * x_46_im) + Float64(x_46_im * x_46_re)) * x_46_im)) end
function tmp = code(x_46_re, x_46_im) tmp = (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im); end
code[x$46$re_, x$46$im_] := N[(N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] - N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision] * x$46$re), $MachinePrecision] - N[(N[(N[(x$46$re * x$46$im), $MachinePrecision] + N[(x$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] * x$46$im), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\end{array}
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(let* ((t_0 (* x.re_m (- x.re_m x.im_m))))
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
1e+177)
(fma (* x.re_m (+ x.im_m x.im_m)) (- x.im_m) (* (+ x.re_m x.im_m) t_0))
(fma
(+ x.im_m x.im_m)
(- x.im_m)
(* t_0 (fma (/ 1.0 x.im_m) (* x.re_m x.im_m) x.im_m)))))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double t_0 = x_46_re_m * (x_46_re_m - x_46_im_m);
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= 1e+177) {
tmp = fma((x_46_re_m * (x_46_im_m + x_46_im_m)), -x_46_im_m, ((x_46_re_m + x_46_im_m) * t_0));
} else {
tmp = fma((x_46_im_m + x_46_im_m), -x_46_im_m, (t_0 * fma((1.0 / x_46_im_m), (x_46_re_m * x_46_im_m), x_46_im_m)));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) t_0 = Float64(x_46_re_m * Float64(x_46_re_m - x_46_im_m)) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= 1e+177) tmp = fma(Float64(x_46_re_m * Float64(x_46_im_m + x_46_im_m)), Float64(-x_46_im_m), Float64(Float64(x_46_re_m + x_46_im_m) * t_0)); else tmp = fma(Float64(x_46_im_m + x_46_im_m), Float64(-x_46_im_m), Float64(t_0 * fma(Float64(1.0 / x_46_im_m), Float64(x_46_re_m * x_46_im_m), x_46_im_m))); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := Block[{t$95$0 = N[(x$46$re$95$m * N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision]), $MachinePrecision]}, N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e+177], N[(N[(x$46$re$95$m * N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] * (-x$46$im$95$m) + N[(N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision] * (-x$46$im$95$m) + N[(t$95$0 * N[(N[(1.0 / x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
\begin{array}{l}
t_0 := x.re\_m \cdot \left(x.re\_m - x.im\_m\right)\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq 10^{+177}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m \cdot \left(x.im\_m + x.im\_m\right), -x.im\_m, \left(x.re\_m + x.im\_m\right) \cdot t\_0\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.im\_m + x.im\_m, -x.im\_m, t\_0 \cdot \mathsf{fma}\left(\frac{1}{x.im\_m}, x.re\_m \cdot x.im\_m, x.im\_m\right)\right)\\
\end{array}
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < 1e177Initial program 92.0%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8
Applied egg-rr99.8%
if 1e177 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 61.6%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6480.4
Applied egg-rr80.4%
Taylor expanded in x.im around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6477.3
Simplified77.3%
distribute-lft-inN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
+-lowering-+.f6484.3
Applied egg-rr84.3%
associate-*r/N/A
div-invN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6478.1
Applied egg-rr78.1%
Final simplification92.5%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(let* ((t_0
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))))
(*
x.re_s
(if (<= t_0 -2e-161)
(* -3.0 (* x.im_m (* x.re_m x.im_m)))
(if (<= t_0 1e-99)
(* x.re_m (fma x.re_m x.re_m (* (* x.im_m x.im_m) -3.0)))
(fma
(- x.re_m x.im_m)
(* x.re_m (+ x.re_m x.im_m))
(+ x.im_m x.im_m)))))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double t_0 = (x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)));
double tmp;
if (t_0 <= -2e-161) {
tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m));
} else if (t_0 <= 1e-99) {
tmp = x_46_re_m * fma(x_46_re_m, x_46_re_m, ((x_46_im_m * x_46_im_m) * -3.0));
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) t_0 = Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) tmp = 0.0 if (t_0 <= -2e-161) tmp = Float64(-3.0 * Float64(x_46_im_m * Float64(x_46_re_m * x_46_im_m))); elseif (t_0 <= 1e-99) tmp = Float64(x_46_re_m * fma(x_46_re_m, x_46_re_m, Float64(Float64(x_46_im_m * x_46_im_m) * -3.0))); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := Block[{t$95$0 = N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(x$46$re$95$s * If[LessEqual[t$95$0, -2e-161], N[(-3.0 * N[(x$46$im$95$m * N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e-99], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m + N[(N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
\begin{array}{l}
t_0 := x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right)\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-161}:\\
\;\;\;\;-3 \cdot \left(x.im\_m \cdot \left(x.re\_m \cdot x.im\_m\right)\right)\\
\mathbf{elif}\;t\_0 \leq 10^{-99}:\\
\;\;\;\;x.re\_m \cdot \mathsf{fma}\left(x.re\_m, x.re\_m, \left(x.im\_m \cdot x.im\_m\right) \cdot -3\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.00000000000000006e-161Initial program 86.2%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval42.5
Simplified42.5%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6456.1
Applied egg-rr56.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6456.1
Applied egg-rr56.1%
if -2.00000000000000006e-161 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < 1e-99Initial program 99.8%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
unpow2N/A
accelerator-lowering-fma.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval99.8
Simplified99.8%
if 1e-99 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 69.0%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6484.1
Applied egg-rr84.1%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr82.4%
Final simplification75.9%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(let* ((t_0
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))))
(*
x.re_s
(if (<= t_0 -2e-316)
(* -3.0 (* x.im_m (* x.re_m x.im_m)))
(if (<= t_0 1e-99)
(* x.re_m (* x.re_m x.re_m))
(fma
(- x.re_m x.im_m)
(* x.re_m (+ x.re_m x.im_m))
(+ x.im_m x.im_m)))))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double t_0 = (x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)));
double tmp;
if (t_0 <= -2e-316) {
tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m));
} else if (t_0 <= 1e-99) {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) t_0 = Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) tmp = 0.0 if (t_0 <= -2e-316) tmp = Float64(-3.0 * Float64(x_46_im_m * Float64(x_46_re_m * x_46_im_m))); elseif (t_0 <= 1e-99) tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := Block[{t$95$0 = N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(x$46$re$95$s * If[LessEqual[t$95$0, -2e-316], N[(-3.0 * N[(x$46$im$95$m * N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e-99], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
\begin{array}{l}
t_0 := x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right)\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-316}:\\
\;\;\;\;-3 \cdot \left(x.im\_m \cdot \left(x.re\_m \cdot x.im\_m\right)\right)\\
\mathbf{elif}\;t\_0 \leq 10^{-99}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.000000017e-316Initial program 87.4%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval44.2
Simplified44.2%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
if -2.000000017e-316 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < 1e-99Initial program 99.9%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6479.0
Simplified79.0%
if 1e-99 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 69.0%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6484.1
Applied egg-rr84.1%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr82.4%
Final simplification71.2%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
1e-99)
(fma
(* x.re_m (+ x.im_m x.im_m))
(- x.im_m)
(* (+ x.re_m x.im_m) (* x.re_m (- x.re_m x.im_m))))
(fma (- x.re_m x.im_m) (* x.re_m (+ x.re_m x.im_m)) (+ x.im_m x.im_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= 1e-99) {
tmp = fma((x_46_re_m * (x_46_im_m + x_46_im_m)), -x_46_im_m, ((x_46_re_m + x_46_im_m) * (x_46_re_m * (x_46_re_m - x_46_im_m))));
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= 1e-99) tmp = fma(Float64(x_46_re_m * Float64(x_46_im_m + x_46_im_m)), Float64(-x_46_im_m), Float64(Float64(x_46_re_m + x_46_im_m) * Float64(x_46_re_m * Float64(x_46_re_m - x_46_im_m)))); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e-99], N[(N[(x$46$re$95$m * N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] * (-x$46$im$95$m) + N[(N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq 10^{-99}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m \cdot \left(x.im\_m + x.im\_m\right), -x.im\_m, \left(x.re\_m + x.im\_m\right) \cdot \left(x.re\_m \cdot \left(x.re\_m - x.im\_m\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < 1e-99Initial program 90.9%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8
Applied egg-rr99.8%
if 1e-99 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 69.0%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6484.1
Applied egg-rr84.1%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr82.4%
Final simplification92.5%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
1e-99)
(-
(* (+ x.re_m x.im_m) (* x.re_m (- x.re_m x.im_m)))
(* x.im_m (* x.re_m (+ x.im_m x.im_m))))
(fma (- x.re_m x.im_m) (* x.re_m (+ x.re_m x.im_m)) (+ x.im_m x.im_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= 1e-99) {
tmp = ((x_46_re_m + x_46_im_m) * (x_46_re_m * (x_46_re_m - x_46_im_m))) - (x_46_im_m * (x_46_re_m * (x_46_im_m + x_46_im_m)));
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= 1e-99) tmp = Float64(Float64(Float64(x_46_re_m + x_46_im_m) * Float64(x_46_re_m * Float64(x_46_re_m - x_46_im_m))) - Float64(x_46_im_m * Float64(x_46_re_m * Float64(x_46_im_m + x_46_im_m)))); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e-99], N[(N[(N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(x$46$re$95$m * N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq 10^{-99}:\\
\;\;\;\;\left(x.re\_m + x.im\_m\right) \cdot \left(x.re\_m \cdot \left(x.re\_m - x.im\_m\right)\right) - x.im\_m \cdot \left(x.re\_m \cdot \left(x.im\_m + x.im\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < 1e-99Initial program 90.9%
--lowering--.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
if 1e-99 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 69.0%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6484.1
Applied egg-rr84.1%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr82.4%
Final simplification92.5%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-2e-316)
(* -3.0 (* x.im_m (* x.re_m x.im_m)))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m));
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-2d-316)) then
tmp = (-3.0d0) * (x_46im_m * (x_46re_m * x_46im_m))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m));
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316: tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m)) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = Float64(-3.0 * Float64(x_46_im_m * Float64(x_46_re_m * x_46_im_m))); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = -3.0 * (x_46_im_m * (x_46_re_m * x_46_im_m)); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -2e-316], N[(-3.0 * N[(x$46$im$95$m * N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -2 \cdot 10^{-316}:\\
\;\;\;\;-3 \cdot \left(x.im\_m \cdot \left(x.re\_m \cdot x.im\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.000000017e-316Initial program 87.4%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval44.2
Simplified44.2%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
if -2.000000017e-316 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 77.6%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.1
Simplified65.1%
Final simplification61.6%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-2e-316)
(* (* x.re_m x.im_m) (* x.im_m -3.0))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = (x_46_re_m * x_46_im_m) * (x_46_im_m * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-2d-316)) then
tmp = (x_46re_m * x_46im_m) * (x_46im_m * (-3.0d0))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = (x_46_re_m * x_46_im_m) * (x_46_im_m * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316: tmp = (x_46_re_m * x_46_im_m) * (x_46_im_m * -3.0) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = Float64(Float64(x_46_re_m * x_46_im_m) * Float64(x_46_im_m * -3.0)); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = (x_46_re_m * x_46_im_m) * (x_46_im_m * -3.0); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -2e-316], N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] * N[(x$46$im$95$m * -3.0), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -2 \cdot 10^{-316}:\\
\;\;\;\;\left(x.re\_m \cdot x.im\_m\right) \cdot \left(x.im\_m \cdot -3\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.000000017e-316Initial program 87.4%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval44.2
Simplified44.2%
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
if -2.000000017e-316 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 77.6%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.1
Simplified65.1%
Final simplification61.6%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-2e-316)
(* x.im_m (* (* x.re_m x.im_m) -3.0))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = x_46_im_m * ((x_46_re_m * x_46_im_m) * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-2d-316)) then
tmp = x_46im_m * ((x_46re_m * x_46im_m) * (-3.0d0))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = x_46_im_m * ((x_46_re_m * x_46_im_m) * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316: tmp = x_46_im_m * ((x_46_re_m * x_46_im_m) * -3.0) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) * -3.0)); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = x_46_im_m * ((x_46_re_m * x_46_im_m) * -3.0); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -2e-316], N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -2 \cdot 10^{-316}:\\
\;\;\;\;x.im\_m \cdot \left(\left(x.re\_m \cdot x.im\_m\right) \cdot -3\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.000000017e-316Initial program 87.4%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval44.2
Simplified44.2%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6456.6
Applied egg-rr56.6%
if -2.000000017e-316 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 77.6%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.1
Simplified65.1%
Final simplification61.6%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-2e-316)
(* x.re_m (* (* x.im_m x.im_m) -3.0))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = x_46_re_m * ((x_46_im_m * x_46_im_m) * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-2d-316)) then
tmp = x_46re_m * ((x_46im_m * x_46im_m) * (-3.0d0))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) {
tmp = x_46_re_m * ((x_46_im_m * x_46_im_m) * -3.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316: tmp = x_46_re_m * ((x_46_im_m * x_46_im_m) * -3.0) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = Float64(x_46_re_m * Float64(Float64(x_46_im_m * x_46_im_m) * -3.0)); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -2e-316) tmp = x_46_re_m * ((x_46_im_m * x_46_im_m) * -3.0); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -2e-316], N[(x$46$re$95$m * N[(N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -2 \cdot 10^{-316}:\\
\;\;\;\;x.re\_m \cdot \left(\left(x.im\_m \cdot x.im\_m\right) \cdot -3\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -2.000000017e-316Initial program 87.4%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval44.2
Simplified44.2%
if -2.000000017e-316 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 77.6%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.1
Simplified65.1%
Final simplification56.4%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-1e-71)
(* x.im_m (* x.im_m (- -2.0 x.re_m)))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) {
tmp = x_46_im_m * (x_46_im_m * (-2.0 - x_46_re_m));
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-1d-71)) then
tmp = x_46im_m * (x_46im_m * ((-2.0d0) - x_46re_m))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) {
tmp = x_46_im_m * (x_46_im_m * (-2.0 - x_46_re_m));
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71: tmp = x_46_im_m * (x_46_im_m * (-2.0 - x_46_re_m)) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -1e-71) tmp = Float64(x_46_im_m * Float64(x_46_im_m * Float64(-2.0 - x_46_re_m))); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) tmp = x_46_im_m * (x_46_im_m * (-2.0 - x_46_re_m)); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -1e-71], N[(x$46$im$95$m * N[(x$46$im$95$m * N[(-2.0 - x$46$re$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -1 \cdot 10^{-71}:\\
\;\;\;\;x.im\_m \cdot \left(x.im\_m \cdot \left(-2 - x.re\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -9.9999999999999992e-72Initial program 85.5%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8
Applied egg-rr99.8%
Taylor expanded in x.im around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6496.7
Simplified96.7%
distribute-lft-inN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
+-lowering-+.f6460.9
Applied egg-rr60.9%
Taylor expanded in x.im around inf
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6420.1
Simplified20.1%
if -9.9999999999999992e-72 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 79.5%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.0
Simplified63.0%
Final simplification47.6%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<=
(-
(* x.re_m (- (* x.re_m x.re_m) (* x.im_m x.im_m)))
(* x.im_m (+ (* x.re_m x.im_m) (* x.re_m x.im_m))))
-1e-71)
(* x.im_m (* x.im_m -2.0))
(* x.re_m (* x.re_m x.re_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) {
tmp = x_46_im_m * (x_46_im_m * -2.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
real(8) :: tmp
if (((x_46re_m * ((x_46re_m * x_46re_m) - (x_46im_m * x_46im_m))) - (x_46im_m * ((x_46re_m * x_46im_m) + (x_46re_m * x_46im_m)))) <= (-1d-71)) then
tmp = x_46im_m * (x_46im_m * (-2.0d0))
else
tmp = x_46re_m * (x_46re_m * x_46re_m)
end if
code = x_46re_s * tmp
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) {
tmp = x_46_im_m * (x_46_im_m * -2.0);
} else {
tmp = x_46_re_m * (x_46_re_m * x_46_re_m);
}
return x_46_re_s * tmp;
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): tmp = 0 if ((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71: tmp = x_46_im_m * (x_46_im_m * -2.0) else: tmp = x_46_re_m * (x_46_re_m * x_46_re_m) return x_46_re_s * tmp
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (Float64(Float64(x_46_re_m * Float64(Float64(x_46_re_m * x_46_re_m) - Float64(x_46_im_m * x_46_im_m))) - Float64(x_46_im_m * Float64(Float64(x_46_re_m * x_46_im_m) + Float64(x_46_re_m * x_46_im_m)))) <= -1e-71) tmp = Float64(x_46_im_m * Float64(x_46_im_m * -2.0)); else tmp = Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp_2 = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0; if (((x_46_re_m * ((x_46_re_m * x_46_re_m) - (x_46_im_m * x_46_im_m))) - (x_46_im_m * ((x_46_re_m * x_46_im_m) + (x_46_re_m * x_46_im_m)))) <= -1e-71) tmp = x_46_im_m * (x_46_im_m * -2.0); else tmp = x_46_re_m * (x_46_re_m * x_46_re_m); end tmp_2 = x_46_re_s * tmp; end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[N[(N[(x$46$re$95$m * N[(N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision] - N[(x$46$im$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x$46$im$95$m * N[(N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision] + N[(x$46$re$95$m * x$46$im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -1e-71], N[(x$46$im$95$m * N[(x$46$im$95$m * -2.0), $MachinePrecision]), $MachinePrecision], N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m - x.im\_m \cdot x.im\_m\right) - x.im\_m \cdot \left(x.re\_m \cdot x.im\_m + x.re\_m \cdot x.im\_m\right) \leq -1 \cdot 10^{-71}:\\
\;\;\;\;x.im\_m \cdot \left(x.im\_m \cdot -2\right)\\
\mathbf{else}:\\
\;\;\;\;x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\\
\end{array}
\end{array}
if (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) < -9.9999999999999992e-72Initial program 85.5%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8
Applied egg-rr99.8%
Taylor expanded in x.im around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6496.7
Simplified96.7%
distribute-lft-inN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
+-lowering-+.f6460.9
Applied egg-rr60.9%
Taylor expanded in x.re around 0
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6418.6
Simplified18.6%
if -9.9999999999999992e-72 < (-.f64 (*.f64 (-.f64 (*.f64 x.re x.re) (*.f64 x.im x.im)) x.re) (*.f64 (+.f64 (*.f64 x.re x.im) (*.f64 x.im x.re)) x.im)) Initial program 79.5%
Taylor expanded in x.re around inf
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.0
Simplified63.0%
Final simplification47.1%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<= x.re_m 6.5e+20)
(fma x.im_m (* x.im_m (* x.re_m -3.0)) (* x.re_m (* x.re_m x.re_m)))
(fma (- x.re_m x.im_m) (* x.re_m (+ x.re_m x.im_m)) (+ x.im_m x.im_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (x_46_re_m <= 6.5e+20) {
tmp = fma(x_46_im_m, (x_46_im_m * (x_46_re_m * -3.0)), (x_46_re_m * (x_46_re_m * x_46_re_m)));
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (x_46_re_m <= 6.5e+20) tmp = fma(x_46_im_m, Float64(x_46_im_m * Float64(x_46_re_m * -3.0)), Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m))); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[x$46$re$95$m, 6.5e+20], N[(x$46$im$95$m * N[(x$46$im$95$m * N[(x$46$re$95$m * -3.0), $MachinePrecision]), $MachinePrecision] + N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \leq 6.5 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(x.im\_m, x.im\_m \cdot \left(x.re\_m \cdot -3\right), x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
if x.re < 6.5e20Initial program 82.3%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
unpow2N/A
accelerator-lowering-fma.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval85.3
Simplified85.3%
+-commutativeN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6491.7
Applied egg-rr91.7%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6491.7
Applied egg-rr91.7%
if 6.5e20 < x.re Initial program 79.5%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6488.0
Applied egg-rr88.0%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr99.8%
Final simplification93.6%
x.im_m = (fabs.f64 x.im)
x.re\_m = (fabs.f64 x.re)
x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re)
(FPCore (x.re_s x.re_m x.im_m)
:precision binary64
(*
x.re_s
(if (<= x.re_m 6.5e+20)
(fma x.im_m (* x.re_m (* x.im_m -3.0)) (* x.re_m (* x.re_m x.re_m)))
(fma (- x.re_m x.im_m) (* x.re_m (+ x.re_m x.im_m)) (+ x.im_m x.im_m)))))x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
double tmp;
if (x_46_re_m <= 6.5e+20) {
tmp = fma(x_46_im_m, (x_46_re_m * (x_46_im_m * -3.0)), (x_46_re_m * (x_46_re_m * x_46_re_m)));
} else {
tmp = fma((x_46_re_m - x_46_im_m), (x_46_re_m * (x_46_re_m + x_46_im_m)), (x_46_im_m + x_46_im_m));
}
return x_46_re_s * tmp;
}
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = 0.0 if (x_46_re_m <= 6.5e+20) tmp = fma(x_46_im_m, Float64(x_46_re_m * Float64(x_46_im_m * -3.0)), Float64(x_46_re_m * Float64(x_46_re_m * x_46_re_m))); else tmp = fma(Float64(x_46_re_m - x_46_im_m), Float64(x_46_re_m * Float64(x_46_re_m + x_46_im_m)), Float64(x_46_im_m + x_46_im_m)); end return Float64(x_46_re_s * tmp) end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * If[LessEqual[x$46$re$95$m, 6.5e+20], N[(x$46$im$95$m * N[(x$46$re$95$m * N[(x$46$im$95$m * -3.0), $MachinePrecision]), $MachinePrecision] + N[(x$46$re$95$m * N[(x$46$re$95$m * x$46$re$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re$95$m - x$46$im$95$m), $MachinePrecision] * N[(x$46$re$95$m * N[(x$46$re$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision] + N[(x$46$im$95$m + x$46$im$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \begin{array}{l}
\mathbf{if}\;x.re\_m \leq 6.5 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(x.im\_m, x.re\_m \cdot \left(x.im\_m \cdot -3\right), x.re\_m \cdot \left(x.re\_m \cdot x.re\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x.re\_m - x.im\_m, x.re\_m \cdot \left(x.re\_m + x.im\_m\right), x.im\_m + x.im\_m\right)\\
\end{array}
\end{array}
if x.re < 6.5e20Initial program 82.3%
Taylor expanded in x.re around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
unpow2N/A
accelerator-lowering-fma.f64N/A
distribute-rgt-out--N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval85.3
Simplified85.3%
+-commutativeN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6491.7
Applied egg-rr91.7%
if 6.5e20 < x.re Initial program 79.5%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6488.0
Applied egg-rr88.0%
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-outN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
associate-*r/N/A
+-inversesN/A
distribute-lft-out--N/A
+-inversesN/A
distribute-neg-frac2N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
flip-+N/A
distribute-lft-inN/A
Applied egg-rr99.8%
Final simplification93.5%
x.im_m = (fabs.f64 x.im) x.re\_m = (fabs.f64 x.re) x.re\_s = (copysign.f64 #s(literal 1 binary64) x.re) (FPCore (x.re_s x.re_m x.im_m) :precision binary64 (* x.re_s (* x.im_m (* x.im_m -2.0))))
x.im_m = fabs(x_46_im);
x.re\_m = fabs(x_46_re);
x.re\_s = copysign(1.0, x_46_re);
double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
return x_46_re_s * (x_46_im_m * (x_46_im_m * -2.0));
}
x.im_m = abs(x_46im)
x.re\_m = abs(x_46re)
x.re\_s = copysign(1.0d0, x_46re)
real(8) function code(x_46re_s, x_46re_m, x_46im_m)
real(8), intent (in) :: x_46re_s
real(8), intent (in) :: x_46re_m
real(8), intent (in) :: x_46im_m
code = x_46re_s * (x_46im_m * (x_46im_m * (-2.0d0)))
end function
x.im_m = Math.abs(x_46_im);
x.re\_m = Math.abs(x_46_re);
x.re\_s = Math.copySign(1.0, x_46_re);
public static double code(double x_46_re_s, double x_46_re_m, double x_46_im_m) {
return x_46_re_s * (x_46_im_m * (x_46_im_m * -2.0));
}
x.im_m = math.fabs(x_46_im) x.re\_m = math.fabs(x_46_re) x.re\_s = math.copysign(1.0, x_46_re) def code(x_46_re_s, x_46_re_m, x_46_im_m): return x_46_re_s * (x_46_im_m * (x_46_im_m * -2.0))
x.im_m = abs(x_46_im) x.re\_m = abs(x_46_re) x.re\_s = copysign(1.0, x_46_re) function code(x_46_re_s, x_46_re_m, x_46_im_m) return Float64(x_46_re_s * Float64(x_46_im_m * Float64(x_46_im_m * -2.0))) end
x.im_m = abs(x_46_im); x.re\_m = abs(x_46_re); x.re\_s = sign(x_46_re) * abs(1.0); function tmp = code(x_46_re_s, x_46_re_m, x_46_im_m) tmp = x_46_re_s * (x_46_im_m * (x_46_im_m * -2.0)); end
x.im_m = N[Abs[x$46$im], $MachinePrecision]
x.re\_m = N[Abs[x$46$re], $MachinePrecision]
x.re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[x$46$re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[x$46$re$95$s_, x$46$re$95$m_, x$46$im$95$m_] := N[(x$46$re$95$s * N[(x$46$im$95$m * N[(x$46$im$95$m * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x.im_m = \left|x.im\right|
\\
x.re\_m = \left|x.re\right|
\\
x.re\_s = \mathsf{copysign}\left(1, x.re\right)
\\
x.re\_s \cdot \left(x.im\_m \cdot \left(x.im\_m \cdot -2\right)\right)
\end{array}
Initial program 81.7%
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6493.2
Applied egg-rr93.2%
Taylor expanded in x.im around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6491.0
Simplified91.0%
distribute-lft-inN/A
flip-+N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
+-lowering-+.f6464.1
Applied egg-rr64.1%
Taylor expanded in x.re around 0
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6417.3
Simplified17.3%
(FPCore (x.re x.im) :precision binary64 (+ (* (* x.re x.re) (- x.re x.im)) (* (* x.re x.im) (- x.re (* 3.0 x.im)))))
double code(double x_46_re, double x_46_im) {
return ((x_46_re * x_46_re) * (x_46_re - x_46_im)) + ((x_46_re * x_46_im) * (x_46_re - (3.0 * x_46_im)));
}
real(8) function code(x_46re, x_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
code = ((x_46re * x_46re) * (x_46re - x_46im)) + ((x_46re * x_46im) * (x_46re - (3.0d0 * x_46im)))
end function
public static double code(double x_46_re, double x_46_im) {
return ((x_46_re * x_46_re) * (x_46_re - x_46_im)) + ((x_46_re * x_46_im) * (x_46_re - (3.0 * x_46_im)));
}
def code(x_46_re, x_46_im): return ((x_46_re * x_46_re) * (x_46_re - x_46_im)) + ((x_46_re * x_46_im) * (x_46_re - (3.0 * x_46_im)))
function code(x_46_re, x_46_im) return Float64(Float64(Float64(x_46_re * x_46_re) * Float64(x_46_re - x_46_im)) + Float64(Float64(x_46_re * x_46_im) * Float64(x_46_re - Float64(3.0 * x_46_im)))) end
function tmp = code(x_46_re, x_46_im) tmp = ((x_46_re * x_46_re) * (x_46_re - x_46_im)) + ((x_46_re * x_46_im) * (x_46_re - (3.0 * x_46_im))); end
code[x$46$re_, x$46$im_] := N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] * N[(x$46$re - x$46$im), $MachinePrecision]), $MachinePrecision] + N[(N[(x$46$re * x$46$im), $MachinePrecision] * N[(x$46$re - N[(3.0 * x$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x.re \cdot x.re\right) \cdot \left(x.re - x.im\right) + \left(x.re \cdot x.im\right) \cdot \left(x.re - 3 \cdot x.im\right)
\end{array}
herbie shell --seed 2024205
(FPCore (x.re x.im)
:name "math.cube on complex, real part"
:precision binary64
:alt
(! :herbie-platform default (+ (* (* x.re x.re) (- x.re x.im)) (* (* x.re x.im) (- x.re (* 3 x.im)))))
(- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))