
(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 23 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
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -0.1)
t_1
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.9999999995)
t_1
(*
(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 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(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -0.1) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= 0.9999999995) {
tmp = t_1;
} 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)) 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(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_1; elseif (t_0 <= 0.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.9999999995) tmp = t_1; 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]}, 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[(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$1, 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.9999999995], 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]), $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(\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\_1\\
\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.9999999995:\\
\;\;\;\;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(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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6457.6
Applied rewrites57.6%
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.99999999949999996Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.9
Applied rewrites98.9%
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-*.f6474.1
Applied rewrites74.1%
Taylor expanded in im around inf
Applied rewrites74.1%
if 0.99999999949999996 < (*.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%
Taylor expanded in im around inf
Applied rewrites100.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.f6488.2
Applied rewrites88.2%
(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
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -0.1)
t_1
(if (<= t_0 0.0)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(if (<= t_0 0.9999999995)
t_1
(*
(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 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(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -0.1) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} else if (t_0 <= 0.9999999995) {
tmp = t_1;
} 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)) 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(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_1; elseif (t_0 <= 0.0) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; elseif (t_0 <= 0.9999999995) tmp = t_1; 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]}, 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[(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$1, If[LessEqual[t$95$0, 0.0], N[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], If[LessEqual[t$95$0, 0.9999999995], 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]), $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(\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\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{elif}\;t\_0 \leq 0.9999999995:\\
\;\;\;\;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(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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6457.6
Applied rewrites57.6%
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.99999999949999996Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.9
Applied rewrites98.9%
if -0.10000000000000001 < (*.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.5%
Applied rewrites2.5%
Taylor expanded in im around 0
Applied rewrites46.8%
if 0.99999999949999996 < (*.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%
Taylor expanded in im around inf
Applied rewrites100.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.f6488.2
Applied rewrites88.2%
Final simplification82.7%
(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
(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)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(if (<= t_0 0.9999999995)
(cos im)
(*
(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 = fma(fma(0.5, re, 1.0), re, 1.0) * 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 = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} else if (t_0 <= 0.9999999995) {
tmp = cos(im);
} 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(fma(fma(0.5, re, 1.0), re, 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.1) tmp = cos(im); elseif (t_0 <= 0.0) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; elseif (t_0 <= 0.9999999995) tmp = cos(im); 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[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * 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[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], If[LessEqual[t$95$0, 0.9999999995], 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]), $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(\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:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{elif}\;t\_0 \leq 0.9999999995:\\
\;\;\;\;\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(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 re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6457.6
Applied rewrites57.6%
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.99999999949999996Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6498.6
Applied rewrites98.6%
if -0.10000000000000001 < (*.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.5%
Applied rewrites2.5%
Taylor expanded in im around 0
Applied rewrites46.8%
if 0.99999999949999996 < (*.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%
Taylor expanded in im around inf
Applied rewrites100.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.f6488.2
Applied rewrites88.2%
Final simplification82.6%
(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)
(cos im))))
(if (<= t_0 -0.1)
t_1
(if (<= t_0 0.0)
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 0.9999999995)
t_1
(*
(exp re)
(fma (* 0.041666666666666664 (* im im)) (* 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) * cos(im);
double tmp;
if (t_0 <= -0.1) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= 0.9999999995) {
tmp = t_1;
} else {
tmp = exp(re) * fma((0.041666666666666664 * (im * im)), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)) tmp = 0.0 if (t_0 <= -0.1) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= 0.9999999995) tmp = t_1; else tmp = Float64(exp(re) * 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]}, Block[{t$95$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.1], t$95$1, 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.9999999995], t$95$1, N[(N[Exp[re], $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\\
t_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{if}\;t\_0 \leq -0.1:\\
\;\;\;\;t\_1\\
\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.9999999995:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;e^{re} \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)) < -0.10000000000000001 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99999999949999996Initial 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.1
Applied rewrites95.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-*.f6474.1
Applied rewrites74.1%
Taylor expanded in im around inf
Applied rewrites74.1%
if 0.99999999949999996 < (*.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%
Taylor expanded in im around inf
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.1)
(*
(fma (fma 0.5 re 1.0) re 1.0)
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 0.5)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(*
(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 <= -0.1) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 0.5) {
tmp = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} 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 <= -0.1) tmp = Float64(fma(fma(0.5, re, 1.0), re, 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.5) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; 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, -0.1], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * 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.5], N[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $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 -0.1:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \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.5:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-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)) < -0.10000000000000001Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6484.6
Applied rewrites84.6%
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-*.f6438.6
Applied rewrites38.6%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.5Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6412.2
Applied rewrites12.2%
Taylor expanded in im around 0
Applied rewrites2.3%
Applied rewrites2.3%
Taylor expanded in im around 0
Applied rewrites42.7%
if 0.5 < (*.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.f6487.5
Applied rewrites87.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-*.f6471.6
Applied rewrites71.6%
Final simplification57.3%
(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 -0.1)
(* t_1 (fma (* im im) -0.5 1.0))
(if (<= t_0 0.5)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(*
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 <= -0.1) {
tmp = t_1 * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.5) {
tmp = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} 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 <= -0.1) tmp = Float64(t_1 * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.5) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; 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, -0.1], N[(t$95$1 * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.5], N[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $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 -0.1:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.5:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.5Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6412.2
Applied rewrites12.2%
Taylor expanded in im around 0
Applied rewrites2.3%
Applied rewrites2.3%
Taylor expanded in im around 0
Applied rewrites42.7%
if 0.5 < (*.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.f6487.5
Applied rewrites87.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-*.f6471.6
Applied rewrites71.6%
Final simplification57.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)))
(if (<= t_0 -0.1)
(* t_1 (fma (* im im) -0.5 1.0))
(if (<= t_0 0.9999999995)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(* t_1 (fma (* 0.041666666666666664 (* im im)) (* 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 <= -0.1) {
tmp = t_1 * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.9999999995) {
tmp = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} else {
tmp = t_1 * fma((0.041666666666666664 * (im * im)), (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 <= -0.1) tmp = Float64(t_1 * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.9999999995) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; else tmp = Float64(t_1 * 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]}, 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, -0.1], N[(t$95$1 * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.9999999995], N[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(t$95$1 * 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\\
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 -0.1:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.9999999995:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99999999949999996Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6438.7
Applied rewrites38.7%
Taylor expanded in im around 0
Applied rewrites4.0%
Applied rewrites4.0%
Taylor expanded in im around 0
Applied rewrites32.5%
if 0.99999999949999996 < (*.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%
Taylor expanded in im around inf
Applied rewrites100.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.f6488.2
Applied rewrites88.2%
Final simplification57.0%
(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)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.5)
(pow
(fma (fma (fma 0.125 (* im im) 0.25) (* im im) 0.5) (* im im) 1.0)
-1.0)
(*
(fma (fma 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 <= -0.1) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.5) {
tmp = pow(fma(fma(fma(0.125, (im * im), 0.25), (im * im), 0.5), (im * im), 1.0), -1.0);
} else {
tmp = fma(fma(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 <= -0.1) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.5) tmp = fma(fma(fma(0.125, Float64(im * im), 0.25), Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; else tmp = Float64(fma(fma(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, -0.1], 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], If[LessEqual[t$95$0, 0.5], N[Power[N[(N[(N[(0.125 * N[(im * im), $MachinePrecision] + 0.25), $MachinePrecision] * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(0.5 * 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 -0.1:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq 0.5:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.125, im \cdot im, 0.25\right), im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, 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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.5Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6412.2
Applied rewrites12.2%
Taylor expanded in im around 0
Applied rewrites2.3%
Applied rewrites2.3%
Taylor expanded in im around 0
Applied rewrites42.7%
if 0.5 < (*.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.f6483.7
Applied rewrites83.7%
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.3
Applied rewrites69.3%
Final simplification55.8%
(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)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.99998)
(pow (fma (fma 0.25 (* im im) 0.5) (* im im) 1.0) -1.0)
(*
(fma (fma 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 <= -0.1) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.99998) {
tmp = pow(fma(fma(0.25, (im * im), 0.5), (im * im), 1.0), -1.0);
} else {
tmp = fma(fma(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 <= -0.1) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.99998) tmp = fma(fma(0.25, Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; else tmp = Float64(fma(fma(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, -0.1], 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], If[LessEqual[t$95$0, 0.99998], N[Power[N[(N[(0.25 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(0.5 * 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 -0.1:\\
\;\;\;\;\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)\\
\mathbf{elif}\;t\_0 \leq 0.99998:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, 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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99997999999999998Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6437.3
Applied rewrites37.3%
Taylor expanded in im around 0
Applied rewrites1.9%
Applied rewrites1.9%
Taylor expanded in im around 0
Applied rewrites26.9%
if 0.99997999999999998 < (*.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.f6480.4
Applied rewrites80.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-*.f6485.5
Applied rewrites85.5%
Final simplification54.3%
(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)
(fma (* im im) -0.5 1.0))
(if (<= t_0 2.0)
(pow (fma (fma 0.25 (* im im) 0.5) (* im im) 1.0) -1.0)
(*
(* (fma 0.5 re 1.0) 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) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 2.0) {
tmp = pow(fma(fma(0.25, (im * im), 0.5), (im * im), 1.0), -1.0);
} else {
tmp = (fma(0.5, re, 1.0) * 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) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 2.0) tmp = fma(fma(0.25, Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; else tmp = Float64(Float64(fma(0.5, re, 1.0) * 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[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2.0], N[Power[N[(N[(0.25 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * 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 \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.5, re, 1\right) \cdot re\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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6463.0
Applied rewrites63.0%
Taylor expanded in im around 0
Applied rewrites42.7%
Applied rewrites42.7%
Taylor expanded in im around 0
Applied rewrites57.1%
if 2 < (*.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.f6450.4
Applied rewrites50.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-*.f6463.8
Applied rewrites63.8%
Taylor expanded in re around inf
Applied rewrites63.8%
Final simplification54.0%
(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)
(fma (* im im) -0.5 1.0))
(if (<= t_0 2.0)
(pow (fma (fma 0.25 (* im im) 0.5) (* im im) 1.0) -1.0)
(*
(* (* re re) 0.5)
(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) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 2.0) {
tmp = pow(fma(fma(0.25, (im * im), 0.5), (im * im), 1.0), -1.0);
} else {
tmp = ((re * re) * 0.5) * 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) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 2.0) tmp = fma(fma(0.25, Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -1.0; else tmp = Float64(Float64(Float64(re * re) * 0.5) * 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[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2.0], N[Power[N[(N[(0.25 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $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 \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\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)) < -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.f6493.1
Applied rewrites93.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.4
Applied rewrites37.4%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6463.0
Applied rewrites63.0%
Taylor expanded in im around 0
Applied rewrites42.7%
Applied rewrites42.7%
Taylor expanded in im around 0
Applied rewrites57.1%
if 2 < (*.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.f6450.4
Applied rewrites50.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-*.f6463.8
Applied rewrites63.8%
Taylor expanded in re around inf
Applied rewrites63.8%
Final simplification54.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.1)
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.9999999995)
(pow (fma (* im im) 0.5 1.0) -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 <= -0.1) {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.9999999995) {
tmp = pow(fma((im * im), 0.5, 1.0), -1.0);
} else {
tmp = 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 <= -0.1) tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.9999999995) tmp = fma(Float64(im * im), 0.5, 1.0) ^ -1.0; else tmp = 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, -0.1], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.9999999995], N[Power[N[(N[(im * im), $MachinePrecision] * 0.5 + 1.0), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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 -0.1:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.9999999995:\\
\;\;\;\;{\left(\mathsf{fma}\left(im \cdot im, 0.5, 1\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\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)) < -0.10000000000000001Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6465.6
Applied rewrites65.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6432.2
Applied rewrites32.2%
if -0.10000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99999999949999996Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6438.7
Applied rewrites38.7%
Taylor expanded in im around 0
Applied rewrites4.0%
Applied rewrites4.0%
Taylor expanded in im around 0
Applied rewrites21.4%
if 0.99999999949999996 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6460.7
Applied rewrites60.7%
Taylor expanded in im around 0
Applied rewrites71.5%
Taylor expanded in im around inf
Applied rewrites71.5%
Final simplification44.8%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.715)
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(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 <= -0.715) {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else {
tmp = 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 <= -0.715) tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = 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, -0.715], 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 * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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 -0.715:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\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)) < -0.714999999999999969Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6447.5
Applied rewrites47.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6446.9
Applied rewrites46.9%
if -0.714999999999999969 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6426.8
Applied rewrites26.8%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites17.4%
if 0.0 < (*.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 rewrites56.1%
Taylor expanded in im around inf
Applied rewrites55.6%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 0.0)
(pow (fma (fma 0.25 (* im im) 0.5) (* im im) 1.0) -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 (exp(re) <= 0.0) {
tmp = pow(fma(fma(0.25, (im * im), 0.5), (im * im), 1.0), -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 (exp(re) <= 0.0) tmp = fma(fma(0.25, Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -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[N[Exp[re], $MachinePrecision], 0.0], N[Power[N[(N[(0.25 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], -1.0], $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:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-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(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.5%
Applied rewrites2.5%
Taylor expanded in im around 0
Applied rewrites40.4%
if 0.0 < (exp.f64 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.f6489.4
Applied rewrites89.4%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6453.5
Applied rewrites53.5%
Final simplification50.6%
(FPCore (re im) :precision binary64 (if (<= (exp re) 0.0) (pow (fma (fma 0.25 (* im im) 0.5) (* im im) 1.0) -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 (exp(re) <= 0.0) {
tmp = pow(fma(fma(0.25, (im * im), 0.5), (im * im), 1.0), -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 (exp(re) <= 0.0) tmp = fma(fma(0.25, Float64(im * im), 0.5), Float64(im * im), 1.0) ^ -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[N[Exp[re], $MachinePrecision], 0.0], N[Power[N[(N[(0.25 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $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}\;e^{re} \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, im \cdot im, 0.5\right), im \cdot im, 1\right)\right)}^{-1}\\
\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.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.5%
Applied rewrites2.5%
Taylor expanded in im around 0
Applied rewrites40.4%
if 0.0 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6484.5
Applied rewrites84.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6451.5
Applied rewrites51.5%
Final simplification48.9%
(FPCore (re im) :precision binary64 (if (<= (exp re) 1.0) (pow (fma (* im im) 0.5 1.0) -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 (exp(re) <= 1.0) {
tmp = pow(fma((im * im), 0.5, 1.0), -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 (exp(re) <= 1.0) tmp = fma(Float64(im * im), 0.5, 1.0) ^ -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[N[Exp[re], $MachinePrecision], 1.0], N[Power[N[(N[(im * im), $MachinePrecision] * 0.5 + 1.0), $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}\;e^{re} \leq 1:\\
\;\;\;\;{\left(\mathsf{fma}\left(im \cdot im, 0.5, 1\right)\right)}^{-1}\\
\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) < 1Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6470.3
Applied rewrites70.3%
Taylor expanded in im around 0
Applied rewrites36.0%
Applied rewrites36.0%
Taylor expanded in im around 0
Applied rewrites44.2%
if 1 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6452.7
Applied rewrites52.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6452.5
Applied rewrites52.5%
Final simplification46.3%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 0.000104)
(* (+ 1.0 re) (cos im))
(if (<= re 4.8e+99)
(* (exp re) (fma (* im im) -0.5 1.0))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -1.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 0.000104) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 4.8e+99) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * 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 <= 0.000104) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 4.8e+99) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * 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, 0.000104], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 4.8e+99], 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 + 1.0), $MachinePrecision] * re + 1.0), $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 0.000104:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 4.8 \cdot 10^{+99}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 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 \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-*.f6474.1
Applied rewrites74.1%
Taylor expanded in im around inf
Applied rewrites74.1%
if -1 < re < 1.03999999999999994e-4Initial program 100.0%
Taylor expanded in re around 0
lower-+.f64100.0
Applied rewrites100.0%
if 1.03999999999999994e-4 < re < 4.8000000000000002e99Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.3
Applied rewrites78.3%
if 4.8000000000000002e99 < 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.f6497.9
Applied rewrites97.9%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 0.000104)
(* (+ 1.0 re) (cos im))
(if (<= re 1.45e+153)
(* (exp re) (fma (* im im) -0.5 1.0))
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))))))
double code(double re, double im) {
double tmp;
if (re <= -1.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 0.000104) {
tmp = (1.0 + re) * cos(im);
} else if (re <= 1.45e+153) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * 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 <= 0.000104) tmp = Float64(Float64(1.0 + re) * cos(im)); elseif (re <= 1.45e+153) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * 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, 0.000104], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[re, 1.45e+153], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $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 0.000104:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{elif}\;re \leq 1.45 \cdot 10^{+153}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\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-*.f6474.1
Applied rewrites74.1%
Taylor expanded in im around inf
Applied rewrites74.1%
if -1 < re < 1.03999999999999994e-4Initial program 100.0%
Taylor expanded in re around 0
lower-+.f64100.0
Applied rewrites100.0%
if 1.03999999999999994e-4 < re < 1.45000000000000001e153Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.8
Applied rewrites78.8%
if 1.45000000000000001e153 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6497.3
Applied rewrites97.3%
(FPCore (re im)
:precision binary64
(if (<= re -1.0)
(* (exp re) (* (* im im) -0.5))
(if (<= re 0.000104)
(* (+ 1.0 re) (cos im))
(* (exp re) (fma (* im im) -0.5 1.0)))))
double code(double re, double im) {
double tmp;
if (re <= -1.0) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (re <= 0.000104) {
tmp = (1.0 + re) * cos(im);
} else {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
}
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 <= 0.000104) tmp = Float64(Float64(1.0 + re) * cos(im)); else tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); 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, 0.000104], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $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 0.000104:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\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-*.f6474.1
Applied rewrites74.1%
Taylor expanded in im around inf
Applied rewrites74.1%
if -1 < re < 1.03999999999999994e-4Initial program 100.0%
Taylor expanded in re around 0
lower-+.f64100.0
Applied rewrites100.0%
if 1.03999999999999994e-4 < re Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.5
Applied rewrites78.5%
(FPCore (re im) :precision binary64 (if (<= re -1.45) (* (* im im) -0.5) (* (+ 1.0 re) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -1.45) {
tmp = (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 (re <= -1.45) tmp = 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[re, -1.45], N[(N[(im * im), $MachinePrecision] * -0.5), $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}\;re \leq -1.45:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -1.44999999999999996Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites21.7%
if -1.44999999999999996 < re Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6469.4
Applied rewrites69.4%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6443.6
Applied rewrites43.6%
(FPCore (re im) :precision binary64 (if (<= re -430.0) (* (* im im) -0.5) (fma (* im im) -0.5 1.0)))
double code(double re, double im) {
double tmp;
if (re <= -430.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -430.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = fma(Float64(im * im), -0.5, 1.0); end return tmp end
code[re_, im_] := If[LessEqual[re, -430.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -430:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -430Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.1
Applied rewrites3.1%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites21.7%
if -430 < re Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6468.3
Applied rewrites68.3%
Taylor expanded in im around 0
Applied rewrites40.1%
(FPCore (re im) :precision binary64 (* (* im im) -0.5))
double code(double re, double im) {
return (im * im) * -0.5;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (im * im) * (-0.5d0)
end function
public static double code(double re, double im) {
return (im * im) * -0.5;
}
def code(re, im): return (im * im) * -0.5
function code(re, im) return Float64(Float64(im * im) * -0.5) end
function tmp = code(re, im) tmp = (im * im) * -0.5; end
code[re_, im_] := N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(im \cdot im\right) \cdot -0.5
\end{array}
Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6453.5
Applied rewrites53.5%
Taylor expanded in im around 0
Applied rewrites31.6%
Taylor expanded in im around inf
Applied rewrites10.6%
herbie shell --seed 2024322
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))