
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 28 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Initial program 100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (fma im (* im -0.5) 1.0))
(if (<= t_0 -0.01)
(* (cos im) (fma re (fma re 0.5 1.0) 1.0))
(if (<= t_0 2e-26)
(exp re)
(if (<= t_0 0.9999999)
(* (cos im) (fma (* re re) 0.5 (+ re 1.0)))
(exp re)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * fma(im, (im * -0.5), 1.0);
} else if (t_0 <= -0.01) {
tmp = cos(im) * fma(re, fma(re, 0.5, 1.0), 1.0);
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.9999999) {
tmp = cos(im) * fma((re * re), 0.5, (re + 1.0));
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * fma(im, Float64(im * -0.5), 1.0)); elseif (t_0 <= -0.01) tmp = Float64(cos(im) * fma(re, fma(re, 0.5, 1.0), 1.0)); elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.9999999) tmp = Float64(cos(im) * fma(Float64(re * re), 0.5, Float64(re + 1.0))); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(im * N[(im * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.01], N[(N[Cos[im], $MachinePrecision] * N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.9999999], N[(N[Cos[im], $MachinePrecision] * N[(N[(re * re), $MachinePrecision] * 0.5 + N[(re + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Exp[re], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im, im \cdot -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq -0.01:\\
\;\;\;\;\cos im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.9999999:\\
\;\;\;\;\cos im \cdot \mathsf{fma}\left(re \cdot re, 0.5, re + 1\right)\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0
Simplified100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6498.0
Simplified98.0%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.999999900000000053 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.7
Simplified99.7%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.999999900000000053Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6496.3
Simplified96.3%
distribute-lft-inN/A
*-rgt-identityN/A
associate-+l+N/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6496.3
Applied egg-rr96.3%
Final simplification99.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (cos im) (fma re (fma re 0.5 1.0) 1.0)))
(t_1 (* (exp re) (cos im))))
(if (<= t_1 (- INFINITY))
(* (exp re) (fma im (* im -0.5) 1.0))
(if (<= t_1 -0.01)
t_0
(if (<= t_1 2e-26) (exp re) (if (<= t_1 0.9999999) t_0 (exp re)))))))
double code(double re, double im) {
double t_0 = cos(im) * fma(re, fma(re, 0.5, 1.0), 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = exp(re) * fma(im, (im * -0.5), 1.0);
} else if (t_1 <= -0.01) {
tmp = t_0;
} else if (t_1 <= 2e-26) {
tmp = exp(re);
} else if (t_1 <= 0.9999999) {
tmp = t_0;
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(cos(im) * fma(re, fma(re, 0.5, 1.0), 1.0)) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(exp(re) * fma(im, Float64(im * -0.5), 1.0)); elseif (t_1 <= -0.01) tmp = t_0; elseif (t_1 <= 2e-26) tmp = exp(re); elseif (t_1 <= 0.9999999) tmp = t_0; else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Cos[im], $MachinePrecision] * N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(im * N[(im * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -0.01], t$95$0, If[LessEqual[t$95$1, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$1, 0.9999999], t$95$0, N[Exp[re], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im, im \cdot -0.5, 1\right)\\
\mathbf{elif}\;t\_1 \leq -0.01:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_1 \leq 0.9999999:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0
Simplified100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.999999900000000053Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6497.3
Simplified97.3%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.999999900000000053 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.7
Simplified99.7%
Final simplification99.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (cos im) (+ re 1.0))) (t_1 (* (exp re) (cos im))))
(if (<= t_1 (- INFINITY))
(* (exp re) (fma im (* im -0.5) 1.0))
(if (<= t_1 -0.01)
t_0
(if (<= t_1 2e-26) (exp re) (if (<= t_1 0.999) t_0 (exp re)))))))
double code(double re, double im) {
double t_0 = cos(im) * (re + 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = exp(re) * fma(im, (im * -0.5), 1.0);
} else if (t_1 <= -0.01) {
tmp = t_0;
} else if (t_1 <= 2e-26) {
tmp = exp(re);
} else if (t_1 <= 0.999) {
tmp = t_0;
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(cos(im) * Float64(re + 1.0)) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(exp(re) * fma(im, Float64(im * -0.5), 1.0)); elseif (t_1 <= -0.01) tmp = t_0; elseif (t_1 <= 2e-26) tmp = exp(re); elseif (t_1 <= 0.999) tmp = t_0; else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Cos[im], $MachinePrecision] * N[(re + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(im * N[(im * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -0.01], t$95$0, If[LessEqual[t$95$1, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$1, 0.999], t$95$0, N[Exp[re], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos im \cdot \left(re + 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im, im \cdot -0.5, 1\right)\\
\mathbf{elif}\;t\_1 \leq -0.01:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_1 \leq 0.999:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0
Simplified100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998999999999999999Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6497.0
Simplified97.0%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.998999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.4
Simplified99.4%
Final simplification98.8%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))) (t_1 (* (cos im) (+ re 1.0))))
(if (<= t_0 (- INFINITY))
(fma
(* im im)
(fma
re
(* im (* im (fma (* im im) -0.001388888888888889 0.041666666666666664)))
(* re -0.5))
re)
(if (<= t_0 -0.01)
t_1
(if (<= t_0 2e-26) (exp re) (if (<= t_0 0.999) t_1 (exp re)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = cos(im) * (re + 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma((im * im), fma(re, (im * (im * fma((im * im), -0.001388888888888889, 0.041666666666666664))), (re * -0.5)), re);
} else if (t_0 <= -0.01) {
tmp = t_1;
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.999) {
tmp = t_1;
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = Float64(cos(im) * Float64(re + 1.0)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = fma(Float64(im * im), fma(re, Float64(im * Float64(im * fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664))), Float64(re * -0.5)), re); elseif (t_0 <= -0.01) tmp = t_1; elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.999) tmp = t_1; else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[im], $MachinePrecision] * N[(re + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(im * im), $MachinePrecision] * N[(re * N[(im * N[(im * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(re * -0.5), $MachinePrecision]), $MachinePrecision] + re), $MachinePrecision], If[LessEqual[t$95$0, -0.01], t$95$1, If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.999], t$95$1, N[Exp[re], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \cos im \cdot \left(re + 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(re, im \cdot \left(im \cdot \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right)\right), re \cdot -0.5\right), re\right)\\
\mathbf{elif}\;t\_0 \leq -0.01:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.999:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f645.3
Simplified5.3%
Taylor expanded in re around inf
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f645.3
Simplified5.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified94.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998999999999999999Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6497.0
Simplified97.0%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.998999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.4
Simplified99.4%
Final simplification98.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(fma
(* im im)
(fma
re
(* im (* im (fma (* im im) -0.001388888888888889 0.041666666666666664)))
(* re -0.5))
re)
(if (<= t_0 -0.01)
(cos im)
(if (<= t_0 2e-26) (exp re) (if (<= t_0 0.999) (cos im) (exp re)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma((im * im), fma(re, (im * (im * fma((im * im), -0.001388888888888889, 0.041666666666666664))), (re * -0.5)), re);
} else if (t_0 <= -0.01) {
tmp = cos(im);
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.999) {
tmp = cos(im);
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = fma(Float64(im * im), fma(re, Float64(im * Float64(im * fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664))), Float64(re * -0.5)), re); elseif (t_0 <= -0.01) tmp = cos(im); elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.999) tmp = cos(im); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(im * im), $MachinePrecision] * N[(re * N[(im * N[(im * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(re * -0.5), $MachinePrecision]), $MachinePrecision] + re), $MachinePrecision], If[LessEqual[t$95$0, -0.01], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.999], N[Cos[im], $MachinePrecision], N[Exp[re], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(re, im \cdot \left(im \cdot \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right)\right), re \cdot -0.5\right), re\right)\\
\mathbf{elif}\;t\_0 \leq -0.01:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.999:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f645.3
Simplified5.3%
Taylor expanded in re around inf
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f645.3
Simplified5.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified94.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998999999999999999Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6495.8
Simplified95.8%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.998999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.4
Simplified99.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(fma
(* im im)
(fma
re
(* im (* im (fma (* im im) -0.001388888888888889 0.041666666666666664)))
(* re -0.5))
re)
(if (<= t_0 -0.01)
(cos im)
(if (<= t_0 0.0)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(if (<= t_0 0.999)
(cos im)
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma
im
(* im (fma im (* im 0.041666666666666664) -0.5))
1.0))))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma((im * im), fma(re, (im * (im * fma((im * im), -0.001388888888888889, 0.041666666666666664))), (re * -0.5)), re);
} else if (t_0 <= -0.01) {
tmp = cos(im);
} else if (t_0 <= 0.0) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else if (t_0 <= 0.999) {
tmp = cos(im);
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, (im * fma(im, (im * 0.041666666666666664), -0.5)), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = fma(Float64(im * im), fma(re, Float64(im * Float64(im * fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664))), Float64(re * -0.5)), re); elseif (t_0 <= -0.01) tmp = cos(im); elseif (t_0 <= 0.0) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); elseif (t_0 <= 0.999) tmp = cos(im); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, Float64(im * fma(im, Float64(im * 0.041666666666666664), -0.5)), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(im * im), $MachinePrecision] * N[(re * N[(im * N[(im * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(re * -0.5), $MachinePrecision]), $MachinePrecision] + re), $MachinePrecision], If[LessEqual[t$95$0, -0.01], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.999], N[Cos[im], $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(im * N[(im * N[(im * 0.041666666666666664), $MachinePrecision] + -0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(re, im \cdot \left(im \cdot \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right)\right), re \cdot -0.5\right), re\right)\\
\mathbf{elif}\;t\_0 \leq -0.01:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0.999:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im, im \cdot \mathsf{fma}\left(im, im \cdot 0.041666666666666664, -0.5\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f645.3
Simplified5.3%
Taylor expanded in re around inf
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f645.3
Simplified5.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified94.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998999999999999999Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6493.6
Simplified93.6%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.7
Simplified35.7%
if 0.998999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.8
Simplified87.8%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6491.5
Simplified91.5%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(* (cos im) (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))
(if (<= t_0 2e-26)
(exp re)
(if (<= t_0 0.999)
(*
(cos im)
(+ re (fma re (* re (fma re 0.16666666666666666 0.5)) 1.0)))
(*
(exp re)
(fma im (* im (fma (* im im) 0.041666666666666664 -0.5)) 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.999) {
tmp = cos(im) * (re + fma(re, (re * fma(re, 0.16666666666666666, 0.5)), 1.0));
} else {
tmp = exp(re) * fma(im, (im * fma((im * im), 0.041666666666666664, -0.5)), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = Float64(cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.999) tmp = Float64(cos(im) * Float64(re + fma(re, Float64(re * fma(re, 0.16666666666666666, 0.5)), 1.0))); else tmp = Float64(exp(re) * fma(im, Float64(im * fma(Float64(im * im), 0.041666666666666664, -0.5)), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[Cos[im], $MachinePrecision] * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.999], N[(N[Cos[im], $MachinePrecision] * N[(re + N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Exp[re], $MachinePrecision] * N[(im * N[(im * N[(N[(im * im), $MachinePrecision] * 0.041666666666666664 + -0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\cos im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.999:\\
\;\;\;\;\cos im \cdot \left(re + \mathsf{fma}\left(re, re \cdot \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im, im \cdot \mathsf{fma}\left(im \cdot im, 0.041666666666666664, -0.5\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6492.4
Simplified92.4%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.3
Simplified99.3%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998999999999999999Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6497.8
Simplified97.8%
distribute-lft-inN/A
*-rgt-identityN/A
associate-+l+N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f6497.8
Applied egg-rr97.8%
if 0.998999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
associate-+l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-lft1-inN/A
distribute-rgt-outN/A
associate-+l+N/A
Simplified100.0%
Final simplification97.7%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(* (cos im) (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))
(if (<= t_0 2e-26)
(exp re)
(if (<= t_0 0.9999999)
(*
(cos im)
(+ re (fma re (* re (fma re 0.16666666666666666 0.5)) 1.0)))
(exp re))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.9999999) {
tmp = cos(im) * (re + fma(re, (re * fma(re, 0.16666666666666666, 0.5)), 1.0));
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = Float64(cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.9999999) tmp = Float64(cos(im) * Float64(re + fma(re, Float64(re * fma(re, 0.16666666666666666, 0.5)), 1.0))); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[Cos[im], $MachinePrecision] * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.9999999], N[(N[Cos[im], $MachinePrecision] * N[(re + N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Exp[re], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\cos im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.9999999:\\
\;\;\;\;\cos im \cdot \left(re + \mathsf{fma}\left(re, re \cdot \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6492.4
Simplified92.4%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.999999900000000053 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.7
Simplified99.7%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.999999900000000053Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6497.6
Simplified97.6%
distribute-lft-inN/A
*-rgt-identityN/A
associate-+l+N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f6497.6
Applied egg-rr97.6%
Final simplification97.7%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im)))
(t_1
(*
(cos im)
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
(if (<= t_0 -0.01)
t_1
(if (<= t_0 2e-26) (exp re) (if (<= t_0 0.9999999) t_1 (exp re))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
double tmp;
if (t_0 <= -0.01) {
tmp = t_1;
} else if (t_0 <= 2e-26) {
tmp = exp(re);
} else if (t_0 <= 0.9999999) {
tmp = t_1;
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = Float64(cos(im) * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)) tmp = 0.0 if (t_0 <= -0.01) tmp = t_1; elseif (t_0 <= 2e-26) tmp = exp(re); elseif (t_0 <= 0.9999999) tmp = t_1; else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[im], $MachinePrecision] * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], t$95$1, If[LessEqual[t$95$0, 2e-26], N[Exp[re], $MachinePrecision], If[LessEqual[t$95$0, 0.9999999], t$95$1, N[Exp[re], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \cos im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;t\_0 \leq 0.9999999:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002 or 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.999999900000000053Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6494.3
Simplified94.3%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26 or 0.999999900000000053 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6499.7
Simplified99.7%
Final simplification97.7%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(fma
(* im im)
(fma
(* im im)
(fma (* im im) -0.001388888888888889 0.041666666666666664)
-0.5)
1.0)
(if (<= t_0 0.0)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(if (<= t_0 0.9995)
(fma re (/ 1.0 (fma re (fma re 0.08333333333333333 -0.5) 1.0)) 1.0)
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma im (* im (fma im (* im 0.041666666666666664) -0.5)) 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = fma((im * im), fma((im * im), fma((im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0);
} else if (t_0 <= 0.0) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else if (t_0 <= 0.9995) {
tmp = fma(re, (1.0 / fma(re, fma(re, 0.08333333333333333, -0.5), 1.0)), 1.0);
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, (im * fma(im, (im * 0.041666666666666664), -0.5)), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = fma(Float64(im * im), fma(Float64(im * im), fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0); elseif (t_0 <= 0.0) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); elseif (t_0 <= 0.9995) tmp = fma(re, Float64(1.0 / fma(re, fma(re, 0.08333333333333333, -0.5), 1.0)), 1.0); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, Float64(im * fma(im, Float64(im * 0.041666666666666664), -0.5)), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision] + -0.5), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.9995], N[(re * N[(1.0 / N[(re * N[(re * 0.08333333333333333 + -0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(im * N[(im * N[(im * 0.041666666666666664), $MachinePrecision] + -0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right), -0.5\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0.9995:\\
\;\;\;\;\mathsf{fma}\left(re, \frac{1}{\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.08333333333333333, -0.5\right), 1\right)}, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im, im \cdot \mathsf{fma}\left(im, im \cdot 0.041666666666666664, -0.5\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6470.1
Simplified70.1%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.7
Simplified35.7%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99950000000000006Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6492.2
Simplified92.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6421.7
Simplified21.7%
flip3-+N/A
clear-numN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
accelerator-lowering-fma.f6421.7
Applied egg-rr21.7%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6422.4
Simplified22.4%
if 0.99950000000000006 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.8
Simplified87.8%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6491.6
Simplified91.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(fma
(* im im)
(fma
(* im im)
(fma (* im im) -0.001388888888888889 0.041666666666666664)
-0.5)
1.0)
(if (<= t_0 2e-26)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(fma
(* (* re re) (fma (* re re) 0.027777777777777776 -0.25))
(/ 1.0 (fma re 0.16666666666666666 -0.5))
(+ re 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = fma((im * im), fma((im * im), fma((im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0);
} else if (t_0 <= 2e-26) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else {
tmp = fma(((re * re) * fma((re * re), 0.027777777777777776, -0.25)), (1.0 / fma(re, 0.16666666666666666, -0.5)), (re + 1.0));
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = fma(Float64(im * im), fma(Float64(im * im), fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0); elseif (t_0 <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); else tmp = fma(Float64(Float64(re * re) * fma(Float64(re * re), 0.027777777777777776, -0.25)), Float64(1.0 / fma(re, 0.16666666666666666, -0.5)), Float64(re + 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision] + -0.5), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * N[(N[(re * re), $MachinePrecision] * 0.027777777777777776 + -0.25), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(re * 0.16666666666666666 + -0.5), $MachinePrecision]), $MachinePrecision] + N[(re + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right), -0.5\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(re \cdot re\right) \cdot \mathsf{fma}\left(re \cdot re, 0.027777777777777776, -0.25\right), \frac{1}{\mathsf{fma}\left(re, 0.16666666666666666, -0.5\right)}, re + 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6470.1
Simplified70.1%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6434.6
Simplified34.6%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
distribute-lft-inN/A
*-rgt-identityN/A
associate-+r+N/A
associate-*r*N/A
flip-+N/A
associate-*r/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
Applied egg-rr73.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(fma
(* im im)
(fma
(* im im)
(fma (* im im) -0.001388888888888889 0.041666666666666664)
-0.5)
1.0)
(if (<= t_0 2e-26)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = fma((im * im), fma((im * im), fma((im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0);
} else if (t_0 <= 2e-26) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = fma(Float64(im * im), fma(Float64(im * im), fma(Float64(im * im), -0.001388888888888889, 0.041666666666666664), -0.5), 1.0); elseif (t_0 <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); else tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision] + -0.5), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.001388888888888889, 0.041666666666666664\right), -0.5\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6470.1
Simplified70.1%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6434.6
Simplified34.6%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(* (fma im (* im -0.5) 1.0) (fma re (fma re 0.5 1.0) 1.0))
(if (<= t_0 2e-26)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = fma(im, (im * -0.5), 1.0) * fma(re, fma(re, 0.5, 1.0), 1.0);
} else if (t_0 <= 2e-26) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = Float64(fma(im, Float64(im * -0.5), 1.0) * fma(re, fma(re, 0.5, 1.0), 1.0)); elseif (t_0 <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); else tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(N[(im * N[(im * -0.5), $MachinePrecision] + 1.0), $MachinePrecision] * N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot -0.5, 1\right) \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6430.7
Simplified30.7%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6427.7
Simplified27.7%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6434.6
Simplified34.6%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
Final simplification52.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(* im (* im (fma re -0.5 -0.5)))
(if (<= t_0 2e-26)
(* (+ re 1.0) (* 0.041666666666666664 (* (* im im) (* im im))))
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else if (t_0 <= 2e-26) {
tmp = (re + 1.0) * (0.041666666666666664 * ((im * im) * (im * im)));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); elseif (t_0 <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(0.041666666666666664 * Float64(Float64(im * im) * Float64(im * im)))); else tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(0.041666666666666664 * N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(0.041666666666666664 \cdot \left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6471.9
Simplified71.9%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f640.5
Simplified0.5%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6423.3
Simplified23.3%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6423.3
Simplified23.3%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6434.6
Simplified34.6%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.01)
(* im (* im (fma re -0.5 -0.5)))
(if (<= t_0 2e-26)
(*
(* (* im im) (* im im))
(fma re 0.041666666666666664 0.041666666666666664))
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.01) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else if (t_0 <= 2e-26) {
tmp = ((im * im) * (im * im)) * fma(re, 0.041666666666666664, 0.041666666666666664);
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.01) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); elseif (t_0 <= 2e-26) tmp = Float64(Float64(Float64(im * im) * Float64(im * im)) * fma(re, 0.041666666666666664, 0.041666666666666664)); else tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.01], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-26], N[(N[(N[(im * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] * N[(re * 0.041666666666666664 + 0.041666666666666664), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.01:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(\left(im \cdot im\right) \cdot \left(im \cdot im\right)\right) \cdot \mathsf{fma}\left(re, 0.041666666666666664, 0.041666666666666664\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0100000000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6471.9
Simplified71.9%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f640.5
Simplified0.5%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6423.3
Simplified23.3%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6423.3
Simplified23.3%
if -0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f642.2
Simplified2.2%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.8
Simplified1.8%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6433.0
Simplified33.0%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 5e-127)
(* im (* im (fma re -0.5 -0.5)))
(if (<= t_0 2.0)
(fma re (fma re 0.5 1.0) 1.0)
(* (fma re 0.16666666666666666 0.5) (* re re))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= 5e-127) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else if (t_0 <= 2.0) {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0);
} else {
tmp = fma(re, 0.16666666666666666, 0.5) * (re * re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= 5e-127) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); elseif (t_0 <= 2.0) tmp = fma(re, fma(re, 0.5, 1.0), 1.0); else tmp = Float64(fma(re, 0.16666666666666666, 0.5) * Float64(re * re)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 5e-127], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2.0], N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] * N[(re * re), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq 5 \cdot 10^{-127}:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, 0.16666666666666666, 0.5\right) \cdot \left(re \cdot re\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 4.9999999999999997e-127Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.6
Simplified37.6%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.7
Simplified12.7%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6417.0
Simplified17.0%
if 4.9999999999999997e-127 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6471.4
Simplified71.4%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6470.0
Simplified70.0%
if 2 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in re around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f6471.3
Simplified71.3%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 5e-127)
(* im (* im (fma re -0.5 -0.5)))
(if (<= t_0 2.0)
(fma re (fma re 0.5 1.0) 1.0)
(* 0.16666666666666666 (* re (* re re)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= 5e-127) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else if (t_0 <= 2.0) {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0);
} else {
tmp = 0.16666666666666666 * (re * (re * re));
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= 5e-127) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); elseif (t_0 <= 2.0) tmp = fma(re, fma(re, 0.5, 1.0), 1.0); else tmp = Float64(0.16666666666666666 * Float64(re * Float64(re * re))); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 5e-127], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2.0], N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(0.16666666666666666 * N[(re * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq 5 \cdot 10^{-127}:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;0.16666666666666666 \cdot \left(re \cdot \left(re \cdot re\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 4.9999999999999997e-127Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.6
Simplified37.6%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.7
Simplified12.7%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6417.0
Simplified17.0%
if 4.9999999999999997e-127 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6471.4
Simplified71.4%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6470.0
Simplified70.0%
if 2 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in re around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6471.3
Simplified71.3%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 2e-26)
(fma -0.5 (* im im) 1.0)
(if (<= t_0 2.0)
(fma re (fma re 0.5 1.0) 1.0)
(* 0.16666666666666666 (* re (* re re)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= 2e-26) {
tmp = fma(-0.5, (im * im), 1.0);
} else if (t_0 <= 2.0) {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0);
} else {
tmp = 0.16666666666666666 * (re * (re * re));
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= 2e-26) tmp = fma(-0.5, Float64(im * im), 1.0); elseif (t_0 <= 2.0) tmp = fma(re, fma(re, 0.5, 1.0), 1.0); else tmp = Float64(0.16666666666666666 * Float64(re * Float64(re * re))); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 2e-26], N[(-0.5 * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[t$95$0, 2.0], N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(0.16666666666666666 * N[(re * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\mathsf{fma}\left(-0.5, im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;0.16666666666666666 \cdot \left(re \cdot \left(re \cdot re\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6436.9
Simplified36.9%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f649.7
Simplified9.7%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6471.1
Simplified71.1%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6470.7
Simplified70.7%
if 2 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.3
Simplified71.3%
Taylor expanded in re around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6471.3
Simplified71.3%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (* (+ re 1.0) (* (* im im) -0.5)) (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = (re + 1.0) * ((im * im) * -0.5);
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(Float64(im * im) * -0.5)); else tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.3
Simplified37.3%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.6
Simplified12.6%
Taylor expanded in im around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6422.5
Simplified22.5%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
Final simplification48.1%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (* (+ re 1.0) (* (* im im) -0.5)) (fma re (fma re (* re 0.16666666666666666) 1.0) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = (re + 1.0) * ((im * im) * -0.5);
} else {
tmp = fma(re, fma(re, (re * 0.16666666666666666), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = Float64(Float64(re + 1.0) * Float64(Float64(im * im) * -0.5)); else tmp = fma(re, fma(re, Float64(re * 0.16666666666666666), 1.0), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(N[(re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\left(re + 1\right) \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, re \cdot 0.16666666666666666, 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.3
Simplified37.3%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.6
Simplified12.6%
Taylor expanded in im around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6422.5
Simplified22.5%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
Taylor expanded in re around inf
*-commutativeN/A
*-lowering-*.f6470.3
Simplified70.3%
Final simplification47.7%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (* im (* im (fma re -0.5 -0.5))) (fma re (fma re (* re 0.16666666666666666) 1.0) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else {
tmp = fma(re, fma(re, (re * 0.16666666666666666), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); else tmp = fma(re, fma(re, Float64(re * 0.16666666666666666), 1.0), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, re \cdot 0.16666666666666666, 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.3
Simplified37.3%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.6
Simplified12.6%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6416.9
Simplified16.9%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
Taylor expanded in re around inf
*-commutativeN/A
*-lowering-*.f6470.3
Simplified70.3%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (* im (* im (fma re -0.5 -0.5))) (fma re (* re (* re 0.16666666666666666)) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = im * (im * fma(re, -0.5, -0.5));
} else {
tmp = fma(re, (re * (re * 0.16666666666666666)), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = Float64(im * Float64(im * fma(re, -0.5, -0.5))); else tmp = fma(re, Float64(re * Float64(re * 0.16666666666666666)), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(im * N[(im * N[(re * -0.5 + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(re * N[(re * N[(re * 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;im \cdot \left(im \cdot \mathsf{fma}\left(re, -0.5, -0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, re \cdot \left(re \cdot 0.16666666666666666\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6437.3
Simplified37.3%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f641.2
Simplified1.2%
Taylor expanded in im around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6412.6
Simplified12.6%
Taylor expanded in im around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
metadata-evalN/A
accelerator-lowering-fma.f6416.9
Simplified16.9%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6490.2
Simplified90.2%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6471.1
Simplified71.1%
Taylor expanded in re around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6469.3
Simplified69.3%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (fma -0.5 (* im im) 1.0) (fma re (fma re 0.5 1.0) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = fma(-0.5, (im * im), 1.0);
} else {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = fma(-0.5, Float64(im * im), 1.0); else tmp = fma(re, fma(re, 0.5, 1.0), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(-0.5 * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\mathsf{fma}\left(-0.5, im \cdot im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6436.9
Simplified36.9%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f649.7
Simplified9.7%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6480.3
Simplified80.3%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6467.4
Simplified67.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 2e-26) (fma -0.5 (* im im) 1.0) (+ re 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 2e-26) {
tmp = fma(-0.5, (im * im), 1.0);
} else {
tmp = re + 1.0;
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 2e-26) tmp = fma(-0.5, Float64(im * im), 1.0); else tmp = Float64(re + 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 2e-26], N[(-0.5 * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], N[(re + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 2 \cdot 10^{-26}:\\
\;\;\;\;\mathsf{fma}\left(-0.5, im \cdot im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;re + 1\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < 2.0000000000000001e-26Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6436.9
Simplified36.9%
Taylor expanded in im around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f649.7
Simplified9.7%
if 2.0000000000000001e-26 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6480.3
Simplified80.3%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6449.3
Simplified49.3%
(FPCore (re im)
:precision binary64
(if (<= re -0.00096)
(exp re)
(if (<= re 11.0)
(* (cos im) (fma (* re re) 0.5 (+ re 1.0)))
(if (<= re 1e+103)
(exp re)
(* (cos im) (* re (* 0.16666666666666666 (* re re))))))))
double code(double re, double im) {
double tmp;
if (re <= -0.00096) {
tmp = exp(re);
} else if (re <= 11.0) {
tmp = cos(im) * fma((re * re), 0.5, (re + 1.0));
} else if (re <= 1e+103) {
tmp = exp(re);
} else {
tmp = cos(im) * (re * (0.16666666666666666 * (re * re)));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -0.00096) tmp = exp(re); elseif (re <= 11.0) tmp = Float64(cos(im) * fma(Float64(re * re), 0.5, Float64(re + 1.0))); elseif (re <= 1e+103) tmp = exp(re); else tmp = Float64(cos(im) * Float64(re * Float64(0.16666666666666666 * Float64(re * re)))); end return tmp end
code[re_, im_] := If[LessEqual[re, -0.00096], N[Exp[re], $MachinePrecision], If[LessEqual[re, 11.0], N[(N[Cos[im], $MachinePrecision] * N[(N[(re * re), $MachinePrecision] * 0.5 + N[(re + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1e+103], N[Exp[re], $MachinePrecision], N[(N[Cos[im], $MachinePrecision] * N[(re * N[(0.16666666666666666 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -0.00096:\\
\;\;\;\;e^{re}\\
\mathbf{elif}\;re \leq 11:\\
\;\;\;\;\cos im \cdot \mathsf{fma}\left(re \cdot re, 0.5, re + 1\right)\\
\mathbf{elif}\;re \leq 10^{+103}:\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\cos im \cdot \left(re \cdot \left(0.16666666666666666 \cdot \left(re \cdot re\right)\right)\right)\\
\end{array}
\end{array}
if re < -9.60000000000000024e-4 or 11 < re < 1e103Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6493.4
Simplified93.4%
if -9.60000000000000024e-4 < re < 11Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6499.0
Simplified99.0%
distribute-lft-inN/A
*-rgt-identityN/A
associate-+l+N/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6499.0
Applied egg-rr99.0%
if 1e103 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64100.0
Simplified100.0%
Taylor expanded in re around inf
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0
Simplified100.0%
Final simplification97.4%
(FPCore (re im) :precision binary64 (+ re 1.0))
double code(double re, double im) {
return re + 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = re + 1.0d0
end function
public static double code(double re, double im) {
return re + 1.0;
}
def code(re, im): return re + 1.0
function code(re, im) return Float64(re + 1.0) end
function tmp = code(re, im) tmp = re + 1.0; end
code[re_, im_] := N[(re + 1.0), $MachinePrecision]
\begin{array}{l}
\\
re + 1
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
exp-lowering-exp.f6465.9
Simplified65.9%
Taylor expanded in re around 0
+-commutativeN/A
+-lowering-+.f6426.8
Simplified26.8%
(FPCore (re im) :precision binary64 1.0)
double code(double re, double im) {
return 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = 1.0d0
end function
public static double code(double re, double im) {
return 1.0;
}
def code(re, im): return 1.0
function code(re, im) return 1.0 end
function tmp = code(re, im) tmp = 1.0; end
code[re_, im_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
Taylor expanded in re around 0
cos-lowering-cos.f6452.4
Simplified52.4%
Taylor expanded in im around 0
Simplified26.2%
herbie shell --seed 2024204
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))