
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
Initial program 100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))) (t_1 (* (exp re) im)))
(if (<= t_0 (- INFINITY))
(*
(* 0.16666666666666666 (* re (* re re)))
(fma
(fma
(* im im)
(fma (* im im) -0.0001984126984126984 0.008333333333333333)
-0.16666666666666666)
(* im (* im im))
im))
(if (<= t_0 -0.02)
(sin im)
(if (<= t_0 5e-103) t_1 (if (<= t_0 1.0) (sin im) t_1))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double t_1 = exp(re) * im;
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (0.16666666666666666 * (re * (re * re))) * fma(fma((im * im), fma((im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), (im * (im * im)), im);
} else if (t_0 <= -0.02) {
tmp = sin(im);
} else if (t_0 <= 5e-103) {
tmp = t_1;
} else if (t_0 <= 1.0) {
tmp = sin(im);
} else {
tmp = t_1;
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) t_1 = Float64(exp(re) * im) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(0.16666666666666666 * Float64(re * Float64(re * re))) * fma(fma(Float64(im * im), fma(Float64(im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), Float64(im * Float64(im * im)), im)); elseif (t_0 <= -0.02) tmp = sin(im); elseif (t_0 <= 5e-103) tmp = t_1; elseif (t_0 <= 1.0) tmp = sin(im); else tmp = t_1; end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * im), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(0.16666666666666666 * N[(re * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision] * N[(im * N[(im * im), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.02], N[Sin[im], $MachinePrecision], If[LessEqual[t$95$0, 5e-103], t$95$1, If[LessEqual[t$95$0, 1.0], N[Sin[im], $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
t_1 := e^{re} \cdot im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\left(0.16666666666666666 \cdot \left(re \cdot \left(re \cdot re\right)\right)\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), im \cdot \left(im \cdot im\right), im\right)\\
\mathbf{elif}\;t\_0 \leq -0.02:\\
\;\;\;\;\sin im\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{-103}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\sin im\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6483.7
Applied rewrites83.7%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6461.6
Applied rewrites61.6%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
*-rgt-identityN/A
lower-fma.f64N/A
Applied rewrites65.7%
Taylor expanded in re around inf
Applied rewrites65.7%
if -inf.0 < (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004 or 4.99999999999999966e-103 < (*.f64 (exp.f64 re) (sin.f64 im)) < 1Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f64100.0
Applied rewrites100.0%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 4.99999999999999966e-103 or 1 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6495.1
Applied rewrites95.1%
Final simplification91.5%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))))
(if (<= t_0 (- INFINITY))
(*
(* 0.16666666666666666 (* re (* re re)))
(fma
(fma
(* im im)
(fma (* im im) -0.0001984126984126984 0.008333333333333333)
-0.16666666666666666)
(* im (* im im))
im))
(if (<= t_0 -0.02)
(sin im)
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(if (<= t_0 1.0)
(sin im)
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (0.16666666666666666 * (re * (re * re))) * fma(fma((im * im), fma((im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), (im * (im * im)), im);
} else if (t_0 <= -0.02) {
tmp = sin(im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else if (t_0 <= 1.0) {
tmp = sin(im);
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(0.16666666666666666 * Float64(re * Float64(re * re))) * fma(fma(Float64(im * im), fma(Float64(im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), Float64(im * Float64(im * im)), im)); elseif (t_0 <= -0.02) tmp = sin(im); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); elseif (t_0 <= 1.0) tmp = sin(im); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(0.16666666666666666 * N[(re * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision] * N[(im * N[(im * im), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -0.02], N[Sin[im], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1.0], N[Sin[im], $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\left(0.16666666666666666 \cdot \left(re \cdot \left(re \cdot re\right)\right)\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), im \cdot \left(im \cdot im\right), im\right)\\
\mathbf{elif}\;t\_0 \leq -0.02:\\
\;\;\;\;\sin im\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\sin im\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6483.7
Applied rewrites83.7%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6461.6
Applied rewrites61.6%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
*-rgt-identityN/A
lower-fma.f64N/A
Applied rewrites65.7%
Taylor expanded in re around inf
Applied rewrites65.7%
if -inf.0 < (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004 or 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) < 1Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6499.6
Applied rewrites99.6%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 1 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6473.7
Applied rewrites73.7%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6463.5
Applied rewrites63.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6452.9
Applied rewrites52.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))
(t_1 (* (exp re) (sin im))))
(if (<= t_1 -0.02)
(*
(fma
(fma
(* im im)
(fma (* im im) -0.0001984126984126984 0.008333333333333333)
-0.16666666666666666)
(* im (* im im))
im)
t_0)
(if (<= t_1 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
t_0
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))
double code(double re, double im) {
double t_0 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
double t_1 = exp(re) * sin(im);
double tmp;
if (t_1 <= -0.02) {
tmp = fma(fma((im * im), fma((im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), (im * (im * im)), im) * t_0;
} else if (t_1 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = t_0 * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) t_1 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_1 <= -0.02) tmp = Float64(fma(fma(Float64(im * im), fma(Float64(im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), Float64(im * Float64(im * im)), im) * t_0); elseif (t_1 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(t_0 * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.02], N[(N[(N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision] * N[(im * N[(im * im), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
t_1 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_1 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), im \cdot \left(im \cdot im\right), im\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6440.4
Applied rewrites40.4%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
*-rgt-identityN/A
lower-fma.f64N/A
Applied rewrites43.0%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.8
Applied rewrites53.8%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6451.5
Applied rewrites51.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
Final simplification42.3%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))))
(if (<= t_0 -0.04)
(*
(* 0.16666666666666666 (* re (* re re)))
(fma
(fma
(* im im)
(fma (* im im) -0.0001984126984126984 0.008333333333333333)
-0.16666666666666666)
(* im (* im im))
im))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double tmp;
if (t_0 <= -0.04) {
tmp = (0.16666666666666666 * (re * (re * re))) * fma(fma((im * im), fma((im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), (im * (im * im)), im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_0 <= -0.04) tmp = Float64(Float64(0.16666666666666666 * Float64(re * Float64(re * re))) * fma(fma(Float64(im * im), fma(Float64(im * im), -0.0001984126984126984, 0.008333333333333333), -0.16666666666666666), Float64(im * Float64(im * im)), im)); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.04], N[(N[(0.16666666666666666 * N[(re * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision] * N[(im * N[(im * im), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_0 \leq -0.04:\\
\;\;\;\;\left(0.16666666666666666 \cdot \left(re \cdot \left(re \cdot re\right)\right)\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), im \cdot \left(im \cdot im\right), im\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0400000000000000008Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6440.4
Applied rewrites40.4%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
*-rgt-identityN/A
lower-fma.f64N/A
Applied rewrites43.0%
Taylor expanded in re around inf
Applied rewrites42.5%
if -0.0400000000000000008 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.8
Applied rewrites53.8%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6451.5
Applied rewrites51.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))))
(if (<= t_0 -0.02)
(*
(fma re (fma re 0.5 1.0) 1.0)
(*
im
(fma
(* im im)
(fma
(* im im)
(fma im (* im -0.0001984126984126984) 0.008333333333333333)
-0.16666666666666666)
1.0)))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double tmp;
if (t_0 <= -0.02) {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0) * (im * fma((im * im), fma((im * im), fma(im, (im * -0.0001984126984126984), 0.008333333333333333), -0.16666666666666666), 1.0));
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_0 <= -0.02) tmp = Float64(fma(re, fma(re, 0.5, 1.0), 1.0) * Float64(im * fma(Float64(im * im), fma(Float64(im * im), fma(im, Float64(im * -0.0001984126984126984), 0.008333333333333333), -0.16666666666666666), 1.0))); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.02], N[(N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(im * -0.0001984126984126984), $MachinePrecision] + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_0 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right) \cdot \left(im \cdot \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im, im \cdot -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), 1\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6458.7
Applied rewrites58.7%
Taylor expanded in im around 0
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6435.8
Applied rewrites35.8%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.8
Applied rewrites53.8%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6451.5
Applied rewrites51.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))))
(if (<= t_0 -0.02)
(*
(fma re (* re 0.5) re)
(fma
im
(*
(* im im)
(fma
(* im im)
(fma im (* im -0.0001984126984126984) 0.008333333333333333)
-0.16666666666666666))
im))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double tmp;
if (t_0 <= -0.02) {
tmp = fma(re, (re * 0.5), re) * fma(im, ((im * im) * fma((im * im), fma(im, (im * -0.0001984126984126984), 0.008333333333333333), -0.16666666666666666)), im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_0 <= -0.02) tmp = Float64(fma(re, Float64(re * 0.5), re) * fma(im, Float64(Float64(im * im) * fma(Float64(im * im), fma(im, Float64(im * -0.0001984126984126984), 0.008333333333333333), -0.16666666666666666)), im)); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.02], N[(N[(re * N[(re * 0.5), $MachinePrecision] + re), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(im * -0.0001984126984126984), $MachinePrecision] + 0.008333333333333333), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_0 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(re, re \cdot 0.5, re\right) \cdot \mathsf{fma}\left(im, \left(im \cdot im\right) \cdot \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im, im \cdot -0.0001984126984126984, 0.008333333333333333\right), -0.16666666666666666\right), im\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6458.7
Applied rewrites58.7%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6434.4
Applied rewrites34.4%
Taylor expanded in re around inf
Applied rewrites34.5%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
Applied rewrites35.9%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.8
Applied rewrites53.8%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6451.5
Applied rewrites51.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im)))
(t_1 (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))
(if (<= t_0 -0.02)
(* t_1 (fma im (* (* im im) -0.16666666666666666) im))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
t_1
(fma
(* im im)
(* im (fma (* im im) 0.008333333333333333 -0.16666666666666666))
im))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double t_1 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
double tmp;
if (t_0 <= -0.02) {
tmp = t_1 * fma(im, ((im * im) * -0.16666666666666666), im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = t_1 * fma((im * im), (im * fma((im * im), 0.008333333333333333, -0.16666666666666666)), im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) t_1 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) tmp = 0.0 if (t_0 <= -0.02) tmp = Float64(t_1 * fma(im, Float64(Float64(im * im) * -0.16666666666666666), im)); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(t_1 * fma(Float64(im * im), Float64(im * fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666)), im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, -0.02], N[(t$95$1 * N[(im * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[(N[(im * im), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
t_1 := \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{if}\;t\_0 \leq -0.02:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im, \left(im \cdot im\right) \cdot -0.16666666666666666, im\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im \cdot im, im \cdot \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), im\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6440.4
Applied rewrites40.4%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.8
Applied rewrites53.8%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6451.5
Applied rewrites51.5%
Taylor expanded in im around 0
distribute-rgt-inN/A
*-lft-identityN/A
+-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
Final simplification41.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im))))
(if (<= t_0 -0.02)
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma im (* (* im im) -0.16666666666666666) im))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(*
(fma re (fma re 0.5 1.0) 1.0)
(*
im
(fma
(* im im)
(fma (* im im) 0.008333333333333333 -0.16666666666666666)
1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double tmp;
if (t_0 <= -0.02) {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, ((im * im) * -0.16666666666666666), im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0) * (im * fma((im * im), fma((im * im), 0.008333333333333333, -0.16666666666666666), 1.0));
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) tmp = 0.0 if (t_0 <= -0.02) tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, Float64(Float64(im * im) * -0.16666666666666666), im)); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(fma(re, fma(re, 0.5, 1.0), 1.0) * Float64(im * fma(Float64(im * im), fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666), 1.0))); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.02], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
\mathbf{if}\;t\_0 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im, \left(im \cdot im\right) \cdot -0.16666666666666666, im\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right) \cdot \left(im \cdot \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), 1\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6440.4
Applied rewrites40.4%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6489.9
Applied rewrites89.9%
Taylor expanded in im around 0
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6449.2
Applied rewrites49.2%
Final simplification41.6%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (sin im)))
(t_1 (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)))
(if (<= t_0 -0.02)
(* t_1 (fma im (* (* im im) -0.16666666666666666) im))
(if (<= t_0 0.0)
(/ (- (* im im) (* re (* re (* im im)))) (- im (* re im)))
(* im t_1)))))
double code(double re, double im) {
double t_0 = exp(re) * sin(im);
double t_1 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
double tmp;
if (t_0 <= -0.02) {
tmp = t_1 * fma(im, ((im * im) * -0.16666666666666666), im);
} else if (t_0 <= 0.0) {
tmp = ((im * im) - (re * (re * (im * im)))) / (im - (re * im));
} else {
tmp = im * t_1;
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * sin(im)) t_1 = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) tmp = 0.0 if (t_0 <= -0.02) tmp = Float64(t_1 * fma(im, Float64(Float64(im * im) * -0.16666666666666666), im)); elseif (t_0 <= 0.0) tmp = Float64(Float64(Float64(im * im) - Float64(re * Float64(re * Float64(im * im)))) / Float64(im - Float64(re * im))); else tmp = Float64(im * t_1); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[t$95$0, -0.02], N[(t$95$1 * N[(im * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(N[(im * im), $MachinePrecision] - N[(re * N[(re * N[(im * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(im - N[(re * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im * t$95$1), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \sin im\\
t_1 := \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{if}\;t\_0 \leq -0.02:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im, \left(im \cdot im\right) \cdot -0.16666666666666666, im\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{im \cdot im - re \cdot \left(re \cdot \left(im \cdot im\right)\right)}{im - re \cdot im}\\
\mathbf{else}:\\
\;\;\;\;im \cdot t\_1\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < -0.0200000000000000004Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6454.6
Applied rewrites54.6%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6440.4
Applied rewrites40.4%
if -0.0200000000000000004 < (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
Applied rewrites33.3%
Applied rewrites36.1%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
Final simplification41.0%
(FPCore (re im)
:precision binary64
(if (<= (* (exp re) (sin im)) 0.0)
(*
(fma re (fma re 0.5 1.0) 1.0)
(fma im (* im (* im -0.16666666666666666)) im))
(* im (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(re, fma(re, 0.5, 1.0), 1.0) * fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = Float64(fma(re, fma(re, 0.5, 1.0), 1.0) * fma(im, Float64(im * Float64(im * -0.16666666666666666)), im)); else tmp = Float64(im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right) \cdot \mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6443.2
Applied rewrites43.2%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6433.6
Applied rewrites33.6%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.0) (* (fma re (* re 0.5) re) (fma im (* im (* im -0.16666666666666666)) im)) (* im (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(re, (re * 0.5), re) * fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = Float64(fma(re, Float64(re * 0.5), re) * fma(im, Float64(im * Float64(im * -0.16666666666666666)), im)); else tmp = Float64(im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(re * N[(re * 0.5), $MachinePrecision] + re), $MachinePrecision] * N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(re, re \cdot 0.5, re\right) \cdot \mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6443.2
Applied rewrites43.2%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6433.6
Applied rewrites33.6%
Taylor expanded in re around inf
Applied rewrites15.4%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.0) (* (fma im (* im (* im -0.16666666666666666)) im) (* 0.5 (* re re))) (* im (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im) * (0.5 * (re * re));
} else {
tmp = im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = Float64(fma(im, Float64(im * Float64(im * -0.16666666666666666)), im) * Float64(0.5 * Float64(re * re))); else tmp = Float64(im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision] * N[(0.5 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right) \cdot \left(0.5 \cdot \left(re \cdot re\right)\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6443.2
Applied rewrites43.2%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6433.6
Applied rewrites33.6%
Taylor expanded in re around inf
Applied rewrites15.3%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
Final simplification26.0%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.0) (* (fma im (* (* im im) -0.16666666666666666) im) (+ re 1.0)) (* im (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(im, ((im * im) * -0.16666666666666666), im) * (re + 1.0);
} else {
tmp = im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = Float64(fma(im, Float64(Float64(im * im) * -0.16666666666666666), im) * Float64(re + 1.0)); else tmp = Float64(im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(im * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] + im), $MachinePrecision] * N[(re + 1.0), $MachinePrecision]), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(im, \left(im \cdot im\right) \cdot -0.16666666666666666, im\right) \cdot \left(re + 1\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.1
Applied rewrites71.1%
Taylor expanded in re around 0
+-commutativeN/A
lower-+.f6428.7
Applied rewrites28.7%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
Final simplification34.9%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.0) (fma im (* im (* im -0.16666666666666666)) im) (* im (fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6435.6
Applied rewrites35.6%
Taylor expanded in im around 0
Applied rewrites28.4%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites47.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 4e-10) (fma im (* im (* im -0.16666666666666666)) im) (* im (fma (fma re 0.16666666666666666 0.5) (* re re) re))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 4e-10) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * fma(fma(re, 0.16666666666666666, 0.5), (re * re), re);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 4e-10) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * fma(fma(re, 0.16666666666666666, 0.5), Float64(re * re), re)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 4e-10], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] * N[(re * re), $MachinePrecision] + re), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 4 \cdot 10^{-10}:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(\mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), re \cdot re, re\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 4.00000000000000015e-10Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6445.3
Applied rewrites45.3%
Taylor expanded in im around 0
Applied rewrites39.2%
if 4.00000000000000015e-10 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6424.3
Applied rewrites24.3%
Taylor expanded in re around 0
Applied rewrites17.3%
Taylor expanded in re around inf
Applied rewrites17.6%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 4e-10) (fma im (* im (* im -0.16666666666666666)) im) (* im (* re (* re (fma re 0.16666666666666666 0.5))))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 4e-10) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * (re * (re * fma(re, 0.16666666666666666, 0.5)));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 4e-10) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * Float64(re * Float64(re * fma(re, 0.16666666666666666, 0.5)))); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 4e-10], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 4 \cdot 10^{-10}:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \left(re \cdot \left(re \cdot \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 4.00000000000000015e-10Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6445.3
Applied rewrites45.3%
Taylor expanded in im around 0
Applied rewrites39.2%
if 4.00000000000000015e-10 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6424.3
Applied rewrites24.3%
Taylor expanded in re around 0
Applied rewrites17.3%
Taylor expanded in re around inf
Applied rewrites17.2%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 4e-10) (fma im (* im (* im -0.16666666666666666)) im) (* im (* re (* re (* re 0.16666666666666666))))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 4e-10) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * (re * (re * (re * 0.16666666666666666)));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 4e-10) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * Float64(re * Float64(re * Float64(re * 0.16666666666666666)))); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 4e-10], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(re * N[(re * N[(re * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 4 \cdot 10^{-10}:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \left(re \cdot \left(re \cdot \left(re \cdot 0.16666666666666666\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 4.00000000000000015e-10Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6445.3
Applied rewrites45.3%
Taylor expanded in im around 0
Applied rewrites39.2%
if 4.00000000000000015e-10 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6424.3
Applied rewrites24.3%
Taylor expanded in re around 0
Applied rewrites17.3%
Taylor expanded in re around inf
Applied rewrites17.3%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.0) (fma im (* im (* im -0.16666666666666666)) im) (* im (fma re (fma re 0.5 1.0) 1.0))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.0) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * fma(re, fma(re, 0.5, 1.0), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.0) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * fma(re, fma(re, 0.5, 1.0), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.0], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(re * N[(re * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.5, 1\right), 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6435.6
Applied rewrites35.6%
Taylor expanded in im around 0
Applied rewrites28.4%
if 0.0 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6451.9
Applied rewrites51.9%
Taylor expanded in re around 0
Applied rewrites46.4%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 4e-10) (fma im (* im (* im -0.16666666666666666)) im) (* im (* 0.5 (* re re)))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 4e-10) {
tmp = fma(im, (im * (im * -0.16666666666666666)), im);
} else {
tmp = im * (0.5 * (re * re));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 4e-10) tmp = fma(im, Float64(im * Float64(im * -0.16666666666666666)), im); else tmp = Float64(im * Float64(0.5 * Float64(re * re))); end return tmp end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 4e-10], N[(im * N[(im * N[(im * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + im), $MachinePrecision], N[(im * N[(0.5 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 4 \cdot 10^{-10}:\\
\;\;\;\;\mathsf{fma}\left(im, im \cdot \left(im \cdot -0.16666666666666666\right), im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \left(0.5 \cdot \left(re \cdot re\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 4.00000000000000015e-10Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6445.3
Applied rewrites45.3%
Taylor expanded in im around 0
Applied rewrites39.2%
if 4.00000000000000015e-10 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6424.3
Applied rewrites24.3%
Taylor expanded in re around 0
Applied rewrites10.4%
Taylor expanded in re around inf
Applied rewrites15.6%
(FPCore (re im) :precision binary64 (if (<= (* (exp re) (sin im)) 0.96) (* im 1.0) (* im (* 0.5 (* re re)))))
double code(double re, double im) {
double tmp;
if ((exp(re) * sin(im)) <= 0.96) {
tmp = im * 1.0;
} else {
tmp = im * (0.5 * (re * re));
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((exp(re) * sin(im)) <= 0.96d0) then
tmp = im * 1.0d0
else
tmp = im * (0.5d0 * (re * re))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((Math.exp(re) * Math.sin(im)) <= 0.96) {
tmp = im * 1.0;
} else {
tmp = im * (0.5 * (re * re));
}
return tmp;
}
def code(re, im): tmp = 0 if (math.exp(re) * math.sin(im)) <= 0.96: tmp = im * 1.0 else: tmp = im * (0.5 * (re * re)) return tmp
function code(re, im) tmp = 0.0 if (Float64(exp(re) * sin(im)) <= 0.96) tmp = Float64(im * 1.0); else tmp = Float64(im * Float64(0.5 * Float64(re * re))); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((exp(re) * sin(im)) <= 0.96) tmp = im * 1.0; else tmp = im * (0.5 * (re * re)); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision], 0.96], N[(im * 1.0), $MachinePrecision], N[(im * N[(0.5 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \cdot \sin im \leq 0.96:\\
\;\;\;\;im \cdot 1\\
\mathbf{else}:\\
\;\;\;\;im \cdot \left(0.5 \cdot \left(re \cdot re\right)\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (sin.f64 im)) < 0.95999999999999996Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6469.6
Applied rewrites69.6%
Taylor expanded in re around 0
Applied rewrites29.4%
if 0.95999999999999996 < (*.f64 (exp.f64 re) (sin.f64 im)) Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6454.8
Applied rewrites54.8%
Taylor expanded in re around 0
Applied rewrites20.7%
Taylor expanded in re around inf
Applied rewrites33.6%
(FPCore (re im)
:precision binary64
(if (<= (sin im) 0.01)
(*
(fma re (fma re (fma re 0.16666666666666666 0.5) 1.0) 1.0)
(fma im (* (* im im) -0.16666666666666666) im))
(*
im
(fma
(* im im)
(fma (* im im) 0.008333333333333333 -0.16666666666666666)
1.0))))
double code(double re, double im) {
double tmp;
if (sin(im) <= 0.01) {
tmp = fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, ((im * im) * -0.16666666666666666), im);
} else {
tmp = im * fma((im * im), fma((im * im), 0.008333333333333333, -0.16666666666666666), 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (sin(im) <= 0.01) tmp = Float64(fma(re, fma(re, fma(re, 0.16666666666666666, 0.5), 1.0), 1.0) * fma(im, Float64(Float64(im * im) * -0.16666666666666666), im)); else tmp = Float64(im * fma(Float64(im * im), fma(Float64(im * im), 0.008333333333333333, -0.16666666666666666), 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[N[Sin[im], $MachinePrecision], 0.01], N[(N[(re * N[(re * N[(re * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] + im), $MachinePrecision]), $MachinePrecision], N[(im * N[(N[(im * im), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * 0.008333333333333333 + -0.16666666666666666), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\sin im \leq 0.01:\\
\;\;\;\;\mathsf{fma}\left(re, \mathsf{fma}\left(re, \mathsf{fma}\left(re, 0.16666666666666666, 0.5\right), 1\right), 1\right) \cdot \mathsf{fma}\left(im, \left(im \cdot im\right) \cdot -0.16666666666666666, im\right)\\
\mathbf{else}:\\
\;\;\;\;im \cdot \mathsf{fma}\left(im \cdot im, \mathsf{fma}\left(im \cdot im, 0.008333333333333333, -0.16666666666666666\right), 1\right)\\
\end{array}
\end{array}
if (sin.f64 im) < 0.0100000000000000002Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6478.9
Applied rewrites78.9%
Taylor expanded in re around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6449.6
Applied rewrites49.6%
if 0.0100000000000000002 < (sin.f64 im) Initial program 100.0%
Taylor expanded in re around 0
lower-sin.f6459.6
Applied rewrites59.6%
Taylor expanded in im around 0
Applied rewrites8.5%
Final simplification40.0%
(FPCore (re im) :precision binary64 (fma im re im))
double code(double re, double im) {
return fma(im, re, im);
}
function code(re, im) return fma(im, re, im) end
code[re_, im_] := N[(im * re + im), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(im, re, im\right)
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6468.3
Applied rewrites68.3%
Taylor expanded in re around 0
Applied rewrites27.7%
(FPCore (re im) :precision binary64 (* im 1.0))
double code(double re, double im) {
return im * 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = im * 1.0d0
end function
public static double code(double re, double im) {
return im * 1.0;
}
def code(re, im): return im * 1.0
function code(re, im) return Float64(im * 1.0) end
function tmp = code(re, im) tmp = im * 1.0; end
code[re_, im_] := N[(im * 1.0), $MachinePrecision]
\begin{array}{l}
\\
im \cdot 1
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6468.3
Applied rewrites68.3%
Taylor expanded in re around 0
Applied rewrites27.0%
(FPCore (re im) :precision binary64 (* re im))
double code(double re, double im) {
return re * im;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = re * im
end function
public static double code(double re, double im) {
return re * im;
}
def code(re, im): return re * im
function code(re, im) return Float64(re * im) end
function tmp = code(re, im) tmp = re * im; end
code[re_, im_] := N[(re * im), $MachinePrecision]
\begin{array}{l}
\\
re \cdot im
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
lower-*.f64N/A
lower-exp.f6468.3
Applied rewrites68.3%
Taylor expanded in re around 0
Applied rewrites27.7%
Taylor expanded in re around inf
Applied rewrites4.5%
Final simplification4.5%
herbie shell --seed 2024237
(FPCore (re im)
:name "math.exp on complex, imaginary part"
:precision binary64
(* (exp re) (sin im)))