
(FPCore (x) :precision binary64 (/ 2.0 (+ (exp x) (exp (- x)))))
double code(double x) {
return 2.0 / (exp(x) + exp(-x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (exp(x) + exp(-x))
end function
public static double code(double x) {
return 2.0 / (Math.exp(x) + Math.exp(-x));
}
def code(x): return 2.0 / (math.exp(x) + math.exp(-x))
function code(x) return Float64(2.0 / Float64(exp(x) + exp(Float64(-x)))) end
function tmp = code(x) tmp = 2.0 / (exp(x) + exp(-x)); end
code[x_] := N[(2.0 / N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{e^{x} + e^{-x}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ 2.0 (+ (exp x) (exp (- x)))))
double code(double x) {
return 2.0 / (exp(x) + exp(-x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (exp(x) + exp(-x))
end function
public static double code(double x) {
return 2.0 / (Math.exp(x) + Math.exp(-x));
}
def code(x): return 2.0 / (math.exp(x) + math.exp(-x))
function code(x) return Float64(2.0 / Float64(exp(x) + exp(Float64(-x)))) end
function tmp = code(x) tmp = 2.0 / (exp(x) + exp(-x)); end
code[x_] := N[(2.0 / N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{e^{x} + e^{-x}}
\end{array}
(FPCore (x) :precision binary64 (/ 1.0 (cosh x)))
double code(double x) {
return 1.0 / cosh(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / cosh(x)
end function
public static double code(double x) {
return 1.0 / Math.cosh(x);
}
def code(x): return 1.0 / math.cosh(x)
function code(x) return Float64(1.0 / cosh(x)) end
function tmp = code(x) tmp = 1.0 / cosh(x); end
code[x_] := N[(1.0 / N[Cosh[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\cosh x}
\end{array}
Initial program 100.0%
cosh-undefN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f64N/A
lower-cosh.f64100.0
Applied egg-rr100.0%
(FPCore (x) :precision binary64 (if (<= (+ (exp x) (exp (- x))) 4.0) (fma (* x x) (fma x (* x 0.20833333333333334) -0.5) 1.0) (/ 360.0 (* x (* x (* x x))))))
double code(double x) {
double tmp;
if ((exp(x) + exp(-x)) <= 4.0) {
tmp = fma((x * x), fma(x, (x * 0.20833333333333334), -0.5), 1.0);
} else {
tmp = 360.0 / (x * (x * (x * x)));
}
return tmp;
}
function code(x) tmp = 0.0 if (Float64(exp(x) + exp(Float64(-x))) <= 4.0) tmp = fma(Float64(x * x), fma(x, Float64(x * 0.20833333333333334), -0.5), 1.0); else tmp = Float64(360.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
code[x_] := If[LessEqual[N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision], 4.0], N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * 0.20833333333333334), $MachinePrecision] + -0.5), $MachinePrecision] + 1.0), $MachinePrecision], N[(360.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{x} + e^{-x} \leq 4:\\
\;\;\;\;\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot 0.20833333333333334, -0.5\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{360}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if (+.f64 (exp.f64 x) (exp.f64 (neg.f64 x))) < 4Initial program 100.0%
cosh-undefN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f64N/A
lower-cosh.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
sub-negN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6499.7
Simplified99.7%
if 4 < (+.f64 (exp.f64 x) (exp.f64 (neg.f64 x))) Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6482.7
Simplified82.7%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6462.4
Simplified62.4%
Taylor expanded in x around inf
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in x around inf
lower-/.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
(FPCore (x) :precision binary64 (if (<= (+ (exp x) (exp (- x))) 2e+50) (/ 2.0 (fma x x 2.0)) (/ 360.0 (* x (* x x)))))
double code(double x) {
double tmp;
if ((exp(x) + exp(-x)) <= 2e+50) {
tmp = 2.0 / fma(x, x, 2.0);
} else {
tmp = 360.0 / (x * (x * x));
}
return tmp;
}
function code(x) tmp = 0.0 if (Float64(exp(x) + exp(Float64(-x))) <= 2e+50) tmp = Float64(2.0 / fma(x, x, 2.0)); else tmp = Float64(360.0 / Float64(x * Float64(x * x))); end return tmp end
code[x_] := If[LessEqual[N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision], 2e+50], N[(2.0 / N[(x * x + 2.0), $MachinePrecision]), $MachinePrecision], N[(360.0 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{x} + e^{-x} \leq 2 \cdot 10^{+50}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(x, x, 2\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{360}{x \cdot \left(x \cdot x\right)}\\
\end{array}
\end{array}
if (+.f64 (exp.f64 x) (exp.f64 (neg.f64 x))) < 2.0000000000000002e50Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6498.2
Simplified98.2%
if 2.0000000000000002e50 < (+.f64 (exp.f64 x) (exp.f64 (neg.f64 x))) Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6483.8
Simplified83.8%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6463.3
Simplified63.3%
Taylor expanded in x around -inf
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.3
Simplified63.3%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x x))))
(/
2.0
(fma
(fma
(* x x)
(fma 0.002777777777777778 (* t_0 t_0) 0.08333333333333333)
1.0)
(* x x)
2.0))))
double code(double x) {
double t_0 = x * (x * x);
return 2.0 / fma(fma((x * x), fma(0.002777777777777778, (t_0 * t_0), 0.08333333333333333), 1.0), (x * x), 2.0);
}
function code(x) t_0 = Float64(x * Float64(x * x)) return Float64(2.0 / fma(fma(Float64(x * x), fma(0.002777777777777778, Float64(t_0 * t_0), 0.08333333333333333), 1.0), Float64(x * x), 2.0)) end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, N[(2.0 / N[(N[(N[(x * x), $MachinePrecision] * N[(0.002777777777777778 * N[(t$95$0 * t$95$0), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] + 1.0), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
\frac{2}{\mathsf{fma}\left(\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(0.002777777777777778, t\_0 \cdot t\_0, 0.08333333333333333\right), 1\right), x \cdot x, 2\right)}
\end{array}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
pow-sqrN/A
metadata-evalN/A
lft-mult-inverseN/A
metadata-evalN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
Simplified92.1%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
lower-fma.f6492.1
Applied egg-rr92.1%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6493.2
Simplified93.2%
(FPCore (x) :precision binary64 (/ 2.0 (fma (* x x) (* (* x (* x (* x x))) (* x (* x (* x (- x))))) 2.0)))
double code(double x) {
return 2.0 / fma((x * x), ((x * (x * (x * x))) * (x * (x * (x * -x)))), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(x * Float64(x * Float64(x * Float64(-x))))), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * (-x)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(-x\right)\right)\right)\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f6490.9
Simplified90.9%
Taylor expanded in x around inf
mul-1-negN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-neg-outN/A
mul-1-negN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
Simplified92.9%
(FPCore (x) :precision binary64 (let* ((t_0 (* x (* x x)))) (/ 2.0 (fma (* 0.08333333333333333 (* x (* t_0 t_0))) (* x x) 2.0))))
double code(double x) {
double t_0 = x * (x * x);
return 2.0 / fma((0.08333333333333333 * (x * (t_0 * t_0))), (x * x), 2.0);
}
function code(x) t_0 = Float64(x * Float64(x * x)) return Float64(2.0 / fma(Float64(0.08333333333333333 * Float64(x * Float64(t_0 * t_0))), Float64(x * x), 2.0)) end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, N[(2.0 / N[(N[(0.08333333333333333 * N[(x * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
\frac{2}{\mathsf{fma}\left(0.08333333333333333 \cdot \left(x \cdot \left(t\_0 \cdot t\_0\right)\right), x \cdot x, 2\right)}
\end{array}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
pow-sqrN/A
metadata-evalN/A
lft-mult-inverseN/A
metadata-evalN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
Simplified92.1%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
lower-fma.f6492.1
Applied egg-rr92.1%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6492.2
Simplified92.2%
Final simplification92.2%
(FPCore (x)
:precision binary64
(/
2.0
(fma
(* x x)
(fma
(* x x)
(fma
(* x x)
(fma (* x x) 9.259259259259259e-5 0.002777777777777778)
0.08333333333333333)
1.0)
2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma((x * x), fma((x * x), fma((x * x), 9.259259259259259e-5, 0.002777777777777778), 0.08333333333333333), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(Float64(x * x), fma(Float64(x * x), fma(Float64(x * x), 9.259259259259259e-5, 0.002777777777777778), 0.08333333333333333), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * 9.259259259259259e-5 + 0.002777777777777778), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, 9.259259259259259 \cdot 10^{-5}, 0.002777777777777778\right), 0.08333333333333333\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
lift-*.f64N/A
flip-+N/A
lower-/.f64N/A
Applied egg-rr66.1%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6492.2
Simplified92.2%
(FPCore (x)
:precision binary64
(/
2.0
(fma
(* x x)
(fma
x
(* x (fma x (* 0.002777777777777778 (* x (* x x))) 0.08333333333333333))
1.0)
2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma(x, (x * fma(x, (0.002777777777777778 * (x * (x * x))), 0.08333333333333333)), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(x, Float64(x * fma(x, Float64(0.002777777777777778 * Float64(x * Float64(x * x))), 0.08333333333333333)), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(x * N[(0.002777777777777778 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.08333333333333333), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, 0.002777777777777778 \cdot \left(x \cdot \left(x \cdot x\right)\right), 0.08333333333333333\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
pow-sqrN/A
metadata-evalN/A
lft-mult-inverseN/A
metadata-evalN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
Simplified92.1%
Final simplification92.1%
(FPCore (x) :precision binary64 (/ 2.0 (fma (fma (* x x) (fma (* x x) (* x x) 0.08333333333333333) 1.0) (* x x) 2.0)))
double code(double x) {
return 2.0 / fma(fma((x * x), fma((x * x), (x * x), 0.08333333333333333), 1.0), (x * x), 2.0);
}
function code(x) return Float64(2.0 / fma(fma(Float64(x * x), fma(Float64(x * x), Float64(x * x), 0.08333333333333333), 1.0), Float64(x * x), 2.0)) end
code[x_] := N[(2.0 / N[(N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] + 1.0), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, x \cdot x, 0.08333333333333333\right), 1\right), x \cdot x, 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
pow-sqrN/A
metadata-evalN/A
lft-mult-inverseN/A
metadata-evalN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
Simplified92.1%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
lower-fma.f6492.1
Applied egg-rr92.1%
Taylor expanded in x around inf
unpow2N/A
lower-*.f6492.1
Simplified92.1%
(FPCore (x)
:precision binary64
(/
2.0
(fma
(* x x)
(fma
x
(* x (fma (* x (* x x)) 0.002777777777777778 0.08333333333333333))
1.0)
2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma(x, (x * fma((x * (x * x)), 0.002777777777777778, 0.08333333333333333)), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(x, Float64(x * fma(Float64(x * Float64(x * x)), 0.002777777777777778, 0.08333333333333333)), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * 0.002777777777777778 + 0.08333333333333333), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x \cdot \left(x \cdot x\right), 0.002777777777777778, 0.08333333333333333\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6492.0
Simplified92.0%
(FPCore (x) :precision binary64 (/ 2.0 (fma (* x x) (fma x (* x (fma x (* x x) 0.08333333333333333)) 1.0) 2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma(x, (x * fma(x, (x * x), 0.08333333333333333)), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(x, Float64(x * fma(x, Float64(x * x), 0.08333333333333333)), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(x * N[(x * x), $MachinePrecision] + 0.08333333333333333), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, x \cdot x, 0.08333333333333333\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
pow-sqrN/A
metadata-evalN/A
lft-mult-inverseN/A
metadata-evalN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
Simplified92.1%
Taylor expanded in x around inf
unpow2N/A
lower-*.f6492.0
Simplified92.0%
(FPCore (x) :precision binary64 (/ 2.0 (fma (* x x) (fma x (* x (* x (* x x))) 1.0) 2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma(x, (x * (x * (x * x))), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(x, Float64(x * Float64(x * Float64(x * x))), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \left(x \cdot \left(x \cdot x\right)\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f6490.9
Simplified90.9%
Taylor expanded in x around -inf
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6491.9
Simplified91.9%
(FPCore (x) :precision binary64 (/ 2.0 (fma (* x x) (fma x (* x (fma x (- x) 0.08333333333333333)) 1.0) 2.0)))
double code(double x) {
return 2.0 / fma((x * x), fma(x, (x * fma(x, -x, 0.08333333333333333)), 1.0), 2.0);
}
function code(x) return Float64(2.0 / fma(Float64(x * x), fma(x, Float64(x * fma(x, Float64(-x), 0.08333333333333333)), 1.0), 2.0)) end
code[x_] := N[(2.0 / N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(x * (-x) + 0.08333333333333333), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, -x, 0.08333333333333333\right), 1\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f6490.9
Simplified90.9%
(FPCore (x) :precision binary64 (/ 1.0 (fma (* x x) (fma (* x x) 0.041666666666666664 0.5) 1.0)))
double code(double x) {
return 1.0 / fma((x * x), fma((x * x), 0.041666666666666664, 0.5), 1.0);
}
function code(x) return Float64(1.0 / fma(Float64(x * x), fma(Float64(x * x), 0.041666666666666664, 0.5), 1.0)) end
code[x_] := N[(1.0 / N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * 0.041666666666666664 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, 0.041666666666666664, 0.5\right), 1\right)}
\end{array}
Initial program 100.0%
cosh-undefN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f64N/A
lower-cosh.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.1
Simplified85.1%
(FPCore (x) :precision binary64 (/ 2.0 (fma x (* x (* x x)) 2.0)))
double code(double x) {
return 2.0 / fma(x, (x * (x * x)), 2.0);
}
function code(x) return Float64(2.0 / fma(x, Float64(x * Float64(x * x)), 2.0)) end
code[x_] := N[(2.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x, x \cdot \left(x \cdot x\right), 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6480.0
Simplified80.0%
Taylor expanded in x around inf
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.9
Simplified84.9%
Taylor expanded in x around inf
unpow2N/A
lower-*.f6484.9
Simplified84.9%
(FPCore (x) :precision binary64 (/ 2.0 (fma x (* x x) 2.0)))
double code(double x) {
return 2.0 / fma(x, (x * x), 2.0);
}
function code(x) return Float64(2.0 / fma(x, Float64(x * x), 2.0)) end
code[x_] := N[(2.0 / N[(x * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x, x \cdot x, 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6491.1
Simplified91.1%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6480.0
Simplified80.0%
Taylor expanded in x around inf
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.9
Simplified84.9%
Taylor expanded in x around inf
unpow2N/A
lower-*.f6480.4
Simplified80.4%
(FPCore (x) :precision binary64 (/ 2.0 (fma x x 2.0)))
double code(double x) {
return 2.0 / fma(x, x, 2.0);
}
function code(x) return Float64(2.0 / fma(x, x, 2.0)) end
code[x_] := N[(2.0 / N[(x * x + 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\mathsf{fma}\left(x, x, 2\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.9
Simplified74.9%
(FPCore (x) :precision binary64 1.0)
double code(double x) {
return 1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double x) {
return 1.0;
}
def code(x): return 1.0
function code(x) return 1.0 end
function tmp = code(x) tmp = 1.0; end
code[x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
cosh-undefN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f64N/A
lower-cosh.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
Simplified50.1%
(FPCore (x) :precision binary64 0.5)
double code(double x) {
return 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.5d0
end function
public static double code(double x) {
return 0.5;
}
def code(x): return 0.5
function code(x) return 0.5 end
function tmp = code(x) tmp = 0.5; end
code[x_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
Simplified10.8%
herbie shell --seed 2024214
(FPCore (x)
:name "Hyperbolic secant"
:precision binary64
(/ 2.0 (+ (exp x) (exp (- x)))))