
(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 21 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)))
(t_1 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0))
(t_2 (* (+ 1.0 re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -0.1)
t_2
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.99)
t_2
(*
t_1
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double t_2 = (1.0 + re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -0.1) {
tmp = t_2;
} else if (t_0 <= 0.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= 0.99) {
tmp = t_2;
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) t_2 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= -0.1) tmp = t_2; elseif (t_0 <= 0.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.99) tmp = t_2; else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); 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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.1], t$95$2, If[LessEqual[t$95$0, 0.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.99], t$95$2, N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
t_2 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq -0.1:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.99:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites83.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.10000000000000001 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.98999999999999999Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.1
Applied rewrites98.1%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in im around inf
Applied rewrites79.7%
if 0.98999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.6
Applied rewrites84.6%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites79.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6494.4
Applied rewrites94.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im)))
(t_1 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0))
(t_2 (* (+ 1.0 re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -0.1)
t_2
(if (<= t_0 0.0)
(* 1.0 (* (* im im) -0.5))
(if (<= t_0 0.99)
t_2
(*
t_1
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double t_2 = (1.0 + re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -0.1) {
tmp = t_2;
} else if (t_0 <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else if (t_0 <= 0.99) {
tmp = t_2;
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) t_2 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= -0.1) tmp = t_2; elseif (t_0 <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.99) tmp = t_2; else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); 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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.1], t$95$2, If[LessEqual[t$95$0, 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.99], t$95$2, N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
t_2 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq -0.1:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.99:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites83.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.10000000000000001 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.98999999999999999Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.1
Applied rewrites98.1%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.98999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.6
Applied rewrites84.6%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites79.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6494.4
Applied rewrites94.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im)))
(t_1 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -0.1)
(cos im)
(if (<= t_0 0.0)
(* 1.0 (* (* im im) -0.5))
(if (<= t_0 0.99)
(cos im)
(*
t_1
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -0.1) {
tmp = cos(im);
} else if (t_0 <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else if (t_0 <= 0.99) {
tmp = cos(im);
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= -0.1) tmp = cos(im); elseif (t_0 <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.99) tmp = cos(im); else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); 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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.1], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.99], N[Cos[im], $MachinePrecision], N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq -0.1:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.99:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites83.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.10000000000000001 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.98999999999999999Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6497.3
Applied rewrites97.3%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.98999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.6
Applied rewrites84.6%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites79.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6494.4
Applied rewrites94.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.1)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= t_0 0.0)
(* (exp re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.99)
(*
(fma (fma (/ -0.25 (fma 0.16666666666666666 re -0.5)) re 1.0) re 1.0)
(cos im))
(*
(exp re)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.1) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
} else if (t_0 <= 0.0) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.99) {
tmp = fma(fma((-0.25 / fma(0.16666666666666666, re, -0.5)), re, 1.0), re, 1.0) * cos(im);
} else {
tmp = exp(re) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.1) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)); elseif (t_0 <= 0.0) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.99) tmp = Float64(fma(fma(Float64(-0.25 / fma(0.16666666666666666, re, -0.5)), re, 1.0), re, 1.0) * cos(im)); else tmp = Float64(exp(re) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 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.1], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.99], N[(N[(N[(N[(-0.25 / N[(0.16666666666666666 * re + -0.5), $MachinePrecision]), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], N[(N[Exp[re], $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.1:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.99:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\frac{-0.25}{\mathsf{fma}\left(0.16666666666666666, re, -0.5\right)}, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.10000000000000001Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.6
Applied rewrites90.6%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.98999999999999999Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6495.9
Applied rewrites95.9%
Applied rewrites95.9%
Taylor expanded in re around 0
Applied rewrites96.4%
if 0.98999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6499.2
Applied rewrites99.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.1)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= t_0 0.0)
(* (exp re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.99)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(*
(exp re)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.1) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
} else if (t_0 <= 0.0) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.99) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else {
tmp = exp(re) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= -0.1) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)); elseif (t_0 <= 0.0) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.99) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); else tmp = Float64(exp(re) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 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.1], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.99], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], N[(N[Exp[re], $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.1:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.99:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.10000000000000001Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.6
Applied rewrites90.6%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.98999999999999999Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6496.4
Applied rewrites96.4%
if 0.98999999999999999 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6499.2
Applied rewrites99.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im)))
(t_1 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 0.01)
(* 1.0 (* (* im im) -0.5))
(*
t_1
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 0.01) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 0.01) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); 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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.01], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.01:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites83.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6457.9
Applied rewrites57.9%
Taylor expanded in re around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites15.5%
if 0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6471.0
Applied rewrites71.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites67.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6479.5
Applied rewrites79.5%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
1.0
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 0.01)
(* 1.0 (* (* im im) -0.5))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = 1.0 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 0.01) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(1.0 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 0.01) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 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[(1.0 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.01], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.01:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites48.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6491.3
Applied rewrites91.3%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6457.9
Applied rewrites57.9%
Taylor expanded in re around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites15.5%
if 0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6471.0
Applied rewrites71.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites67.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6479.5
Applied rewrites79.5%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
1.0
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 0.01)
(* 1.0 (* (* im im) -0.5))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* 0.041666666666666664 (* im im)) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = 1.0 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 0.01) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((0.041666666666666664 * (im * im)), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(1.0 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 0.01) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(0.041666666666666664 * Float64(im * im)), Float64(im * im), 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[(1.0 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.01], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.01:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites48.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6491.3
Applied rewrites91.3%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6457.9
Applied rewrites57.9%
Taylor expanded in re around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites15.5%
if 0.0100000000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6471.0
Applied rewrites71.0%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites67.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6479.5
Applied rewrites79.5%
Taylor expanded in im around inf
Applied rewrites79.1%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 0.0) (* 1.0 (* (* im im) -0.5)) (* 1.0 (fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = 1.0 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(1.0 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(1.0 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.4
Applied rewrites62.4%
Taylor expanded in re around 0
Applied rewrites7.4%
Taylor expanded in im around inf
Applied rewrites18.8%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6470.5
Applied rewrites70.5%
Taylor expanded in re around 0
Applied rewrites48.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6457.1
Applied rewrites57.1%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 0.0)
(* 1.0 (* (* im im) -0.5))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.7
Applied rewrites61.7%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites57.6%
(FPCore (re im)
:precision binary64
(if (<= re -5.2)
(* (exp re) (* (* im im) -0.5))
(if (<= re 420000.0)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= re 1e+103)
(* (exp re) (fma (* im im) -0.5 1.0))
(* (* (* (fma 0.16666666666666666 re 0.5) re) re) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -5.2) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 420000.0) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1e+103) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -5.2) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 420000.0) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1e+103) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * cos(im)); end return tmp end
code[re_, im_] := If[LessEqual[re, -5.2], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 420000.0], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1e+103], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -5.2:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 420000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 10^{+103}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot re\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -5.20000000000000018Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in im around inf
Applied rewrites79.7%
if -5.20000000000000018 < re < 4.2e5Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6497.7
Applied rewrites97.7%
if 4.2e5 < re < 1e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6476.9
Applied rewrites76.9%
if 1e103 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in re around inf
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (fma (* im im) -0.5 1.0))))
(if (<= re -0.00011)
t_0
(if (<= re 420000.0)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(if (<= re 1e+103)
t_0
(* (* (* (fma 0.16666666666666666 re 0.5) re) re) (cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * fma((im * im), -0.5, 1.0);
double tmp;
if (re <= -0.00011) {
tmp = t_0;
} else if (re <= 420000.0) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1e+103) {
tmp = t_0;
} else {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)) tmp = 0.0 if (re <= -0.00011) tmp = t_0; elseif (re <= 420000.0) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1e+103) tmp = t_0; else tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * cos(im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[re, -0.00011], t$95$0, If[LessEqual[re, 420000.0], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1e+103], t$95$0, N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{if}\;re \leq -0.00011:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 420000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 10^{+103}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot re\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1.10000000000000004e-4 or 4.2e5 < re < 1e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.7
Applied rewrites78.7%
if -1.10000000000000004e-4 < re < 4.2e5Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6498.5
Applied rewrites98.5%
if 1e103 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in re around inf
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 0.0)
(* 1.0 (* (* im im) -0.5))
(*
(fma (* (fma 0.16666666666666666 re 0.5) re) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.7
Applied rewrites61.7%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites57.6%
Taylor expanded in re around inf
Applied rewrites57.0%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 0.0)
(* 1.0 (* (* im im) -0.5))
(*
(fma (* (* 0.16666666666666666 re) re) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma(((0.16666666666666666 * re) * re), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(Float64(Float64(0.16666666666666666 * re) * re), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re), $MachinePrecision] * re), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(0.16666666666666666 \cdot re\right) \cdot re, re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.7
Applied rewrites61.7%
Taylor expanded in re around 0
+-commutativeN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
rgt-mult-inverseN/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites57.6%
Taylor expanded in re around inf
Applied rewrites57.0%
Applied rewrites57.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (fma (* im im) -0.5 1.0))))
(if (<= re -0.00011)
t_0
(if (<= re 420000.0)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(if (<= re 1.35e+154) t_0 (* (* (* re re) 0.5) (cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * fma((im * im), -0.5, 1.0);
double tmp;
if (re <= -0.00011) {
tmp = t_0;
} else if (re <= 420000.0) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1.35e+154) {
tmp = t_0;
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)) tmp = 0.0 if (re <= -0.00011) tmp = t_0; elseif (re <= 420000.0) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1.35e+154) tmp = t_0; else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[re, -0.00011], t$95$0, If[LessEqual[re, 420000.0], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.35e+154], t$95$0, N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{if}\;re \leq -0.00011:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 420000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1.10000000000000004e-4 or 4.2e5 < re < 1.35000000000000003e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6476.5
Applied rewrites76.5%
if -1.10000000000000004e-4 < re < 4.2e5Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6498.5
Applied rewrites98.5%
if 1.35000000000000003e154 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in re around inf
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (fma (* im im) -0.5 1.0))))
(if (<= re -0.0001)
t_0
(if (<= re 420000.0)
(* (+ 1.0 re) (cos im))
(if (<= re 1.35e+154) t_0 (* (* (* re re) 0.5) (cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * fma((im * im), -0.5, 1.0);
double tmp;
if (re <= -0.0001) {
tmp = t_0;
} else if (re <= 420000.0) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 1.35e+154) {
tmp = t_0;
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)) tmp = 0.0 if (re <= -0.0001) tmp = t_0; elseif (re <= 420000.0) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 1.35e+154) tmp = t_0; else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[re, -0.0001], t$95$0, If[LessEqual[re, 420000.0], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.35e+154], t$95$0, N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{if}\;re \leq -0.0001:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 420000:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1.00000000000000005e-4 or 4.2e5 < re < 1.35000000000000003e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6476.5
Applied rewrites76.5%
if -1.00000000000000005e-4 < re < 4.2e5Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.5
Applied rewrites98.5%
if 1.35000000000000003e154 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in re around inf
Applied rewrites100.0%
(FPCore (re im) :precision binary64 (if (<= (exp re) 0.01) (* 1.0 (* (* im im) -0.5)) (* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.01) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 0.01) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.01], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0.01:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.7
Applied rewrites78.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites19.9%
if 0.0100000000000000002 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.0
Applied rewrites62.0%
Taylor expanded in re around 0
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
cube-multN/A
rgt-mult-inverseN/A
associate-*l*N/A
unpow2N/A
cube-multN/A
distribute-lft-inN/A
*-commutativeN/A
Applied rewrites56.2%
(FPCore (re im) :precision binary64 (if (<= (exp re) 0.0) (* 1.0 (* (* im im) -0.5)) (* (+ 1.0 re) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 0.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.7
Applied rewrites79.7%
Taylor expanded in re around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites20.2%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.7
Applied rewrites61.7%
Taylor expanded in re around 0
lower-+.f6443.9
Applied rewrites43.9%
(FPCore (re im) :precision binary64 (if (<= re -11200000.0) (* 1.0 (* (* im im) -0.5)) (* 1.0 (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -11200000.0) {
tmp = 1.0 * ((im * im) * -0.5);
} else {
tmp = 1.0 * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -11200000.0) tmp = Float64(1.0 * Float64(Float64(im * im) * -0.5)); else tmp = Float64(1.0 * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -11200000.0], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -11200000:\\
\;\;\;\;1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;1 \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -1.12e7Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6480.6
Applied rewrites80.6%
Taylor expanded in re around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites20.7%
if -1.12e7 < re Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6461.6
Applied rewrites61.6%
Taylor expanded in re around 0
Applied rewrites40.6%
(FPCore (re im) :precision binary64 (* 1.0 (* (* im im) -0.5)))
double code(double re, double im) {
return 1.0 * ((im * im) * -0.5);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = 1.0d0 * ((im * im) * (-0.5d0))
end function
public static double code(double re, double im) {
return 1.0 * ((im * im) * -0.5);
}
def code(re, im): return 1.0 * ((im * im) * -0.5)
function code(re, im) return Float64(1.0 * Float64(Float64(im * im) * -0.5)) end
function tmp = code(re, im) tmp = 1.0 * ((im * im) * -0.5); end
code[re_, im_] := N[(1.0 * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6466.9
Applied rewrites66.9%
Taylor expanded in re around 0
Applied rewrites29.9%
Taylor expanded in im around inf
Applied rewrites9.4%
herbie shell --seed 2024315
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))