
(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 56.7%
sinh-defN/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
(FPCore (x)
:precision binary64
(let* ((t_0
(*
(* x x)
(+
0.16666666666666666
(*
x
(*
x
(+ 0.008333333333333333 (* (* x x) 0.0001984126984126984))))))))
(if (<= x 5e+44)
(/ (* x (- 1.0 (* t_0 t_0))) (- 1.0 t_0))
(* x (* x (* x (* 0.0001984126984126984 (* x (* x (* x x))))))))))
double code(double x) {
double t_0 = (x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * 0.0001984126984126984)))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0);
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = (x * x) * (0.16666666666666666d0 + (x * (x * (0.008333333333333333d0 + ((x * x) * 0.0001984126984126984d0)))))
if (x <= 5d+44) then
tmp = (x * (1.0d0 - (t_0 * t_0))) / (1.0d0 - t_0)
else
tmp = x * (x * (x * (0.0001984126984126984d0 * (x * (x * (x * x))))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * 0.0001984126984126984)))));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0);
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
def code(x): t_0 = (x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * 0.0001984126984126984))))) tmp = 0 if x <= 5e+44: tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0) else: tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * 0.0001984126984126984)))))) tmp = 0.0 if (x <= 5e+44) tmp = Float64(Float64(x * Float64(1.0 - Float64(t_0 * t_0))) / Float64(1.0 - t_0)); else tmp = Float64(x * Float64(x * Float64(x * Float64(0.0001984126984126984 * Float64(x * Float64(x * Float64(x * x))))))); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * 0.0001984126984126984))))); tmp = 0.0; if (x <= 5e+44) tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0); else tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5e+44], N[(N[(x * N[(1.0 - N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$0), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(0.0001984126984126984 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)\\
\mathbf{if}\;x \leq 5 \cdot 10^{+44}:\\
\;\;\;\;\frac{x \cdot \left(1 - t\_0 \cdot t\_0\right)}{1 - t\_0}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(0.0001984126984126984 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 4.9999999999999996e44Initial program 44.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr67.6%
if 4.9999999999999996e44 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
metadata-evalN/A
metadata-evalN/A
pow-plusN/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
*-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.8%
(FPCore (x)
:precision binary64
(let* ((t_0
(*
(* x x)
(+
0.16666666666666666
(* (* x x) (* (* x x) 0.0001984126984126984))))))
(if (<= x 5e+44)
(/ (* x (- 1.0 (* t_0 t_0))) (- 1.0 t_0))
(* x (* x (* x (* 0.0001984126984126984 (* x (* x (* x x))))))))))
double code(double x) {
double t_0 = (x * x) * (0.16666666666666666 + ((x * x) * ((x * x) * 0.0001984126984126984)));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0);
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = (x * x) * (0.16666666666666666d0 + ((x * x) * ((x * x) * 0.0001984126984126984d0)))
if (x <= 5d+44) then
tmp = (x * (1.0d0 - (t_0 * t_0))) / (1.0d0 - t_0)
else
tmp = x * (x * (x * (0.0001984126984126984d0 * (x * (x * (x * x))))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (0.16666666666666666 + ((x * x) * ((x * x) * 0.0001984126984126984)));
double tmp;
if (x <= 5e+44) {
tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0);
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
def code(x): t_0 = (x * x) * (0.16666666666666666 + ((x * x) * ((x * x) * 0.0001984126984126984))) tmp = 0 if x <= 5e+44: tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0) else: tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(Float64(x * x) * Float64(Float64(x * x) * 0.0001984126984126984)))) tmp = 0.0 if (x <= 5e+44) tmp = Float64(Float64(x * Float64(1.0 - Float64(t_0 * t_0))) / Float64(1.0 - t_0)); else tmp = Float64(x * Float64(x * Float64(x * Float64(0.0001984126984126984 * Float64(x * Float64(x * Float64(x * x))))))); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (0.16666666666666666 + ((x * x) * ((x * x) * 0.0001984126984126984))); tmp = 0.0; if (x <= 5e+44) tmp = (x * (1.0 - (t_0 * t_0))) / (1.0 - t_0); else tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5e+44], N[(N[(x * N[(1.0 - N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$0), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(0.0001984126984126984 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\\
\mathbf{if}\;x \leq 5 \cdot 10^{+44}:\\
\;\;\;\;\frac{x \cdot \left(1 - t\_0 \cdot t\_0\right)}{1 - t\_0}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(0.0001984126984126984 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 4.9999999999999996e44Initial program 44.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.4%
Simplified87.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr67.4%
if 4.9999999999999996e44 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
metadata-evalN/A
metadata-evalN/A
pow-plusN/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
*-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.7%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.008333333333333333 (* (* x x) 0.0001984126984126984)))
(t_1 (* x (* x t_0))))
(if (<= x 1.2e+62)
(*
x
(+
1.0
(/
(* (* x x) (- 0.027777777777777776 (* t_0 (* x (* x t_1)))))
(- 0.16666666666666666 t_1))))
(* x (+ 1.0 (* 0.008333333333333333 (* x (* x (* x x)))))))))
double code(double x) {
double t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984);
double t_1 = x * (x * t_0);
double tmp;
if (x <= 1.2e+62) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * t_1))))) / (0.16666666666666666 - t_1)));
} else {
tmp = x * (1.0 + (0.008333333333333333 * (x * (x * (x * x)))));
}
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 = 0.008333333333333333d0 + ((x * x) * 0.0001984126984126984d0)
t_1 = x * (x * t_0)
if (x <= 1.2d+62) then
tmp = x * (1.0d0 + (((x * x) * (0.027777777777777776d0 - (t_0 * (x * (x * t_1))))) / (0.16666666666666666d0 - t_1)))
else
tmp = x * (1.0d0 + (0.008333333333333333d0 * (x * (x * (x * x)))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984);
double t_1 = x * (x * t_0);
double tmp;
if (x <= 1.2e+62) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * t_1))))) / (0.16666666666666666 - t_1)));
} else {
tmp = x * (1.0 + (0.008333333333333333 * (x * (x * (x * x)))));
}
return tmp;
}
def code(x): t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984) t_1 = x * (x * t_0) tmp = 0 if x <= 1.2e+62: tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * t_1))))) / (0.16666666666666666 - t_1))) else: tmp = x * (1.0 + (0.008333333333333333 * (x * (x * (x * x))))) return tmp
function code(x) t_0 = Float64(0.008333333333333333 + Float64(Float64(x * x) * 0.0001984126984126984)) t_1 = Float64(x * Float64(x * t_0)) tmp = 0.0 if (x <= 1.2e+62) tmp = Float64(x * Float64(1.0 + Float64(Float64(Float64(x * x) * Float64(0.027777777777777776 - Float64(t_0 * Float64(x * Float64(x * t_1))))) / Float64(0.16666666666666666 - t_1)))); else tmp = Float64(x * Float64(1.0 + Float64(0.008333333333333333 * Float64(x * Float64(x * Float64(x * x)))))); end return tmp end
function tmp_2 = code(x) t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984); t_1 = x * (x * t_0); tmp = 0.0; if (x <= 1.2e+62) tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * t_1))))) / (0.16666666666666666 - t_1))); else tmp = x * (1.0 + (0.008333333333333333 * (x * (x * (x * x))))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.2e+62], N[(x * N[(1.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(0.027777777777777776 - N[(t$95$0 * N[(x * N[(x * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.16666666666666666 - t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(1.0 + N[(0.008333333333333333 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.008333333333333333 + \left(x \cdot x\right) \cdot 0.0001984126984126984\\
t_1 := x \cdot \left(x \cdot t\_0\right)\\
\mathbf{if}\;x \leq 1.2 \cdot 10^{+62}:\\
\;\;\;\;x \cdot \left(1 + \frac{\left(x \cdot x\right) \cdot \left(0.027777777777777776 - t\_0 \cdot \left(x \cdot \left(x \cdot t\_1\right)\right)\right)}{0.16666666666666666 - t\_1}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(1 + 0.008333333333333333 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 1.2e62Initial program 46.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.0%
Simplified88.0%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr67.2%
if 1.2e62 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
fma-defineN/A
*-commutativeN/A
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/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 simplification73.8%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.008333333333333333 (* (* x x) 0.0001984126984126984))))
(if (<= x 5e+101)
(*
x
(+
1.0
(/
(*
(* x x)
(- 0.027777777777777776 (* t_0 (* x (* x (* x (* x t_0)))))))
(- 0.16666666666666666 (* (* x x) 0.008333333333333333)))))
(* 0.16666666666666666 (* x (* x x))))))
double code(double x) {
double t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984);
double tmp;
if (x <= 5e+101) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * (x * (x * t_0))))))) / (0.16666666666666666 - ((x * x) * 0.008333333333333333))));
} else {
tmp = 0.16666666666666666 * (x * (x * x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = 0.008333333333333333d0 + ((x * x) * 0.0001984126984126984d0)
if (x <= 5d+101) then
tmp = x * (1.0d0 + (((x * x) * (0.027777777777777776d0 - (t_0 * (x * (x * (x * (x * t_0))))))) / (0.16666666666666666d0 - ((x * x) * 0.008333333333333333d0))))
else
tmp = 0.16666666666666666d0 * (x * (x * x))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984);
double tmp;
if (x <= 5e+101) {
tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * (x * (x * t_0))))))) / (0.16666666666666666 - ((x * x) * 0.008333333333333333))));
} else {
tmp = 0.16666666666666666 * (x * (x * x));
}
return tmp;
}
def code(x): t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984) tmp = 0 if x <= 5e+101: tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * (x * (x * t_0))))))) / (0.16666666666666666 - ((x * x) * 0.008333333333333333)))) else: tmp = 0.16666666666666666 * (x * (x * x)) return tmp
function code(x) t_0 = Float64(0.008333333333333333 + Float64(Float64(x * x) * 0.0001984126984126984)) tmp = 0.0 if (x <= 5e+101) tmp = Float64(x * Float64(1.0 + Float64(Float64(Float64(x * x) * Float64(0.027777777777777776 - Float64(t_0 * Float64(x * Float64(x * Float64(x * Float64(x * t_0))))))) / Float64(0.16666666666666666 - Float64(Float64(x * x) * 0.008333333333333333))))); else tmp = Float64(0.16666666666666666 * Float64(x * Float64(x * x))); end return tmp end
function tmp_2 = code(x) t_0 = 0.008333333333333333 + ((x * x) * 0.0001984126984126984); tmp = 0.0; if (x <= 5e+101) tmp = x * (1.0 + (((x * x) * (0.027777777777777776 - (t_0 * (x * (x * (x * (x * t_0))))))) / (0.16666666666666666 - ((x * x) * 0.008333333333333333)))); else tmp = 0.16666666666666666 * (x * (x * x)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5e+101], N[(x * N[(1.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(0.027777777777777776 - N[(t$95$0 * N[(x * N[(x * N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.16666666666666666 - N[(N[(x * x), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.16666666666666666 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.008333333333333333 + \left(x \cdot x\right) \cdot 0.0001984126984126984\\
\mathbf{if}\;x \leq 5 \cdot 10^{+101}:\\
\;\;\;\;x \cdot \left(1 + \frac{\left(x \cdot x\right) \cdot \left(0.027777777777777776 - t\_0 \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot t\_0\right)\right)\right)\right)\right)}{0.16666666666666666 - \left(x \cdot x\right) \cdot 0.008333333333333333}\right)\\
\mathbf{else}:\\
\;\;\;\;0.16666666666666666 \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < 4.99999999999999989e101Initial program 48.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.6%
Simplified88.6%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr65.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.6%
Simplified78.6%
if 4.99999999999999989e101 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification81.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x (* x x)))))
(if (<= x 1.2e+62)
(*
x
(+
1.0
(/
(*
(* x x)
(- 0.027777777777777776 (* t_0 (* t_0 3.936759889140842e-8))))
(-
0.16666666666666666
(* (* x x) (* (* x x) 0.0001984126984126984))))))
(* x (+ 1.0 (* 0.008333333333333333 t_0))))))
double code(double x) {
double t_0 = x * (x * (x * x));
double tmp;
if (x <= 1.2e+62) {
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 * (1.0 + (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 <= 1.2d+62) 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 * (1.0d0 + (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 <= 1.2e+62) {
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 * (1.0 + (0.008333333333333333 * t_0));
}
return tmp;
}
def code(x): t_0 = x * (x * (x * x)) tmp = 0 if x <= 1.2e+62: 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 * (1.0 + (0.008333333333333333 * t_0)) return tmp
function code(x) t_0 = Float64(x * Float64(x * Float64(x * x))) tmp = 0.0 if (x <= 1.2e+62) 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(Float64(x * x) * Float64(Float64(x * x) * 0.0001984126984126984)))))); else tmp = Float64(x * Float64(1.0 + 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 <= 1.2e+62) 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 * (1.0 + (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, 1.2e+62], 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[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(1.0 + N[(0.008333333333333333 * t$95$0), $MachinePrecision]), $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 1.2 \cdot 10^{+62}:\\
\;\;\;\;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 - \left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot 0.0001984126984126984\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(1 + 0.008333333333333333 \cdot t\_0\right)\\
\end{array}
\end{array}
if x < 1.2e62Initial program 46.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.0%
Simplified88.0%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.8%
Simplified87.8%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr67.0%
if 1.2e62 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
fma-defineN/A
*-commutativeN/A
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/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 simplification73.6%
(FPCore (x)
:precision binary64
(if (<= x 5.5)
(*
x
(+
1.0
(* x (* x (+ 0.16666666666666666 (* x (* x 0.008333333333333333)))))))
(*
x
(*
(+ 0.008333333333333333 (* (* x x) 0.0001984126984126984))
(* x (* x (* x x)))))))
double code(double x) {
double tmp;
if (x <= 5.5) {
tmp = x * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333))))));
} else {
tmp = x * ((0.008333333333333333 + ((x * x) * 0.0001984126984126984)) * (x * (x * (x * x))));
}
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 + (x * (x * 0.008333333333333333d0))))))
else
tmp = x * ((0.008333333333333333d0 + ((x * x) * 0.0001984126984126984d0)) * (x * (x * (x * x))))
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 + (x * (x * 0.008333333333333333))))));
} else {
tmp = x * ((0.008333333333333333 + ((x * x) * 0.0001984126984126984)) * (x * (x * (x * x))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 5.5: tmp = x * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333)))))) else: tmp = x * ((0.008333333333333333 + ((x * x) * 0.0001984126984126984)) * (x * (x * (x * x)))) return tmp
function code(x) tmp = 0.0 if (x <= 5.5) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(x * Float64(x * 0.008333333333333333))))))); else tmp = Float64(x * Float64(Float64(0.008333333333333333 + Float64(Float64(x * x) * 0.0001984126984126984)) * Float64(x * Float64(x * Float64(x * x))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 5.5) tmp = x * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * 0.008333333333333333)))))); else tmp = x * ((0.008333333333333333 + ((x * x) * 0.0001984126984126984)) * (x * (x * (x * x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 5.5], N[(x * N[(1.0 + N[(x * N[(x * N[(0.16666666666666666 + N[(x * N[(x * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * x), $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 \left(0.16666666666666666 + x \cdot \left(x \cdot 0.008333333333333333\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(0.008333333333333333 + \left(x \cdot x\right) \cdot 0.0001984126984126984\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 5.5Initial program 40.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
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6492.7%
Simplified92.7%
if 5.5 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.3%
Simplified81.3%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
Taylor expanded in x around inf
Simplified81.3%
(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(Float64(x * x) * Float64(0.008333333333333333 + Float64(Float64(x * 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[(N[(x * x), $MachinePrecision] * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot 0.0001984126984126984\right)\right)\right)
\end{array}
Initial program 56.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
(FPCore (x)
:precision binary64
(if (<= x 7.5)
(*
x
(+
1.0
(* x (* x (+ 0.16666666666666666 (* x (* x 0.008333333333333333)))))))
(* x (* x (* x (* 0.0001984126984126984 (* x (* x (* x x)))))))))
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 * (0.0001984126984126984 * (x * (x * (x * x))))));
}
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 * (0.0001984126984126984d0 * (x * (x * (x * x))))))
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 * (0.0001984126984126984 * (x * (x * (x * x))))));
}
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 * (0.0001984126984126984 * (x * (x * (x * x)))))) 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(x * Float64(x * 0.008333333333333333))))))); else tmp = Float64(x * Float64(x * Float64(x * Float64(0.0001984126984126984 * Float64(x * Float64(x * Float64(x * x))))))); 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 * (0.0001984126984126984 * (x * (x * (x * x)))))); 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[(x * N[(x * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * N[(0.0001984126984126984 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $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 + x \cdot \left(x \cdot 0.008333333333333333\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(0.0001984126984126984 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 7.5Initial program 40.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
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6492.7%
Simplified92.7%
if 7.5 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.3%
Simplified81.3%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
Taylor expanded in x around inf
metadata-evalN/A
metadata-evalN/A
pow-plusN/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
*-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-*.f6481.3%
Simplified81.3%
(FPCore (x) :precision binary64 (if (<= x 5.5) (/ (+ (* x (* (* x x) 0.3333333333333333)) (* x 2.0)) 2.0) (* x (* x (* x (* 0.0001984126984126984 (* x (* x (* x x)))))))))
double code(double x) {
double tmp;
if (x <= 5.5) {
tmp = ((x * ((x * x) * 0.3333333333333333)) + (x * 2.0)) / 2.0;
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 5.5d0) then
tmp = ((x * ((x * x) * 0.3333333333333333d0)) + (x * 2.0d0)) / 2.0d0
else
tmp = x * (x * (x * (0.0001984126984126984d0 * (x * (x * (x * x))))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 5.5) {
tmp = ((x * ((x * x) * 0.3333333333333333)) + (x * 2.0)) / 2.0;
} else {
tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x))))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 5.5: tmp = ((x * ((x * x) * 0.3333333333333333)) + (x * 2.0)) / 2.0 else: tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))) return tmp
function code(x) tmp = 0.0 if (x <= 5.5) tmp = Float64(Float64(Float64(x * Float64(Float64(x * x) * 0.3333333333333333)) + Float64(x * 2.0)) / 2.0); else tmp = Float64(x * Float64(x * Float64(x * Float64(0.0001984126984126984 * Float64(x * Float64(x * Float64(x * x))))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 5.5) tmp = ((x * ((x * x) * 0.3333333333333333)) + (x * 2.0)) / 2.0; else tmp = x * (x * (x * (0.0001984126984126984 * (x * (x * (x * x)))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 5.5], N[(N[(N[(x * N[(N[(x * x), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]), $MachinePrecision] + N[(x * 2.0), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(x * N[(x * N[(x * N[(0.0001984126984126984 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.5:\\
\;\;\;\;\frac{x \cdot \left(\left(x \cdot x\right) \cdot 0.3333333333333333\right) + x \cdot 2}{2}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot \left(0.0001984126984126984 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < 5.5Initial program 40.1%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.0%
Simplified87.0%
+-commutativeN/A
distribute-rgt-inN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6487.0%
Applied egg-rr87.0%
if 5.5 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.3%
Simplified81.3%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6481.3%
Applied egg-rr81.3%
Taylor expanded in x around inf
metadata-evalN/A
metadata-evalN/A
pow-plusN/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
*-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-*.f6481.3%
Simplified81.3%
(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 56.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.2%
Simplified90.2%
(FPCore (x)
:precision binary64
(*
x
(+
1.0
(*
(* x x)
(* x (* x (+ 0.008333333333333333 (* x (* x 0.0001984126984126984)))))))))
double code(double x) {
return x * (1.0 + ((x * x) * (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984)))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + ((x * x) * (x * (x * (0.008333333333333333d0 + (x * (x * 0.0001984126984126984d0)))))))
end function
public static double code(double x) {
return x * (1.0 + ((x * x) * (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984)))))));
}
def code(x): return x * (1.0 + ((x * x) * (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984)))))))
function code(x) return Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(x * Float64(x * Float64(0.008333333333333333 + Float64(x * Float64(x * 0.0001984126984126984)))))))) end
function tmp = code(x) tmp = x * (1.0 + ((x * x) * (x * (x * (0.008333333333333333 + (x * (x * 0.0001984126984126984))))))); end
code[x_] := N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(0.008333333333333333 + N[(x * N[(x * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(x \cdot \left(x \cdot \left(0.008333333333333333 + x \cdot \left(x \cdot 0.0001984126984126984\right)\right)\right)\right)\right)
\end{array}
Initial program 56.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Applied egg-rr90.4%
Taylor expanded in x around inf
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
fma-defineN/A
*-commutativeN/A
Simplified90.0%
(FPCore (x) :precision binary64 (* x (+ 1.0 (* (* x x) (* x (* x (* x (* x 0.0001984126984126984))))))))
double code(double x) {
return x * (1.0 + ((x * x) * (x * (x * (x * (x * 0.0001984126984126984))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + ((x * x) * (x * (x * (x * (x * 0.0001984126984126984d0))))))
end function
public static double code(double x) {
return x * (1.0 + ((x * x) * (x * (x * (x * (x * 0.0001984126984126984))))));
}
def code(x): return x * (1.0 + ((x * x) * (x * (x * (x * (x * 0.0001984126984126984))))))
function code(x) return Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(x * Float64(x * Float64(x * Float64(x * 0.0001984126984126984))))))) end
function tmp = code(x) tmp = x * (1.0 + ((x * x) * (x * (x * (x * (x * 0.0001984126984126984)))))); end
code[x_] := N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * N[(x * N[(x * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot 0.0001984126984126984\right)\right)\right)\right)\right)
\end{array}
Initial program 56.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Applied egg-rr90.4%
Taylor expanded in x around inf
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
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6490.0%
Simplified90.0%
(FPCore (x) :precision binary64 (* x (+ 1.0 (* 0.008333333333333333 (* x (* x (* x x)))))))
double code(double x) {
return x * (1.0 + (0.008333333333333333 * (x * (x * (x * x)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (1.0d0 + (0.008333333333333333d0 * (x * (x * (x * x)))))
end function
public static double code(double x) {
return x * (1.0 + (0.008333333333333333 * (x * (x * (x * x)))));
}
def code(x): return x * (1.0 + (0.008333333333333333 * (x * (x * (x * x)))))
function code(x) return Float64(x * Float64(1.0 + Float64(0.008333333333333333 * Float64(x * Float64(x * Float64(x * x)))))) end
function tmp = code(x) tmp = x * (1.0 + (0.008333333333333333 * (x * (x * (x * x))))); end
code[x_] := N[(x * N[(1.0 + N[(0.008333333333333333 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + 0.008333333333333333 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)
\end{array}
Initial program 56.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.4%
Simplified90.4%
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.4%
Applied egg-rr90.4%
Taylor expanded in x around inf
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
fma-defineN/A
*-commutativeN/A
Simplified90.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/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-*.f6487.0%
Simplified87.0%
(FPCore (x) :precision binary64 (if (<= x 2.5) x (* 0.16666666666666666 (* x (* x x)))))
double code(double x) {
double tmp;
if (x <= 2.5) {
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.5d0) 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.5) {
tmp = x;
} else {
tmp = 0.16666666666666666 * (x * (x * x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.5: tmp = x else: tmp = 0.16666666666666666 * (x * (x * x)) return tmp
function code(x) tmp = 0.0 if (x <= 2.5) 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.5) tmp = x; else tmp = 0.16666666666666666 * (x * (x * x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.5], x, N[(0.16666666666666666 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.5:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;0.16666666666666666 \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < 2.5Initial program 40.1%
Taylor expanded in x around 0
Simplified66.6%
if 2.5 < x Initial program 100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.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(Float64(x * 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[(N[(x * x), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(x \cdot x\right) \cdot 0.16666666666666666\right)
\end{array}
Initial program 56.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-*.f6479.2%
Simplified79.2%
Final simplification79.2%
(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 56.7%
Taylor expanded in x around 0
Simplified49.6%
herbie shell --seed 2024150
(FPCore (x)
:name "Hyperbolic sine"
:precision binary64
(/ (- (exp x) (exp (- x))) 2.0))