
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
(FPCore (re im) :precision binary64 (* (exp re) (cos im)))
double code(double re, double im) {
return exp(re) * cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * cos(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.cos(im);
}
def code(re, im): return math.exp(re) * math.cos(im)
function code(re, im) return Float64(exp(re) * cos(im)) end
function tmp = code(re, im) tmp = exp(re) * cos(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Cos[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \cos im
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (re im) :precision binary64 (if (or (<= re -5.5e-10) (not (<= re 2150000.0))) (exp re) (/ (- (cos im)) (+ re -1.0))))
double code(double re, double im) {
double tmp;
if ((re <= -5.5e-10) || !(re <= 2150000.0)) {
tmp = exp(re);
} else {
tmp = -cos(im) / (re + -1.0);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((re <= (-5.5d-10)) .or. (.not. (re <= 2150000.0d0))) then
tmp = exp(re)
else
tmp = -cos(im) / (re + (-1.0d0))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((re <= -5.5e-10) || !(re <= 2150000.0)) {
tmp = Math.exp(re);
} else {
tmp = -Math.cos(im) / (re + -1.0);
}
return tmp;
}
def code(re, im): tmp = 0 if (re <= -5.5e-10) or not (re <= 2150000.0): tmp = math.exp(re) else: tmp = -math.cos(im) / (re + -1.0) return tmp
function code(re, im) tmp = 0.0 if ((re <= -5.5e-10) || !(re <= 2150000.0)) tmp = exp(re); else tmp = Float64(Float64(-cos(im)) / Float64(re + -1.0)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((re <= -5.5e-10) || ~((re <= 2150000.0))) tmp = exp(re); else tmp = -cos(im) / (re + -1.0); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[re, -5.5e-10], N[Not[LessEqual[re, 2150000.0]], $MachinePrecision]], N[Exp[re], $MachinePrecision], N[((-N[Cos[im], $MachinePrecision]) / N[(re + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -5.5 \cdot 10^{-10} \lor \neg \left(re \leq 2150000\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\cos im}{re + -1}\\
\end{array}
\end{array}
if re < -5.4999999999999996e-10 or 2.15e6 < re Initial program 100.0%
Taylor expanded in im around 0 81.0%
if -5.4999999999999996e-10 < re < 2.15e6Initial program 100.0%
Taylor expanded in re around 0 99.0%
distribute-rgt1-in99.0%
Simplified99.0%
flip-+99.1%
associate-*l/99.1%
metadata-eval99.1%
fma-neg99.1%
metadata-eval99.1%
sub-neg99.1%
metadata-eval99.1%
Applied egg-rr99.1%
Taylor expanded in re around 0 99.0%
neg-mul-199.0%
Simplified99.0%
Final simplification89.9%
(FPCore (re im) :precision binary64 (if (or (<= re -4.8e-6) (not (<= re 2150000.0))) (exp re) (* (cos im) (+ re 1.0))))
double code(double re, double im) {
double tmp;
if ((re <= -4.8e-6) || !(re <= 2150000.0)) {
tmp = exp(re);
} else {
tmp = cos(im) * (re + 1.0);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((re <= (-4.8d-6)) .or. (.not. (re <= 2150000.0d0))) then
tmp = exp(re)
else
tmp = cos(im) * (re + 1.0d0)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((re <= -4.8e-6) || !(re <= 2150000.0)) {
tmp = Math.exp(re);
} else {
tmp = Math.cos(im) * (re + 1.0);
}
return tmp;
}
def code(re, im): tmp = 0 if (re <= -4.8e-6) or not (re <= 2150000.0): tmp = math.exp(re) else: tmp = math.cos(im) * (re + 1.0) return tmp
function code(re, im) tmp = 0.0 if ((re <= -4.8e-6) || !(re <= 2150000.0)) tmp = exp(re); else tmp = Float64(cos(im) * Float64(re + 1.0)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((re <= -4.8e-6) || ~((re <= 2150000.0))) tmp = exp(re); else tmp = cos(im) * (re + 1.0); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[re, -4.8e-6], N[Not[LessEqual[re, 2150000.0]], $MachinePrecision]], N[Exp[re], $MachinePrecision], N[(N[Cos[im], $MachinePrecision] * N[(re + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -4.8 \cdot 10^{-6} \lor \neg \left(re \leq 2150000\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\cos im \cdot \left(re + 1\right)\\
\end{array}
\end{array}
if re < -4.7999999999999998e-6 or 2.15e6 < re Initial program 100.0%
Taylor expanded in im around 0 80.8%
if -4.7999999999999998e-6 < re < 2.15e6Initial program 100.0%
Taylor expanded in re around 0 99.1%
distribute-rgt1-in99.0%
Simplified99.0%
Final simplification89.9%
(FPCore (re im) :precision binary64 (if (or (<= re -2.7e-6) (not (<= re 2150000.0))) (exp re) (+ re (cos im))))
double code(double re, double im) {
double tmp;
if ((re <= -2.7e-6) || !(re <= 2150000.0)) {
tmp = exp(re);
} else {
tmp = re + cos(im);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((re <= (-2.7d-6)) .or. (.not. (re <= 2150000.0d0))) then
tmp = exp(re)
else
tmp = re + cos(im)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((re <= -2.7e-6) || !(re <= 2150000.0)) {
tmp = Math.exp(re);
} else {
tmp = re + Math.cos(im);
}
return tmp;
}
def code(re, im): tmp = 0 if (re <= -2.7e-6) or not (re <= 2150000.0): tmp = math.exp(re) else: tmp = re + math.cos(im) return tmp
function code(re, im) tmp = 0.0 if ((re <= -2.7e-6) || !(re <= 2150000.0)) tmp = exp(re); else tmp = Float64(re + cos(im)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((re <= -2.7e-6) || ~((re <= 2150000.0))) tmp = exp(re); else tmp = re + cos(im); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[re, -2.7e-6], N[Not[LessEqual[re, 2150000.0]], $MachinePrecision]], N[Exp[re], $MachinePrecision], N[(re + N[Cos[im], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -2.7 \cdot 10^{-6} \lor \neg \left(re \leq 2150000\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;re + \cos im\\
\end{array}
\end{array}
if re < -2.69999999999999998e-6 or 2.15e6 < re Initial program 100.0%
Taylor expanded in im around 0 80.8%
if -2.69999999999999998e-6 < re < 2.15e6Initial program 100.0%
Taylor expanded in re around 0 99.1%
Taylor expanded in im around 0 98.5%
Final simplification89.6%
(FPCore (re im) :precision binary64 (if (or (<= re -1.14e-10) (not (<= re 2150000.0))) (exp re) (cos im)))
double code(double re, double im) {
double tmp;
if ((re <= -1.14e-10) || !(re <= 2150000.0)) {
tmp = exp(re);
} else {
tmp = cos(im);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((re <= (-1.14d-10)) .or. (.not. (re <= 2150000.0d0))) then
tmp = exp(re)
else
tmp = cos(im)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((re <= -1.14e-10) || !(re <= 2150000.0)) {
tmp = Math.exp(re);
} else {
tmp = Math.cos(im);
}
return tmp;
}
def code(re, im): tmp = 0 if (re <= -1.14e-10) or not (re <= 2150000.0): tmp = math.exp(re) else: tmp = math.cos(im) return tmp
function code(re, im) tmp = 0.0 if ((re <= -1.14e-10) || !(re <= 2150000.0)) tmp = exp(re); else tmp = cos(im); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((re <= -1.14e-10) || ~((re <= 2150000.0))) tmp = exp(re); else tmp = cos(im); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[re, -1.14e-10], N[Not[LessEqual[re, 2150000.0]], $MachinePrecision]], N[Exp[re], $MachinePrecision], N[Cos[im], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1.14 \cdot 10^{-10} \lor \neg \left(re \leq 2150000\right):\\
\;\;\;\;e^{re}\\
\mathbf{else}:\\
\;\;\;\;\cos im\\
\end{array}
\end{array}
if re < -1.1399999999999999e-10 or 2.15e6 < re Initial program 100.0%
Taylor expanded in im around 0 81.0%
if -1.1399999999999999e-10 < re < 2.15e6Initial program 100.0%
Taylor expanded in re around 0 98.2%
Final simplification89.4%
(FPCore (re im) :precision binary64 (cos im))
double code(double re, double im) {
return cos(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = cos(im)
end function
public static double code(double re, double im) {
return Math.cos(im);
}
def code(re, im): return math.cos(im)
function code(re, im) return cos(im) end
function tmp = code(re, im) tmp = cos(im); end
code[re_, im_] := N[Cos[im], $MachinePrecision]
\begin{array}{l}
\\
\cos im
\end{array}
Initial program 100.0%
Taylor expanded in re around 0 50.7%
Final simplification50.7%
(FPCore (re im) :precision binary64 (+ re 1.0))
double code(double re, double im) {
return re + 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = re + 1.0d0
end function
public static double code(double re, double im) {
return re + 1.0;
}
def code(re, im): return re + 1.0
function code(re, im) return Float64(re + 1.0) end
function tmp = code(re, im) tmp = re + 1.0; end
code[re_, im_] := N[(re + 1.0), $MachinePrecision]
\begin{array}{l}
\\
re + 1
\end{array}
Initial program 100.0%
Taylor expanded in re around 0 51.8%
distribute-rgt1-in51.8%
Simplified51.8%
Taylor expanded in im around 0 29.7%
Final simplification29.7%
(FPCore (re im) :precision binary64 1.0)
double code(double re, double im) {
return 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = 1.0d0
end function
public static double code(double re, double im) {
return 1.0;
}
def code(re, im): return 1.0
function code(re, im) return 1.0 end
function tmp = code(re, im) tmp = 1.0; end
code[re_, im_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
Taylor expanded in re around 0 51.8%
distribute-rgt1-in51.8%
Simplified51.8%
Taylor expanded in im around 0 35.2%
associate-+r+35.2%
+-commutative35.2%
associate-*r*35.2%
+-commutative35.2%
distribute-rgt1-in35.2%
+-commutative35.2%
Simplified35.2%
Taylor expanded in re around 0 31.9%
Taylor expanded in im around 0 29.1%
Final simplification29.1%
herbie shell --seed 2024044
(FPCore (re im)
:name "math.exp on complex, real part"
:precision binary64
(* (exp re) (cos im)))