
(FPCore (x) :precision binary64 (* (fmod (exp x) (sqrt (cos x))) (exp (- x))))
double code(double x) {
return fmod(exp(x), sqrt(cos(x))) * exp(-x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x))) * exp(-x)
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x))) * math.exp(-x)
function code(x) return Float64(rem(exp(x), sqrt(cos(x))) * exp(Float64(-x))) end
code[x_] := N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] * N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right) \cdot e^{-x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (* (fmod (exp x) (sqrt (cos x))) (exp (- x))))
double code(double x) {
return fmod(exp(x), sqrt(cos(x))) * exp(-x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x))) * exp(-x)
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x))) * math.exp(-x)
function code(x) return Float64(rem(exp(x), sqrt(cos(x))) * exp(Float64(-x))) end
code[x_] := N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] * N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right) \cdot e^{-x}
\end{array}
(FPCore (x) :precision binary64 (+ (+ 1.0 (/ (fmod (exp x) (sqrt (cos x))) (exp x))) -1.0))
double code(double x) {
return (1.0 + (fmod(exp(x), sqrt(cos(x))) / exp(x))) + -1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (1.0d0 + (mod(exp(x), sqrt(cos(x))) / exp(x))) + (-1.0d0)
end function
def code(x): return (1.0 + (math.fmod(math.exp(x), math.sqrt(math.cos(x))) / math.exp(x))) + -1.0
function code(x) return Float64(Float64(1.0 + Float64(rem(exp(x), sqrt(cos(x))) / exp(x))) + -1.0) end
code[x_] := N[(N[(1.0 + N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] / N[Exp[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]
\begin{array}{l}
\\
\left(1 + \frac{\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right)}{e^{x}}\right) + -1
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
expm1-log1p-u7.8%
expm1-udef7.8%
log1p-udef7.8%
add-exp-log7.8%
Applied egg-rr7.8%
Final simplification7.8%
(FPCore (x) :precision binary64 (/ 1.0 (/ (exp x) (fmod (exp x) (sqrt (cos x))))))
double code(double x) {
return 1.0 / (exp(x) / fmod(exp(x), sqrt(cos(x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (exp(x) / mod(exp(x), sqrt(cos(x))))
end function
def code(x): return 1.0 / (math.exp(x) / math.fmod(math.exp(x), math.sqrt(math.cos(x))))
function code(x) return Float64(1.0 / Float64(exp(x) / rem(exp(x), sqrt(cos(x))))) end
code[x_] := N[(1.0 / N[(N[Exp[x], $MachinePrecision] / N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{e^{x}}{\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right)}}
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
expm1-log1p-u7.8%
expm1-udef7.8%
log1p-udef7.8%
add-exp-log7.8%
Applied egg-rr7.8%
add-exp-log7.8%
log1p-udef7.8%
expm1-udef7.8%
expm1-log1p-u7.8%
clear-num7.8%
Applied egg-rr7.8%
Final simplification7.8%
(FPCore (x) :precision binary64 (/ (fmod (exp x) (sqrt (cos x))) (exp x)))
double code(double x) {
return fmod(exp(x), sqrt(cos(x))) / exp(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x))) / exp(x)
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x))) / math.exp(x)
function code(x) return Float64(rem(exp(x), sqrt(cos(x))) / exp(x)) end
code[x_] := N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] / N[Exp[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right)}{e^{x}}
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
Final simplification7.8%
(FPCore (x) :precision binary64 (/ (fmod (exp x) (sqrt (cos x))) (/ (+ 1.0 x) (- 1.0 (* x x)))))
double code(double x) {
return fmod(exp(x), sqrt(cos(x))) / ((1.0 + x) / (1.0 - (x * x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x))) / ((1.0d0 + x) / (1.0d0 - (x * x)))
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x))) / ((1.0 + x) / (1.0 - (x * x)))
function code(x) return Float64(rem(exp(x), sqrt(cos(x))) / Float64(Float64(1.0 + x) / Float64(1.0 - Float64(x * x)))) end
code[x_] := N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] / N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right)}{\frac{1 + x}{1 - x \cdot x}}
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
Taylor expanded in x around 0 7.0%
mul-1-neg7.0%
unsub-neg7.0%
*-rgt-identity7.0%
distribute-lft-out--7.0%
Simplified7.0%
flip--7.0%
associate-*r/7.0%
metadata-eval7.0%
Applied egg-rr7.0%
associate-/l*7.0%
+-commutative7.0%
Simplified7.0%
Final simplification7.0%
(FPCore (x) :precision binary64 (* (fmod (exp x) (sqrt (cos x))) (- 1.0 x)))
double code(double x) {
return fmod(exp(x), sqrt(cos(x))) * (1.0 - x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x))) * (1.0d0 - x)
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x))) * (1.0 - x)
function code(x) return Float64(rem(exp(x), sqrt(cos(x))) * Float64(1.0 - x)) end
code[x_] := N[(N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision] * N[(1.0 - x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right) \cdot \left(1 - x\right)
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
Taylor expanded in x around 0 7.0%
mul-1-neg7.0%
unsub-neg7.0%
*-rgt-identity7.0%
distribute-lft-out--7.0%
Simplified7.0%
Final simplification7.0%
(FPCore (x) :precision binary64 (fmod (exp x) (sqrt (cos x))))
double code(double x) {
return fmod(exp(x), sqrt(cos(x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = mod(exp(x), sqrt(cos(x)))
end function
def code(x): return math.fmod(math.exp(x), math.sqrt(math.cos(x)))
function code(x) return rem(exp(x), sqrt(cos(x))) end
code[x_] := N[With[{TMP1 = N[Exp[x], $MachinePrecision], TMP2 = N[Sqrt[N[Cos[x], $MachinePrecision]], $MachinePrecision]}, Mod[Abs[TMP1], Abs[TMP2]] * Sign[TMP1]], $MachinePrecision]
\begin{array}{l}
\\
\left(\left(e^{x}\right) \bmod \left(\sqrt{\cos x}\right)\right)
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
Taylor expanded in x around 0 6.8%
Final simplification6.8%
(FPCore (x) :precision binary64 (exp (- x)))
double code(double x) {
return exp(-x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp(-x)
end function
public static double code(double x) {
return Math.exp(-x);
}
def code(x): return math.exp(-x)
function code(x) return exp(Float64(-x)) end
function tmp = code(x) tmp = exp(-x); end
code[x_] := N[Exp[(-x)], $MachinePrecision]
\begin{array}{l}
\\
e^{-x}
\end{array}
Initial program 7.8%
exp-neg7.8%
associate-*r/7.8%
*-rgt-identity7.8%
Simplified7.8%
add-exp-log7.8%
div-exp7.8%
Applied egg-rr7.8%
Taylor expanded in x around inf 5.2%
neg-mul-15.2%
Simplified5.2%
Final simplification5.2%
herbie shell --seed 2023278
(FPCore (x)
:name "expfmod (used to be hard to sample)"
:precision binary64
(* (fmod (exp x) (sqrt (cos x))) (exp (- x))))