
(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 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Initial program 100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 -2e-5)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im))
(if (<= t_0 2e-11)
(* (exp re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.9999999999999927)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(cos im))
(*
(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 tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= -2e-5) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
} else if (t_0 <= 2e-11) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.9999999999999927) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
} 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)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= -2e-5) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); elseif (t_0 <= 2e-11) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.9999999999999927) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im)); 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]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -2e-5], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-11], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.9999999999999927], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $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]), $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:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -2 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-11}:\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0.9999999999999927:\\
\;\;\;\;\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{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)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -2.00000000000000016e-5Initial program 99.9%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6495.0
Applied rewrites95.0%
if -2.00000000000000016e-5 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6466.7
Applied rewrites66.7%
if 1.99999999999999988e-11 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99999999999999267Initial 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.3
Applied rewrites95.3%
if 0.99999999999999267 < (*.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 (* (+ 1.0 re) (cos im)))
(t_1 (* (exp re) (cos im)))
(t_2 (* (exp re) (* (* im im) -0.5))))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 -0.001)
t_0
(if (<= t_1 2e-11)
t_2
(if (<= t_1 2.0)
t_0
(*
(* (* (fma 0.16666666666666666 re 0.5) re) re)
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = (1.0 + re) * cos(im);
double t_1 = exp(re) * cos(im);
double t_2 = exp(re) * ((im * im) * -0.5);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= -0.001) {
tmp = t_0;
} else if (t_1 <= 2e-11) {
tmp = t_2;
} else if (t_1 <= 2.0) {
tmp = t_0;
} else {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(Float64(1.0 + re) * cos(im)) t_1 = Float64(exp(re) * cos(im)) t_2 = Float64(exp(re) * Float64(Float64(im * im) * -0.5)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= -0.001) tmp = t_0; elseif (t_1 <= 2e-11) tmp = t_2; elseif (t_1 <= 2.0) tmp = t_0; else tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * 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[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, -0.001], t$95$0, If[LessEqual[t$95$1, 2e-11], t$95$2, If[LessEqual[t$95$1, 2.0], t$95$0, N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $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 := \left(1 + re\right) \cdot \cos im\\
t_1 := e^{re} \cdot \cos im\\
t_2 := e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq -0.001:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-11}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 2:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\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)) < -inf.0 or -1e-3 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6471.1
Applied rewrites71.1%
Taylor expanded in im around inf
Applied rewrites68.3%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -1e-3 or 1.99999999999999988e-11 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6497.7
Applied rewrites97.7%
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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6461.6
Applied rewrites61.6%
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-*.f6474.6
Applied rewrites74.6%
Taylor expanded in re around inf
Applied rewrites74.6%
(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 (fma 0.16666666666666666 re 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.001)
t_1
(if (<= t_0 2e-11)
(* (* im im) -0.5)
(if (<= t_0 2.0)
t_1
(*
(* (* (fma 0.16666666666666666 re 0.5) re) 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 t_1 = (1.0 + re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(fma(0.16666666666666666, re, 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.001) {
tmp = t_1;
} else if (t_0 <= 2e-11) {
tmp = (im * im) * -0.5;
} else if (t_0 <= 2.0) {
tmp = t_1;
} else {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * 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)) t_1 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 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.001) tmp = t_1; elseif (t_0 <= 2e-11) tmp = Float64(Float64(im * im) * -0.5); elseif (t_0 <= 2.0) tmp = t_1; else tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * 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]}, Block[{t$95$1 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * 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.001], t$95$1, If[LessEqual[t$95$0, 2e-11], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[t$95$0, 2.0], t$95$1, N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $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\\
t_1 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\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(\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.001:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-11}:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;t\_0 \leq 2:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\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)) < -inf.0Initial 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.f6463.2
Applied rewrites63.2%
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-*.f6492.8
Applied rewrites92.8%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -1e-3 or 1.99999999999999988e-11 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6497.7
Applied rewrites97.7%
if -1e-3 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.8
Applied rewrites3.8%
Taylor expanded in im around 0
Applied rewrites2.7%
Taylor expanded in im around inf
Applied rewrites21.8%
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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6461.6
Applied rewrites61.6%
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-*.f6474.6
Applied rewrites74.6%
Taylor expanded in re around inf
Applied rewrites74.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)))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 -2e-5)
(cos im)
(if (<= t_0 0.0)
(* (* im im) -0.5)
(if (<= t_0 0.999998)
(cos im)
(*
t_1
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= -2e-5) {
tmp = cos(im);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else if (t_0 <= 0.999998) {
tmp = cos(im);
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= -2e-5) tmp = cos(im); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); elseif (t_0 <= 0.999998) tmp = cos(im); else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -2e-5], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[t$95$0, 0.999998], N[Cos[im], $MachinePrecision], N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq -2 \cdot 10^{-5}:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;t\_0 \leq 0.999998:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6463.2
Applied rewrites63.2%
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-*.f6492.8
Applied rewrites92.8%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -2.00000000000000016e-5 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.999998000000000054Initial program 99.9%
Taylor expanded in re around 0
lower-cos.f6490.2
Applied rewrites90.2%
if -2.00000000000000016e-5 < (*.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 rewrites23.7%
if 0.999998000000000054 < (*.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.2
Applied rewrites86.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-*.f6490.9
Applied rewrites90.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))
(t_1 (* (exp re) (cos im))))
(if (<= t_1 (- INFINITY))
(*
(+ 1.0 re)
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_1 0.0)
(* (* im im) -0.5)
(if (<= t_1 2.0)
(* (fma (fma 0.5 re 1.0) re 1.0) t_0)
(* (* (* (fma 0.16666666666666666 re 0.5) re) re) t_0))))))
double code(double re, double im) {
double t_0 = fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = (1.0 + re) * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_1 <= 0.0) {
tmp = (im * im) * -0.5;
} else if (t_1 <= 2.0) {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * t_0;
} else {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * t_0;
}
return tmp;
}
function code(re, im) t_0 = fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(Float64(1.0 + re) * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_1 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); elseif (t_1 <= 2.0) tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * t_0); else tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * t_0); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(N[(1.0 + re), $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$1, 0.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[t$95$1, 2.0], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * t$95$0), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * t$95$0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\left(1 + re\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\_1 \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;t\_1 \leq 2:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot re\right) \cdot t\_0\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
lower-+.f644.9
Applied rewrites4.9%
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-*.f6479.1
Applied rewrites79.1%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6447.4
Applied rewrites47.4%
Taylor expanded in im around 0
Applied rewrites3.6%
Taylor expanded in im around inf
Applied rewrites14.5%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6495.4
Applied rewrites95.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-*.f6464.2
Applied rewrites64.2%
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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6461.6
Applied rewrites61.6%
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-*.f6474.6
Applied rewrites74.6%
Taylor expanded in re around inf
Applied rewrites74.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)))
(if (<= t_0 (- INFINITY))
(*
t_1
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 2e-11)
(* (* im im) -0.5)
(*
t_1
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = t_1 * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 2e-11) {
tmp = (im * im) * -0.5;
} else {
tmp = t_1 * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(t_1 * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 2e-11) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(t_1 * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(t$95$1 * N[(N[(N[(-0.001388888888888889 * N[(im * im), $MachinePrecision] + 0.041666666666666664), $MachinePrecision] * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e-11], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(t$95$1 * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001388888888888889, im \cdot im, 0.041666666666666664\right), im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-11}:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6463.2
Applied rewrites63.2%
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-*.f6492.8
Applied rewrites92.8%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6446.3
Applied rewrites46.3%
Taylor expanded in im around 0
Applied rewrites3.7%
Taylor expanded in im around inf
Applied rewrites14.2%
if 1.99999999999999988e-11 < (*.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.f6488.6
Applied rewrites88.6%
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-*.f6468.1
Applied rewrites68.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(fma (* re re) 0.5 re)
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 2e-11)
(* (* im im) -0.5)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma((re * re), 0.5, re) * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 2e-11) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(Float64(re * re), 0.5, re) * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 2e-11) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(re * re), $MachinePrecision] * 0.5 + re), $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, 2e-11], N[(N[(im * im), $MachinePrecision] * -0.5), $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(re \cdot re, 0.5, re\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 2 \cdot 10^{-11}:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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.f6442.0
Applied rewrites42.0%
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-*.f6492.7
Applied rewrites92.7%
Taylor expanded in re around inf
Applied rewrites92.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6446.3
Applied rewrites46.3%
Taylor expanded in im around 0
Applied rewrites3.7%
Taylor expanded in im around inf
Applied rewrites14.2%
if 1.99999999999999988e-11 < (*.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.f6488.6
Applied rewrites88.6%
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-*.f6468.1
Applied rewrites68.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(+ 1.0 re)
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 2e-11)
(* (* im im) -0.5)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (1.0 + re) * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 2e-11) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(1.0 + re) * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 2e-11) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(1.0 + re), $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, 2e-11], N[(N[(im * im), $MachinePrecision] * -0.5), $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:\\
\;\;\;\;\left(1 + re\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 2 \cdot 10^{-11}:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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
lower-+.f644.9
Applied rewrites4.9%
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-*.f6479.1
Applied rewrites79.1%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 1.99999999999999988e-11Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6446.3
Applied rewrites46.3%
Taylor expanded in im around 0
Applied rewrites3.7%
Taylor expanded in im around inf
Applied rewrites14.2%
if 1.99999999999999988e-11 < (*.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.f6488.6
Applied rewrites88.6%
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-*.f6468.1
Applied rewrites68.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(+ 1.0 re)
(fma
(fma
(fma -0.001388888888888889 (* im im) 0.041666666666666664)
(* im im)
-0.5)
(* im im)
1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(*
(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 <= -((double) INFINITY)) {
tmp = (1.0 + re) * fma(fma(fma(-0.001388888888888889, (im * im), 0.041666666666666664), (im * im), -0.5), (im * im), 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} 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 <= Float64(-Inf)) tmp = Float64(Float64(1.0 + re) * fma(fma(fma(-0.001388888888888889, Float64(im * im), 0.041666666666666664), Float64(im * im), -0.5), Float64(im * im), 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); 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, (-Infinity)], N[(N[(1.0 + re), $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.0], N[(N[(im * im), $MachinePrecision] * -0.5), $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 -\infty:\\
\;\;\;\;\left(1 + re\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:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
lower-+.f644.9
Applied rewrites4.9%
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-*.f6479.1
Applied rewrites79.1%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6447.4
Applied rewrites47.4%
Taylor expanded in im around 0
Applied rewrites3.6%
Taylor expanded in im around inf
Applied rewrites14.5%
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.f6483.4
Applied rewrites83.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.6
Applied rewrites63.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.995)
(*
(* (* (fma 0.16666666666666666 re 0.5) re) re)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(*
(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.995) {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} 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.995) tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); 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.995], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * 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[(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.995:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot 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(\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.994999999999999996Initial 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.f6468.1
Applied rewrites68.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in re around inf
Applied rewrites68.6%
if -0.994999999999999996 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6446.3
Applied rewrites46.3%
Taylor expanded in im around 0
Applied rewrites3.6%
Taylor expanded in im around inf
Applied rewrites14.8%
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.f6483.4
Applied rewrites83.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.6
Applied rewrites63.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.995)
(*
(* (* (fma 0.16666666666666666 re 0.5) re) re)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(*
(fma (fma 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 <= -0.995) {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(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 <= -0.995) tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(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, -0.995], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * 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[(N[(0.5 * 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 -0.995:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot 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(\mathsf{fma}\left(0.5, 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)) < -0.994999999999999996Initial 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.f6468.1
Applied rewrites68.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in re around inf
Applied rewrites68.6%
if -0.994999999999999996 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6446.3
Applied rewrites46.3%
Taylor expanded in im around 0
Applied rewrites3.6%
Taylor expanded in im around inf
Applied rewrites14.8%
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.f6486.7
Applied rewrites86.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-*.f6466.7
Applied rewrites66.7%
Taylor expanded in re around 0
Applied rewrites63.6%
Taylor expanded in im around inf
Applied rewrites63.0%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 0.0) (* (* im im) -0.5) (fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0)))
double code(double re, double im) {
double tmp;
if ((exp(re) * cos(im)) <= 0.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * cos(im)) <= 0.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(im * im), $MachinePrecision] * -0.5), $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}
\mathbf{if}\;e^{re} \cdot \cos im \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6442.2
Applied rewrites42.2%
Taylor expanded in im around 0
Applied rewrites8.0%
Taylor expanded in im around inf
Applied rewrites17.7%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6470.8
Applied rewrites70.8%
Taylor expanded in im around 0
Applied rewrites54.3%
(FPCore (re im) :precision binary64 (if (or (<= re -0.0125) (not (or (<= re 9e+15) (not (<= re 1.05e+103))))) (* (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 <= -0.0125) || !((re <= 9e+15) || !(re <= 1.05e+103))) {
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 <= -0.0125) || !((re <= 9e+15) || !(re <= 1.05e+103))) 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[Or[LessEqual[re, -0.0125], N[Not[Or[LessEqual[re, 9e+15], N[Not[LessEqual[re, 1.05e+103]], $MachinePrecision]]], $MachinePrecision]], 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 -0.0125 \lor \neg \left(re \leq 9 \cdot 10^{+15} \lor \neg \left(re \leq 1.05 \cdot 10^{+103}\right)\right):\\
\;\;\;\;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 < -0.012500000000000001 or 9e15 < re < 1.0500000000000001e103Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6465.1
Applied rewrites65.1%
if -0.012500000000000001 < re < 9e15 or 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.f6498.7
Applied rewrites98.7%
Final simplification88.6%
(FPCore (re im) :precision binary64 (if (or (<= re -0.012) (not (or (<= re 9e+15) (not (<= re 1.9e+154))))) (* (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 <= -0.012) || !((re <= 9e+15) || !(re <= 1.9e+154))) {
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 <= -0.012) || !((re <= 9e+15) || !(re <= 1.9e+154))) 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[Or[LessEqual[re, -0.012], N[Not[Or[LessEqual[re, 9e+15], N[Not[LessEqual[re, 1.9e+154]], $MachinePrecision]]], $MachinePrecision]], 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 -0.012 \lor \neg \left(re \leq 9 \cdot 10^{+15} \lor \neg \left(re \leq 1.9 \cdot 10^{+154}\right)\right):\\
\;\;\;\;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 < -0.012 or 9e15 < re < 1.8999999999999999e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6467.2
Applied rewrites67.2%
if -0.012 < re < 9e15 or 1.8999999999999999e154 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6498.4
Applied rewrites98.4%
Final simplification88.1%
(FPCore (re im) :precision binary64 (if (or (<= re -0.00034) (not (<= re 9e+15))) (* (exp re) (fma (* im im) -0.5 1.0)) (* (+ 1.0 re) (cos im))))
double code(double re, double im) {
double tmp;
if ((re <= -0.00034) || !(re <= 9e+15)) {
tmp = exp(re) * fma((im * im), -0.5, 1.0);
} else {
tmp = (1.0 + re) * cos(im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if ((re <= -0.00034) || !(re <= 9e+15)) tmp = Float64(exp(re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(Float64(1.0 + re) * cos(im)); end return tmp end
code[re_, im_] := If[Or[LessEqual[re, -0.00034], N[Not[LessEqual[re, 9e+15]], $MachinePrecision]], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -0.00034 \lor \neg \left(re \leq 9 \cdot 10^{+15}\right):\\
\;\;\;\;e^{re} \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\end{array}
\end{array}
if re < -3.4e-4 or 9e15 < re Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6469.3
Applied rewrites69.3%
if -3.4e-4 < re < 9e15Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6498.0
Applied rewrites98.0%
Final simplification86.0%
(FPCore (re im) :precision binary64 (if (<= (exp re) 0.2) (* (* im im) -0.5) (* (+ 1.0 re) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.2) {
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 (exp(re) <= 0.2) 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[N[Exp[re], $MachinePrecision], 0.2], 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}\;e^{re} \leq 0.2:\\
\;\;\;\;\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 (exp.f64 re) < 0.20000000000000001Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f644.0
Applied rewrites4.0%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites21.5%
if 0.20000000000000001 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6474.9
Applied rewrites74.9%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6439.6
Applied rewrites39.6%
(FPCore (re im)
:precision binary64
(if (<= re -0.15)
(* (* im im) -0.5)
(if (<= re 5e+30)
(*
(+ 1.0 re)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0))
(* (* (* (* 0.16666666666666666 re) re) re) (fma (* im im) -0.5 1.0)))))
double code(double re, double im) {
double tmp;
if (re <= -0.15) {
tmp = (im * im) * -0.5;
} else if (re <= 5e+30) {
tmp = (1.0 + re) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
} else {
tmp = (((0.16666666666666666 * re) * re) * re) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -0.15) tmp = Float64(Float64(im * im) * -0.5); elseif (re <= 5e+30) tmp = Float64(Float64(1.0 + re) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); else tmp = Float64(Float64(Float64(Float64(0.16666666666666666 * re) * re) * re) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -0.15], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[re, 5e+30], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -0.15:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;re \leq 5 \cdot 10^{+30}:\\
\;\;\;\;\left(1 + re\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(\left(0.16666666666666666 \cdot re\right) \cdot re\right) \cdot re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -0.149999999999999994Initial program 99.9%
Taylor expanded in re around 0
lower-cos.f644.3
Applied rewrites4.3%
Taylor expanded in im around 0
Applied rewrites3.0%
Taylor expanded in im around inf
Applied rewrites21.2%
if -0.149999999999999994 < re < 4.9999999999999998e30Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6493.9
Applied rewrites93.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-*.f6445.4
Applied rewrites45.4%
if 4.9999999999999998e30 < 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.f6473.1
Applied rewrites73.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.9
Applied rewrites62.9%
Taylor expanded in re around inf
Applied rewrites62.9%
Taylor expanded in re around inf
Applied rewrites62.9%
(FPCore (re im)
:precision binary64
(if (<= re -1.6)
(* (* im im) -0.5)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -1.6) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.6) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.6], N[(N[(im * im), $MachinePrecision] * -0.5), $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.6:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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 < -1.6000000000000001Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f644.0
Applied rewrites4.0%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites21.5%
if -1.6000000000000001 < 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.7
Applied rewrites89.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6448.2
Applied rewrites48.2%
(FPCore (re im)
:precision binary64
(if (<= re -430.0)
(* (* im im) -0.5)
(if (<= re 5e+30)
(fma (fma 0.041666666666666664 (* im im) -0.5) (* im im) 1.0)
(* (* (* 0.5 re) re) (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 if (re <= 5e+30) {
tmp = fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
} else {
tmp = ((0.5 * re) * re) * 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); elseif (re <= 5e+30) tmp = fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0); else tmp = Float64(Float64(Float64(0.5 * re) * re) * 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], If[LessEqual[re, 5e+30], N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(N[(0.5 * re), $MachinePrecision] * re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -430:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;re \leq 5 \cdot 10^{+30}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(0.5 \cdot re\right) \cdot re\right) \cdot \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.4%
Taylor expanded in im around inf
Applied rewrites23.7%
if -430 < re < 4.9999999999999998e30Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6489.8
Applied rewrites89.8%
Taylor expanded in im around 0
Applied rewrites43.5%
if 4.9999999999999998e30 < 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.f6473.1
Applied rewrites73.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.9
Applied rewrites62.9%
Taylor expanded in re around inf
Applied rewrites62.9%
Taylor expanded in re around 0
Applied rewrites53.9%
(FPCore (re im)
:precision binary64
(if (<= re -1.8)
(* (* im im) -0.5)
(*
(fma (* (* re re) 0.16666666666666666) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -1.8) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(((re * re) * 0.16666666666666666), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.8) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(Float64(Float64(re * re) * 0.16666666666666666), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.8], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(N[(N[(re * re), $MachinePrecision] * 0.16666666666666666), $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.8:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(re \cdot re\right) \cdot 0.16666666666666666, re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -1.80000000000000004Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f644.0
Applied rewrites4.0%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites21.5%
if -1.80000000000000004 < 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.7
Applied rewrites89.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6448.2
Applied rewrites48.2%
Taylor expanded in re around inf
Applied rewrites48.0%
(FPCore (re im) :precision binary64 (if (<= re -102.0) (* (* im im) -0.5) (* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -102.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -102.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -102.0], N[(N[(im * im), $MachinePrecision] * -0.5), $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 -102:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\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 < -102Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites23.2%
if -102 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.8
Applied rewrites83.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6445.4
Applied rewrites45.4%
(FPCore (re im) :precision binary64 (if (<= re -102.0) (* (* im im) -0.5) (* (fma (* 0.5 re) re 1.0) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -102.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma((0.5 * re), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -102.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(Float64(0.5 * re), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -102.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(N[(0.5 * re), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -102:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(0.5 \cdot re, re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -102Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.4%
Taylor expanded in im around inf
Applied rewrites23.2%
if -102 < re Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.8
Applied rewrites83.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6445.4
Applied rewrites45.4%
Taylor expanded in re around inf
Applied rewrites45.2%
(FPCore (re im) :precision binary64 (if (<= re -420.0) (* (* im im) -0.5) (fma (* im im) -0.5 1.0)))
double code(double re, double im) {
double tmp;
if (re <= -420.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 <= -420.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, -420.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 -420:\\
\;\;\;\;\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 < -420Initial 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 rewrites23.7%
if -420 < re Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6472.0
Applied rewrites72.0%
Taylor expanded in im around 0
Applied rewrites36.8%
(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.f6458.3
Applied rewrites58.3%
Taylor expanded in im around 0
Applied rewrites30.0%
Taylor expanded in im around inf
Applied rewrites8.8%
herbie shell --seed 2024303
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))