
(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 22 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 (* (+ 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.004)
t_0
(if (<= t_1 0.0)
t_2
(if (<= t_1 0.995)
t_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 = (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.004) {
tmp = t_0;
} else if (t_1 <= 0.0) {
tmp = t_2;
} else if (t_1 <= 0.995) {
tmp = t_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(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.004) tmp = t_0; elseif (t_1 <= 0.0) tmp = t_2; elseif (t_1 <= 0.995) tmp = t_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[(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.004], t$95$0, If[LessEqual[t$95$1, 0.0], t$95$2, If[LessEqual[t$95$1, 0.995], t$95$0, 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 := \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.004:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0.995:\\
\;\;\;\;t\_0\\
\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.0 or -0.0040000000000000001 < (*.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-*.f6481.4
Applied rewrites81.4%
Taylor expanded in im around inf
Applied rewrites81.4%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0040000000000000001 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.994999999999999996Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.1
Applied rewrites99.1%
if 0.994999999999999996 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6489.2
Applied rewrites89.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-*.f6492.1
Applied rewrites92.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0))
(t_1 (* (exp re) (cos im)))
(t_2 (* (+ 1.0 re) (cos im))))
(if (<= t_1 (- INFINITY))
(* (* (* (fma 0.16666666666666666 re 0.5) re) re) t_0)
(if (<= t_1 -0.004)
t_2
(if (<= t_1 0.0)
(*
(pow
(fma (fma (fma -0.16666666666666666 re 0.5) re -1.0) re 1.0)
-1.0)
t_0)
(if (<= t_1 0.995)
t_2
(*
(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 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double t_2 = (1.0 + re) * cos(im);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * t_0;
} else if (t_1 <= -0.004) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = pow(fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0), -1.0) * t_0;
} else if (t_1 <= 0.995) {
tmp = t_2;
} 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 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) t_2 = Float64(Float64(1.0 + re) * cos(im)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * t_0); elseif (t_1 <= -0.004) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64((fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0) ^ -1.0) * t_0); elseif (t_1 <= 0.995) tmp = t_2; 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[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, -0.004], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[Power[N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.995], t$95$2, 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 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
t_2 := \left(1 + re\right) \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq -0.004:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right), re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0.995:\\
\;\;\;\;t\_2\\
\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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6477.5
Applied rewrites77.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6495.5
Applied rewrites95.5%
Taylor expanded in re around inf
Applied rewrites95.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0040000000000000001 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.994999999999999996Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.1
Applied rewrites99.1%
if -0.0040000000000000001 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.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.f641.7
Applied rewrites1.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f641.6
Applied rewrites1.6%
Applied rewrites1.6%
Taylor expanded in re around 0
Applied rewrites57.0%
if 0.994999999999999996 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6489.2
Applied rewrites89.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-*.f6492.1
Applied rewrites92.1%
Final simplification84.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0)) (t_1 (* (exp re) (cos im))))
(if (<= t_1 (- INFINITY))
(* (* (* (fma 0.16666666666666666 re 0.5) re) re) t_0)
(if (<= t_1 -5e-108)
(cos im)
(if (<= t_1 0.0)
(*
(pow
(fma (fma (fma -0.16666666666666666 re 0.5) re -1.0) re 1.0)
-1.0)
t_0)
(if (<= t_1 0.995)
(cos im)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma
(fma 0.041666666666666664 (* im im) -0.5)
(* im im)
1.0))))))))
double code(double re, double im) {
double t_0 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = ((fma(0.16666666666666666, re, 0.5) * re) * re) * t_0;
} else if (t_1 <= -5e-108) {
tmp = cos(im);
} else if (t_1 <= 0.0) {
tmp = pow(fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0), -1.0) * t_0;
} else if (t_1 <= 0.995) {
tmp = cos(im);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, (im * im), -0.5), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(Float64(Float64(fma(0.16666666666666666, re, 0.5) * re) * re) * t_0); elseif (t_1 <= -5e-108) tmp = cos(im); elseif (t_1 <= 0.0) tmp = Float64((fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0) ^ -1.0) * t_0); elseif (t_1 <= 0.995) tmp = cos(im); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(fma(0.041666666666666664, Float64(im * im), -0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[(im * im), $MachinePrecision] * -0.5 + 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[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, -5e-108], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[Power[N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.995], N[Cos[im], $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re\right) \cdot re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-108}:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right), re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0.995:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right), im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6477.5
Applied rewrites77.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6495.5
Applied rewrites95.5%
Taylor expanded in re around inf
Applied rewrites95.5%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -5e-108 or -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.994999999999999996Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6496.0
Applied rewrites96.0%
if -5e-108 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.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.f641.7
Applied rewrites1.7%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f641.6
Applied rewrites1.6%
Applied rewrites1.6%
Taylor expanded in re around 0
Applied rewrites57.7%
if 0.994999999999999996 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6489.2
Applied rewrites89.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-*.f6492.1
Applied rewrites92.1%
Final simplification84.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.982)
(*
(fma (* (fma 0.16666666666666666 re 0.5) re) re re)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(if (<= t_0 2.0)
(*
(+ 1.0 re)
(fma (* (fma 0.041666666666666664 (* im im) -0.5) im) im 1.0))
(*
(* (fma 0.5 re 1.0) 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 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else if (t_0 <= 2.0) {
tmp = (1.0 + re) * fma((fma(0.041666666666666664, (im * im), -0.5) * im), im, 1.0);
} else {
tmp = (fma(0.5, re, 1.0) * 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 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); elseif (t_0 <= 2.0) tmp = Float64(Float64(1.0 + re) * fma(Float64(fma(0.041666666666666664, Float64(im * im), -0.5) * im), im, 1.0)); else tmp = Float64(Float64(fma(0.5, re, 1.0) * 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, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + 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], If[LessEqual[t$95$0, 2.0], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re), $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.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, 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{elif}\;t\_0 \leq 2:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right) \cdot im, im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.5, re, 1\right) \cdot re\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.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6431.0
Applied rewrites31.0%
Taylor expanded in im around 0
Applied rewrites3.2%
Taylor expanded in im around inf
Applied rewrites20.9%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.0
Applied rewrites99.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-*.f6464.8
Applied rewrites64.8%
Applied rewrites64.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
lower-fma.f6461.9
Applied rewrites61.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-*.f6472.9
Applied rewrites72.9%
Taylor expanded in re around inf
Applied rewrites72.9%
Taylor expanded in im around inf
Applied rewrites72.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im)))
(t_1 (fma 0.041666666666666664 (* im im) -0.5)))
(if (<= t_0 -0.982)
(*
(fma (* (fma 0.16666666666666666 re 0.5) re) re re)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(if (<= t_0 2.0)
(* (+ 1.0 re) (fma (* t_1 im) im 1.0))
(* (* (* re re) 0.5) (fma t_1 (* im im) 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double t_1 = fma(0.041666666666666664, (im * im), -0.5);
double tmp;
if (t_0 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else if (t_0 <= 2.0) {
tmp = (1.0 + re) * fma((t_1 * im), im, 1.0);
} else {
tmp = ((re * re) * 0.5) * fma(t_1, (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) t_1 = fma(0.041666666666666664, Float64(im * im), -0.5) tmp = 0.0 if (t_0 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); elseif (t_0 <= 2.0) tmp = Float64(Float64(1.0 + re) * fma(Float64(t_1 * im), im, 1.0)); else tmp = Float64(Float64(Float64(re * re) * 0.5) * fma(t_1, 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[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision]}, If[LessEqual[t$95$0, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + 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], If[LessEqual[t$95$0, 2.0], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(t$95$1 * im), $MachinePrecision] * im + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[(t$95$1 * 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(0.041666666666666664, im \cdot im, -0.5\right)\\
\mathbf{if}\;t\_0 \leq -0.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, 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{elif}\;t\_0 \leq 2:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(t\_1 \cdot im, im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \mathsf{fma}\left(t\_1, im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6431.0
Applied rewrites31.0%
Taylor expanded in im around 0
Applied rewrites3.2%
Taylor expanded in im around inf
Applied rewrites20.9%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 2Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.0
Applied rewrites99.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-*.f6464.8
Applied rewrites64.8%
Applied rewrites64.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
lower-fma.f6461.9
Applied rewrites61.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-*.f6472.9
Applied rewrites72.9%
Taylor expanded in re around inf
Applied rewrites72.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0)) (t_1 (* (exp re) (cos im))))
(if (<= t_1 -0.982)
(* (fma (* (fma 0.16666666666666666 re 0.5) re) re re) t_0)
(if (<= t_1 0.0)
(*
(pow (fma (fma (fma -0.16666666666666666 re 0.5) re -1.0) re 1.0) -1.0)
t_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 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0;
} else if (t_1 <= 0.0) {
tmp = pow(fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0), -1.0) * t_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 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0); elseif (t_1 <= 0.0) tmp = Float64((fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0) ^ -1.0) * t_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[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[Power[N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$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 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -0.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right), re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\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.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.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.f6431.2
Applied rewrites31.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f642.5
Applied rewrites2.5%
Applied rewrites2.5%
Taylor expanded in re around 0
Applied rewrites41.2%
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.f6492.0
Applied rewrites92.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-*.f6468.7
Applied rewrites68.7%
Final simplification59.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0)) (t_1 (* (exp re) (cos im))))
(if (<= t_1 -0.982)
(* (fma (* (fma 0.16666666666666666 re 0.5) re) re re) t_0)
(if (<= t_1 0.0)
(*
(pow (fma (fma (fma -0.16666666666666666 re 0.5) re -1.0) re 1.0) -1.0)
t_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 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0;
} else if (t_1 <= 0.0) {
tmp = pow(fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0), -1.0) * t_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 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0); elseif (t_1 <= 0.0) tmp = Float64((fma(fma(fma(-0.16666666666666666, re, 0.5), re, -1.0), re, 1.0) ^ -1.0) * t_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[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[Power[N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$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 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -0.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right), re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\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.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.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.f6431.2
Applied rewrites31.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f642.5
Applied rewrites2.5%
Applied rewrites2.5%
Taylor expanded in re around 0
Applied rewrites41.2%
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.f6489.5
Applied rewrites89.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-*.f6467.1
Applied rewrites67.1%
Final simplification58.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0)) (t_1 (* (exp re) (cos im))))
(if (<= t_1 -0.982)
(* (fma (* (fma 0.16666666666666666 re 0.5) re) re re) t_0)
(if (<= t_1 0.0)
(* (pow (fma (fma 0.5 re -1.0) re 1.0) -1.0) t_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 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0;
} else if (t_1 <= 0.0) {
tmp = pow(fma(fma(0.5, re, -1.0), re, 1.0), -1.0) * t_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 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0); elseif (t_1 <= 0.0) tmp = Float64((fma(fma(0.5, re, -1.0), re, 1.0) ^ -1.0) * t_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[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[Power[N[(N[(0.5 * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$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 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -0.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\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.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.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.f6431.2
Applied rewrites31.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f642.5
Applied rewrites2.5%
Applied rewrites2.5%
Taylor expanded in re around 0
Applied rewrites31.7%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6489.5
Applied rewrites89.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-*.f6467.1
Applied rewrites67.1%
Final simplification54.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0)) (t_1 (* (exp re) (cos im))))
(if (<= t_1 -0.982)
(* (fma (* (fma 0.16666666666666666 re 0.5) re) re re) t_0)
(if (<= t_1 2.0)
(* (pow (fma (fma 0.5 re -1.0) re 1.0) -1.0) t_0)
(*
(* (fma 0.5 re 1.0) re)
(fma (* 0.041666666666666664 (* im im)) (* im im) 1.0))))))
double code(double re, double im) {
double t_0 = fma((im * im), -0.5, 1.0);
double t_1 = exp(re) * cos(im);
double tmp;
if (t_1 <= -0.982) {
tmp = fma((fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0;
} else if (t_1 <= 2.0) {
tmp = pow(fma(fma(0.5, re, -1.0), re, 1.0), -1.0) * t_0;
} else {
tmp = (fma(0.5, re, 1.0) * re) * fma((0.041666666666666664 * (im * im)), (im * im), 1.0);
}
return tmp;
}
function code(re, im) t_0 = fma(Float64(im * im), -0.5, 1.0) t_1 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_1 <= -0.982) tmp = Float64(fma(Float64(fma(0.16666666666666666, re, 0.5) * re), re, re) * t_0); elseif (t_1 <= 2.0) tmp = Float64((fma(fma(0.5, re, -1.0), re, 1.0) ^ -1.0) * t_0); else tmp = Float64(Float64(fma(0.5, re, 1.0) * 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[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.982], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] * re + re), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 2.0], N[(N[Power[N[(N[(0.5 * re + -1.0), $MachinePrecision] * re + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * t$95$0), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * 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 := \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
t_1 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_1 \leq -0.982:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right) \cdot re, re, re\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 2:\\
\;\;\;\;{\left(\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, -1\right), re, 1\right)\right)}^{-1} \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.5, re, 1\right) \cdot re\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.98199999999999998Initial 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.f6480.3
Applied rewrites80.3%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6484.1
Applied rewrites84.1%
Taylor expanded in re around inf
Applied rewrites84.2%
if -0.98199999999999998 < (*.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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6465.0
Applied rewrites65.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6432.9
Applied rewrites32.9%
Applied rewrites32.9%
Taylor expanded in re around 0
Applied rewrites47.6%
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.f6461.9
Applied rewrites61.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-*.f6472.9
Applied rewrites72.9%
Taylor expanded in re around inf
Applied rewrites72.9%
Taylor expanded in im around inf
Applied rewrites72.9%
Final simplification54.5%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.55)
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(if (<= t_0 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 t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -0.55) {
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((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.55) 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(fma(0.041666666666666664, Float64(im * im), -0.5) * 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.55], 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[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + -0.5), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.55:\\
\;\;\;\;\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(\mathsf{fma}\left(0.041666666666666664, im \cdot im, -0.5\right) \cdot im, im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.55000000000000004Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6447.6
Applied rewrites47.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6452.4
Applied rewrites52.4%
if -0.55000000000000004 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6420.8
Applied rewrites20.8%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites23.5%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6473.2
Applied rewrites73.2%
Taylor expanded in im around 0
Applied rewrites53.4%
Applied rewrites53.4%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.55)
(* (+ 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.55) {
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.55) 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.55], 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.55:\\
\;\;\;\;\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.55000000000000004Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6447.6
Applied rewrites47.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6452.4
Applied rewrites52.4%
if -0.55000000000000004 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6420.8
Applied rewrites20.8%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites23.5%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6473.2
Applied rewrites73.2%
Taylor expanded in im around 0
Applied rewrites53.4%
Taylor expanded in im around inf
Applied rewrites53.3%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 1e-107)
(* (* im im) -0.5)
(if (<= (exp re) 1.001)
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(* (* (fma 0.5 re 1.0) re) (fma -0.5 (* im im) 1.0)))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 1e-107) {
tmp = (im * im) * -0.5;
} else if (exp(re) <= 1.001) {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
} else {
tmp = (fma(0.5, re, 1.0) * re) * fma(-0.5, (im * im), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (exp(re) <= 1e-107) tmp = Float64(Float64(im * im) * -0.5); elseif (exp(re) <= 1.001) tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); else tmp = Float64(Float64(fma(0.5, re, 1.0) * re) * fma(-0.5, Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 1e-107], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[N[Exp[re], $MachinePrecision], 1.001], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(-0.5 * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 10^{-107}:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;e^{re} \leq 1.001:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.5, re, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.5, im \cdot im, 1\right)\\
\end{array}
\end{array}
if (exp.f64 re) < 1e-107Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites27.9%
if 1e-107 < (exp.f64 re) < 1.0009999999999999Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.6
Applied rewrites99.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.1
Applied rewrites49.1%
if 1.0009999999999999 < (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.f6462.9
Applied rewrites62.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-*.f6444.8
Applied rewrites44.8%
Taylor expanded in re around inf
Applied rewrites44.8%
Taylor expanded in im around 0
Applied rewrites60.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (cos im)) 0.0) (* (* im im) -0.5) (fma (* im im) -0.5 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((im * im), -0.5, 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(Float64(im * im), -0.5, 1.0); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], 0.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}\;e^{re} \cdot \cos im \leq 0:\\
\;\;\;\;\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 (*.f64 (exp.f64 re) (cos.f64 im)) < -0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6428.0
Applied rewrites28.0%
Taylor expanded in im around 0
Applied rewrites12.7%
Taylor expanded in im around inf
Applied rewrites27.0%
if -0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6473.2
Applied rewrites73.2%
Taylor expanded in im around 0
Applied rewrites47.2%
(FPCore (re im) :precision binary64 (if (or (<= re -0.00024) (not (or (<= re 0.295) (not (<= re 4e+95))))) (* (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.00024) || !((re <= 0.295) || !(re <= 4e+95))) {
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.00024) || !((re <= 0.295) || !(re <= 4e+95))) 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.00024], N[Not[Or[LessEqual[re, 0.295], N[Not[LessEqual[re, 4e+95]], $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.00024 \lor \neg \left(re \leq 0.295 \lor \neg \left(re \leq 4 \cdot 10^{+95}\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 < -2.40000000000000006e-4 or 0.294999999999999984 < re < 4.00000000000000008e95Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6478.0
Applied rewrites78.0%
if -2.40000000000000006e-4 < re < 0.294999999999999984 or 4.00000000000000008e95 < 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.f6499.4
Applied rewrites99.4%
Final simplification92.2%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 1e-107)
(* (* im im) -0.5)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 1e-107) {
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 (exp(re) <= 1e-107) 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[N[Exp[re], $MachinePrecision], 1e-107], 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}\;e^{re} \leq 10^{-107}:\\
\;\;\;\;\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 (exp.f64 re) < 1e-107Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites27.9%
if 1e-107 < (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.f6491.8
Applied rewrites91.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6453.3
Applied rewrites53.3%
(FPCore (re im)
:precision binary64
(if (<= (exp re) 1e-107)
(* (* 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 (exp(re) <= 1e-107) {
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 (exp(re) <= 1e-107) 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[N[Exp[re], $MachinePrecision], 1e-107], 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}\;e^{re} \leq 10^{-107}:\\
\;\;\;\;\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 (exp.f64 re) < 1e-107Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites27.9%
if 1e-107 < (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.f6491.8
Applied rewrites91.8%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6453.3
Applied rewrites53.3%
Taylor expanded in re around inf
Applied rewrites52.9%
(FPCore (re im) :precision binary64 (if (or (<= re -0.00024) (not (or (<= re 0.295) (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.00024) || !((re <= 0.295) || !(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.00024) || !((re <= 0.295) || !(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.00024], N[Not[Or[LessEqual[re, 0.295], 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.00024 \lor \neg \left(re \leq 0.295 \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 < -2.40000000000000006e-4 or 0.294999999999999984 < re < 1.8999999999999999e154Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6476.4
Applied rewrites76.4%
if -2.40000000000000006e-4 < re < 0.294999999999999984 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.f6499.8
Applied rewrites99.8%
Final simplification91.3%
(FPCore (re im) :precision binary64 (if (<= (exp re) 1e-107) (* (* im im) -0.5) (* (fma (fma 0.5 re 1.0) re 1.0) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 1e-107) {
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 (exp(re) <= 1e-107) 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[N[Exp[re], $MachinePrecision], 1e-107], 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}\;e^{re} \leq 10^{-107}:\\
\;\;\;\;\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 (exp.f64 re) < 1e-107Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites27.9%
if 1e-107 < (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.f6488.5
Applied rewrites88.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6452.7
Applied rewrites52.7%
(FPCore (re im) :precision binary64 (if (or (<= re -9.2e-5) (not (<= re 0.0047))) (* (exp re) (fma (* im im) -0.5 1.0)) (* (+ 1.0 re) (cos im))))
double code(double re, double im) {
double tmp;
if ((re <= -9.2e-5) || !(re <= 0.0047)) {
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 <= -9.2e-5) || !(re <= 0.0047)) 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, -9.2e-5], N[Not[LessEqual[re, 0.0047]], $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 -9.2 \cdot 10^{-5} \lor \neg \left(re \leq 0.0047\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 < -9.20000000000000001e-5 or 0.00470000000000000018 < re Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6477.2
Applied rewrites77.2%
if -9.20000000000000001e-5 < re < 0.00470000000000000018Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6499.5
Applied rewrites99.5%
Final simplification88.5%
(FPCore (re im) :precision binary64 (if (<= (exp re) 1e-107) (* (* 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) <= 1e-107) {
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) <= 1e-107) 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], 1e-107], 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 10^{-107}:\\
\;\;\;\;\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) < 1e-107Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.6%
Taylor expanded in im around inf
Applied rewrites27.9%
if 1e-107 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6471.0
Applied rewrites71.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6445.0
Applied rewrites45.0%
(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.f6451.3
Applied rewrites51.3%
Taylor expanded in im around 0
Applied rewrites30.5%
Taylor expanded in im around inf
Applied rewrites14.0%
herbie shell --seed 2024321
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))