
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (+ (exp (- im)) (exp im))))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp(-im) + exp(im));
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * cos(re)) * (exp(-im) + exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.cos(re)) * (Math.exp(-im) + Math.exp(im));
}
def code(re, im): return (0.5 * math.cos(re)) * (math.exp(-im) + math.exp(im))
function code(re, im) return Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(-im)) + exp(im))) end
function tmp = code(re, im) tmp = (0.5 * cos(re)) * (exp(-im) + exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \cos re\right) \cdot \left(e^{-im} + e^{im}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (+ (exp (- im)) (exp im))))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp(-im) + exp(im));
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * cos(re)) * (exp(-im) + exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.cos(re)) * (Math.exp(-im) + Math.exp(im));
}
def code(re, im): return (0.5 * math.cos(re)) * (math.exp(-im) + math.exp(im))
function code(re, im) return Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(-im)) + exp(im))) end
function tmp = code(re, im) tmp = (0.5 * cos(re)) * (exp(-im) + exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \cos re\right) \cdot \left(e^{-im} + e^{im}\right)
\end{array}
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (+ (exp (- im)) (exp im))))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp(-im) + exp(im));
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * cos(re)) * (exp(-im) + exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.cos(re)) * (Math.exp(-im) + Math.exp(im));
}
def code(re, im): return (0.5 * math.cos(re)) * (math.exp(-im) + math.exp(im))
function code(re, im) return Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(-im)) + exp(im))) end
function tmp = code(re, im) tmp = (0.5 * cos(re)) * (exp(-im) + exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \cos re\right) \cdot \left(e^{-im} + e^{im}\right)
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (re im) :precision binary64 (if (or (<= im 1120.0) (not (<= im 1.32e+154))) (* (cos re) (+ (* 0.5 (* im im)) 1.0)) (* 0.5 (+ (exp (- im)) (exp im)))))
double code(double re, double im) {
double tmp;
if ((im <= 1120.0) || !(im <= 1.32e+154)) {
tmp = cos(re) * ((0.5 * (im * im)) + 1.0);
} else {
tmp = 0.5 * (exp(-im) + exp(im));
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((im <= 1120.0d0) .or. (.not. (im <= 1.32d+154))) then
tmp = cos(re) * ((0.5d0 * (im * im)) + 1.0d0)
else
tmp = 0.5d0 * (exp(-im) + exp(im))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((im <= 1120.0) || !(im <= 1.32e+154)) {
tmp = Math.cos(re) * ((0.5 * (im * im)) + 1.0);
} else {
tmp = 0.5 * (Math.exp(-im) + Math.exp(im));
}
return tmp;
}
def code(re, im): tmp = 0 if (im <= 1120.0) or not (im <= 1.32e+154): tmp = math.cos(re) * ((0.5 * (im * im)) + 1.0) else: tmp = 0.5 * (math.exp(-im) + math.exp(im)) return tmp
function code(re, im) tmp = 0.0 if ((im <= 1120.0) || !(im <= 1.32e+154)) tmp = Float64(cos(re) * Float64(Float64(0.5 * Float64(im * im)) + 1.0)); else tmp = Float64(0.5 * Float64(exp(Float64(-im)) + exp(im))); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((im <= 1120.0) || ~((im <= 1.32e+154))) tmp = cos(re) * ((0.5 * (im * im)) + 1.0); else tmp = 0.5 * (exp(-im) + exp(im)); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[im, 1120.0], N[Not[LessEqual[im, 1.32e+154]], $MachinePrecision]], N[(N[Cos[re], $MachinePrecision] * N[(N[(0.5 * N[(im * im), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(0.5 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 1120 \lor \neg \left(im \leq 1.32 \cdot 10^{+154}\right):\\
\;\;\;\;\cos re \cdot \left(0.5 \cdot \left(im \cdot im\right) + 1\right)\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(e^{-im} + e^{im}\right)\\
\end{array}
\end{array}
if im < 1120 or 1.31999999999999998e154 < im Initial program 100.0%
Taylor expanded in im around 0 85.0%
Simplified85.0%
unpow251.8%
Applied egg-rr85.0%
if 1120 < im < 1.31999999999999998e154Initial program 100.0%
Taylor expanded in re around 0 82.4%
Final simplification84.6%
(FPCore (re im)
:precision binary64
(if (<= im 3.5e+18)
(cos re)
(if (<= im 3.9e+152)
(+ 0.25 (* 0.25 (pow re 2.0)))
(+ (* 0.5 (* im im)) 1.0))))
double code(double re, double im) {
double tmp;
if (im <= 3.5e+18) {
tmp = cos(re);
} else if (im <= 3.9e+152) {
tmp = 0.25 + (0.25 * pow(re, 2.0));
} else {
tmp = (0.5 * (im * im)) + 1.0;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (im <= 3.5d+18) then
tmp = cos(re)
else if (im <= 3.9d+152) then
tmp = 0.25d0 + (0.25d0 * (re ** 2.0d0))
else
tmp = (0.5d0 * (im * im)) + 1.0d0
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (im <= 3.5e+18) {
tmp = Math.cos(re);
} else if (im <= 3.9e+152) {
tmp = 0.25 + (0.25 * Math.pow(re, 2.0));
} else {
tmp = (0.5 * (im * im)) + 1.0;
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 3.5e+18: tmp = math.cos(re) elif im <= 3.9e+152: tmp = 0.25 + (0.25 * math.pow(re, 2.0)) else: tmp = (0.5 * (im * im)) + 1.0 return tmp
function code(re, im) tmp = 0.0 if (im <= 3.5e+18) tmp = cos(re); elseif (im <= 3.9e+152) tmp = Float64(0.25 + Float64(0.25 * (re ^ 2.0))); else tmp = Float64(Float64(0.5 * Float64(im * im)) + 1.0); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (im <= 3.5e+18) tmp = cos(re); elseif (im <= 3.9e+152) tmp = 0.25 + (0.25 * (re ^ 2.0)); else tmp = (0.5 * (im * im)) + 1.0; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 3.5e+18], N[Cos[re], $MachinePrecision], If[LessEqual[im, 3.9e+152], N[(0.25 + N[(0.25 * N[Power[re, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * N[(im * im), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 3.5 \cdot 10^{+18}:\\
\;\;\;\;\cos re\\
\mathbf{elif}\;im \leq 3.9 \cdot 10^{+152}:\\
\;\;\;\;0.25 + 0.25 \cdot {re}^{2}\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(im \cdot im\right) + 1\\
\end{array}
\end{array}
if im < 3.5e18Initial program 100.0%
Taylor expanded in im around 0 64.7%
if 3.5e18 < im < 3.90000000000000011e152Initial program 100.0%
Applied egg-rr2.7%
Taylor expanded in re around 0 20.9%
*-commutative20.9%
Simplified20.9%
if 3.90000000000000011e152 < im Initial program 100.0%
Taylor expanded in im around 0 100.0%
Simplified100.0%
Taylor expanded in re around 0 75.9%
unpow275.9%
Applied egg-rr75.9%
Final simplification60.5%
(FPCore (re im) :precision binary64 (* (cos re) (+ (* 0.5 (* im im)) 1.0)))
double code(double re, double im) {
return cos(re) * ((0.5 * (im * im)) + 1.0);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = cos(re) * ((0.5d0 * (im * im)) + 1.0d0)
end function
public static double code(double re, double im) {
return Math.cos(re) * ((0.5 * (im * im)) + 1.0);
}
def code(re, im): return math.cos(re) * ((0.5 * (im * im)) + 1.0)
function code(re, im) return Float64(cos(re) * Float64(Float64(0.5 * Float64(im * im)) + 1.0)) end
function tmp = code(re, im) tmp = cos(re) * ((0.5 * (im * im)) + 1.0); end
code[re_, im_] := N[(N[Cos[re], $MachinePrecision] * N[(N[(0.5 * N[(im * im), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos re \cdot \left(0.5 \cdot \left(im \cdot im\right) + 1\right)
\end{array}
Initial program 100.0%
Taylor expanded in im around 0 74.4%
Simplified74.4%
unpow245.5%
Applied egg-rr74.4%
Final simplification74.4%
(FPCore (re im) :precision binary64 (if (<= im 1.05e+37) (cos re) (+ (* 0.5 (* im im)) 1.0)))
double code(double re, double im) {
double tmp;
if (im <= 1.05e+37) {
tmp = cos(re);
} else {
tmp = (0.5 * (im * im)) + 1.0;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (im <= 1.05d+37) then
tmp = cos(re)
else
tmp = (0.5d0 * (im * im)) + 1.0d0
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (im <= 1.05e+37) {
tmp = Math.cos(re);
} else {
tmp = (0.5 * (im * im)) + 1.0;
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 1.05e+37: tmp = math.cos(re) else: tmp = (0.5 * (im * im)) + 1.0 return tmp
function code(re, im) tmp = 0.0 if (im <= 1.05e+37) tmp = cos(re); else tmp = Float64(Float64(0.5 * Float64(im * im)) + 1.0); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (im <= 1.05e+37) tmp = cos(re); else tmp = (0.5 * (im * im)) + 1.0; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 1.05e+37], N[Cos[re], $MachinePrecision], N[(N[(0.5 * N[(im * im), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 1.05 \cdot 10^{+37}:\\
\;\;\;\;\cos re\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(im \cdot im\right) + 1\\
\end{array}
\end{array}
if im < 1.0500000000000001e37Initial program 100.0%
Taylor expanded in im around 0 64.1%
if 1.0500000000000001e37 < im Initial program 100.0%
Taylor expanded in im around 0 51.8%
Simplified51.8%
Taylor expanded in re around 0 39.6%
unpow239.6%
Applied egg-rr39.6%
Final simplification58.5%
(FPCore (re im) :precision binary64 (+ (* 0.5 (* im im)) 1.0))
double code(double re, double im) {
return (0.5 * (im * im)) + 1.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * (im * im)) + 1.0d0
end function
public static double code(double re, double im) {
return (0.5 * (im * im)) + 1.0;
}
def code(re, im): return (0.5 * (im * im)) + 1.0
function code(re, im) return Float64(Float64(0.5 * Float64(im * im)) + 1.0) end
function tmp = code(re, im) tmp = (0.5 * (im * im)) + 1.0; end
code[re_, im_] := N[(N[(0.5 * N[(im * im), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \left(im \cdot im\right) + 1
\end{array}
Initial program 100.0%
Taylor expanded in im around 0 74.4%
Simplified74.4%
Taylor expanded in re around 0 45.5%
unpow245.5%
Applied egg-rr45.5%
Final simplification45.5%
(FPCore (re im) :precision binary64 0.25)
double code(double re, double im) {
return 0.25;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = 0.25d0
end function
public static double code(double re, double im) {
return 0.25;
}
def code(re, im): return 0.25
function code(re, im) return 0.25 end
function tmp = code(re, im) tmp = 0.25; end
code[re_, im_] := 0.25
\begin{array}{l}
\\
0.25
\end{array}
Initial program 100.0%
Applied egg-rr8.1%
Taylor expanded in re around 0 8.2%
Final simplification8.2%
herbie shell --seed 2023321
(FPCore (re im)
:name "math.cos on complex, real part"
:precision binary64
(* (* 0.5 (cos re)) (+ (exp (- im)) (exp im))))