
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
(FPCore (re im) :precision binary64 (* (exp re) (sin im)))
double code(double re, double im) {
return exp(re) * sin(im);
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = exp(re) * sin(im)
end function
public static double code(double re, double im) {
return Math.exp(re) * Math.sin(im);
}
def code(re, im): return math.exp(re) * math.sin(im)
function code(re, im) return Float64(exp(re) * sin(im)) end
function tmp = code(re, im) tmp = exp(re) * sin(im); end
code[re_, im_] := N[(N[Exp[re], $MachinePrecision] * N[Sin[im], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{re} \cdot \sin im
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (re im) :precision binary64 (if (<= (exp re) 0.0) 0.0 (if (<= (exp re) 1.000000000005) (* (sin im) (+ re 1.0)) (* (exp re) im))))
double code(double re, double im) {
double tmp;
if (exp(re) <= 0.0) {
tmp = 0.0;
} else if (exp(re) <= 1.000000000005) {
tmp = sin(im) * (re + 1.0);
} else {
tmp = exp(re) * im;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (exp(re) <= 0.0d0) then
tmp = 0.0d0
else if (exp(re) <= 1.000000000005d0) then
tmp = sin(im) * (re + 1.0d0)
else
tmp = exp(re) * im
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (Math.exp(re) <= 0.0) {
tmp = 0.0;
} else if (Math.exp(re) <= 1.000000000005) {
tmp = Math.sin(im) * (re + 1.0);
} else {
tmp = Math.exp(re) * im;
}
return tmp;
}
def code(re, im): tmp = 0 if math.exp(re) <= 0.0: tmp = 0.0 elif math.exp(re) <= 1.000000000005: tmp = math.sin(im) * (re + 1.0) else: tmp = math.exp(re) * im return tmp
function code(re, im) tmp = 0.0 if (exp(re) <= 0.0) tmp = 0.0; elseif (exp(re) <= 1.000000000005) tmp = Float64(sin(im) * Float64(re + 1.0)); else tmp = Float64(exp(re) * im); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (exp(re) <= 0.0) tmp = 0.0; elseif (exp(re) <= 1.000000000005) tmp = sin(im) * (re + 1.0); else tmp = exp(re) * im; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[N[Exp[re], $MachinePrecision], 0.0], 0.0, If[LessEqual[N[Exp[re], $MachinePrecision], 1.000000000005], N[(N[Sin[im], $MachinePrecision] * N[(re + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[Exp[re], $MachinePrecision] * im), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 0:\\
\;\;\;\;0\\
\mathbf{elif}\;e^{re} \leq 1.000000000005:\\
\;\;\;\;\sin im \cdot \left(re + 1\right)\\
\mathbf{else}:\\
\;\;\;\;e^{re} \cdot im\\
\end{array}
\end{array}
if (exp.f64 re) < 0.0Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if 0.0 < (exp.f64 re) < 1.000000000005Initial program 99.9%
Taylor expanded in re around 0 99.5%
distribute-rgt1-in99.5%
Simplified99.5%
if 1.000000000005 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0 79.6%
Final simplification94.7%
(FPCore (re im) :precision binary64 (if (or (<= (exp re) 1.0) (not (<= (exp re) 1.000000000005))) (* (exp re) im) (sin im)))
double code(double re, double im) {
double tmp;
if ((exp(re) <= 1.0) || !(exp(re) <= 1.000000000005)) {
tmp = exp(re) * im;
} else {
tmp = sin(im);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if ((exp(re) <= 1.0d0) .or. (.not. (exp(re) <= 1.000000000005d0))) then
tmp = exp(re) * im
else
tmp = sin(im)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((Math.exp(re) <= 1.0) || !(Math.exp(re) <= 1.000000000005)) {
tmp = Math.exp(re) * im;
} else {
tmp = Math.sin(im);
}
return tmp;
}
def code(re, im): tmp = 0 if (math.exp(re) <= 1.0) or not (math.exp(re) <= 1.000000000005): tmp = math.exp(re) * im else: tmp = math.sin(im) return tmp
function code(re, im) tmp = 0.0 if ((exp(re) <= 1.0) || !(exp(re) <= 1.000000000005)) tmp = Float64(exp(re) * im); else tmp = sin(im); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((exp(re) <= 1.0) || ~((exp(re) <= 1.000000000005))) tmp = exp(re) * im; else tmp = sin(im); end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[N[Exp[re], $MachinePrecision], 1.0], N[Not[LessEqual[N[Exp[re], $MachinePrecision], 1.000000000005]], $MachinePrecision]], N[(N[Exp[re], $MachinePrecision] * im), $MachinePrecision], N[Sin[im], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{re} \leq 1 \lor \neg \left(e^{re} \leq 1.000000000005\right):\\
\;\;\;\;e^{re} \cdot im\\
\mathbf{else}:\\
\;\;\;\;\sin im\\
\end{array}
\end{array}
if (exp.f64 re) < 1 or 1.000000000005 < (exp.f64 re) Initial program 100.0%
Taylor expanded in im around 0 66.8%
if 1 < (exp.f64 re) < 1.000000000005Initial program 100.0%
Taylor expanded in re around 0 95.3%
Final simplification66.9%
(FPCore (re im) :precision binary64 (if (<= re -160.0) 0.0 (if (<= re 70.0) (sin im) (pow E re))))
double code(double re, double im) {
double tmp;
if (re <= -160.0) {
tmp = 0.0;
} else if (re <= 70.0) {
tmp = sin(im);
} else {
tmp = pow(((double) M_E), re);
}
return tmp;
}
public static double code(double re, double im) {
double tmp;
if (re <= -160.0) {
tmp = 0.0;
} else if (re <= 70.0) {
tmp = Math.sin(im);
} else {
tmp = Math.pow(Math.E, re);
}
return tmp;
}
def code(re, im): tmp = 0 if re <= -160.0: tmp = 0.0 elif re <= 70.0: tmp = math.sin(im) else: tmp = math.pow(math.e, re) return tmp
function code(re, im) tmp = 0.0 if (re <= -160.0) tmp = 0.0; elseif (re <= 70.0) tmp = sin(im); else tmp = exp(1) ^ re; end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= -160.0) tmp = 0.0; elseif (re <= 70.0) tmp = sin(im); else tmp = 2.71828182845904523536 ^ re; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, -160.0], 0.0, If[LessEqual[re, 70.0], N[Sin[im], $MachinePrecision], N[Power[E, re], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -160:\\
\;\;\;\;0\\
\mathbf{elif}\;re \leq 70:\\
\;\;\;\;\sin im\\
\mathbf{else}:\\
\;\;\;\;{e}^{re}\\
\end{array}
\end{array}
if re < -160Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if -160 < re < 70Initial program 99.9%
Taylor expanded in re around 0 97.4%
if 70 < re Initial program 100.0%
add-exp-log50.0%
*-un-lft-identity50.0%
exp-prod50.0%
exp-1-e50.0%
log-prod50.0%
add-log-exp50.0%
Applied egg-rr50.0%
Taylor expanded in re around inf 50.0%
Final simplification86.6%
(FPCore (re im) :precision binary64 (if (<= re -130.0) 0.0 (if (<= re 6.8e-12) (sin im) (+ im (* re im)))))
double code(double re, double im) {
double tmp;
if (re <= -130.0) {
tmp = 0.0;
} else if (re <= 6.8e-12) {
tmp = sin(im);
} else {
tmp = im + (re * im);
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (re <= (-130.0d0)) then
tmp = 0.0d0
else if (re <= 6.8d-12) then
tmp = sin(im)
else
tmp = im + (re * im)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (re <= -130.0) {
tmp = 0.0;
} else if (re <= 6.8e-12) {
tmp = Math.sin(im);
} else {
tmp = im + (re * im);
}
return tmp;
}
def code(re, im): tmp = 0 if re <= -130.0: tmp = 0.0 elif re <= 6.8e-12: tmp = math.sin(im) else: tmp = im + (re * im) return tmp
function code(re, im) tmp = 0.0 if (re <= -130.0) tmp = 0.0; elseif (re <= 6.8e-12) tmp = sin(im); else tmp = Float64(im + Float64(re * im)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= -130.0) tmp = 0.0; elseif (re <= 6.8e-12) tmp = sin(im); else tmp = im + (re * im); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, -130.0], 0.0, If[LessEqual[re, 6.8e-12], N[Sin[im], $MachinePrecision], N[(im + N[(re * im), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -130:\\
\;\;\;\;0\\
\mathbf{elif}\;re \leq 6.8 \cdot 10^{-12}:\\
\;\;\;\;\sin im\\
\mathbf{else}:\\
\;\;\;\;im + re \cdot im\\
\end{array}
\end{array}
if re < -130Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if -130 < re < 6.8000000000000001e-12Initial program 99.9%
Taylor expanded in re around 0 97.9%
if 6.8000000000000001e-12 < re Initial program 100.0%
Taylor expanded in im around 0 79.6%
Taylor expanded in re around 0 22.9%
Final simplification79.7%
(FPCore (re im) :precision binary64 (if (<= re -90.0) 0.0 (if (<= re 1.0) im (* re im))))
double code(double re, double im) {
double tmp;
if (re <= -90.0) {
tmp = 0.0;
} else if (re <= 1.0) {
tmp = im;
} else {
tmp = re * im;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (re <= (-90.0d0)) then
tmp = 0.0d0
else if (re <= 1.0d0) then
tmp = im
else
tmp = re * im
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (re <= -90.0) {
tmp = 0.0;
} else if (re <= 1.0) {
tmp = im;
} else {
tmp = re * im;
}
return tmp;
}
def code(re, im): tmp = 0 if re <= -90.0: tmp = 0.0 elif re <= 1.0: tmp = im else: tmp = re * im return tmp
function code(re, im) tmp = 0.0 if (re <= -90.0) tmp = 0.0; elseif (re <= 1.0) tmp = im; else tmp = Float64(re * im); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= -90.0) tmp = 0.0; elseif (re <= 1.0) tmp = im; else tmp = re * im; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, -90.0], 0.0, If[LessEqual[re, 1.0], im, N[(re * im), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -90:\\
\;\;\;\;0\\
\mathbf{elif}\;re \leq 1:\\
\;\;\;\;im\\
\mathbf{else}:\\
\;\;\;\;re \cdot im\\
\end{array}
\end{array}
if re < -90Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if -90 < re < 1Initial program 99.9%
Taylor expanded in im around 0 42.4%
Taylor expanded in re around 0 40.7%
if 1 < re Initial program 100.0%
Taylor expanded in re around 0 4.7%
distribute-rgt1-in4.7%
Simplified4.7%
Taylor expanded in re around inf 4.7%
*-commutative4.7%
Simplified4.7%
Taylor expanded in im around 0 20.6%
Final simplification51.6%
(FPCore (re im) :precision binary64 (if (<= re -1.15) 0.0 (+ im (* re im))))
double code(double re, double im) {
double tmp;
if (re <= -1.15) {
tmp = 0.0;
} else {
tmp = im + (re * 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.15d0)) then
tmp = 0.0d0
else
tmp = im + (re * im)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (re <= -1.15) {
tmp = 0.0;
} else {
tmp = im + (re * im);
}
return tmp;
}
def code(re, im): tmp = 0 if re <= -1.15: tmp = 0.0 else: tmp = im + (re * im) return tmp
function code(re, im) tmp = 0.0 if (re <= -1.15) tmp = 0.0; else tmp = Float64(im + Float64(re * im)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= -1.15) tmp = 0.0; else tmp = im + (re * im); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, -1.15], 0.0, N[(im + N[(re * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1.15:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;im + re \cdot im\\
\end{array}
\end{array}
if re < -1.1499999999999999Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if -1.1499999999999999 < re Initial program 100.0%
Taylor expanded in im around 0 54.5%
Taylor expanded in re around 0 35.2%
Final simplification52.4%
(FPCore (re im) :precision binary64 (if (<= re -62.0) 0.0 im))
double code(double re, double im) {
double tmp;
if (re <= -62.0) {
tmp = 0.0;
} else {
tmp = im;
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (re <= (-62.0d0)) then
tmp = 0.0d0
else
tmp = im
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (re <= -62.0) {
tmp = 0.0;
} else {
tmp = im;
}
return tmp;
}
def code(re, im): tmp = 0 if re <= -62.0: tmp = 0.0 else: tmp = im return tmp
function code(re, im) tmp = 0.0 if (re <= -62.0) tmp = 0.0; else tmp = im; end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= -62.0) tmp = 0.0; else tmp = im; end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, -62.0], 0.0, im]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -62:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;im\\
\end{array}
\end{array}
if re < -62Initial program 100.0%
Taylor expanded in re around 0 2.8%
distribute-rgt1-in2.8%
Simplified2.8%
expm1-log1p-u2.8%
expm1-undefine55.4%
log1p-undefine55.4%
rem-exp-log55.4%
Applied egg-rr55.4%
sub-neg55.4%
distribute-rgt-in55.4%
+-commutative55.4%
metadata-eval55.4%
Applied egg-rr55.4%
Taylor expanded in im around 0 8.8%
associate-+r+100.0%
distribute-rgt1-in100.0%
metadata-eval100.0%
mul0-lft100.0%
Simplified100.0%
if -62 < re Initial program 100.0%
Taylor expanded in im around 0 54.5%
Taylor expanded in re around 0 28.3%
Final simplification47.3%
(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%
Taylor expanded in re around 0 50.8%
distribute-rgt1-in50.8%
Simplified50.8%
expm1-log1p-u50.6%
expm1-undefine46.0%
log1p-undefine46.0%
rem-exp-log46.0%
Applied egg-rr46.0%
sub-neg46.0%
distribute-rgt-in46.0%
+-commutative46.0%
metadata-eval46.0%
Applied egg-rr46.0%
Taylor expanded in im around 0 4.8%
associate-+r+29.0%
distribute-rgt1-in29.0%
metadata-eval29.0%
mul0-lft29.0%
Simplified29.0%
Final simplification29.0%
herbie shell --seed 2024062
(FPCore (re im)
:name "math.exp on complex, imaginary part"
:precision binary64
(* (exp re) (sin im)))