
(FPCore (x) :precision binary64 (/ (- (exp x) (exp (- x))) 2.0))
double code(double x) {
return (exp(x) - exp(-x)) / 2.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (exp(x) - exp(-x)) / 2.0d0
end function
public static double code(double x) {
return (Math.exp(x) - Math.exp(-x)) / 2.0;
}
def code(x): return (math.exp(x) - math.exp(-x)) / 2.0
function code(x) return Float64(Float64(exp(x) - exp(Float64(-x))) / 2.0) end
function tmp = code(x) tmp = (exp(x) - exp(-x)) / 2.0; end
code[x_] := N[(N[(N[Exp[x], $MachinePrecision] - N[Exp[(-x)], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x} - e^{-x}}{2}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ (- (exp x) (exp (- x))) 2.0))
double code(double x) {
return (exp(x) - exp(-x)) / 2.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (exp(x) - exp(-x)) / 2.0d0
end function
public static double code(double x) {
return (Math.exp(x) - Math.exp(-x)) / 2.0;
}
def code(x): return (math.exp(x) - math.exp(-x)) / 2.0
function code(x) return Float64(Float64(exp(x) - exp(Float64(-x))) / 2.0) end
function tmp = code(x) tmp = (exp(x) - exp(-x)) / 2.0; end
code[x_] := N[(N[(N[Exp[x], $MachinePrecision] - N[Exp[(-x)], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x} - e^{-x}}{2}
\end{array}
(FPCore (x) :precision binary64 (sinh x))
double code(double x) {
return sinh(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = sinh(x)
end function
public static double code(double x) {
return Math.sinh(x);
}
def code(x): return math.sinh(x)
function code(x) return sinh(x) end
function tmp = code(x) tmp = sinh(x); end
code[x_] := N[Sinh[x], $MachinePrecision]
\begin{array}{l}
\\
\sinh x
\end{array}
Initial program 54.8%
lift--.f64N/A
lift-exp.f64N/A
lift-exp.f64N/A
lift-neg.f64N/A
sinh-undefN/A
*-commutativeN/A
lower-*.f64N/A
lower-sinh.f64100.0
Applied rewrites100.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
metadata-evalN/A
*-rgt-identity100.0
Applied rewrites100.0%
(FPCore (x)
:precision binary64
(*
0.5
(*
(fma
(fma
(fma 0.0003968253968253968 (* x x) 0.016666666666666666)
(* x x)
0.3333333333333333)
(* x x)
2.0)
x)))
double code(double x) {
return 0.5 * (fma(fma(fma(0.0003968253968253968, (x * x), 0.016666666666666666), (x * x), 0.3333333333333333), (x * x), 2.0) * x);
}
function code(x) return Float64(0.5 * Float64(fma(fma(fma(0.0003968253968253968, Float64(x * x), 0.016666666666666666), Float64(x * x), 0.3333333333333333), Float64(x * x), 2.0) * x)) end
code[x_] := N[(0.5 * N[(N[(N[(N[(0.0003968253968253968 * N[(x * x), $MachinePrecision] + 0.016666666666666666), $MachinePrecision] * N[(x * x), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.0003968253968253968, x \cdot x, 0.016666666666666666\right), x \cdot x, 0.3333333333333333\right), x \cdot x, 2\right) \cdot x\right)
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6489.7
Applied rewrites89.7%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval89.7
Applied rewrites89.7%
Final simplification89.7%
(FPCore (x)
:precision binary64
(*
(*
(fma
(fma (* (* x x) 0.0003968253968253968) (* x x) 0.3333333333333333)
(* x x)
2.0)
x)
0.5))
double code(double x) {
return (fma(fma(((x * x) * 0.0003968253968253968), (x * x), 0.3333333333333333), (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(fma(Float64(Float64(x * x) * 0.0003968253968253968), Float64(x * x), 0.3333333333333333), Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(N[(N[(N[(x * x), $MachinePrecision] * 0.0003968253968253968), $MachinePrecision] * N[(x * x), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(\mathsf{fma}\left(\left(x \cdot x\right) \cdot 0.0003968253968253968, x \cdot x, 0.3333333333333333\right), x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6489.7
Applied rewrites89.7%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval89.7
Applied rewrites89.7%
Taylor expanded in x around inf
Applied rewrites89.6%
Final simplification89.6%
(FPCore (x)
:precision binary64
(*
(*
(fma
(* (* (fma 0.0003968253968253968 (* x x) 0.016666666666666666) x) x)
(* x x)
2.0)
x)
0.5))
double code(double x) {
return (fma(((fma(0.0003968253968253968, (x * x), 0.016666666666666666) * x) * x), (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(Float64(Float64(fma(0.0003968253968253968, Float64(x * x), 0.016666666666666666) * x) * x), Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(N[(N[(N[(0.0003968253968253968 * N[(x * x), $MachinePrecision] + 0.016666666666666666), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(\left(\mathsf{fma}\left(0.0003968253968253968, x \cdot x, 0.016666666666666666\right) \cdot x\right) \cdot x, x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6489.7
Applied rewrites89.7%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval89.7
Applied rewrites89.7%
Taylor expanded in x around inf
Applied rewrites89.3%
(FPCore (x) :precision binary64 (* (* (fma (* (* (* (* x x) 0.0003968253968253968) x) x) (* x x) 2.0) x) 0.5))
double code(double x) {
return (fma(((((x * x) * 0.0003968253968253968) * x) * x), (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(Float64(Float64(Float64(Float64(x * x) * 0.0003968253968253968) * x) * x), Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(N[(N[(N[(N[(x * x), $MachinePrecision] * 0.0003968253968253968), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(\left(\left(\left(x \cdot x\right) \cdot 0.0003968253968253968\right) \cdot x\right) \cdot x, x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6489.7
Applied rewrites89.7%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval89.7
Applied rewrites89.7%
Taylor expanded in x around inf
Applied rewrites89.3%
Taylor expanded in x around inf
Applied rewrites89.3%
(FPCore (x)
:precision binary64
(if (<= x 3.3)
(* (* (fma 0.3333333333333333 (* x x) 2.0) x) 0.5)
(*
(* (* (* (fma 0.016666666666666666 (* x x) 0.3333333333333333) x) x) x)
0.5)))
double code(double x) {
double tmp;
if (x <= 3.3) {
tmp = (fma(0.3333333333333333, (x * x), 2.0) * x) * 0.5;
} else {
tmp = (((fma(0.016666666666666666, (x * x), 0.3333333333333333) * x) * x) * x) * 0.5;
}
return tmp;
}
function code(x) tmp = 0.0 if (x <= 3.3) tmp = Float64(Float64(fma(0.3333333333333333, Float64(x * x), 2.0) * x) * 0.5); else tmp = Float64(Float64(Float64(Float64(fma(0.016666666666666666, Float64(x * x), 0.3333333333333333) * x) * x) * x) * 0.5); end return tmp end
code[x_] := If[LessEqual[x, 3.3], N[(N[(N[(0.3333333333333333 * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(N[(N[(0.016666666666666666 * N[(x * x), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.3:\\
\;\;\;\;\left(\mathsf{fma}\left(0.3333333333333333, x \cdot x, 2\right) \cdot x\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\mathsf{fma}\left(0.016666666666666666, x \cdot x, 0.3333333333333333\right) \cdot x\right) \cdot x\right) \cdot x\right) \cdot 0.5\\
\end{array}
\end{array}
if x < 3.2999999999999998Initial program 43.0%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.1
Applied rewrites87.1%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-evalN/A
Applied rewrites87.1%
if 3.2999999999999998 < x Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.8
Applied rewrites67.8%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval67.8
Applied rewrites67.8%
Taylor expanded in x around inf
Applied rewrites67.8%
(FPCore (x) :precision binary64 (* (* (fma (fma 0.016666666666666666 (* x x) 0.3333333333333333) (* x x) 2.0) x) 0.5))
double code(double x) {
return (fma(fma(0.016666666666666666, (x * x), 0.3333333333333333), (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(fma(0.016666666666666666, Float64(x * x), 0.3333333333333333), Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(N[(0.016666666666666666 * N[(x * x), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(\mathsf{fma}\left(0.016666666666666666, x \cdot x, 0.3333333333333333\right), x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6486.2
Applied rewrites86.2%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval86.2
Applied rewrites86.2%
(FPCore (x) :precision binary64 (* (* (fma (* 0.016666666666666666 (* x x)) (* x x) 2.0) x) 0.5))
double code(double x) {
return (fma((0.016666666666666666 * (x * x)), (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(Float64(0.016666666666666666 * Float64(x * x)), Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(N[(0.016666666666666666 * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(0.016666666666666666 \cdot \left(x \cdot x\right), x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6486.2
Applied rewrites86.2%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval86.2
Applied rewrites86.2%
Taylor expanded in x around inf
Applied rewrites85.8%
(FPCore (x) :precision binary64 (* (* (fma 0.3333333333333333 (* x x) 2.0) x) 0.5))
double code(double x) {
return (fma(0.3333333333333333, (x * x), 2.0) * x) * 0.5;
}
function code(x) return Float64(Float64(fma(0.3333333333333333, Float64(x * x), 2.0) * x) * 0.5) end
code[x_] := N[(N[(N[(0.3333333333333333 * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\mathsf{fma}\left(0.3333333333333333, x \cdot x, 2\right) \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6481.2
Applied rewrites81.2%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-evalN/A
Applied rewrites81.2%
(FPCore (x) :precision binary64 (* (* 2.0 x) 0.5))
double code(double x) {
return (2.0 * x) * 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (2.0d0 * x) * 0.5d0
end function
public static double code(double x) {
return (2.0 * x) * 0.5;
}
def code(x): return (2.0 * x) * 0.5
function code(x) return Float64(Float64(2.0 * x) * 0.5) end
function tmp = code(x) tmp = (2.0 * x) * 0.5; end
code[x_] := N[(N[(2.0 * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(2 \cdot x\right) \cdot 0.5
\end{array}
Initial program 54.8%
Taylor expanded in x around 0
lower-*.f6451.8
Applied rewrites51.8%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval51.8
Applied rewrites51.8%
herbie shell --seed 2024283
(FPCore (x)
:name "Hyperbolic sine"
:precision binary64
(/ (- (exp x) (exp (- x))) 2.0))