
(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 17 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 58.7%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) 0.0001984126984126984))
(t_1
(*
x
(*
x
(+ 0.16666666666666666 (* (* x x) (+ 0.008333333333333333 t_0)))))))
(if (<= x 5e+44)
(/ (* x (- 1.0 (* t_1 t_1))) (- 1.0 t_1))
(* x (* x (* x (* x (* x t_0))))))))
double code(double x) {
double t_0 = (x * x) * 0.0001984126984126984;
double t_1 = x * (x * (0.16666666666666666 + ((x * x) * (0.008333333333333333 + t_0))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_1 * t_1))) / (1.0 - t_1);
} else {
tmp = x * (x * (x * (x * (x * t_0))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x * x) * 0.0001984126984126984d0
t_1 = x * (x * (0.16666666666666666d0 + ((x * x) * (0.008333333333333333d0 + t_0))))
if (x <= 5d+44) then
tmp = (x * (1.0d0 - (t_1 * t_1))) / (1.0d0 - t_1)
else
tmp = x * (x * (x * (x * (x * t_0))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * 0.0001984126984126984;
double t_1 = x * (x * (0.16666666666666666 + ((x * x) * (0.008333333333333333 + t_0))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_1 * t_1))) / (1.0 - t_1);
} else {
tmp = x * (x * (x * (x * (x * t_0))));
}
return tmp;
}
def code(x): t_0 = (x * x) * 0.0001984126984126984 t_1 = x * (x * (0.16666666666666666 + ((x * x) * (0.008333333333333333 + t_0)))) tmp = 0 if x <= 5e+44: tmp = (x * (1.0 - (t_1 * t_1))) / (1.0 - t_1) else: tmp = x * (x * (x * (x * (x * t_0)))) return tmp
function code(x) t_0 = Float64(Float64(x * x) * 0.0001984126984126984) t_1 = Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * Float64(0.008333333333333333 + t_0))))) tmp = 0.0 if (x <= 5e+44) tmp = Float64(Float64(x * Float64(1.0 - Float64(t_1 * t_1))) / Float64(1.0 - t_1)); else tmp = Float64(x * Float64(x * Float64(x * Float64(x * Float64(x * t_0))))); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * 0.0001984126984126984; t_1 = x * (x * (0.16666666666666666 + ((x * x) * (0.008333333333333333 + t_0)))); tmp = 0.0; if (x <= 5e+44) tmp = (x * (1.0 - (t_1 * t_1))) / (1.0 - t_1); else tmp = x * (x * (x * (x * (x * t_0)))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.008333333333333333 + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5e+44], N[(N[(x * N[(1.0 - N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$1), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot 0.0001984126984126984\\
t_1 := x \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.008333333333333333 + t\_0\right)\right)\right)\\
\mathbf{if}\;x \leq 5 \cdot 10^{+44}:\\
\;\;\;\;\frac{x \cdot \left(1 - t\_1 \cdot t\_1\right)}{1 - t\_1}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot t\_0\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 4.9999999999999996e44Initial program 46.1%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6491.4%
Simplified91.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr66.2%
if 4.9999999999999996e44 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification74.1%
(FPCore (x)
:precision binary64
(let* ((t_0
(+
0.16666666666666666
(* x (* x (* x (* x 0.0001984126984126984)))))))
(if (<= x 5e+44)
(/
(* x (- 1.0 (* (* x (* x (* x x))) (* t_0 t_0))))
(- 1.0 (* (* x x) t_0)))
(* x (* x (* x (* x (* x (* (* x x) 0.0001984126984126984)))))))))
double code(double x) {
double t_0 = 0.16666666666666666 + (x * (x * (x * (x * 0.0001984126984126984))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - ((x * (x * (x * x))) * (t_0 * t_0)))) / (1.0 - ((x * x) * t_0));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = 0.16666666666666666d0 + (x * (x * (x * (x * 0.0001984126984126984d0))))
if (x <= 5d+44) then
tmp = (x * (1.0d0 - ((x * (x * (x * x))) * (t_0 * t_0)))) / (1.0d0 - ((x * x) * t_0))
else
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984d0)))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.16666666666666666 + (x * (x * (x * (x * 0.0001984126984126984))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - ((x * (x * (x * x))) * (t_0 * t_0)))) / (1.0 - ((x * x) * t_0));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
def code(x): t_0 = 0.16666666666666666 + (x * (x * (x * (x * 0.0001984126984126984)))) tmp = 0 if x <= 5e+44: tmp = (x * (1.0 - ((x * (x * (x * x))) * (t_0 * t_0)))) / (1.0 - ((x * x) * t_0)) else: tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))) return tmp
function code(x) t_0 = Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(x * Float64(x * 0.0001984126984126984))))) tmp = 0.0 if (x <= 5e+44) tmp = Float64(Float64(x * Float64(1.0 - Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(t_0 * t_0)))) / Float64(1.0 - Float64(Float64(x * x) * t_0))); else tmp = Float64(x * Float64(x * Float64(x * Float64(x * Float64(x * Float64(Float64(x * x) * 0.0001984126984126984)))))); end return tmp end
function tmp_2 = code(x) t_0 = 0.16666666666666666 + (x * (x * (x * (x * 0.0001984126984126984)))); tmp = 0.0; if (x <= 5e+44) tmp = (x * (1.0 - ((x * (x * (x * x))) * (t_0 * t_0)))) / (1.0 - ((x * x) * t_0)); else tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.16666666666666666 + N[(x * N[(x * N[(x * N[(x * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5e+44], N[(N[(x * N[(1.0 - N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[(x * x), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.16666666666666666 + x \cdot \left(x \cdot \left(x \cdot \left(x \cdot 0.0001984126984126984\right)\right)\right)\\
\mathbf{if}\;x \leq 5 \cdot 10^{+44}:\\
\;\;\;\;\frac{x \cdot \left(1 - \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(t\_0 \cdot t\_0\right)\right)}{1 - \left(x \cdot x\right) \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 4.9999999999999996e44Initial program 46.1%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6491.4%
Simplified91.4%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.4%
Simplified91.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr66.1%
if 4.9999999999999996e44 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification74.1%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x (* x x)))))
(if (<= x 4e+61)
(*
x
(+
1.0
(/
(*
(* x x)
(- 0.027777777777777776 (* t_0 (* t_0 3.936759889140842e-8))))
(-
0.16666666666666666
(* x (* x (* x (* x 0.0001984126984126984))))))))
(* x (* 0.008333333333333333 t_0)))))
double code(double x) {
double t_0 = x * (x * (x * x));
double tmp;
if (x <= 4e+61) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (t_0 * 3.936759889140842e-8)))) / (0.16666666666666666 - (x * (x * (x * (x * 0.0001984126984126984)))))));
} else {
tmp = x * (0.008333333333333333 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x * (x * (x * x))
if (x <= 4d+61) then
tmp = x * (1.0d0 + (((x * x) * (0.027777777777777776d0 - (t_0 * (t_0 * 3.936759889140842d-8)))) / (0.16666666666666666d0 - (x * (x * (x * (x * 0.0001984126984126984d0)))))))
else
tmp = x * (0.008333333333333333d0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * (x * x));
double tmp;
if (x <= 4e+61) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (t_0 * 3.936759889140842e-8)))) / (0.16666666666666666 - (x * (x * (x * (x * 0.0001984126984126984)))))));
} else {
tmp = x * (0.008333333333333333 * t_0);
}
return tmp;
}
def code(x): t_0 = x * (x * (x * x)) tmp = 0 if x <= 4e+61: tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (t_0 * 3.936759889140842e-8)))) / (0.16666666666666666 - (x * (x * (x * (x * 0.0001984126984126984))))))) else: tmp = x * (0.008333333333333333 * t_0) return tmp
function code(x) t_0 = Float64(x * Float64(x * Float64(x * x))) tmp = 0.0 if (x <= 4e+61) tmp = Float64(x * Float64(1.0 + Float64(Float64(Float64(x * x) * Float64(0.027777777777777776 - Float64(t_0 * Float64(t_0 * 3.936759889140842e-8)))) / Float64(0.16666666666666666 - Float64(x * Float64(x * Float64(x * Float64(x * 0.0001984126984126984)))))))); else tmp = Float64(x * Float64(0.008333333333333333 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * (x * x)); tmp = 0.0; if (x <= 4e+61) tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (t_0 * 3.936759889140842e-8)))) / (0.16666666666666666 - (x * (x * (x * (x * 0.0001984126984126984))))))); else tmp = x * (0.008333333333333333 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 4e+61], N[(x * N[(1.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(0.027777777777777776 - N[(t$95$0 * N[(t$95$0 * 3.936759889140842e-8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.16666666666666666 - N[(x * N[(x * N[(x * N[(x * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.008333333333333333 * t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\mathbf{if}\;x \leq 4 \cdot 10^{+61}:\\
\;\;\;\;x \cdot \left(1 + \frac{\left(x \cdot x\right) \cdot \left(0.027777777777777776 - t\_0 \cdot \left(t\_0 \cdot 3.936759889140842 \cdot 10^{-8}\right)\right)}{0.16666666666666666 - x \cdot \left(x \cdot \left(x \cdot \left(x \cdot 0.0001984126984126984\right)\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0.008333333333333333 \cdot t\_0\right)\\
\end{array}
\end{array}
if x < 3.9999999999999998e61Initial program 46.9%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6491.6%
Simplified91.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.5%
Simplified91.5%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr67.2%
if 3.9999999999999998e61 < x Initial program 100.0%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
associate-*l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
pow-sqrN/A
metadata-evalN/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-*.f64100.0%
Simplified100.0%
Final simplification74.5%
(FPCore (x)
:precision binary64
(if (<= x 6.6)
(/
x
(/
1.0
(+
1.0
(* x (* x (+ 0.16666666666666666 (* x (* x 0.008333333333333333))))))))
(*
x
(*
(* x x)
(+ 0.16666666666666666 (* x (* x (* (* x x) 0.0001984126984126984))))))))
double code(double x) {
double tmp;
if (x <= 6.6) {
tmp = x / (1.0 / (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333)))))));
} else {
tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 6.6d0) then
tmp = x / (1.0d0 / (1.0d0 + (x * (x * (0.16666666666666666d0 + (x * (x * 0.008333333333333333d0)))))))
else
tmp = x * ((x * x) * (0.16666666666666666d0 + (x * (x * ((x * x) * 0.0001984126984126984d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 6.6) {
tmp = x / (1.0 / (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333)))))));
} else {
tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 6.6: tmp = x / (1.0 / (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333))))))) else: tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984))))) return tmp
function code(x) tmp = 0.0 if (x <= 6.6) tmp = Float64(x / Float64(1.0 / Float64(1.0 + Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(x * Float64(x * 0.008333333333333333)))))))); else tmp = Float64(x * Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(Float64(x * x) * 0.0001984126984126984)))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 6.6) tmp = x / (1.0 / (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333))))))); else tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 6.6], N[(x / N[(1.0 / N[(1.0 + N[(x * N[(x * N[(0.16666666666666666 + N[(x * N[(x * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6.6:\\
\;\;\;\;\frac{x}{\frac{1}{1 + x \cdot \left(x \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot 0.008333333333333333\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 6.5999999999999996Initial program 43.5%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.5%
Simplified92.5%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr92.5%
if 6.5999999999999996 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Simplified87.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
Taylor expanded in x around inf
distribute-lft-inN/A
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
associate-*l/N/A
Simplified87.6%
(FPCore (x)
:precision binary64
(if (<= x 6.6)
(*
x
(+
1.0
(* x (* x (+ 0.16666666666666666 (* (* x x) 0.008333333333333333))))))
(*
x
(*
(* x x)
(+ 0.16666666666666666 (* x (* x (* (* x x) 0.0001984126984126984))))))))
double code(double x) {
double tmp;
if (x <= 6.6) {
tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))));
} else {
tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 6.6d0) then
tmp = x * (1.0d0 + (x * (x * (0.16666666666666666d0 + ((x * x) * 0.008333333333333333d0)))))
else
tmp = x * ((x * x) * (0.16666666666666666d0 + (x * (x * ((x * x) * 0.0001984126984126984d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 6.6) {
tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))));
} else {
tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 6.6: tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))))) else: tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984))))) return tmp
function code(x) tmp = 0.0 if (x <= 6.6) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * 0.008333333333333333)))))); else tmp = Float64(x * Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(Float64(x * x) * 0.0001984126984126984)))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 6.6) tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))))); else tmp = x * ((x * x) * (0.16666666666666666 + (x * (x * ((x * x) * 0.0001984126984126984))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 6.6], N[(x * N[(1.0 + N[(x * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6.6:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot 0.008333333333333333\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 6.5999999999999996Initial program 43.5%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.5%
Simplified92.5%
if 6.5999999999999996 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Simplified87.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
Taylor expanded in x around inf
distribute-lft-inN/A
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
associate-*l/N/A
Simplified87.6%
(FPCore (x)
:precision binary64
(*
x
(+
1.0
(*
(* x x)
(+
0.16666666666666666
(*
x
(* x (+ 0.008333333333333333 (* x (* x 0.0001984126984126984))))))))))
double code(double x) {
return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984))))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + ((x * x) * (0.16666666666666666d0 + (x * (x * (0.008333333333333333d0 + (x * (x * 0.0001984126984126984d0))))))))
end function
public static double code(double x) {
return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984))))))));
}
def code(x): return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984))))))))
function code(x) return Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(0.008333333333333333 + Float64(x * Float64(x * 0.0001984126984126984))))))))) end
function tmp = code(x) tmp = x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984)))))))); end
code[x_] := N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * N[(x * N[(0.008333333333333333 + N[(x * N[(x * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot \left(0.008333333333333333 + x \cdot \left(x \cdot 0.0001984126984126984\right)\right)\right)\right)\right)
\end{array}
Initial program 58.7%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.5%
Simplified93.5%
(FPCore (x)
:precision binary64
(if (<= x 7.5)
(*
x
(+
1.0
(* x (* x (+ 0.16666666666666666 (* (* x x) 0.008333333333333333))))))
(* x (* x (* x (* x (* x (* (* x x) 0.0001984126984126984))))))))
double code(double x) {
double tmp;
if (x <= 7.5) {
tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 7.5d0) then
tmp = x * (1.0d0 + (x * (x * (0.16666666666666666d0 + ((x * x) * 0.008333333333333333d0)))))
else
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 7.5) {
tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 7.5: tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))))) else: tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))) return tmp
function code(x) tmp = 0.0 if (x <= 7.5) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * 0.008333333333333333)))))); else tmp = Float64(x * Float64(x * Float64(x * Float64(x * Float64(x * Float64(Float64(x * x) * 0.0001984126984126984)))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 7.5) tmp = x * (1.0 + (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))))); else tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 7.5], N[(x * N[(1.0 + N[(x * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 7.5:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot 0.008333333333333333\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 7.5Initial program 43.5%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.5%
Simplified92.5%
if 7.5 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Simplified87.6%
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Applied egg-rr87.6%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
(FPCore (x) :precision binary64 (if (<= x 5.5) (* x (+ 1.0 (* x (* x 0.16666666666666666)))) (* x (* x (* x (* x (* x (* (* x x) 0.0001984126984126984))))))))
double code(double x) {
double tmp;
if (x <= 5.5) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 5.5d0) then
tmp = x * (1.0d0 + (x * (x * 0.16666666666666666d0)))
else
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 5.5) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984)))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 5.5: tmp = x * (1.0 + (x * (x * 0.16666666666666666))) else: tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))) return tmp
function code(x) tmp = 0.0 if (x <= 5.5) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * 0.16666666666666666)))); else tmp = Float64(x * Float64(x * Float64(x * Float64(x * Float64(x * Float64(Float64(x * x) * 0.0001984126984126984)))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 5.5) tmp = x * (1.0 + (x * (x * 0.16666666666666666))); else tmp = x * (x * (x * (x * (x * ((x * x) * 0.0001984126984126984))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 5.5], N[(x * N[(1.0 + N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.5:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot 0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 5.5Initial program 43.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6485.6%
Simplified85.6%
if 5.5 < x Initial program 100.0%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Simplified87.6%
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.6%
Applied egg-rr87.6%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
(FPCore (x)
:precision binary64
(*
x
(+
1.0
(*
(* x x)
(+ 0.16666666666666666 (* x (* 0.0001984126984126984 (* x (* x x)))))))))
double code(double x) {
return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (0.0001984126984126984 * (x * (x * x)))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + ((x * x) * (0.16666666666666666d0 + (x * (0.0001984126984126984d0 * (x * (x * x)))))))
end function
public static double code(double x) {
return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (0.0001984126984126984 * (x * (x * x)))))));
}
def code(x): return x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (0.0001984126984126984 * (x * (x * x)))))))
function code(x) return Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * Float64(0.0001984126984126984 * Float64(x * Float64(x * x)))))))) end
function tmp = code(x) tmp = x * (1.0 + ((x * x) * (0.16666666666666666 + (x * (0.0001984126984126984 * (x * (x * x))))))); end
code[x_] := N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * N[(0.0001984126984126984 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot \left(0.0001984126984126984 \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\right)
\end{array}
Initial program 58.7%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-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
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.5%
Simplified93.5%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.4%
Simplified93.4%
(FPCore (x) :precision binary64 (if (<= x 3.3) (* x (+ 1.0 (* x (* x 0.16666666666666666)))) (* x (* x (* x (+ 0.16666666666666666 (* (* x x) 0.008333333333333333)))))))
double code(double x) {
double tmp;
if (x <= 3.3) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 3.3d0) then
tmp = x * (1.0d0 + (x * (x * 0.16666666666666666d0)))
else
tmp = x * (x * (x * (0.16666666666666666d0 + ((x * x) * 0.008333333333333333d0))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 3.3) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 3.3: tmp = x * (1.0 + (x * (x * 0.16666666666666666))) else: tmp = x * (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))) return tmp
function code(x) tmp = 0.0 if (x <= 3.3) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * 0.16666666666666666)))); else tmp = Float64(x * Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * 0.008333333333333333))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 3.3) tmp = x * (1.0 + (x * (x * 0.16666666666666666))); else tmp = x * (x * (x * (0.16666666666666666 + ((x * x) * 0.008333333333333333)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 3.3], N[(x * N[(1.0 + N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.3:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot 0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot 0.008333333333333333\right)\right)\right)\\
\end{array}
\end{array}
if x < 3.2999999999999998Initial program 43.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6485.6%
Simplified85.6%
if 3.2999999999999998 < x Initial program 100.0%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.5%
Simplified83.5%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
associate-*r*N/A
cube-multN/A
unpow2N/A
associate-*l*N/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
associate-*l*N/A
*-lowering-*.f64N/A
Simplified83.5%
(FPCore (x) :precision binary64 (if (<= x 5.0) (* x (+ 1.0 (* x (* x 0.16666666666666666)))) (* x (* 0.008333333333333333 (* x (* x (* x x)))))))
double code(double x) {
double tmp;
if (x <= 5.0) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (0.008333333333333333 * (x * (x * (x * x))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 5.0d0) then
tmp = x * (1.0d0 + (x * (x * 0.16666666666666666d0)))
else
tmp = x * (0.008333333333333333d0 * (x * (x * (x * x))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 5.0) {
tmp = x * (1.0 + (x * (x * 0.16666666666666666)));
} else {
tmp = x * (0.008333333333333333 * (x * (x * (x * x))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 5.0: tmp = x * (1.0 + (x * (x * 0.16666666666666666))) else: tmp = x * (0.008333333333333333 * (x * (x * (x * x)))) return tmp
function code(x) tmp = 0.0 if (x <= 5.0) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * 0.16666666666666666)))); else tmp = Float64(x * Float64(0.008333333333333333 * Float64(x * Float64(x * Float64(x * x))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 5.0) tmp = x * (1.0 + (x * (x * 0.16666666666666666))); else tmp = x * (0.008333333333333333 * (x * (x * (x * x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 5.0], N[(x * N[(1.0 + N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.008333333333333333 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot 0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0.008333333333333333 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 5Initial program 43.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6485.6%
Simplified85.6%
if 5 < x Initial program 100.0%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.5%
Simplified83.5%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
associate-*l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
pow-sqrN/A
metadata-evalN/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-*.f6483.5%
Simplified83.5%
(FPCore (x) :precision binary64 (* x (+ 1.0 (* x (* 0.008333333333333333 (* x (* x x)))))))
double code(double x) {
return x * (1.0 + (x * (0.008333333333333333 * (x * (x * x)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + (x * (0.008333333333333333d0 * (x * (x * x)))))
end function
public static double code(double x) {
return x * (1.0 + (x * (0.008333333333333333 * (x * (x * x)))));
}
def code(x): return x * (1.0 + (x * (0.008333333333333333 * (x * (x * x)))))
function code(x) return Float64(x * Float64(1.0 + Float64(x * Float64(0.008333333333333333 * Float64(x * Float64(x * x)))))) end
function tmp = code(x) tmp = x * (1.0 + (x * (0.008333333333333333 * (x * (x * x))))); end
code[x_] := N[(x * N[(1.0 + N[(x * N[(0.008333333333333333 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + x \cdot \left(0.008333333333333333 \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)
\end{array}
Initial program 58.7%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6489.9%
Simplified89.9%
(FPCore (x) :precision binary64 (if (<= x 2.4) x (* 0.16666666666666666 (* x (* x x)))))
double code(double x) {
double tmp;
if (x <= 2.4) {
tmp = x;
} else {
tmp = 0.16666666666666666 * (x * (x * x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 2.4d0) then
tmp = x
else
tmp = 0.16666666666666666d0 * (x * (x * x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 2.4) {
tmp = x;
} else {
tmp = 0.16666666666666666 * (x * (x * x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.4: tmp = x else: tmp = 0.16666666666666666 * (x * (x * x)) return tmp
function code(x) tmp = 0.0 if (x <= 2.4) tmp = x; else tmp = Float64(0.16666666666666666 * Float64(x * Float64(x * x))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.4) tmp = x; else tmp = 0.16666666666666666 * (x * (x * x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.4], x, N[(0.16666666666666666 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.4:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;0.16666666666666666 \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < 2.39999999999999991Initial program 43.5%
Taylor expanded in x around 0
Simplified64.4%
if 2.39999999999999991 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
(FPCore (x) :precision binary64 (+ x (* 0.16666666666666666 (* x (* x x)))))
double code(double x) {
return x + (0.16666666666666666 * (x * (x * x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x + (0.16666666666666666d0 * (x * (x * x)))
end function
public static double code(double x) {
return x + (0.16666666666666666 * (x * (x * x)));
}
def code(x): return x + (0.16666666666666666 * (x * (x * x)))
function code(x) return Float64(x + Float64(0.16666666666666666 * Float64(x * Float64(x * x)))) end
function tmp = code(x) tmp = x + (0.16666666666666666 * (x * (x * x))); end
code[x_] := N[(x + N[(0.16666666666666666 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + 0.16666666666666666 \cdot \left(x \cdot \left(x \cdot x\right)\right)
\end{array}
Initial program 58.7%
Taylor expanded in x around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.1%
Simplified90.1%
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.1%
Applied egg-rr90.1%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.7%
Simplified81.7%
Final simplification81.7%
(FPCore (x) :precision binary64 (* x (+ 1.0 (* x (* x 0.16666666666666666)))))
double code(double x) {
return x * (1.0 + (x * (x * 0.16666666666666666)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + (x * (x * 0.16666666666666666d0)))
end function
public static double code(double x) {
return x * (1.0 + (x * (x * 0.16666666666666666)));
}
def code(x): return x * (1.0 + (x * (x * 0.16666666666666666)))
function code(x) return Float64(x * Float64(1.0 + Float64(x * Float64(x * 0.16666666666666666)))) end
function tmp = code(x) tmp = x * (1.0 + (x * (x * 0.16666666666666666))); end
code[x_] := N[(x * N[(1.0 + N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + x \cdot \left(x \cdot 0.16666666666666666\right)\right)
\end{array}
Initial program 58.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6481.3%
Simplified81.3%
(FPCore (x) :precision binary64 x)
double code(double x) {
return x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x
end function
public static double code(double x) {
return x;
}
def code(x): return x
function code(x) return x end
function tmp = code(x) tmp = x; end
code[x_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 58.7%
Taylor expanded in x around 0
Simplified48.5%
herbie shell --seed 2024139
(FPCore (x)
:name "Hyperbolic sine"
:precision binary64
(/ (- (exp x) (exp (- x))) 2.0))