
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
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 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
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);
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Initial program 100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 -5e-50)
(/ (cos im) (fma (- (* (* -0.16666666666666666 re) re) 1.0) re 1.0))
(if (or (<= t_0 0.0) (not (<= t_0 0.9998)))
(exp re)
(/
(cos im)
(fma (- (* (fma -0.16666666666666666 re 0.5) re) 1.0) re 1.0)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= -5e-50) {
tmp = cos(im) / fma((((-0.16666666666666666 * re) * re) - 1.0), re, 1.0);
} else if ((t_0 <= 0.0) || !(t_0 <= 0.9998)) {
tmp = exp(re);
} else {
tmp = cos(im) / fma(((fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= -5e-50) tmp = Float64(cos(im) / fma(Float64(Float64(Float64(-0.16666666666666666 * re) * re) - 1.0), re, 1.0)); elseif ((t_0 <= 0.0) || !(t_0 <= 0.9998)) tmp = exp(re); else tmp = Float64(cos(im) / fma(Float64(Float64(fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -5e-50], N[(N[Cos[im], $MachinePrecision] / N[(N[(N[(N[(-0.16666666666666666 * re), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]], N[Exp[re], $MachinePrecision], N[(N[Cos[im], $MachinePrecision] / N[(N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50}:\\
\;\;\;\;\frac{\cos im}{\mathsf{fma}\left(\left(-0.16666666666666666 \cdot re\right) \cdot re - 1, re, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos im}{\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right) \cdot re - 1, re, 1\right)}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f6499.9
Applied rewrites99.9%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6496.4
Applied rewrites96.4%
Taylor expanded in re around inf
Applied rewrites96.4%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6497.4
Applied rewrites97.4%
Final simplification99.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 -5e-50)
(/ (cos im) (fma (- (* (* -0.16666666666666666 re) re) 1.0) re 1.0))
(if (or (<= t_0 0.0) (not (<= t_0 0.9998)))
(exp re)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= -5e-50) {
tmp = cos(im) / fma((((-0.16666666666666666 * re) * re) - 1.0), re, 1.0);
} else if ((t_0 <= 0.0) || !(t_0 <= 0.9998)) {
tmp = exp(re);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= -5e-50) tmp = Float64(cos(im) / fma(Float64(Float64(Float64(-0.16666666666666666 * re) * re) - 1.0), re, 1.0)); elseif ((t_0 <= 0.0) || !(t_0 <= 0.9998)) tmp = exp(re); 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_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -5e-50], N[(N[Cos[im], $MachinePrecision] / N[(N[(N[(N[(-0.16666666666666666 * re), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]], N[Exp[re], $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}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50}:\\
\;\;\;\;\frac{\cos im}{\mathsf{fma}\left(\left(-0.16666666666666666 \cdot re\right) \cdot re - 1, re, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right):\\
\;\;\;\;e^{re}\\
\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 (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f6499.9
Applied rewrites99.9%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6496.4
Applied rewrites96.4%
Taylor expanded in re around inf
Applied rewrites96.4%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6497.4
Applied rewrites97.4%
Final simplification99.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 -5e-50)
(/ (cos im) (- 1.0 re))
(if (or (<= t_0 0.0) (not (<= t_0 0.9998)))
(exp re)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= -5e-50) {
tmp = cos(im) / (1.0 - re);
} else if ((t_0 <= 0.0) || !(t_0 <= 0.9998)) {
tmp = exp(re);
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= -5e-50) tmp = Float64(cos(im) / Float64(1.0 - re)); elseif ((t_0 <= 0.0) || !(t_0 <= 0.9998)) tmp = exp(re); 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_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -5e-50], N[(N[Cos[im], $MachinePrecision] / N[(1.0 - re), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]], N[Exp[re], $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}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50}:\\
\;\;\;\;\frac{\cos im}{1 - re}\\
\mathbf{elif}\;t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right):\\
\;\;\;\;e^{re}\\
\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 (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f6499.9
Applied rewrites99.9%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f6496.4
Applied rewrites96.4%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6497.4
Applied rewrites97.4%
Final simplification99.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (<= t_0 -5e-50)
(/ (cos im) (- 1.0 re))
(if (or (<= t_0 0.0) (not (<= t_0 0.9998)))
(exp re)
(* (fma (fma 0.5 re 1.0) re 1.0) (cos im)))))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if (t_0 <= -5e-50) {
tmp = cos(im) / (1.0 - re);
} else if ((t_0 <= 0.0) || !(t_0 <= 0.9998)) {
tmp = exp(re);
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * cos(im);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif (t_0 <= -5e-50) tmp = Float64(cos(im) / Float64(1.0 - re)); elseif ((t_0 <= 0.0) || !(t_0 <= 0.9998)) tmp = exp(re); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * cos(im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, -5e-50], N[(N[Cos[im], $MachinePrecision] / N[(1.0 - re), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]], N[Exp[re], $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}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50}:\\
\;\;\;\;\frac{\cos im}{1 - re}\\
\mathbf{elif}\;t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \cos im\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f6499.9
Applied rewrites99.9%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f6496.4
Applied rewrites96.4%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6496.5
Applied rewrites96.5%
Final simplification99.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (exp re) (* (* im im) -0.5))
(if (or (<= t_0 -5e-50) (not (or (<= t_0 0.0) (not (<= t_0 0.9998)))))
(/ (cos im) (- 1.0 re))
(exp re)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = exp(re) * ((im * im) * -0.5);
} else if ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) {
tmp = cos(im) / (1.0 - re);
} else {
tmp = exp(re);
}
return tmp;
}
public static double code(double re, double im) {
double t_0 = Math.exp(re) * Math.cos(im);
double tmp;
if (t_0 <= -Double.POSITIVE_INFINITY) {
tmp = Math.exp(re) * ((im * im) * -0.5);
} else if ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) {
tmp = Math.cos(im) / (1.0 - re);
} else {
tmp = Math.exp(re);
}
return tmp;
}
def code(re, im): t_0 = math.exp(re) * math.cos(im) tmp = 0 if t_0 <= -math.inf: tmp = math.exp(re) * ((im * im) * -0.5) elif (t_0 <= -5e-50) or not ((t_0 <= 0.0) or not (t_0 <= 0.9998)): tmp = math.cos(im) / (1.0 - re) else: tmp = math.exp(re) return tmp
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(exp(re) * Float64(Float64(im * im) * -0.5)); elseif ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) tmp = Float64(cos(im) / Float64(1.0 - re)); else tmp = exp(re); end return tmp end
function tmp_2 = code(re, im) t_0 = exp(re) * cos(im); tmp = 0.0; if (t_0 <= -Inf) tmp = exp(re) * ((im * im) * -0.5); elseif ((t_0 <= -5e-50) || ~(((t_0 <= 0.0) || ~((t_0 <= 0.9998))))) tmp = cos(im) / (1.0 - re); else tmp = exp(re); end tmp_2 = tmp; end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[Exp[re], $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, -5e-50], N[Not[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]]], $MachinePrecision]], N[(N[Cos[im], $MachinePrecision] / N[(1.0 - re), $MachinePrecision]), $MachinePrecision], N[Exp[re], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;e^{re} \cdot \left(\left(im \cdot im\right) \cdot -0.5\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50} \lor \neg \left(t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right)\right):\\
\;\;\;\;\frac{\cos im}{1 - re}\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in im around inf
Applied rewrites100.0%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f6495.9
Applied rewrites95.9%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
Final simplification98.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))
(if (or (<= t_0 -5e-50) (not (or (<= t_0 0.0) (not (<= t_0 0.9998)))))
(/ (cos im) (- 1.0 re))
(exp re)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) {
tmp = cos(im) / (1.0 - re);
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) tmp = Float64(cos(im) / Float64(1.0 - re)); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, -5e-50], N[Not[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]]], $MachinePrecision]], N[(N[Cos[im], $MachinePrecision] / N[(1.0 - re), $MachinePrecision]), $MachinePrecision], N[Exp[re], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50} \lor \neg \left(t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right)\right):\\
\;\;\;\;\frac{\cos im}{1 - re}\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f6495.9
Applied rewrites95.9%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
Final simplification98.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))
(if (or (<= t_0 -0.05) (not (or (<= t_0 0.0) (not (<= t_0 0.9998)))))
(* (+ 1.0 re) (cos im))
(exp re)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if ((t_0 <= -0.05) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) {
tmp = (1.0 + re) * cos(im);
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif ((t_0 <= -0.05) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) tmp = Float64(Float64(1.0 + re) * cos(im)); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, -0.05], N[Not[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]]], $MachinePrecision]], N[(N[(1.0 + re), $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision], N[Exp[re], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq -0.05 \lor \neg \left(t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right)\right):\\
\;\;\;\;\left(1 + re\right) \cdot \cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -0.050000000000000003 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6497.2
Applied rewrites97.2%
if -0.050000000000000003 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f6499.4
Applied rewrites99.4%
Final simplification98.2%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))
(if (or (<= t_0 -5e-50) (not (or (<= t_0 0.0) (not (<= t_0 0.9998)))))
(cos im)
(exp re)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) {
tmp = cos(im);
} else {
tmp = exp(re);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif ((t_0 <= -5e-50) || !((t_0 <= 0.0) || !(t_0 <= 0.9998))) tmp = cos(im); else tmp = exp(re); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$0, -5e-50], N[Not[Or[LessEqual[t$95$0, 0.0], N[Not[LessEqual[t$95$0, 0.9998]], $MachinePrecision]]], $MachinePrecision]], N[Cos[im], $MachinePrecision], N[Exp[re], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq -5 \cdot 10^{-50} \lor \neg \left(t\_0 \leq 0 \lor \neg \left(t\_0 \leq 0.9998\right)\right):\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;e^{re}\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6494.8
Applied rewrites94.8%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0 or 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
exp-negN/A
remove-double-divN/A
lower-exp.f64100.0
Applied rewrites100.0%
Final simplification97.9%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0))
(t_1 (* (exp re) (cos im)))
(t_2 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)))
(if (<= t_1 (- INFINITY))
(* t_2 t_0)
(if (<= t_1 -5e-50)
(cos im)
(if (<= t_1 0.0)
(/ t_0 (fma (- (* (fma -0.16666666666666666 re 0.5) re) 1.0) re 1.0))
(if (<= t_1 0.9998)
(cos im)
(*
t_2
(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 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2 * t_0;
} else if (t_1 <= -5e-50) {
tmp = cos(im);
} else if (t_1 <= 0.0) {
tmp = t_0 / fma(((fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0);
} else if (t_1 <= 0.9998) {
tmp = cos(im);
} else {
tmp = t_2 * 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 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(t_2 * t_0); elseif (t_1 <= -5e-50) tmp = cos(im); elseif (t_1 <= 0.0) tmp = Float64(t_0 / fma(Float64(Float64(fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0)); elseif (t_1 <= 0.9998) tmp = cos(im); else tmp = Float64(t_2 * fma(Float64(Float64(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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(t$95$2 * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, -5e-50], N[Cos[im], $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(t$95$0 / N[(N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.9998], N[Cos[im], $MachinePrecision], N[(t$95$2 * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2 \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-50}:\\
\;\;\;\;\cos im\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right) \cdot re - 1, re, 1\right)}\\
\mathbf{elif}\;t\_1 \leq 0.9998:\\
\;\;\;\;\cos im\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5, 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
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < -4.99999999999999968e-50 or 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.99980000000000002Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6494.8
Applied rewrites94.8%
if -4.99999999999999968e-50 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6454.6
Applied rewrites54.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.1
Applied rewrites37.1%
if 0.99980000000000002 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6478.5
Applied rewrites78.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.4
Applied rewrites84.4%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6493.4
Applied rewrites93.4%
Final simplification81.7%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* im im) -0.5 1.0))
(t_1 (* (exp re) (cos im)))
(t_2 (fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)))
(if (<= t_1 (- INFINITY))
(* t_2 t_0)
(if (<= t_1 0.0)
(/ t_0 (fma (- (* (fma -0.16666666666666666 re 0.5) re) 1.0) re 1.0))
(*
t_2
(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 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2 * t_0;
} else if (t_1 <= 0.0) {
tmp = t_0 / fma(((fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0);
} else {
tmp = t_2 * 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 = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(t_2 * t_0); elseif (t_1 <= 0.0) tmp = Float64(t_0 / fma(Float64(Float64(fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0)); else tmp = Float64(t_2 * fma(Float64(Float64(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[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(t$95$2 * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(t$95$0 / N[(N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2 \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right) \cdot re - 1, re, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5, 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
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6468.2
Applied rewrites68.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6426.3
Applied rewrites26.3%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.4
Applied rewrites83.4%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.1
Applied rewrites64.1%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6470.8
Applied rewrites70.8%
(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 (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0) t_0)
(if (<= t_1 0.0)
(/ t_0 (fma (- (* (fma -0.16666666666666666 re 0.5) re) 1.0) re 1.0))
(*
(fma (fma 0.5 re 1.0) re 1.0)
(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(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * t_0;
} else if (t_1 <= 0.0) {
tmp = t_0 / fma(((fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0);
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * 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(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * t_0); elseif (t_1 <= 0.0) tmp = Float64(t_0 / fma(Float64(Float64(fma(-0.16666666666666666, re, 0.5) * re) - 1.0), re, 1.0)); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(Float64(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 + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(t$95$0 / N[(N[(N[(N[(-0.16666666666666666 * re + 0.5), $MachinePrecision] * re), $MachinePrecision] - 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot t\_0\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, re, 0.5\right) \cdot re - 1, re, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5, 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
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6447.9
Applied rewrites47.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.0
Applied rewrites86.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6490.7
Applied rewrites90.7%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
lift-*.f64N/A
lift-exp.f64N/A
remove-double-negN/A
rec-expN/A
associate-*l/N/A
*-commutativeN/A
lift-cos.f64N/A
sin-PI/2N/A
lower-/.f64N/A
lift-cos.f64N/A
sin-PI/2N/A
*-commutativeN/A
*-lft-identityN/A
lower-exp.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6468.2
Applied rewrites68.2%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6426.3
Applied rewrites26.3%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.4
Applied rewrites83.4%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.1
Applied rewrites64.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 -0.9995)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(*
(fma (fma 0.5 re 1.0) re 1.0)
(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.9995) {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(0.5, re, 1.0), re, 1.0) * 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.9995) tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(0.5, re, 1.0), re, 1.0) * fma(Float64(Float64(0.041666666666666664 * Float64(im * im)) - 0.5), Float64(im * im), 1.0)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.9995], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(N[(0.5 * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.5), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{re} \cdot \cos im\\
\mathbf{if}\;t\_0 \leq -0.9995:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(0.5, re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5, im \cdot im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -0.99950000000000006Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6450.4
Applied rewrites50.4%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6482.2
Applied rewrites82.2%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.7
Applied rewrites86.7%
if -0.99950000000000006 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6432.1
Applied rewrites32.1%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites19.2%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.4
Applied rewrites83.4%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.1
Applied rewrites64.1%
Final simplification52.1%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* (exp re) (cos im))))
(if (<= t_0 (- INFINITY))
(* (+ 1.0 re) (fma (* im im) -0.5 1.0))
(if (<= t_0 0.0)
(* (* im im) -0.5)
(fma (* (- (* 0.041666666666666664 (* im im)) 0.5) im) im 1.0)))))
double code(double re, double im) {
double t_0 = exp(re) * cos(im);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
} else if (t_0 <= 0.0) {
tmp = (im * im) * -0.5;
} else {
tmp = fma((((0.041666666666666664 * (im * im)) - 0.5) * im), im, 1.0);
}
return tmp;
}
function code(re, im) t_0 = Float64(exp(re) * cos(im)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); elseif (t_0 <= 0.0) tmp = Float64(Float64(im * im) * -0.5); else tmp = fma(Float64(Float64(Float64(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, (-Infinity)], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $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 -\infty:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5\right) \cdot im, im, 1\right)\\
\end{array}
\end{array}
if (*.f64 (exp.f64 re) (cos.f64 im)) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
lower-+.f645.4
Applied rewrites5.4%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6467.4
Applied rewrites67.4%
if -inf.0 < (*.f64 (exp.f64 re) (cos.f64 im)) < 0.0Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6432.9
Applied rewrites32.9%
Taylor expanded in im around 0
Applied rewrites2.9%
Taylor expanded in im around inf
Applied rewrites19.1%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6462.2
Applied rewrites62.2%
Taylor expanded in im around 0
Applied rewrites44.5%
Applied rewrites44.5%
Final simplification38.3%
(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.f6427.0
Applied rewrites27.0%
Taylor expanded in im around 0
Applied rewrites9.2%
Taylor expanded in im around inf
Applied rewrites22.1%
if 0.0 < (*.f64 (exp.f64 re) (cos.f64 im)) Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6462.2
Applied rewrites62.2%
Taylor expanded in im around 0
Applied rewrites37.6%
Final simplification31.6%
(FPCore (re im) :precision binary64 (if (<= (exp re) 2e-49) (* (* 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) <= 2e-49) {
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) <= 2e-49) 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], 2e-49], 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 2 \cdot 10^{-49}:\\
\;\;\;\;\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) < 1.99999999999999987e-49Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites26.0%
if 1.99999999999999987e-49 < (exp.f64 re) Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6481.9
Applied rewrites81.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6451.5
Applied rewrites51.5%
Final simplification46.0%
(FPCore (re im)
:precision binary64
(if (<= re -1.6)
(* (* im im) -0.5)
(*
(fma (fma (fma 0.16666666666666666 re 0.5) re 1.0) re 1.0)
(fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -1.6) {
tmp = (im * im) * -0.5;
} else {
tmp = fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.6) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(fma(fma(fma(0.16666666666666666, re, 0.5), re, 1.0), re, 1.0) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.6], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * re + 0.5), $MachinePrecision] * re + 1.0), $MachinePrecision] * re + 1.0), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1.6:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, re, 0.5\right), re, 1\right), re, 1\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -1.6000000000000001Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites26.0%
if -1.6000000000000001 < re Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6481.9
Applied rewrites81.9%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6451.5
Applied rewrites51.5%
Taylor expanded in re around 0
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6456.3
Applied rewrites56.3%
Final simplification49.8%
(FPCore (re im)
:precision binary64
(if (<= re -110.0)
(* (* im im) -0.5)
(if (<= re 60000000.0)
(fma (* (- (* 0.041666666666666664 (* im im)) 0.5) im) im 1.0)
(* (* (* re re) 0.5) (fma (* im im) -0.5 1.0)))))
double code(double re, double im) {
double tmp;
if (re <= -110.0) {
tmp = (im * im) * -0.5;
} else if (re <= 60000000.0) {
tmp = fma((((0.041666666666666664 * (im * im)) - 0.5) * im), im, 1.0);
} else {
tmp = ((re * re) * 0.5) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -110.0) tmp = Float64(Float64(im * im) * -0.5); elseif (re <= 60000000.0) tmp = fma(Float64(Float64(Float64(0.041666666666666664 * Float64(im * im)) - 0.5) * im), im, 1.0); else tmp = Float64(Float64(Float64(re * re) * 0.5) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -110.0], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], If[LessEqual[re, 60000000.0], N[(N[(N[(N[(0.041666666666666664 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.5), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision], N[(N[(N[(re * re), $MachinePrecision] * 0.5), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -110:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{elif}\;re \leq 60000000:\\
\;\;\;\;\mathsf{fma}\left(\left(0.041666666666666664 \cdot \left(im \cdot im\right) - 0.5\right) \cdot im, im, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(re \cdot re\right) \cdot 0.5\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -110Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites26.0%
if -110 < re < 6e7Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f6495.9
Applied rewrites95.9%
Taylor expanded in im around 0
Applied rewrites48.0%
Applied rewrites48.0%
if 6e7 < re Initial program 100.0%
Taylor expanded in re around 0
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6455.4
Applied rewrites55.4%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6456.9
Applied rewrites56.9%
Taylor expanded in re around inf
Applied rewrites56.9%
Final simplification45.9%
(FPCore (re im) :precision binary64 (if (<= re -1.4) (* (* im im) -0.5) (* (+ 1.0 re) (fma (* im im) -0.5 1.0))))
double code(double re, double im) {
double tmp;
if (re <= -1.4) {
tmp = (im * im) * -0.5;
} else {
tmp = (1.0 + re) * fma((im * im), -0.5, 1.0);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= -1.4) tmp = Float64(Float64(im * im) * -0.5); else tmp = Float64(Float64(1.0 + re) * fma(Float64(im * im), -0.5, 1.0)); end return tmp end
code[re_, im_] := If[LessEqual[re, -1.4], N[(N[(im * im), $MachinePrecision] * -0.5), $MachinePrecision], N[(N[(1.0 + re), $MachinePrecision] * N[(N[(im * im), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1.4:\\
\;\;\;\;\left(im \cdot im\right) \cdot -0.5\\
\mathbf{else}:\\
\;\;\;\;\left(1 + re\right) \cdot \mathsf{fma}\left(im \cdot im, -0.5, 1\right)\\
\end{array}
\end{array}
if re < -1.3999999999999999Initial program 100.0%
Taylor expanded in re around 0
lower-cos.f643.2
Applied rewrites3.2%
Taylor expanded in im around 0
Applied rewrites2.5%
Taylor expanded in im around inf
Applied rewrites26.0%
if -1.3999999999999999 < re Initial program 100.0%
Taylor expanded in re around 0
lower-+.f6462.7
Applied rewrites62.7%
Taylor expanded in im around 0
+-commutativeN/A
/-rgt-identityN/A
/-rgt-identityN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6437.6
Applied rewrites37.6%
Final simplification35.1%
(FPCore (re im) :precision binary64 (* (* im im) -0.5))
double code(double re, double im) {
return (im * im) * -0.5;
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
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.f6448.5
Applied rewrites48.5%
Taylor expanded in im around 0
Applied rewrites26.5%
Taylor expanded in im around inf
Applied rewrites9.6%
Final simplification9.6%
herbie shell --seed 2024359
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))