
(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 24 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 (* (+ 1.0 re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))
(if (<= t_0 -0.05)
t_1
(if (<= t_0 0.0)
(* (+ 1.0 re) (* (pow im 4.0) 0.041666666666666664))
(if (<= t_0 0.998)
t_1
(*
(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 t_1 = (1.0 + re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= -0.05) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (1.0 + re) * (pow(im, 4.0) * 0.041666666666666664);
} else if (t_0 <= 0.998) {
tmp = t_1;
} 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)) t_1 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= -0.05) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(Float64(1.0 + re) * Float64((im ^ 4.0) * 0.041666666666666664)); elseif (t_0 <= 0.998) tmp = t_1; 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]}, Block[{t$95$1 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.05], t$95$1, If[LessEqual[t$95$0, 0.0], N[(N[(1.0 + re), $MachinePrecision] * N[(N[Power[im, 4.0], $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.998], t$95$1, 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\\
t_1 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(1 + re\right) \cdot \left({im}^{4} \cdot 0.041666666666666664\right)\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;t\_1\\
\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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6480.5
Applied rewrites80.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.050000000000000003 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
Taylor expanded in re around 0
lower-+.f6496.6
Applied rewrites96.6%
if -0.050000000000000003 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-+.f642.2
Applied rewrites2.2%
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-*.f641.8
Applied rewrites1.8%
Taylor expanded in im around inf
Applied rewrites29.5%
if 0.998 < (*.f64 (exp.f64 re) (cos.f64 im)) 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.f6486.8
Applied rewrites86.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-*.f6490.1
Applied rewrites90.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))) (t_1 (* (+ 1.0 re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))
(if (<= t_0 -0.05)
t_1
(if (<= t_0 0.0)
(* (* im -0.5) im)
(if (<= t_0 0.998)
t_1
(*
(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 t_1 = (1.0 + re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= -0.05) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (im * -0.5) * im;
} else if (t_0 <= 0.998) {
tmp = t_1;
} 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)) t_1 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= -0.05) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(Float64(im * -0.5) * im); elseif (t_0 <= 0.998) tmp = t_1; 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]}, Block[{t$95$1 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.05], t$95$1, If[LessEqual[t$95$0, 0.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[t$95$0, 0.998], t$95$1, 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\\
t_1 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;t\_1\\
\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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6480.5
Applied rewrites80.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.050000000000000003 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
Taylor expanded in re around 0
lower-+.f6496.6
Applied rewrites96.6%
if -0.050000000000000003 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if 0.998 < (*.f64 (exp.f64 re) (cos.f64 im)) 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.f6486.8
Applied rewrites86.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-*.f6490.1
Applied rewrites90.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))
(if (<= t_0 -0.05)
(cos im)
(if (<= t_0 0.0)
(* (* im -0.5) im)
(if (<= t_0 0.998)
(cos im)
(*
(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 = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= -0.05) {
tmp = cos(im);
} else if (t_0 <= 0.0) {
tmp = (im * -0.5) * im;
} else if (t_0 <= 0.998) {
tmp = cos(im);
} 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(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= -0.05) tmp = cos(im); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * -0.5) * im); elseif (t_0 <= 0.998) tmp = cos(im); 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[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.05], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[t$95$0, 0.998], N[Cos[im], $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:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq -0.05:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;\cos im\\
\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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6480.5
Applied rewrites80.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.050000000000000003 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
Taylor expanded in re around 0
lower-cos.f6495.8
Applied rewrites95.8%
if -0.050000000000000003 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if 0.998 < (*.f64 (exp.f64 re) (cos.f64 im)) 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.f6486.8
Applied rewrites86.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-*.f6490.1
Applied rewrites90.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.05)
(*
(fma
(fma
(/
(* (* re re) 0.027777777777777776)
(fma 0.16666666666666666 re -0.5))
re
1.0)
re
1.0)
(cos im))
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.998)
(*
(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.05) {
tmp = fma(fma((((re * re) * 0.027777777777777776) / fma(0.16666666666666666, re, -0.5)), re, 1.0), re, 1.0) * cos(im);
} else if (t_0 <= 0.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= 0.998) {
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.05) tmp = Float64(fma(fma(Float64(Float64(Float64(re * re) * 0.027777777777777776) / fma(0.16666666666666666, re, -0.5)), re, 1.0), re, 1.0) * cos(im)); elseif (t_0 <= 0.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.998) 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.05], N[(N[(N[(N[(N[(N[(re * re), $MachinePrecision] * 0.027777777777777776), $MachinePrecision] / N[(0.16666666666666666 * re + -0.5), $MachinePrecision]), $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), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.998], 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.05:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\frac{\left(re \cdot re\right) \cdot 0.027777777777777776}{\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 \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;\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.050000000000000003Initial 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.8
Applied rewrites97.8%
Applied rewrites97.8%
Taylor expanded in re around inf
Applied rewrites97.8%
if -0.050000000000000003 < (*.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-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
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.9
Applied rewrites97.9%
Applied rewrites97.9%
Taylor expanded in re around 0
Applied rewrites98.0%
if 0.998 < (*.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-*.f64100.0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.05)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.998)
(*
(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.05) {
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) * ((im * im) * -0.5);
} else if (t_0 <= 0.998) {
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.05) 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) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.998) 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.05], 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), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.998], 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.05:\\
\;\;\;\;\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 \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;\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.050000000000000003Initial 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.8
Applied rewrites97.8%
if -0.050000000000000003 < (*.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-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
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.9
Applied rewrites97.9%
Applied rewrites97.9%
Taylor expanded in re around 0
Applied rewrites98.0%
if 0.998 < (*.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-*.f64100.0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.05)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.998)
(*
(+ re (fma (* re re) (fma 0.16666666666666666 re 0.5) 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.05) {
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) * ((im * im) * -0.5);
} else if (t_0 <= 0.998) {
tmp = (re + fma((re * re), fma(0.16666666666666666, re, 0.5), 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.05) 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) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.998) tmp = Float64(Float64(re + fma(Float64(re * re), fma(0.16666666666666666, re, 0.5), 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.05], 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), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.998], N[(N[(re + N[(N[(re * re), $MachinePrecision] * N[(0.16666666666666666 * re + 0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $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.05:\\
\;\;\;\;\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 \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq 0.998:\\
\;\;\;\;\left(re + \mathsf{fma}\left(re \cdot re, \mathsf{fma}\left(0.16666666666666666, re, 0.5\right), 1\right)\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.050000000000000003Initial 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.8
Applied rewrites97.8%
if -0.050000000000000003 < (*.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-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.998Initial program 99.9%
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.9
Applied rewrites97.9%
Applied rewrites97.9%
if 0.998 < (*.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-*.f64100.0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im -0.5) im)
(*
(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 = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * -0.5) * im;
} 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(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * -0.5) * im); 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[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * -0.5), $MachinePrecision] * im), $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:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6480.5
Applied rewrites80.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6432.3
Applied rewrites32.3%
Taylor expanded in im around 0
Applied rewrites3.0%
Taylor expanded in im around inf
Applied rewrites16.7%
Applied rewrites16.7%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) 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.f6489.4
Applied rewrites89.4%
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-*.f6469.8
Applied rewrites69.8%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))) (t_1 (fma (fma 0.5 re 1.0) re 1.0)))
(if (<= t_0 (- INFINITY))
(* t_1 (fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im -0.5) 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(0.5, re, 1.0), re, 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * -0.5) * 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(0.5, re, 1.0), re, 1.0) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * -0.5) * 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[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * -0.5), $MachinePrecision] * 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(0.5, re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot 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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6480.5
Applied rewrites80.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6432.3
Applied rewrites32.3%
Taylor expanded in im around 0
Applied rewrites3.0%
Taylor expanded in im around inf
Applied rewrites16.7%
Applied rewrites16.7%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6482.9
Applied rewrites82.9%
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-*.f6465.4
Applied rewrites65.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im -0.5) im)
(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 + re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * -0.5) * im;
} else {
tmp = 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(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * -0.5) * im); else tmp = 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[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\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 re around 0
lower-+.f646.0
Applied rewrites6.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6473.5
Applied rewrites73.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6432.3
Applied rewrites32.3%
Taylor expanded in im around 0
Applied rewrites3.0%
Taylor expanded in im around inf
Applied rewrites16.7%
Applied rewrites16.7%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6469.9
Applied rewrites69.9%
Taylor expanded in im around 0
Applied rewrites55.7%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 0.1) (* (* im -0.5) im) (fma (* im im) -0.5 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 0.1) {
tmp = (im * -0.5) * im;
} else {
tmp = fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 0.1) tmp = Float64(Float64(im * -0.5) * im); else tmp = fma(Float64(im * im), -0.5, 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 0.1], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \cos im \leq 0.1:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, -0.5, 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
lower-cos.f6429.2
Applied rewrites29.2%
Taylor expanded in im around 0
Applied rewrites6.8%
Taylor expanded in im around inf
Applied rewrites18.6%
Applied rewrites18.6%
if 0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6469.7
Applied rewrites69.7%
Taylor expanded in im around 0
Applied rewrites47.8%
(FPCore (re im)
:precision binary64
(if (<= re -1.6)
(* (exp re) (* (* im im) -0.5))
(if (<= re 255.0)
(* (+ re (fma (* re re) (fma 0.16666666666666666 re 0.5) 1.0)) (cos im))
(if (<= re 1.05e+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 <= -1.6) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 255.0) {
tmp = (re + fma((re * re), fma(0.16666666666666666, re, 0.5), 1.0)) * cos(im);
} else if (re <= 1.05e+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 <= -1.6) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 255.0) tmp = Float64(Float64(re + fma(Float64(re * re), fma(0.16666666666666666, re, 0.5), 1.0)) * cos(im)); elseif (re <= 1.05e+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, -1.6], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 255.0], N[(N[(re + N[(N[(re * re), $MachinePrecision] * N[(0.16666666666666666 * re + 0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.05e+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 -1.6:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 255:\\
\;\;\;\;\left(re + \mathsf{fma}\left(re \cdot re, \mathsf{fma}\left(0.16666666666666666, re, 0.5\right), 1\right)\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.05 \cdot 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 < -1.6000000000000001Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -1.6000000000000001 < re < 255Initial 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.f6498.7
Applied rewrites98.7%
Applied rewrites98.7%
if 255 < re < 1.0500000000000001e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
if 1.0500000000000001e103 < 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 (<= re -1.6)
(* (exp re) (* (* im im) -0.5))
(if (<= re 255.0)
(* (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) (cos im))
(if (<= re 1.05e+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 <= -1.6) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 255.0) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1.05e+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 <= -1.6) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 255.0) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1.05e+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, -1.6], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 255.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, 1.05e+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 -1.6:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 255:\\
\;\;\;\;\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 1.05 \cdot 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 < -1.6000000000000001Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -1.6000000000000001 < re < 255Initial 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.f6498.7
Applied rewrites98.7%
if 255 < re < 1.0500000000000001e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
if 1.0500000000000001e103 < 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 (<= re -460.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 255.0)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(if (<= re 1.05e+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 <= -460.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 255.0) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1.05e+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 <= -460.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 255.0) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1.05e+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, -460.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 255.0], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.05e+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 -460:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 255:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.05 \cdot 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 < -460Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -460 < re < 255Initial 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 255 < re < 1.0500000000000001e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
if 1.0500000000000001e103 < 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 (<= re -460.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 255.0)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(if (<= re 1.32e+154)
(* (exp re) (fma (* im im) -0.5 1.0))
(* (* (* re re) 0.5) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -460.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 255.0) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else if (re <= 1.32e+154) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -460.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 255.0) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); elseif (re <= 1.32e+154) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := If[LessEqual[re, -460.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 255.0], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.32e+154], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -460:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 255:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -460Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -460 < re < 255Initial 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 255 < re < 1.31999999999999998e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6472.4
Applied rewrites72.4%
if 1.31999999999999998e154 < 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 (<= re -1.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 255.0)
(* (+ 1.0 re) (cos im))
(if (<= re 1.32e+154)
(* (exp re) (fma (* im im) -0.5 1.0))
(* (* (* re re) 0.5) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -1.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 255.0) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 1.32e+154) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 255.0) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 1.32e+154) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 255.0], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.32e+154], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 255:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -1 < re < 255Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.2
Applied rewrites98.2%
if 255 < re < 1.31999999999999998e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6472.4
Applied rewrites72.4%
if 1.31999999999999998e154 < 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 (<= re -1.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 600.0)
(* (+ 1.0 re) (cos im))
(if (<= re 1.32e+154)
(*
(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))
(* (* (* re re) 0.5) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -1.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 600.0) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 1.32e+154) {
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);
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (re <= 600.0) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 1.32e+154) 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)); else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.0], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 600.0], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.32e+154], 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], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;re \leq 600:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;\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)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.8
Applied rewrites68.8%
Taylor expanded in im around inf
Applied rewrites68.8%
if -1 < re < 600Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.2
Applied rewrites98.2%
if 600 < re < 1.31999999999999998e154Initial 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.f6448.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-*.f6450.7
Applied rewrites50.7%
if 1.31999999999999998e154 < 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 (<= re -1500000000.0)
(* (+ 1.0 re) (* (pow im 4.0) 0.041666666666666664))
(if (<= re 600.0)
(* (+ 1.0 re) (cos im))
(if (<= re 1.32e+154)
(*
(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))
(* (* (* re re) 0.5) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -1500000000.0) {
tmp = (1.0 + re) * (pow(im, 4.0) * 0.041666666666666664);
} else if (re <= 600.0) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 1.32e+154) {
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);
} else {
tmp = ((re * re) * 0.5) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1500000000.0) tmp = Float64(Float64(1.0 + re) * Float64((im ^ 4.0) * 0.041666666666666664)); elseif (re <= 600.0) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 1.32e+154) 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)); else tmp = Float64(Float64(Float64(re * re) * 0.5) * cos(im)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1500000000.0], N[(N[(1.0 + re), $MachinePrecision] * N[(N[Power[im, 4.0], $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 600.0], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.32e+154], 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], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1500000000:\\
\;\;\;\;\left(1 + re\right) \cdot \left({im}^{4} \cdot 0.041666666666666664\right)\\
\mathbf{elif}\;re \leq 600:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;\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)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -1.5e9Initial program 100.0%
Taylor expanded in re around 0
lower-+.f642.2
Applied rewrites2.2%
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-*.f641.8
Applied rewrites1.8%
Taylor expanded in im around inf
Applied rewrites29.9%
if -1.5e9 < re < 600Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6497.5
Applied rewrites97.5%
if 600 < re < 1.31999999999999998e154Initial 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.f6448.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-*.f6450.7
Applied rewrites50.7%
if 1.31999999999999998e154 < 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.0) (* (* im -0.5) im) (* (+ 1.0 re) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = (im * -0.5) * im;
} 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(Float64(im * -0.5) * im); 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[(N[(im * -0.5), $MachinePrecision] * im), $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:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\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 re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6470.2
Applied rewrites70.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6443.3
Applied rewrites43.3%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (* im -0.5) im)
(if (<= re 1.6e+111)
(*
(+ 1.0 re)
(fma (* (fma 0.041666666666666664 (* im im) -0.5) im) im 1.0))
(*
(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 (re <= -1.0) {
tmp = (im * -0.5) * im;
} else if (re <= 1.6e+111) {
tmp = (1.0 + re) * fma((fma(0.041666666666666664, (im * im), -0.5) * im), im, 1.0);
} 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 (re <= -1.0) tmp = Float64(Float64(im * -0.5) * im); elseif (re <= 1.6e+111) tmp = Float64(Float64(1.0 + re) * fma(Float64(fma(0.041666666666666664, Float64(im * im), -0.5) * im), im, 1.0)); 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[re, -1.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[re, 1.6e+111], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $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}\;re \leq -1:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;re \leq 1.6 \cdot 10^{+111}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right) \cdot im, im, 1\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 re < -1Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if -1 < re < 1.6e111Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6486.5
Applied rewrites86.5%
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-*.f6450.8
Applied rewrites50.8%
Applied rewrites50.8%
if 1.6e111 < 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 im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6479.5
Applied rewrites79.5%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (* im -0.5) im)
(if (<= re 7.2e+113)
(*
(+ 1.0 re)
(fma (* (fma 0.041666666666666664 (* im im) -0.5) im) im 1.0))
(* (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 (re <= -1.0) {
tmp = (im * -0.5) * im;
} else if (re <= 7.2e+113) {
tmp = (1.0 + re) * fma((fma(0.041666666666666664, (im * im), -0.5) * im), im, 1.0);
} 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 (re <= -1.0) tmp = Float64(Float64(im * -0.5) * im); elseif (re <= 7.2e+113) tmp = Float64(Float64(1.0 + re) * fma(Float64(fma(0.041666666666666664, Float64(im * im), -0.5) * im), im, 1.0)); 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[re, -1.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[re, 7.2e+113], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $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}\;re \leq -1:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;re \leq 7.2 \cdot 10^{+113}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right) \cdot im, im, 1\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 re < -1Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if -1 < re < 7.19999999999999984e113Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6486.5
Applied rewrites86.5%
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-*.f6450.8
Applied rewrites50.8%
Applied rewrites50.8%
if 7.19999999999999984e113 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6476.7
Applied rewrites76.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6465.6
Applied rewrites65.6%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (* im -0.5) im)
(if (<= re 7.2e+113)
(* (+ 1.0 re) (fma (* 0.041666666666666664 (* im im)) (* im im) 1.0))
(* (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 (re <= -1.0) {
tmp = (im * -0.5) * im;
} else if (re <= 7.2e+113) {
tmp = (1.0 + re) * fma((0.041666666666666664 * (im * im)), (im * im), 1.0);
} 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 (re <= -1.0) tmp = Float64(Float64(im * -0.5) * im); elseif (re <= 7.2e+113) tmp = Float64(Float64(1.0 + re) * fma(Float64(0.041666666666666664 * Float64(im * im)), Float64(im * im), 1.0)); 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[re, -1.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[re, 7.2e+113], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $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}\;re \leq -1:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;re \leq 7.2 \cdot 10^{+113}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right), im \cdot im, 1\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 re < -1Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if -1 < re < 7.19999999999999984e113Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6486.5
Applied rewrites86.5%
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-*.f6450.8
Applied rewrites50.8%
Taylor expanded in im around inf
Applied rewrites50.6%
if 7.19999999999999984e113 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6476.7
Applied rewrites76.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6465.6
Applied rewrites65.6%
(FPCore (re im)
:precision binary64
(if (<= re -390.0)
(* (* im -0.5) im)
(if (<= re 7.2e+113)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0)
(* (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 (re <= -390.0) {
tmp = (im * -0.5) * im;
} else if (re <= 7.2e+113) {
tmp = fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
} 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 (re <= -390.0) tmp = Float64(Float64(im * -0.5) * im); elseif (re <= 7.2e+113) tmp = fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0); 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[re, -390.0], N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[re, 7.2e+113], N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $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}\;re \leq -390:\\
\;\;\;\;\left(im \cdot -0.5\right) \cdot im\\
\mathbf{elif}\;re \leq 7.2 \cdot 10^{+113}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\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 re < -390Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites22.4%
Applied rewrites22.4%
if -390 < re < 7.19999999999999984e113Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6485.9
Applied rewrites85.9%
Taylor expanded in im around 0
Applied rewrites50.6%
if 7.19999999999999984e113 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6476.7
Applied rewrites76.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6465.6
Applied rewrites65.6%
(FPCore (re im) :precision binary64 (* (* im -0.5) im))
double code(double re, double im) {
return (im * -0.5) * im;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (im * (-0.5d0)) * im
end function
public static double code(double re, double im) {
return (im * -0.5) * im;
}
def code(re, im): return (im * -0.5) * im
function code(re, im) return Float64(Float64(im * -0.5) * im) end
function tmp = code(re, im) tmp = (im * -0.5) * im; end
code[re_, im_] := N[(N[(im * -0.5), $MachinePrecision] * im), $MachinePrecision]
\begin{array}{l}
\\
\left(im \cdot -0.5\right) \cdot im
\end{array}
Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6452.6
Applied rewrites52.6%
Taylor expanded in im around 0
Applied rewrites30.5%
Taylor expanded in im around inf
Applied rewrites8.8%
Applied rewrites8.8%
herbie shell --seed 2024314
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))