
(FPCore (x) :precision binary64 (/ 2.0 (+ (exp x) (exp (- x)))))
double code(double x) {
return 2.0 / (exp(x) + exp(-x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (exp(x) + exp(-x))
end function
public static double code(double x) {
return 2.0 / (Math.exp(x) + Math.exp(-x));
}
def code(x): return 2.0 / (math.exp(x) + math.exp(-x))
function code(x) return Float64(2.0 / Float64(exp(x) + exp(Float64(-x)))) end
function tmp = code(x) tmp = 2.0 / (exp(x) + exp(-x)); end
code[x_] := N[(2.0 / N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{e^{x} + e^{-x}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ 2.0 (+ (exp x) (exp (- x)))))
double code(double x) {
return 2.0 / (exp(x) + exp(-x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (exp(x) + exp(-x))
end function
public static double code(double x) {
return 2.0 / (Math.exp(x) + Math.exp(-x));
}
def code(x): return 2.0 / (math.exp(x) + math.exp(-x))
function code(x) return Float64(2.0 / Float64(exp(x) + exp(Float64(-x)))) end
function tmp = code(x) tmp = 2.0 / (exp(x) + exp(-x)); end
code[x_] := N[(2.0 / N[(N[Exp[x], $MachinePrecision] + N[Exp[(-x)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{e^{x} + e^{-x}}
\end{array}
(FPCore (x) :precision binary64 (/ 1.0 (cosh x)))
double code(double x) {
return 1.0 / cosh(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / cosh(x)
end function
public static double code(double x) {
return 1.0 / Math.cosh(x);
}
def code(x): return 1.0 / math.cosh(x)
function code(x) return Float64(1.0 / cosh(x)) end
function tmp = code(x) tmp = 1.0 / cosh(x); end
code[x_] := N[(1.0 / N[Cosh[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\cosh x}
\end{array}
Initial program 100.0%
clear-numN/A
/-lowering-/.f64N/A
cosh-defN/A
cosh-lowering-cosh.f64100.0%
Applied egg-rr100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.08333333333333333 (* (* x x) 0.002777777777777778)))
(t_1 (* (* x x) t_0))
(t_2 (* (* x x) (* x x))))
(if (<= x 4e+38)
(/
2.0
(/
(-
(/
(*
(* (* x x) (* (* x x) (+ 1.0 t_1)))
(- 1.0 (* (* x x) (* t_0 t_1))))
(- 1.0 t_1))
4.0)
(- (* (* x x) (+ 1.0 (* x (* x t_0)))) 2.0)))
(/ 21600.0 (* t_2 t_2)))))
double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = (x * x) * t_0;
double t_2 = (x * x) * (x * x);
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / ((((((x * x) * ((x * x) * (1.0 + t_1))) * (1.0 - ((x * x) * (t_0 * t_1)))) / (1.0 - t_1)) - 4.0) / (((x * x) * (1.0 + (x * (x * t_0)))) - 2.0));
} else {
tmp = 21600.0 / (t_2 * t_2);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = 0.08333333333333333d0 + ((x * x) * 0.002777777777777778d0)
t_1 = (x * x) * t_0
t_2 = (x * x) * (x * x)
if (x <= 4d+38) then
tmp = 2.0d0 / ((((((x * x) * ((x * x) * (1.0d0 + t_1))) * (1.0d0 - ((x * x) * (t_0 * t_1)))) / (1.0d0 - t_1)) - 4.0d0) / (((x * x) * (1.0d0 + (x * (x * t_0)))) - 2.0d0))
else
tmp = 21600.0d0 / (t_2 * t_2)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = (x * x) * t_0;
double t_2 = (x * x) * (x * x);
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / ((((((x * x) * ((x * x) * (1.0 + t_1))) * (1.0 - ((x * x) * (t_0 * t_1)))) / (1.0 - t_1)) - 4.0) / (((x * x) * (1.0 + (x * (x * t_0)))) - 2.0));
} else {
tmp = 21600.0 / (t_2 * t_2);
}
return tmp;
}
def code(x): t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778) t_1 = (x * x) * t_0 t_2 = (x * x) * (x * x) tmp = 0 if x <= 4e+38: tmp = 2.0 / ((((((x * x) * ((x * x) * (1.0 + t_1))) * (1.0 - ((x * x) * (t_0 * t_1)))) / (1.0 - t_1)) - 4.0) / (((x * x) * (1.0 + (x * (x * t_0)))) - 2.0)) else: tmp = 21600.0 / (t_2 * t_2) return tmp
function code(x) t_0 = Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778)) t_1 = Float64(Float64(x * x) * t_0) t_2 = Float64(Float64(x * x) * Float64(x * x)) tmp = 0.0 if (x <= 4e+38) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(Float64(x * x) * Float64(Float64(x * x) * Float64(1.0 + t_1))) * Float64(1.0 - Float64(Float64(x * x) * Float64(t_0 * t_1)))) / Float64(1.0 - t_1)) - 4.0) / Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * t_0)))) - 2.0))); else tmp = Float64(21600.0 / Float64(t_2 * t_2)); end return tmp end
function tmp_2 = code(x) t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778); t_1 = (x * x) * t_0; t_2 = (x * x) * (x * x); tmp = 0.0; if (x <= 4e+38) tmp = 2.0 / ((((((x * x) * ((x * x) * (1.0 + t_1))) * (1.0 - ((x * x) * (t_0 * t_1)))) / (1.0 - t_1)) - 4.0) / (((x * x) * (1.0 + (x * (x * t_0)))) - 2.0)); else tmp = 21600.0 / (t_2 * t_2); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x * x), $MachinePrecision] * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 4e+38], N[(2.0 / N[(N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[(N[(x * x), $MachinePrecision] * N[(t$95$0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$1), $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\\
t_1 := \left(x \cdot x\right) \cdot t\_0\\
t_2 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 4 \cdot 10^{+38}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(1 + t\_1\right)\right)\right) \cdot \left(1 - \left(x \cdot x\right) \cdot \left(t\_0 \cdot t\_1\right)\right)}{1 - t\_1} - 4}{\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot t\_0\right)\right) - 2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_2 \cdot t\_2}\\
\end{array}
\end{array}
if x < 3.99999999999999991e38Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.1%
Simplified92.1%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr69.0%
associate-*r*N/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
flip-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr70.0%
if 3.99999999999999991e38 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.2%
Simplified95.2%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr5.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6410.3%
Simplified10.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x)))
(t_1
(+
1.0
(*
x
(* x (+ 0.08333333333333333 (* (* x x) 0.002777777777777778))))))
(t_2 (+ 1.0 (* (* x x) 0.08333333333333333)))
(t_3 (* x t_2)))
(if (<= x 4e+38)
(/
(* 2.0 (+ (* (* x x) (* t_3 t_3)) -4.0))
(* (+ -4.0 (* (* x (* x (* x x))) (* t_1 t_1))) (+ 2.0 (* (* x x) t_2))))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = 1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))));
double t_2 = 1.0 + ((x * x) * 0.08333333333333333);
double t_3 = x * t_2;
double tmp;
if (x <= 4e+38) {
tmp = (2.0 * (((x * x) * (t_3 * t_3)) + -4.0)) / ((-4.0 + ((x * (x * (x * x))) * (t_1 * t_1))) * (2.0 + ((x * x) * t_2)));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_0 = (x * x) * (x * x)
t_1 = 1.0d0 + (x * (x * (0.08333333333333333d0 + ((x * x) * 0.002777777777777778d0))))
t_2 = 1.0d0 + ((x * x) * 0.08333333333333333d0)
t_3 = x * t_2
if (x <= 4d+38) then
tmp = (2.0d0 * (((x * x) * (t_3 * t_3)) + (-4.0d0))) / (((-4.0d0) + ((x * (x * (x * x))) * (t_1 * t_1))) * (2.0d0 + ((x * x) * t_2)))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = 1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))));
double t_2 = 1.0 + ((x * x) * 0.08333333333333333);
double t_3 = x * t_2;
double tmp;
if (x <= 4e+38) {
tmp = (2.0 * (((x * x) * (t_3 * t_3)) + -4.0)) / ((-4.0 + ((x * (x * (x * x))) * (t_1 * t_1))) * (2.0 + ((x * x) * t_2)));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) t_1 = 1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))) t_2 = 1.0 + ((x * x) * 0.08333333333333333) t_3 = x * t_2 tmp = 0 if x <= 4e+38: tmp = (2.0 * (((x * x) * (t_3 * t_3)) + -4.0)) / ((-4.0 + ((x * (x * (x * x))) * (t_1 * t_1))) * (2.0 + ((x * x) * t_2))) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) t_1 = Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778))))) t_2 = Float64(1.0 + Float64(Float64(x * x) * 0.08333333333333333)) t_3 = Float64(x * t_2) tmp = 0.0 if (x <= 4e+38) tmp = Float64(Float64(2.0 * Float64(Float64(Float64(x * x) * Float64(t_3 * t_3)) + -4.0)) / Float64(Float64(-4.0 + Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(t_1 * t_1))) * Float64(2.0 + Float64(Float64(x * x) * t_2)))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); t_1 = 1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))); t_2 = 1.0 + ((x * x) * 0.08333333333333333); t_3 = x * t_2; tmp = 0.0; if (x <= 4e+38) tmp = (2.0 * (((x * x) * (t_3 * t_3)) + -4.0)) / ((-4.0 + ((x * (x * (x * x))) * (t_1 * t_1))) * (2.0 + ((x * x) * t_2))); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(1.0 + N[(N[(x * x), $MachinePrecision] * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x * t$95$2), $MachinePrecision]}, If[LessEqual[x, 4e+38], N[(N[(2.0 * N[(N[(N[(x * x), $MachinePrecision] * N[(t$95$3 * t$95$3), $MachinePrecision]), $MachinePrecision] + -4.0), $MachinePrecision]), $MachinePrecision] / N[(N[(-4.0 + N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(2.0 + N[(N[(x * x), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
t_1 := 1 + x \cdot \left(x \cdot \left(0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\right)\right)\\
t_2 := 1 + \left(x \cdot x\right) \cdot 0.08333333333333333\\
t_3 := x \cdot t\_2\\
\mathbf{if}\;x \leq 4 \cdot 10^{+38}:\\
\;\;\;\;\frac{2 \cdot \left(\left(x \cdot x\right) \cdot \left(t\_3 \cdot t\_3\right) + -4\right)}{\left(-4 + \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(t\_1 \cdot t\_1\right)\right) \cdot \left(2 + \left(x \cdot x\right) \cdot t\_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 3.99999999999999991e38Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.1%
Simplified92.1%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr69.0%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6472.5%
Simplified72.5%
Applied egg-rr68.4%
if 3.99999999999999991e38 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.2%
Simplified95.2%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr5.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6410.3%
Simplified10.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification75.6%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.08333333333333333 (* (* x x) 0.002777777777777778)))
(t_1 (* (* x x) (* x x)))
(t_2 (* x (* x t_0))))
(if (<= x 4e+38)
(/
2.0
(+
2.0
(/
(* (* x x) (+ 1.0 (* t_2 (* t_0 (* (* x x) t_2)))))
(+ 1.0 (* t_2 (+ t_2 -1.0))))))
(/ 21600.0 (* t_1 t_1)))))
double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = (x * x) * (x * x);
double t_2 = x * (x * t_0);
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * (t_0 * ((x * x) * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0)))));
} else {
tmp = 21600.0 / (t_1 * t_1);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = 0.08333333333333333d0 + ((x * x) * 0.002777777777777778d0)
t_1 = (x * x) * (x * x)
t_2 = x * (x * t_0)
if (x <= 4d+38) then
tmp = 2.0d0 / (2.0d0 + (((x * x) * (1.0d0 + (t_2 * (t_0 * ((x * x) * t_2))))) / (1.0d0 + (t_2 * (t_2 + (-1.0d0))))))
else
tmp = 21600.0d0 / (t_1 * t_1)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = (x * x) * (x * x);
double t_2 = x * (x * t_0);
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * (t_0 * ((x * x) * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0)))));
} else {
tmp = 21600.0 / (t_1 * t_1);
}
return tmp;
}
def code(x): t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778) t_1 = (x * x) * (x * x) t_2 = x * (x * t_0) tmp = 0 if x <= 4e+38: tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * (t_0 * ((x * x) * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0))))) else: tmp = 21600.0 / (t_1 * t_1) return tmp
function code(x) t_0 = Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778)) t_1 = Float64(Float64(x * x) * Float64(x * x)) t_2 = Float64(x * Float64(x * t_0)) tmp = 0.0 if (x <= 4e+38) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(t_2 * Float64(t_0 * Float64(Float64(x * x) * t_2))))) / Float64(1.0 + Float64(t_2 * Float64(t_2 + -1.0)))))); else tmp = Float64(21600.0 / Float64(t_1 * t_1)); end return tmp end
function tmp_2 = code(x) t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778); t_1 = (x * x) * (x * x); t_2 = x * (x * t_0); tmp = 0.0; if (x <= 4e+38) tmp = 2.0 / (2.0 + (((x * x) * (1.0 + (t_2 * (t_0 * ((x * x) * t_2))))) / (1.0 + (t_2 * (t_2 + -1.0))))); else tmp = 21600.0 / (t_1 * t_1); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 4e+38], N[(2.0 / N[(2.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(t$95$2 * N[(t$95$0 * N[(N[(x * x), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$2 * N[(t$95$2 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\\
t_1 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
t_2 := x \cdot \left(x \cdot t\_0\right)\\
\mathbf{if}\;x \leq 4 \cdot 10^{+38}:\\
\;\;\;\;\frac{2}{2 + \frac{\left(x \cdot x\right) \cdot \left(1 + t\_2 \cdot \left(t\_0 \cdot \left(\left(x \cdot x\right) \cdot t\_2\right)\right)\right)}{1 + t\_2 \cdot \left(t\_2 + -1\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_1 \cdot t\_1}\\
\end{array}
\end{array}
if x < 3.99999999999999991e38Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.1%
Simplified92.1%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr68.0%
if 3.99999999999999991e38 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.2%
Simplified95.2%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr5.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6410.3%
Simplified10.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification75.2%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x)))
(t_1
(*
(* x x)
(+
1.0
(*
x
(* x (+ 0.08333333333333333 (* (* x x) 0.002777777777777778))))))))
(if (<= x 4e+38)
(/ 2.0 (/ (- (* t_1 t_1) 4.0) (- t_1 2.0)))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))));
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / (((t_1 * t_1) - 4.0) / (t_1 - 2.0));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x * x) * (x * x)
t_1 = (x * x) * (1.0d0 + (x * (x * (0.08333333333333333d0 + ((x * x) * 0.002777777777777778d0)))))
if (x <= 4d+38) then
tmp = 2.0d0 / (((t_1 * t_1) - 4.0d0) / (t_1 - 2.0d0))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double t_1 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))));
double tmp;
if (x <= 4e+38) {
tmp = 2.0 / (((t_1 * t_1) - 4.0) / (t_1 - 2.0));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) t_1 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))))) tmp = 0 if x <= 4e+38: tmp = 2.0 / (((t_1 * t_1) - 4.0) / (t_1 - 2.0)) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) t_1 = Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778)))))) tmp = 0.0 if (x <= 4e+38) tmp = Float64(2.0 / Float64(Float64(Float64(t_1 * t_1) - 4.0) / Float64(t_1 - 2.0))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); t_1 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))))); tmp = 0.0; if (x <= 4e+38) tmp = 2.0 / (((t_1 * t_1) - 4.0) / (t_1 - 2.0)); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 4e+38], N[(2.0 / N[(N[(N[(t$95$1 * t$95$1), $MachinePrecision] - 4.0), $MachinePrecision] / N[(t$95$1 - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
t_1 := \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\right)\right)\right)\\
\mathbf{if}\;x \leq 4 \cdot 10^{+38}:\\
\;\;\;\;\frac{2}{\frac{t\_1 \cdot t\_1 - 4}{t\_1 - 2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 3.99999999999999991e38Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.1%
Simplified92.1%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr69.0%
if 3.99999999999999991e38 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.2%
Simplified95.2%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr5.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6410.3%
Simplified10.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(let* ((t_0
(*
(* x x)
(+
1.0
(*
x
(* x (+ 0.08333333333333333 (* (* x x) 0.002777777777777778))))))))
(if (<= x 1e+77)
(/
2.0
(/
(- (* t_0 t_0) 4.0)
(- (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333)))) 2.0)))
(/ 24.0 (* x (* x (* x x)))))))
double code(double x) {
double t_0 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))));
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((t_0 * t_0) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (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) * (1.0d0 + (x * (x * (0.08333333333333333d0 + ((x * x) * 0.002777777777777778d0)))))
if (x <= 1d+77) then
tmp = 2.0d0 / (((t_0 * t_0) - 4.0d0) / (((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))) - 2.0d0))
else
tmp = 24.0d0 / (x * (x * (x * x)))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778)))));
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((t_0 * t_0) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
def code(x): t_0 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))))) tmp = 0 if x <= 1e+77: tmp = 2.0 / (((t_0 * t_0) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)) else: tmp = 24.0 / (x * (x * (x * x))) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778)))))) tmp = 0.0 if (x <= 1e+77) tmp = Float64(2.0 / Float64(Float64(Float64(t_0 * t_0) - 4.0) / Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))) - 2.0))); else tmp = Float64(24.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (1.0 + (x * (x * (0.08333333333333333 + ((x * x) * 0.002777777777777778))))); tmp = 0.0; if (x <= 1e+77) tmp = 2.0 / (((t_0 * t_0) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)); else tmp = 24.0 / (x * (x * (x * x))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1e+77], N[(2.0 / N[(N[(N[(t$95$0 * t$95$0), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(24.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\right)\right)\right)\\
\mathbf{if}\;x \leq 10^{+77}:\\
\;\;\;\;\frac{2}{\frac{t\_0 \cdot t\_0 - 4}{\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right) - 2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{24}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if x < 9.99999999999999983e76Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr68.5%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6473.3%
Simplified73.3%
if 9.99999999999999983e76 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
/-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%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) 0.002777777777777778)))
(if (<= x 1e+77)
(/
2.0
(/
(-
(*
(* (* x x) (+ 1.0 (* x (* x (+ 0.08333333333333333 t_0)))))
(* (* x x) (+ 1.0 (* x (* x t_0)))))
4.0)
(- (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333)))) 2.0)))
(/ 24.0 (* x (* x (* x x)))))))
double code(double x) {
double t_0 = (x * x) * 0.002777777777777778;
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((((x * x) * (1.0 + (x * (x * (0.08333333333333333 + t_0))))) * ((x * x) * (1.0 + (x * (x * t_0))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (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.002777777777777778d0
if (x <= 1d+77) then
tmp = 2.0d0 / (((((x * x) * (1.0d0 + (x * (x * (0.08333333333333333d0 + t_0))))) * ((x * x) * (1.0d0 + (x * (x * t_0))))) - 4.0d0) / (((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))) - 2.0d0))
else
tmp = 24.0d0 / (x * (x * (x * x)))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * 0.002777777777777778;
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((((x * x) * (1.0 + (x * (x * (0.08333333333333333 + t_0))))) * ((x * x) * (1.0 + (x * (x * t_0))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
def code(x): t_0 = (x * x) * 0.002777777777777778 tmp = 0 if x <= 1e+77: tmp = 2.0 / (((((x * x) * (1.0 + (x * (x * (0.08333333333333333 + t_0))))) * ((x * x) * (1.0 + (x * (x * t_0))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)) else: tmp = 24.0 / (x * (x * (x * x))) return tmp
function code(x) t_0 = Float64(Float64(x * x) * 0.002777777777777778) tmp = 0.0 if (x <= 1e+77) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + t_0))))) * Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * t_0))))) - 4.0) / Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))) - 2.0))); else tmp = Float64(24.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * 0.002777777777777778; tmp = 0.0; if (x <= 1e+77) tmp = 2.0 / (((((x * x) * (1.0 + (x * (x * (0.08333333333333333 + t_0))))) * ((x * x) * (1.0 + (x * (x * t_0))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)); else tmp = 24.0 / (x * (x * (x * x))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]}, If[LessEqual[x, 1e+77], N[(2.0 / N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(24.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot 0.002777777777777778\\
\mathbf{if}\;x \leq 10^{+77}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + t\_0\right)\right)\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot t\_0\right)\right)\right) - 4}{\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right) - 2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{24}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if x < 9.99999999999999983e76Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr68.5%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6473.3%
Simplified73.3%
Taylor expanded in x around inf
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6473.3%
Simplified73.3%
if 9.99999999999999983e76 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
/-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%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.08333333333333333 (* (* x x) 0.002777777777777778)))
(t_1 (* x (* x t_0))))
(if (<= x 1e+77)
(/
2.0
(+ 2.0 (/ (* (* x x) (- 1.0 (* t_0 (* (* x x) t_1)))) (- 1.0 t_1))))
(/ 24.0 (* x (* x (* x x)))))))
double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = x * (x * t_0);
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 - (t_0 * ((x * x) * t_1)))) / (1.0 - t_1)));
} else {
tmp = 24.0 / (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.08333333333333333d0 + ((x * x) * 0.002777777777777778d0)
t_1 = x * (x * t_0)
if (x <= 1d+77) then
tmp = 2.0d0 / (2.0d0 + (((x * x) * (1.0d0 - (t_0 * ((x * x) * t_1)))) / (1.0d0 - t_1)))
else
tmp = 24.0d0 / (x * (x * (x * x)))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778);
double t_1 = x * (x * t_0);
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (2.0 + (((x * x) * (1.0 - (t_0 * ((x * x) * t_1)))) / (1.0 - t_1)));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
def code(x): t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778) t_1 = x * (x * t_0) tmp = 0 if x <= 1e+77: tmp = 2.0 / (2.0 + (((x * x) * (1.0 - (t_0 * ((x * x) * t_1)))) / (1.0 - t_1))) else: tmp = 24.0 / (x * (x * (x * x))) return tmp
function code(x) t_0 = Float64(0.08333333333333333 + Float64(Float64(x * x) * 0.002777777777777778)) t_1 = Float64(x * Float64(x * t_0)) tmp = 0.0 if (x <= 1e+77) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(Float64(x * x) * Float64(1.0 - Float64(t_0 * Float64(Float64(x * x) * t_1)))) / Float64(1.0 - t_1)))); else tmp = Float64(24.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
function tmp_2 = code(x) t_0 = 0.08333333333333333 + ((x * x) * 0.002777777777777778); t_1 = x * (x * t_0); tmp = 0.0; if (x <= 1e+77) tmp = 2.0 / (2.0 + (((x * x) * (1.0 - (t_0 * ((x * x) * t_1)))) / (1.0 - t_1))); else tmp = 24.0 / (x * (x * (x * x))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1e+77], N[(2.0 / N[(2.0 + N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 - N[(t$95$0 * N[(N[(x * x), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(24.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.08333333333333333 + \left(x \cdot x\right) \cdot 0.002777777777777778\\
t_1 := x \cdot \left(x \cdot t\_0\right)\\
\mathbf{if}\;x \leq 10^{+77}:\\
\;\;\;\;\frac{2}{2 + \frac{\left(x \cdot x\right) \cdot \left(1 - t\_0 \cdot \left(\left(x \cdot x\right) \cdot t\_1\right)\right)}{1 - t\_1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{24}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if x < 9.99999999999999983e76Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr72.4%
if 9.99999999999999983e76 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
/-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 simplification78.0%
(FPCore (x)
:precision binary64
(if (<= x 1e+77)
(/
2.0
(/
(-
(*
(* (* x x) (* x x))
(+
1.0
(*
x
(*
x
(+
0.16666666666666666
(* x (* x (+ 0.0125 (* (* x x) 0.000462962962962963)))))))))
4.0)
(- (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333)))) 2.0)))
(/ 24.0 (* x (* x (* x x))))))
double code(double x) {
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((((x * x) * (x * x)) * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * (0.0125 + ((x * x) * 0.000462962962962963))))))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1d+77) then
tmp = 2.0d0 / (((((x * x) * (x * x)) * (1.0d0 + (x * (x * (0.16666666666666666d0 + (x * (x * (0.0125d0 + ((x * x) * 0.000462962962962963d0))))))))) - 4.0d0) / (((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))) - 2.0d0))
else
tmp = 24.0d0 / (x * (x * (x * x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1e+77) {
tmp = 2.0 / (((((x * x) * (x * x)) * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * (0.0125 + ((x * x) * 0.000462962962962963))))))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1e+77: tmp = 2.0 / (((((x * x) * (x * x)) * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * (0.0125 + ((x * x) * 0.000462962962962963))))))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)) else: tmp = 24.0 / (x * (x * (x * x))) return tmp
function code(x) tmp = 0.0 if (x <= 1e+77) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(x * x) * Float64(x * x)) * Float64(1.0 + Float64(x * Float64(x * Float64(0.16666666666666666 + Float64(x * Float64(x * Float64(0.0125 + Float64(Float64(x * x) * 0.000462962962962963))))))))) - 4.0) / Float64(Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))) - 2.0))); else tmp = Float64(24.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1e+77) tmp = 2.0 / (((((x * x) * (x * x)) * (1.0 + (x * (x * (0.16666666666666666 + (x * (x * (0.0125 + ((x * x) * 0.000462962962962963))))))))) - 4.0) / (((x * x) * (1.0 + (x * (x * 0.08333333333333333)))) - 2.0)); else tmp = 24.0 / (x * (x * (x * x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1e+77], N[(2.0 / N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.16666666666666666 + N[(x * N[(x * N[(0.0125 + N[(N[(x * x), $MachinePrecision] * 0.000462962962962963), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(24.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 10^{+77}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.16666666666666666 + x \cdot \left(x \cdot \left(0.0125 + \left(x \cdot x\right) \cdot 0.000462962962962963\right)\right)\right)\right)\right) - 4}{\left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right) - 2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{24}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if x < 9.99999999999999983e76Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr68.5%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6473.3%
Simplified73.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
Simplified72.4%
if 9.99999999999999983e76 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
/-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%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x))))
(if (<= x 5.2)
(/
2.0
(+
2.0
(*
(* x x)
(+
1.0
(/
(* (* x x) 0.006944444444444444)
(+ 0.08333333333333333 (* (* x x) -0.002777777777777778)))))))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 5.2) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (((x * x) * 0.006944444444444444) / (0.08333333333333333 + ((x * x) * -0.002777777777777778))))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = (x * x) * (x * x)
if (x <= 5.2d0) then
tmp = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (((x * x) * 0.006944444444444444d0) / (0.08333333333333333d0 + ((x * x) * (-0.002777777777777778d0)))))))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 5.2) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (((x * x) * 0.006944444444444444) / (0.08333333333333333 + ((x * x) * -0.002777777777777778))))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) tmp = 0 if x <= 5.2: tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (((x * x) * 0.006944444444444444) / (0.08333333333333333 + ((x * x) * -0.002777777777777778)))))) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) tmp = 0.0 if (x <= 5.2) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(Float64(Float64(x * x) * 0.006944444444444444) / Float64(0.08333333333333333 + Float64(Float64(x * x) * -0.002777777777777778))))))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); tmp = 0.0; if (x <= 5.2) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (((x * x) * 0.006944444444444444) / (0.08333333333333333 + ((x * x) * -0.002777777777777778)))))); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5.2], N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(N[(N[(x * x), $MachinePrecision] * 0.006944444444444444), $MachinePrecision] / N[(0.08333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 5.2:\\
\;\;\;\;\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + \frac{\left(x \cdot x\right) \cdot 0.006944444444444444}{0.08333333333333333 + \left(x \cdot x\right) \cdot -0.002777777777777778}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 5.20000000000000018Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.7%
Simplified95.7%
associate-*r*N/A
flip-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr75.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.2%
Simplified67.2%
if 5.20000000000000018 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.3%
Simplified84.3%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr13.8%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6418.3%
Simplified18.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.8%
Simplified88.8%
Final simplification72.8%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x))))
(if (<= x 7.2)
(/
2.0
(+
2.0
(*
(* x x)
(+
1.0
(*
x
(* x (+ 0.08333333333333333 (* x (* x 0.002777777777777778)))))))))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 7.2) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778))))))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = (x * x) * (x * x)
if (x <= 7.2d0) then
tmp = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * (0.08333333333333333d0 + (x * (x * 0.002777777777777778d0))))))))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 7.2) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778))))))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) tmp = 0 if x <= 7.2: tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778)))))))) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) tmp = 0.0 if (x <= 7.2) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * Float64(0.08333333333333333 + Float64(x * Float64(x * 0.002777777777777778))))))))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); tmp = 0.0; if (x <= 7.2) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * (0.08333333333333333 + (x * (x * 0.002777777777777778)))))))); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 7.2], N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * N[(0.08333333333333333 + N[(x * N[(x * 0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 7.2:\\
\;\;\;\;\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot \left(0.08333333333333333 + x \cdot \left(x \cdot 0.002777777777777778\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 7.20000000000000018Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.7%
Simplified95.7%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6495.7%
Applied egg-rr95.7%
if 7.20000000000000018 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.3%
Simplified84.3%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr13.8%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6418.3%
Simplified18.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.8%
Simplified88.8%
Final simplification94.0%
(FPCore (x)
:precision binary64
(/
2.0
(+
2.0
(*
(* x x)
(+
1.0
(*
(* x x)
(+
0.08333333333333333
(* (* (* x x) (* x x)) -0.002777777777777778))))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + ((x * x) * (0.08333333333333333 + (((x * x) * (x * x)) * -0.002777777777777778))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + ((x * x) * (0.08333333333333333d0 + (((x * x) * (x * x)) * (-0.002777777777777778d0)))))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + ((x * x) * (0.08333333333333333 + (((x * x) * (x * x)) * -0.002777777777777778))))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + ((x * x) * (0.08333333333333333 + (((x * x) * (x * x)) * -0.002777777777777778))))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(Float64(x * x) * Float64(0.08333333333333333 + Float64(Float64(Float64(x * x) * Float64(x * x)) * -0.002777777777777778))))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + ((x * x) * (0.08333333333333333 + (((x * x) * (x * x)) * -0.002777777777777778)))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(0.08333333333333333 + N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * -0.002777777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + \left(x \cdot x\right) \cdot \left(0.08333333333333333 + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot -0.002777777777777778\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.8%
Simplified92.8%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr54.6%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6458.4%
Simplified58.4%
Taylor expanded in x around 0
+-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
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.7%
Simplified94.7%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) (* x x))))
(if (<= x 6.0)
(/ 2.0 (+ 2.0 (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333))))))
(/ 21600.0 (* t_0 t_0)))))
double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 6.0) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = (x * x) * (x * x)
if (x <= 6.0d0) then
tmp = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))))
else
tmp = 21600.0d0 / (t_0 * t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * x) * (x * x);
double tmp;
if (x <= 6.0) {
tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
} else {
tmp = 21600.0 / (t_0 * t_0);
}
return tmp;
}
def code(x): t_0 = (x * x) * (x * x) tmp = 0 if x <= 6.0: tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333))))) else: tmp = 21600.0 / (t_0 * t_0) return tmp
function code(x) t_0 = Float64(Float64(x * x) * Float64(x * x)) tmp = 0.0 if (x <= 6.0) tmp = Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))))); else tmp = Float64(21600.0 / Float64(t_0 * t_0)); end return tmp end
function tmp_2 = code(x) t_0 = (x * x) * (x * x); tmp = 0.0; if (x <= 6.0) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333))))); else tmp = 21600.0 / (t_0 * t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 6.0], N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(21600.0 / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 6:\\
\;\;\;\;\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{21600}{t\_0 \cdot t\_0}\\
\end{array}
\end{array}
if x < 6Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.3%
Simplified92.3%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6492.3%
Applied egg-rr92.3%
if 6 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.3%
Simplified84.3%
+-commutativeN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr13.8%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6418.3%
Simplified18.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.8%
Simplified88.8%
Final simplification91.4%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* (* x x) (+ 1.0 (* x (* x 0.08333333333333333)))))))
double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + ((x * x) * (1.0d0 + (x * (x * 0.08333333333333333d0)))))
end function
public static double code(double x) {
return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))));
}
def code(x): return 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333)))))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(Float64(x * x) * Float64(1.0 + Float64(x * Float64(x * 0.08333333333333333)))))) end
function tmp = code(x) tmp = 2.0 / (2.0 + ((x * x) * (1.0 + (x * (x * 0.08333333333333333))))); end
code[x_] := N[(2.0 / N[(2.0 + N[(N[(x * x), $MachinePrecision] * N[(1.0 + N[(x * N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + \left(x \cdot x\right) \cdot \left(1 + x \cdot \left(x \cdot 0.08333333333333333\right)\right)}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6489.1%
Simplified89.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6489.1%
Applied egg-rr89.1%
Final simplification89.1%
(FPCore (x) :precision binary64 (if (<= x 3.75) (/ 2.0 (+ 2.0 (* x x))) (/ 24.0 (* x (* x (* x x))))))
double code(double x) {
double tmp;
if (x <= 3.75) {
tmp = 2.0 / (2.0 + (x * x));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 3.75d0) then
tmp = 2.0d0 / (2.0d0 + (x * x))
else
tmp = 24.0d0 / (x * (x * (x * x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 3.75) {
tmp = 2.0 / (2.0 + (x * x));
} else {
tmp = 24.0 / (x * (x * (x * x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= 3.75: tmp = 2.0 / (2.0 + (x * x)) else: tmp = 24.0 / (x * (x * (x * x))) return tmp
function code(x) tmp = 0.0 if (x <= 3.75) tmp = Float64(2.0 / Float64(2.0 + Float64(x * x))); else tmp = Float64(24.0 / Float64(x * Float64(x * Float64(x * x)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 3.75) tmp = 2.0 / (2.0 + (x * x)); else tmp = 24.0 / (x * (x * (x * x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 3.75], N[(2.0 / N[(2.0 + N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(24.0 / N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.75:\\
\;\;\;\;\frac{2}{2 + x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\frac{24}{x \cdot \left(x \cdot \left(x \cdot x\right)\right)}\\
\end{array}
\end{array}
if x < 3.75Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6487.6%
Simplified87.6%
if 3.75 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6479.9%
Simplified79.9%
Taylor expanded in x around inf
/-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-*.f6479.9%
Simplified79.9%
(FPCore (x) :precision binary64 (if (<= x 1.42) 1.0 (/ 2.0 (* x x))))
double code(double x) {
double tmp;
if (x <= 1.42) {
tmp = 1.0;
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.42d0) then
tmp = 1.0d0
else
tmp = 2.0d0 / (x * x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.42) {
tmp = 1.0;
} else {
tmp = 2.0 / (x * x);
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.42: tmp = 1.0 else: tmp = 2.0 / (x * x) return tmp
function code(x) tmp = 0.0 if (x <= 1.42) tmp = 1.0; else tmp = Float64(2.0 / Float64(x * x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.42) tmp = 1.0; else tmp = 2.0 / (x * x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.42], 1.0, N[(2.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.42:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{x \cdot x}\\
\end{array}
\end{array}
if x < 1.4199999999999999Initial program 100.0%
Taylor expanded in x around 0
Simplified67.3%
if 1.4199999999999999 < x Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.3%
Simplified57.3%
(FPCore (x) :precision binary64 (/ 2.0 (+ 2.0 (* x x))))
double code(double x) {
return 2.0 / (2.0 + (x * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 / (2.0d0 + (x * x))
end function
public static double code(double x) {
return 2.0 / (2.0 + (x * x));
}
def code(x): return 2.0 / (2.0 + (x * x))
function code(x) return Float64(2.0 / Float64(2.0 + Float64(x * x))) end
function tmp = code(x) tmp = 2.0 / (2.0 + (x * x)); end
code[x_] := N[(2.0 / N[(2.0 + N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{2 + x \cdot x}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6479.8%
Simplified79.8%
(FPCore (x) :precision binary64 1.0)
double code(double x) {
return 1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double x) {
return 1.0;
}
def code(x): return 1.0
function code(x) return 1.0 end
function tmp = code(x) tmp = 1.0; end
code[x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
Simplified50.7%
herbie shell --seed 2024158
(FPCore (x)
:name "Hyperbolic secant"
:precision binary64
(/ 2.0 (+ (exp x) (exp (- x)))))