
(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 11 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 (<= im 0.058)
(* (cos re) (+ (* 0.5 (pow im 2.0)) 1.0))
(if (<= im 2.6e+45)
(* 0.5 (+ (exp (- im)) (exp im)))
(* (cos re) (cbrt (* (pow im 6.0) 0.125))))))
double code(double re, double im) {
double tmp;
if (im <= 0.058) {
tmp = cos(re) * ((0.5 * pow(im, 2.0)) + 1.0);
} else if (im <= 2.6e+45) {
tmp = 0.5 * (exp(-im) + exp(im));
} else {
tmp = cos(re) * cbrt((pow(im, 6.0) * 0.125));
}
return tmp;
}
public static double code(double re, double im) {
double tmp;
if (im <= 0.058) {
tmp = Math.cos(re) * ((0.5 * Math.pow(im, 2.0)) + 1.0);
} else if (im <= 2.6e+45) {
tmp = 0.5 * (Math.exp(-im) + Math.exp(im));
} else {
tmp = Math.cos(re) * Math.cbrt((Math.pow(im, 6.0) * 0.125));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (im <= 0.058) tmp = Float64(cos(re) * Float64(Float64(0.5 * (im ^ 2.0)) + 1.0)); elseif (im <= 2.6e+45) tmp = Float64(0.5 * Float64(exp(Float64(-im)) + exp(im))); else tmp = Float64(cos(re) * cbrt(Float64((im ^ 6.0) * 0.125))); end return tmp end
code[re_, im_] := If[LessEqual[im, 0.058], N[(N[Cos[re], $MachinePrecision] * N[(N[(0.5 * N[Power[im, 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 2.6e+45], N[(0.5 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Cos[re], $MachinePrecision] * N[Power[N[(N[Power[im, 6.0], $MachinePrecision] * 0.125), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 0.058:\\
\;\;\;\;\cos re \cdot \left(0.5 \cdot {im}^{2} + 1\right)\\
\mathbf{elif}\;im \leq 2.6 \cdot 10^{+45}:\\
\;\;\;\;0.5 \cdot \left(e^{-im} + e^{im}\right)\\
\mathbf{else}:\\
\;\;\;\;\cos re \cdot \sqrt[3]{{im}^{6} \cdot 0.125}\\
\end{array}
\end{array}
if im < 0.0580000000000000029Initial program 100.0%
Taylor expanded in im around 0 81.4%
associate-*r*81.4%
distribute-rgt1-in81.4%
Simplified81.4%
if 0.0580000000000000029 < im < 2.60000000000000007e45Initial program 100.0%
Taylor expanded in re around 0 83.3%
if 2.60000000000000007e45 < im Initial program 100.0%
Taylor expanded in im around 0 57.6%
associate-*r*57.6%
distribute-rgt1-in57.6%
Simplified57.6%
Taylor expanded in im around inf 57.6%
associate-*r*57.6%
*-commutative57.6%
Simplified57.6%
add-cbrt-cube77.6%
pow1/377.6%
pow377.6%
*-commutative77.6%
unpow-prod-down77.6%
pow-pow77.6%
metadata-eval77.6%
metadata-eval77.6%
Applied egg-rr100.0%
unpow1/377.6%
Simplified100.0%
Final simplification85.1%
(FPCore (re im)
:precision binary64
(if (<= im 0.054)
(cos re)
(if (<= im 1.35e+154)
(* 0.5 (+ (exp (- im)) (exp im)))
(* (cos re) (* 0.5 (pow im 2.0))))))
double code(double re, double im) {
double tmp;
if (im <= 0.054) {
tmp = cos(re);
} else if (im <= 1.35e+154) {
tmp = 0.5 * (exp(-im) + exp(im));
} else {
tmp = cos(re) * (0.5 * pow(im, 2.0));
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (im <= 0.054d0) then
tmp = cos(re)
else if (im <= 1.35d+154) then
tmp = 0.5d0 * (exp(-im) + exp(im))
else
tmp = cos(re) * (0.5d0 * (im ** 2.0d0))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (im <= 0.054) {
tmp = Math.cos(re);
} else if (im <= 1.35e+154) {
tmp = 0.5 * (Math.exp(-im) + Math.exp(im));
} else {
tmp = Math.cos(re) * (0.5 * Math.pow(im, 2.0));
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 0.054: tmp = math.cos(re) elif im <= 1.35e+154: tmp = 0.5 * (math.exp(-im) + math.exp(im)) else: tmp = math.cos(re) * (0.5 * math.pow(im, 2.0)) return tmp
function code(re, im) tmp = 0.0 if (im <= 0.054) tmp = cos(re); elseif (im <= 1.35e+154) tmp = Float64(0.5 * Float64(exp(Float64(-im)) + exp(im))); else tmp = Float64(cos(re) * Float64(0.5 * (im ^ 2.0))); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (im <= 0.054) tmp = cos(re); elseif (im <= 1.35e+154) tmp = 0.5 * (exp(-im) + exp(im)); else tmp = cos(re) * (0.5 * (im ^ 2.0)); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 0.054], N[Cos[re], $MachinePrecision], If[LessEqual[im, 1.35e+154], N[(0.5 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Cos[re], $MachinePrecision] * N[(0.5 * N[Power[im, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 0.054:\\
\;\;\;\;\cos re\\
\mathbf{elif}\;im \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;0.5 \cdot \left(e^{-im} + e^{im}\right)\\
\mathbf{else}:\\
\;\;\;\;\cos re \cdot \left(0.5 \cdot {im}^{2}\right)\\
\end{array}
\end{array}
if im < 0.0539999999999999994Initial program 100.0%
Taylor expanded in im around 0 65.9%
if 0.0539999999999999994 < im < 1.35000000000000003e154Initial program 100.0%
Taylor expanded in re around 0 79.4%
if 1.35000000000000003e154 < im Initial program 100.0%
Taylor expanded in im around 0 100.0%
associate-*r*100.0%
distribute-rgt1-in100.0%
Simplified100.0%
Taylor expanded in im around inf 100.0%
associate-*r*100.0%
*-commutative100.0%
Simplified100.0%
Final simplification71.3%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* 0.5 (pow im 2.0))))
(if (<= im 0.058)
(* (cos re) (+ t_0 1.0))
(if (<= im 1.35e+154)
(* 0.5 (+ (exp (- im)) (exp im)))
(* (cos re) t_0)))))
double code(double re, double im) {
double t_0 = 0.5 * pow(im, 2.0);
double tmp;
if (im <= 0.058) {
tmp = cos(re) * (t_0 + 1.0);
} else if (im <= 1.35e+154) {
tmp = 0.5 * (exp(-im) + exp(im));
} else {
tmp = cos(re) * t_0;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: t_0
real(8) :: tmp
t_0 = 0.5d0 * (im ** 2.0d0)
if (im <= 0.058d0) then
tmp = cos(re) * (t_0 + 1.0d0)
else if (im <= 1.35d+154) then
tmp = 0.5d0 * (exp(-im) + exp(im))
else
tmp = cos(re) * t_0
end if
code = tmp
end function
public static double code(double re, double im) {
double t_0 = 0.5 * Math.pow(im, 2.0);
double tmp;
if (im <= 0.058) {
tmp = Math.cos(re) * (t_0 + 1.0);
} else if (im <= 1.35e+154) {
tmp = 0.5 * (Math.exp(-im) + Math.exp(im));
} else {
tmp = Math.cos(re) * t_0;
}
return tmp;
}
def code(re, im): t_0 = 0.5 * math.pow(im, 2.0) tmp = 0 if im <= 0.058: tmp = math.cos(re) * (t_0 + 1.0) elif im <= 1.35e+154: tmp = 0.5 * (math.exp(-im) + math.exp(im)) else: tmp = math.cos(re) * t_0 return tmp
function code(re, im) t_0 = Float64(0.5 * (im ^ 2.0)) tmp = 0.0 if (im <= 0.058) tmp = Float64(cos(re) * Float64(t_0 + 1.0)); elseif (im <= 1.35e+154) tmp = Float64(0.5 * Float64(exp(Float64(-im)) + exp(im))); else tmp = Float64(cos(re) * t_0); end return tmp end
function tmp_2 = code(re, im) t_0 = 0.5 * (im ^ 2.0); tmp = 0.0; if (im <= 0.058) tmp = cos(re) * (t_0 + 1.0); elseif (im <= 1.35e+154) tmp = 0.5 * (exp(-im) + exp(im)); else tmp = cos(re) * t_0; end tmp_2 = tmp; end
code[re_, im_] := Block[{t$95$0 = N[(0.5 * N[Power[im, 2.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[im, 0.058], N[(N[Cos[re], $MachinePrecision] * N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 1.35e+154], N[(0.5 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot {im}^{2}\\
\mathbf{if}\;im \leq 0.058:\\
\;\;\;\;\cos re \cdot \left(t\_0 + 1\right)\\
\mathbf{elif}\;im \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;0.5 \cdot \left(e^{-im} + e^{im}\right)\\
\mathbf{else}:\\
\;\;\;\;\cos re \cdot t\_0\\
\end{array}
\end{array}
if im < 0.0580000000000000029Initial program 100.0%
Taylor expanded in im around 0 81.4%
associate-*r*81.4%
distribute-rgt1-in81.4%
Simplified81.4%
if 0.0580000000000000029 < im < 1.35000000000000003e154Initial program 100.0%
Taylor expanded in re around 0 79.4%
if 1.35000000000000003e154 < im Initial program 100.0%
Taylor expanded in im around 0 100.0%
associate-*r*100.0%
distribute-rgt1-in100.0%
Simplified100.0%
Taylor expanded in im around inf 100.0%
associate-*r*100.0%
*-commutative100.0%
Simplified100.0%
Final simplification83.1%
(FPCore (re im) :precision binary64 (if (<= im 0.054) (cos re) (* 0.5 (+ (exp (- im)) (exp im)))))
double code(double re, double im) {
double tmp;
if (im <= 0.054) {
tmp = cos(re);
} 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 <= 0.054d0) then
tmp = cos(re)
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 <= 0.054) {
tmp = Math.cos(re);
} else {
tmp = 0.5 * (Math.exp(-im) + Math.exp(im));
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 0.054: tmp = math.cos(re) else: tmp = 0.5 * (math.exp(-im) + math.exp(im)) return tmp
function code(re, im) tmp = 0.0 if (im <= 0.054) tmp = cos(re); 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 <= 0.054) tmp = cos(re); else tmp = 0.5 * (exp(-im) + exp(im)); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 0.054], N[Cos[re], $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 0.054:\\
\;\;\;\;\cos re\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(e^{-im} + e^{im}\right)\\
\end{array}
\end{array}
if im < 0.0539999999999999994Initial program 100.0%
Taylor expanded in im around 0 65.9%
if 0.0539999999999999994 < im Initial program 100.0%
Taylor expanded in re around 0 78.7%
Final simplification69.0%
(FPCore (re im) :precision binary64 (if (<= im 2.9e+52) (cos re) (cbrt (* (pow im 6.0) 0.125))))
double code(double re, double im) {
double tmp;
if (im <= 2.9e+52) {
tmp = cos(re);
} else {
tmp = cbrt((pow(im, 6.0) * 0.125));
}
return tmp;
}
public static double code(double re, double im) {
double tmp;
if (im <= 2.9e+52) {
tmp = Math.cos(re);
} else {
tmp = Math.cbrt((Math.pow(im, 6.0) * 0.125));
}
return tmp;
}
function code(re, im) tmp = 0.0 if (im <= 2.9e+52) tmp = cos(re); else tmp = cbrt(Float64((im ^ 6.0) * 0.125)); end return tmp end
code[re_, im_] := If[LessEqual[im, 2.9e+52], N[Cos[re], $MachinePrecision], N[Power[N[(N[Power[im, 6.0], $MachinePrecision] * 0.125), $MachinePrecision], 1/3], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 2.9 \cdot 10^{+52}:\\
\;\;\;\;\cos re\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{{im}^{6} \cdot 0.125}\\
\end{array}
\end{array}
if im < 2.9e52Initial program 100.0%
Taylor expanded in im around 0 62.0%
if 2.9e52 < im Initial program 100.0%
Taylor expanded in im around 0 58.7%
associate-*r*58.7%
distribute-rgt1-in58.7%
Simplified58.7%
Taylor expanded in im around inf 58.7%
associate-*r*58.7%
*-commutative58.7%
Simplified58.7%
Taylor expanded in re around 0 45.8%
add-cbrt-cube79.2%
pow1/379.2%
pow379.2%
*-commutative79.2%
unpow-prod-down79.2%
pow-pow79.2%
metadata-eval79.2%
metadata-eval79.2%
Applied egg-rr79.2%
unpow1/379.2%
Simplified79.2%
Final simplification65.2%
(FPCore (re im) :precision binary64 (if (<= im 1.32e+67) (cos re) (+ (* 0.5 (pow im 2.0)) 1.0)))
double code(double re, double im) {
double tmp;
if (im <= 1.32e+67) {
tmp = cos(re);
} else {
tmp = (0.5 * pow(im, 2.0)) + 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.32d+67) then
tmp = cos(re)
else
tmp = (0.5d0 * (im ** 2.0d0)) + 1.0d0
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (im <= 1.32e+67) {
tmp = Math.cos(re);
} else {
tmp = (0.5 * Math.pow(im, 2.0)) + 1.0;
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 1.32e+67: tmp = math.cos(re) else: tmp = (0.5 * math.pow(im, 2.0)) + 1.0 return tmp
function code(re, im) tmp = 0.0 if (im <= 1.32e+67) tmp = cos(re); else tmp = Float64(Float64(0.5 * (im ^ 2.0)) + 1.0); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (im <= 1.32e+67) tmp = cos(re); else tmp = (0.5 * (im ^ 2.0)) + 1.0; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 1.32e+67], N[Cos[re], $MachinePrecision], N[(N[(0.5 * N[Power[im, 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 1.32 \cdot 10^{+67}:\\
\;\;\;\;\cos re\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot {im}^{2} + 1\\
\end{array}
\end{array}
if im < 1.3200000000000001e67Initial program 100.0%
Taylor expanded in im around 0 61.2%
if 1.3200000000000001e67 < im Initial program 100.0%
Taylor expanded in im around 0 62.3%
associate-*r*62.3%
distribute-rgt1-in62.3%
Simplified62.3%
Taylor expanded in re around 0 48.6%
Final simplification59.0%
(FPCore (re im) :precision binary64 (if (<= im 2.3e+65) (cos re) (* 0.5 (pow im 2.0))))
double code(double re, double im) {
double tmp;
if (im <= 2.3e+65) {
tmp = cos(re);
} else {
tmp = 0.5 * pow(im, 2.0);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (im <= 2.3d+65) then
tmp = cos(re)
else
tmp = 0.5d0 * (im ** 2.0d0)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (im <= 2.3e+65) {
tmp = Math.cos(re);
} else {
tmp = 0.5 * Math.pow(im, 2.0);
}
return tmp;
}
def code(re, im): tmp = 0 if im <= 2.3e+65: tmp = math.cos(re) else: tmp = 0.5 * math.pow(im, 2.0) return tmp
function code(re, im) tmp = 0.0 if (im <= 2.3e+65) tmp = cos(re); else tmp = Float64(0.5 * (im ^ 2.0)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (im <= 2.3e+65) tmp = cos(re); else tmp = 0.5 * (im ^ 2.0); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[im, 2.3e+65], N[Cos[re], $MachinePrecision], N[(0.5 * N[Power[im, 2.0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 2.3 \cdot 10^{+65}:\\
\;\;\;\;\cos re\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot {im}^{2}\\
\end{array}
\end{array}
if im < 2.3e65Initial program 100.0%
Taylor expanded in im around 0 61.2%
if 2.3e65 < im Initial program 100.0%
Taylor expanded in im around 0 62.3%
associate-*r*62.3%
distribute-rgt1-in62.3%
Simplified62.3%
Taylor expanded in im around inf 62.3%
associate-*r*62.3%
*-commutative62.3%
Simplified62.3%
Taylor expanded in re around 0 48.6%
Final simplification59.0%
(FPCore (re im) :precision binary64 (cos re))
double code(double re, double im) {
return cos(re);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = cos(re)
end function
public static double code(double re, double im) {
return Math.cos(re);
}
def code(re, im): return math.cos(re)
function code(re, im) return cos(re) end
function tmp = code(re, im) tmp = cos(re); end
code[re_, im_] := N[Cos[re], $MachinePrecision]
\begin{array}{l}
\\
\cos re
\end{array}
Initial program 100.0%
Taylor expanded in im around 0 51.0%
Final simplification51.0%
(FPCore (re im) :precision binary64 0.0)
double code(double re, double im) {
return 0.0;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = 0.0d0
end function
public static double code(double re, double im) {
return 0.0;
}
def code(re, im): return 0.0
function code(re, im) return 0.0 end
function tmp = code(re, im) tmp = 0.0; end
code[re_, im_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 100.0%
Applied egg-rr2.4%
pow-base-12.4%
metadata-eval2.4%
Simplified2.4%
Final simplification2.4%
(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%
Applied egg-rr24.2%
+-inverses24.2%
+-rgt-identity24.2%
*-inverses24.2%
Simplified24.2%
Final simplification24.2%
herbie shell --seed 2024082
(FPCore (re im)
:name "math.cos on complex, real part"
:precision binary64
(* (* 0.5 (cos re)) (+ (exp (- im)) (exp im))))