
(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 13 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%
clear-numN/A
/-lowering-/.f64N/A
cosh-defN/A
cosh-lowering-cosh.f64100.0%
Applied egg-rr100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x)))
(t_1 (+ 0.08333333333333333 (* x (* x 0.002777777777777778))))
(t_2 (* (* x x) t_1)))
(if (<= x 5e+38)
(/
2.0
(+
2.0
(/
(* (* x x) (+ 1.0 (* t_2 (* (* x x) (* t_1 t_2)))))
(+ 1.0 (* t_2 (+ t_2 -1.0))))))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = 0.08333333333333333 + (x * (x * 0.002777777777777778));
double t_2 = (x * x) * t_1;
double tmp;
if (x <= 5e+38) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * ((x * x) * (t_1 * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0)))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = (x * x) * (x * x)
t_1 = 0.08333333333333333d0 + (x * (x * 0.002777777777777778d0))
t_2 = (x * x) * t_1
if (x <= 5d+38) then
tmp = 2.0d0 / (2.0d0 + (((x * x) * (1.0d0 + (t_2 * ((x * x) * (t_1 * t_2))))) / (1.0d0 + (t_2 * (t_2 + (-1.0d0))))))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = 0.08333333333333333 + (x * (x * 0.002777777777777778));
double t_2 = (x * x) * t_1;
double tmp;
if (x <= 5e+38) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * ((x * x) * (t_1 * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0)))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) t_1 = 0.08333333333333333 + (x * (x * 0.002777777777777778)) t_2 = (x * x) * t_1 tmp = 0 if x <= 5e+38: tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * ((x * x) * (t_1 * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0))))) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) t_1 = Float64(0.08333333333333333 + Float64(x * Float64(x * 0.002777777777777778))) t_2 = Float64(Float64(x * x) * t_1) tmp = 0.0 if (x <= 5e+38) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(t_2 * Float64(Float64(x * x) * Float64(t_1 * t_2))))) / Float64(1.0 + Float64(t_2 * Float64(t_2 + -1.0)))))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); t_1 = 0.08333333333333333 + (x * (x * 0.002777777777777778)); t_2 = (x * x) * t_1; tmp = 0.0; if (x <= 5e+38) tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * ((x * x) * (t_1 * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0))))); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(0.08333333333333333 + N[(x * N[(x * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * x), $MachinePrecision] * t$95$1), $MachinePrecision]}, If[LessEqual[x, 5e+38], N[(2.0 / N[(2.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(t$95$2 * N[(N[(x * x), $MachinePrecision] * N[(t$95$1 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$2 * N[(t$95$2 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
t_1 := 0.08333333333333333 + x \cdot \left(x \cdot 0.002777777777777778\right)\\
t_2 := \left(x \cdot x\right) \cdot t\_1\\
\mathbf{if}\;x \leq 5 \cdot 10^{+38}:\\
\;\;\;\;\frac{2}{2 + \frac{\left(x \cdot x\right) \cdot \left(1 + t\_2 \cdot \left(\left(x \cdot x\right) \cdot \left(t\_1 \cdot t\_2\right)\right)\right)}{1 + t\_2 \cdot \left(t\_2 + -1\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 4.9999999999999997e38Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6488.8%
Simplified88.8%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr66.0%
if 4.9999999999999997e38 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.5%
Simplified96.5%
distribute-rgt-inN/A
*-lft-identityN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr3.8%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6413.2%
Simplified13.2%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification73.0%
(FPCore (x)
:precision binary64
(/
2.0
(+
2.0
(*
(* x x)
(+
1.0
(*
x
(*
x
(+
0.08333333333333333
(*
(* (* x x) (* x x))
(+ -0.000462962962962963 (* (* x x) -4.6296296296296294e-5)))))))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (((x * x) * (x * x)) * (-0.000462962962962963 + ((x * x) * -4.6296296296296294e-5)))))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * (0.08333333333333333d0 + (((x * x) * (x * x)) * ((-0.000462962962962963d0) + ((x * x) * (-4.6296296296296294d-5))))))))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (((x * x) * (x * x)) * (-0.000462962962962963 + ((x * x) * -4.6296296296296294e-5)))))))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (((x * x) * (x * x)) * (-0.000462962962962963 + ((x * x) * -4.6296296296296294e-5)))))))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(Float64(Float64(x * x) * Float64(x * x)) * Float64(-0.000462962962962963 + Float64(Float64(x * x) * -4.6296296296296294e-5)))))))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (((x * x) * (x * x)) * (-0.000462962962962963 + ((x * x) * -4.6296296296296294e-5))))))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(-0.000462962962962963 + N[(N[(x * x), $MachinePrecision] * -4.6296296296296294e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(-0.000462962962962963 + \left(x \cdot x\right) \cdot -4.6296296296296294 \cdot 10^{-5}\right)\right)\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
distribute-rgt-inN/A
*-lft-identityN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr33.2%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.3%
Simplified36.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified92.5%
(FPCore (x)
:precision binary64
(if (<= x 1.48)
(+ 1.0 (* (* x x) (+ -0.5 (* (* x x) 0.20833333333333334))))
(if (<= x 1.35e+154)
(/ 2.0 (/ (* x (* x (* x x))) (* x x)))
(/ 2.0 (* x x)))))
double code(double x) {
double tmp;
if (x <= 1.48) {
tmp = 1.0 + ((x * x) * (-0.5 + ((x * x) * 0.20833333333333334)));
} else if (x <= 1.35e+154) {
tmp = 2.0 / ((x * (x * (x * x))) / (x * x));
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.48d0) then
tmp = 1.0d0 + ((x * x) * ((-0.5d0) + ((x * x) * 0.20833333333333334d0)))
else if (x <= 1.35d+154) then
tmp = 2.0d0 / ((x * (x * (x * x))) / (x * x))
else
tmp = 2.0d0 / (x * x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.48) {
tmp = 1.0 + ((x * x) * (-0.5 + ((x * x) * 0.20833333333333334)));
} else if (x <= 1.35e+154) {
tmp = 2.0 / ((x * (x * (x * x))) / (x * x));
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.48: tmp = 1.0 + ((x * x) * (-0.5 + ((x * x) * 0.20833333333333334))) elif x <= 1.35e+154: tmp = 2.0 / ((x * (x * (x * x))) / (x * x)) else: tmp = 2.0 / (x * x) return tmp
function code(x) tmp = 0.0 if (x <= 1.48) tmp = Float64(1.0 + Float64(Float64(x * x) * Float64(-0.5 + Float64(Float64(x * x) * 0.20833333333333334)))); elseif (x <= 1.35e+154) tmp = Float64(2.0 / Float64(Float64(x * Float64(x * Float64(x * x))) / Float64(x * x))); else tmp = Float64(2.0 / Float64(x * x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.48) tmp = 1.0 + ((x * x) * (-0.5 + ((x * x) * 0.20833333333333334))); elseif (x <= 1.35e+154) tmp = 2.0 / ((x * (x * (x * x))) / (x * x)); else tmp = 2.0 / (x * x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.48], N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(-0.5 + N[(N[(x * x), $MachinePrecision] * 0.20833333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.35e+154], N[(2.0 / N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.48:\\
\;\;\;\;1 + \left(x \cdot x\right) \cdot \left(-0.5 + \left(x \cdot x\right) \cdot 0.20833333333333334\right)\\
\mathbf{elif}\;x \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;\frac{2}{\frac{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}{x \cdot x}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{x \cdot x}\\
\end{array}
\end{array}
if x < 1.48Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.3%
Simplified65.3%
if 1.48 < x < 1.35000000000000003e154Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f645.7%
Simplified5.7%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f645.7%
Simplified5.7%
pow2N/A
metadata-evalN/A
metadata-evalN/A
pow-divN/A
pow-powN/A
pow2N/A
pow2N/A
pow2N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6452.7%
Applied egg-rr52.7%
if 1.35000000000000003e154 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x) :precision binary64 (let* ((t_0 (* (* x x) (* x x)))) (/ 2.0 (+ 2.0 (* (* t_0 t_0) 9.259259259259259e-5)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
return 2.0 / (2.0 + ((t_0 * t_0) * 9.259259259259259e-5));
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
t_0 = (x * x) * (x * x)
code = 2.0d0 / (2.0d0 + ((t_0 * t_0) * 9.259259259259259d-5))
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
return 2.0 / (2.0 + ((t_0 * t_0) * 9.259259259259259e-5));
}
def code(x): t_0 = (x * x) * (x * x) return 2.0 / (2.0 + ((t_0 * t_0) * 9.259259259259259e-5))
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) return Float64(2.0 / Float64(2.0 + Float64(Float64(t_0 * t_0) * 9.259259259259259e-5))) end
function tmp = code(x) t_0 = (x * x) * (x * x); tmp = 2.0 / (2.0 + ((t_0 * t_0) * 9.259259259259259e-5)); end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, N[(2.0 / N[(2.0 + N[(N[(t$95$0 * t$95$0), $MachinePrecision] * 9.259259259259259e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\frac{2}{2 + \left(t\_0 \cdot t\_0\right) \cdot 9.259259259259259 \cdot 10^{-5}}
\end{array}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
distribute-rgt-inN/A
*-lft-identityN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr33.2%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.3%
Simplified36.3%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.4%
Simplified91.4%
(FPCore (x)
:precision binary64
(/
2.0
(+
2.0
(*
(* x x)
(+
1.0
(* x (* x (+ 0.08333333333333333 (* x (* x 0.002777777777777778))))))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778))))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * (0.08333333333333333d0 + (x * (x * 0.002777777777777778d0))))))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778))))))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778))))))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(x * Float64(x * 0.002777777777777778))))))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778)))))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(x * N[(x * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + x \cdot \left(x \cdot 0.002777777777777778\right)\right)\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* (* x x) (+ 1.0 (* x (* 0.002777777777777778 (* x (* x x)))))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (0.002777777777777778 * (x * (x * x)))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (0.002777777777777778d0 * (x * (x * x)))))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (0.002777777777777778 * (x * (x * x)))))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (0.002777777777777778 * (x * (x * x)))))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(0.002777777777777778 * Float64(x * Float64(x * x)))))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (0.002777777777777778 * (x * (x * x))))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(0.002777777777777778 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(0.002777777777777778 \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* (* x x) (* 0.002777777777777778 (* x (* x (* x x))))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (0.002777777777777778 * (x * (x * (x * x))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (0.002777777777777778d0 * (x * (x * (x * x))))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (0.002777777777777778 * (x * (x * (x * x))))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (0.002777777777777778 * (x * (x * (x * x))))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(0.002777777777777778 * Float64(x * Float64(x * Float64(x * x))))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (0.002777777777777778 * (x * (x * (x * x)))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(0.002777777777777778 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(0.002777777777777778 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
associate-*r*N/A
flip3-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr54.0%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.3%
Simplified90.3%
Final simplification90.3%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333)))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6487.5%
Simplified87.5%
(FPCore (x) :precision binary64 (/ 1.0 (+ 1.0 (* x (* x (+ 0.5 (* (* x x) 0.041666666666666664)))))))
double code(double x) {
return 1.0 / (1.0 + (x * (x * (0.5 + ((x * x) * 0.041666666666666664)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (1.0d0 + (x * (x * (0.5d0 + ((x * x) * 0.041666666666666664d0)))))
end function
public static double code(double x) {
return 1.0 / (1.0 + (x * (x * (0.5 + ((x * x) * 0.041666666666666664)))));
}
def code(x): return 1.0 / (1.0 + (x * (x * (0.5 + ((x * x) * 0.041666666666666664)))))
function code(x) return Float64(1.0 / Float64(1.0 + Float64(x * Float64(x * Float64(0.5 + Float64(Float64(x * x) * 0.041666666666666664)))))) end
function tmp = code(x) tmp = 1.0 / (1.0 + (x * (x * (0.5 + ((x * x) * 0.041666666666666664))))); end
code[x_] := N[(1.0 / N[(1.0 + N[(x * N[(x * N[(0.5 + N[(N[(x * x), $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{1 + x \cdot \left(x \cdot \left(0.5 + \left(x \cdot x\right) \cdot 0.041666666666666664\right)\right)}
\end{array}
Initial program 100.0%
clear-numN/A
/-lowering-/.f64N/A
cosh-defN/A
cosh-lowering-cosh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.2%
Simplified87.2%
(FPCore (x) :precision binary64 (if (<= x 1.4) 1.0 (/ 2.0 (* x x))))
double code(double x) {
double tmp;
if (x <= 1.4) {
tmp = 1.0;
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.4d0) then
tmp = 1.0d0
else
tmp = 2.0d0 / (x * x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.4) {
tmp = 1.0;
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.4: tmp = 1.0 else: tmp = 2.0 / (x * x) return tmp
function code(x) tmp = 0.0 if (x <= 1.4) tmp = 1.0; else tmp = Float64(2.0 / Float64(x * x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.4) tmp = 1.0; else tmp = 2.0 / (x * x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.4], 1.0, N[(2.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.4:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{x \cdot x}\\
\end{array}
\end{array}
if x < 1.3999999999999999Initial program 100.0%
Taylor expanded in x around 0
Simplified65.5%
if 1.3999999999999999 < x Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6453.6%
Simplified53.6%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6453.6%
Simplified53.6%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* x x))))
double code(double x) {
return 2.0 / (2.0 + (x * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + (x * x))
end function
public static double code(double x) {
return 2.0 / (2.0 + (x * x));
}
def code(x): return 2.0 / (2.0 + (x * x))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(x * x))) end
function tmp = code(x) tmp = 2.0 / (2.0 + (x * x)); end
code[x_] := N[(2.0 / N[(2.0 + N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + x \cdot x}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6476.4%
Simplified76.4%
(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%
Taylor expanded in x around 0
Simplified50.8%
herbie shell --seed 2024147
(FPCore (x)
:name "Hyperbolic secant"
:precision binary64
(/ 2.0 (+ (exp x) (exp (- x)))))