
(FPCore (x n) :precision binary64 (- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))
double code(double x, double n) {
return pow((x + 1.0), (1.0 / n)) - pow(x, (1.0 / n));
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
code = ((x + 1.0d0) ** (1.0d0 / n)) - (x ** (1.0d0 / n))
end function
public static double code(double x, double n) {
return Math.pow((x + 1.0), (1.0 / n)) - Math.pow(x, (1.0 / n));
}
def code(x, n): return math.pow((x + 1.0), (1.0 / n)) - math.pow(x, (1.0 / n))
function code(x, n) return Float64((Float64(x + 1.0) ^ Float64(1.0 / n)) - (x ^ Float64(1.0 / n))) end
function tmp = code(x, n) tmp = ((x + 1.0) ^ (1.0 / n)) - (x ^ (1.0 / n)); end
code[x_, n_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x n) :precision binary64 (- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))
double code(double x, double n) {
return pow((x + 1.0), (1.0 / n)) - pow(x, (1.0 / n));
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
code = ((x + 1.0d0) ** (1.0d0 / n)) - (x ** (1.0d0 / n))
end function
public static double code(double x, double n) {
return Math.pow((x + 1.0), (1.0 / n)) - Math.pow(x, (1.0 / n));
}
def code(x, n): return math.pow((x + 1.0), (1.0 / n)) - math.pow(x, (1.0 / n))
function code(x, n) return Float64((Float64(x + 1.0) ^ Float64(1.0 / n)) - (x ^ Float64(1.0 / n))) end
function tmp = code(x, n) tmp = ((x + 1.0) ^ (1.0 / n)) - (x ^ (1.0 / n)); end
code[x_, n_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\end{array}
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 4e-20)
(/ (log (/ x (+ 1.0 x))) (- n))
(- (exp (/ (log1p x) n)) (pow x (/ 1.0 n))))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 4e-20) {
tmp = log((x / (1.0 + x))) / -n;
} else {
tmp = exp((log1p(x) / n)) - pow(x, (1.0 / n));
}
return tmp;
}
public static double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * Math.pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 4e-20) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else {
tmp = Math.exp((Math.log1p(x) / n)) - Math.pow(x, (1.0 / n));
}
return tmp;
}
def code(x, n): tmp = 0 if (1.0 / n) <= -4e-14: tmp = 1.0 / ((n * x) * math.pow(x, (-1.0 / n))) elif (1.0 / n) <= 4e-20: tmp = math.log((x / (1.0 + x))) / -n else: tmp = math.exp((math.log1p(x) / n)) - math.pow(x, (1.0 / n)) return tmp
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 4e-20) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); else tmp = Float64(exp(Float64(log1p(x) / n)) - (x ^ Float64(1.0 / n))); end return tmp end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 4e-20], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], N[(N[Exp[N[(N[Log[1 + x], $MachinePrecision] / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 4 \cdot 10^{-20}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{\mathsf{log1p}\left(x\right)}{n}} - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 3.99999999999999978e-20Initial program 32.0%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6478.1
Simplified78.1%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6478.3
Applied egg-rr78.3%
if 3.99999999999999978e-20 < (/.f64 #s(literal 1 binary64) n) Initial program 58.0%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f6494.7
Applied egg-rr94.7%
Final simplification88.3%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n)))
(t_1 (- (pow (+ 1.0 x) (/ 1.0 n)) t_0))
(t_2 (- 1.0 t_0)))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 0.0) (/ (log (/ x (+ 1.0 x))) (- n)) t_2))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double t_1 = pow((1.0 + x), (1.0 / n)) - t_0;
double t_2 = 1.0 - t_0;
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = log((x / (1.0 + x))) / -n;
} else {
tmp = t_2;
}
return tmp;
}
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double t_1 = Math.pow((1.0 + x), (1.0 / n)) - t_0;
double t_2 = 1.0 - t_0;
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) t_1 = math.pow((1.0 + x), (1.0 / n)) - t_0 t_2 = 1.0 - t_0 tmp = 0 if t_1 <= -math.inf: tmp = t_2 elif t_1 <= 0.0: tmp = math.log((x / (1.0 + x))) / -n else: tmp = t_2 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) t_1 = Float64((Float64(1.0 + x) ^ Float64(1.0 / n)) - t_0) t_2 = Float64(1.0 - t_0) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); t_1 = ((1.0 + x) ^ (1.0 / n)) - t_0; t_2 = 1.0 - t_0; tmp = 0.0; if (t_1 <= -Inf) tmp = t_2; elseif (t_1 <= 0.0) tmp = log((x / (1.0 + x))) / -n; else tmp = t_2; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Power[N[(1.0 + x), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(1.0 - t$95$0), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
t_1 := {\left(1 + x\right)}^{\left(\frac{1}{n}\right)} - t\_0\\
t_2 := 1 - t\_0\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) < -inf.0 or 0.0 < (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) Initial program 81.0%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6479.4
Simplified79.4%
if -inf.0 < (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) < 0.0Initial program 45.7%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6478.6
Simplified78.6%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6478.3
Applied egg-rr78.3%
Final simplification78.7%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n)))
(t_1 (- (pow (+ 1.0 x) (/ 1.0 n)) t_0))
(t_2 (- 1.0 t_0)))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 0.0) (/ (log (/ (+ 1.0 x) x)) n) t_2))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double t_1 = pow((1.0 + x), (1.0 / n)) - t_0;
double t_2 = 1.0 - t_0;
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = log(((1.0 + x) / x)) / n;
} else {
tmp = t_2;
}
return tmp;
}
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double t_1 = Math.pow((1.0 + x), (1.0 / n)) - t_0;
double t_2 = 1.0 - t_0;
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = Math.log(((1.0 + x) / x)) / n;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) t_1 = math.pow((1.0 + x), (1.0 / n)) - t_0 t_2 = 1.0 - t_0 tmp = 0 if t_1 <= -math.inf: tmp = t_2 elif t_1 <= 0.0: tmp = math.log(((1.0 + x) / x)) / n else: tmp = t_2 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) t_1 = Float64((Float64(1.0 + x) ^ Float64(1.0 / n)) - t_0) t_2 = Float64(1.0 - t_0) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(log(Float64(Float64(1.0 + x) / x)) / n); else tmp = t_2; end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); t_1 = ((1.0 + x) ^ (1.0 / n)) - t_0; t_2 = 1.0 - t_0; tmp = 0.0; if (t_1 <= -Inf) tmp = t_2; elseif (t_1 <= 0.0) tmp = log(((1.0 + x) / x)) / n; else tmp = t_2; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Power[N[(1.0 + x), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(1.0 - t$95$0), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[Log[N[(N[(1.0 + x), $MachinePrecision] / x), $MachinePrecision]], $MachinePrecision] / n), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
t_1 := {\left(1 + x\right)}^{\left(\frac{1}{n}\right)} - t\_0\\
t_2 := 1 - t\_0\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{\log \left(\frac{1 + x}{x}\right)}{n}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) < -inf.0 or 0.0 < (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) Initial program 81.0%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6479.4
Simplified79.4%
if -inf.0 < (-.f64 (pow.f64 (+.f64 x #s(literal 1 binary64)) (/.f64 #s(literal 1 binary64) n)) (pow.f64 x (/.f64 #s(literal 1 binary64) n))) < 0.0Initial program 45.7%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6478.6
Simplified78.6%
lift-log1p.f64N/A
lift-log.f64N/A
lift--.f64N/A
lift-/.f6478.6
lift--.f64N/A
lift-log1p.f64N/A
lift-log.f64N/A
diff-logN/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6478.3
Applied egg-rr78.3%
Final simplification78.7%
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(- (exp (/ x n)) (pow x (/ 1.0 n))))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else {
tmp = exp((x / n)) - pow(x, (1.0 / n));
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if ((1.0d0 / n) <= (-4d-14)) then
tmp = 1.0d0 / ((n * x) * (x ** ((-1.0d0) / n)))
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else
tmp = exp((x / n)) - (x ** (1.0d0 / n))
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * Math.pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else {
tmp = Math.exp((x / n)) - Math.pow(x, (1.0 / n));
}
return tmp;
}
def code(x, n): tmp = 0 if (1.0 / n) <= -4e-14: tmp = 1.0 / ((n * x) * math.pow(x, (-1.0 / n))) elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n else: tmp = math.exp((x / n)) - math.pow(x, (1.0 / n)) return tmp
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); else tmp = Float64(exp(Float64(x / n)) - (x ^ Float64(1.0 / n))); end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = 1.0 / ((n * x) * (x ^ (-1.0 / n))); elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; else tmp = exp((x / n)) - (x ^ (1.0 / n)); end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) Initial program 60.5%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f6499.0
Applied egg-rr99.0%
Taylor expanded in x around 0
lower-/.f6498.9
Simplified98.9%
Final simplification88.2%
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 2e+223)
(- (- 1.0 (/ (fma -0.5 (/ (* x x) n) (- x)) n)) (pow x (/ 1.0 n)))
(+ (exp (/ x n)) -1.0)))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 2e+223) {
tmp = (1.0 - (fma(-0.5, ((x * x) / n), -x) / n)) - pow(x, (1.0 / n));
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 2e+223) tmp = Float64(Float64(1.0 - Float64(fma(-0.5, Float64(Float64(x * x) / n), Float64(-x)) / n)) - (x ^ Float64(1.0 / n))); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 2e+223], N[(N[(1.0 - N[(N[(-0.5 * N[(N[(x * x), $MachinePrecision] / n), $MachinePrecision] + (-x)), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 2 \cdot 10^{+223}:\\
\;\;\;\;\left(1 - \frac{\mathsf{fma}\left(-0.5, \frac{x \cdot x}{n}, -x\right)}{n}\right) - {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 2.00000000000000009e223Initial program 75.4%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f6498.7
Applied egg-rr98.7%
Taylor expanded in x around 0
lower-/.f6498.6
Simplified98.6%
Taylor expanded in n around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6476.6
Simplified76.6%
if 2.00000000000000009e223 < (/.f64 #s(literal 1 binary64) n) Initial program 14.3%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified82.3%
Final simplification84.9%
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(-
(fma x (/ (+ 1.0 (fma x -0.5 (* (/ x n) 0.5))) n) 1.0)
(pow x (/ 1.0 n))))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else {
tmp = fma(x, ((1.0 + fma(x, -0.5, ((x / n) * 0.5))) / n), 1.0) - pow(x, (1.0 / n));
}
return tmp;
}
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); else tmp = Float64(fma(x, Float64(Float64(1.0 + fma(x, -0.5, Float64(Float64(x / n) * 0.5))) / n), 1.0) - (x ^ Float64(1.0 / n))); end return tmp end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], N[(N[(x * N[(N[(1.0 + N[(x * -0.5 + N[(N[(x / n), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] + 1.0), $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x, \frac{1 + \mathsf{fma}\left(x, -0.5, \frac{x}{n} \cdot 0.5\right)}{n}, 1\right) - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) Initial program 60.5%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified70.5%
Taylor expanded in n around inf
lower-/.f64N/A
lower-+.f64N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6478.1
Simplified78.1%
Final simplification84.9%
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(- (fma x (fma x (/ 0.5 (* n n)) (/ 1.0 n)) 1.0) (pow x (/ 1.0 n))))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else {
tmp = fma(x, fma(x, (0.5 / (n * n)), (1.0 / n)), 1.0) - pow(x, (1.0 / n));
}
return tmp;
}
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); else tmp = Float64(fma(x, fma(x, Float64(0.5 / Float64(n * n)), Float64(1.0 / n)), 1.0) - (x ^ Float64(1.0 / n))); end return tmp end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], N[(N[(x * N[(x * N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] + N[(1.0 / n), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{0.5}{n \cdot n}, \frac{1}{n}\right), 1\right) - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) Initial program 60.5%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified70.5%
Taylor expanded in n around 0
lower-/.f64N/A
unpow2N/A
lower-*.f6470.5
Simplified70.5%
Final simplification83.7%
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -4e-14)
(/ 1.0 (* (* n x) (pow x (/ -1.0 n))))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 5e+202)
(- (/ (+ n x) n) (pow x (/ 1.0 n)))
(+ (exp (/ x n)) -1.0)))))
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - pow(x, (1.0 / n));
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if ((1.0d0 / n) <= (-4d-14)) then
tmp = 1.0d0 / ((n * x) * (x ** ((-1.0d0) / n)))
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else if ((1.0d0 / n) <= 5d+202) then
tmp = ((n + x) / n) - (x ** (1.0d0 / n))
else
tmp = exp((x / n)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = 1.0 / ((n * x) * Math.pow(x, (-1.0 / n)));
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - Math.pow(x, (1.0 / n));
} else {
tmp = Math.exp((x / n)) + -1.0;
}
return tmp;
}
def code(x, n): tmp = 0 if (1.0 / n) <= -4e-14: tmp = 1.0 / ((n * x) * math.pow(x, (-1.0 / n))) elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n elif (1.0 / n) <= 5e+202: tmp = ((n + x) / n) - math.pow(x, (1.0 / n)) else: tmp = math.exp((x / n)) + -1.0 return tmp
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(1.0 / Float64(Float64(n * x) * (x ^ Float64(-1.0 / n)))); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 5e+202) tmp = Float64(Float64(Float64(n + x) / n) - (x ^ Float64(1.0 / n))); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = 1.0 / ((n * x) * (x ^ (-1.0 / n))); elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; elseif ((1.0 / n) <= 5e+202) tmp = ((n + x) / n) - (x ^ (1.0 / n)); else tmp = exp((x / n)) + -1.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(1.0 / N[(N[(n * x), $MachinePrecision] * N[Power[x, N[(-1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 5e+202], N[(N[(N[(n + x), $MachinePrecision] / n), $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\left(n \cdot x\right) \cdot {x}^{\left(\frac{-1}{n}\right)}}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 5 \cdot 10^{+202}:\\
\;\;\;\;\frac{n + x}{n} - {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
clear-numN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
lift-pow.f64N/A
pow-flipN/A
lower-pow.f64N/A
lower-neg.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 4.9999999999999999e202Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified68.3%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6474.0
Simplified74.0%
Taylor expanded in n around 0
lower-/.f64N/A
+-commutativeN/A
lower-+.f6474.0
Simplified74.0%
if 4.9999999999999999e202 < (/.f64 #s(literal 1 binary64) n) Initial program 19.3%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified77.7%
Final simplification84.4%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n))))
(if (<= (/ 1.0 n) -4e-14)
(/ (/ t_0 n) x)
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 5e+202)
(- (/ (+ n x) n) t_0)
(+ (exp (/ x n)) -1.0))))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = (t_0 / n) / x;
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - t_0;
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = x ** (1.0d0 / n)
if ((1.0d0 / n) <= (-4d-14)) then
tmp = (t_0 / n) / x
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else if ((1.0d0 / n) <= 5d+202) then
tmp = ((n + x) / n) - t_0
else
tmp = exp((x / n)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = (t_0 / n) / x;
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - t_0;
} else {
tmp = Math.exp((x / n)) + -1.0;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) tmp = 0 if (1.0 / n) <= -4e-14: tmp = (t_0 / n) / x elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n elif (1.0 / n) <= 5e+202: tmp = ((n + x) / n) - t_0 else: tmp = math.exp((x / n)) + -1.0 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(Float64(t_0 / n) / x); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 5e+202) tmp = Float64(Float64(Float64(n + x) / n) - t_0); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = (t_0 / n) / x; elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; elseif ((1.0 / n) <= 5e+202) tmp = ((n + x) / n) - t_0; else tmp = exp((x / n)) + -1.0; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(N[(t$95$0 / n), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 5e+202], N[(N[(N[(n + x), $MachinePrecision] / n), $MachinePrecision] - t$95$0), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{\frac{t\_0}{n}}{x}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 5 \cdot 10^{+202}:\\
\;\;\;\;\frac{n + x}{n} - t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
lift-/.f64N/A
lift-pow.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6499.9
Applied egg-rr99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 4.9999999999999999e202Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified68.3%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6474.0
Simplified74.0%
Taylor expanded in n around 0
lower-/.f64N/A
+-commutativeN/A
lower-+.f6474.0
Simplified74.0%
if 4.9999999999999999e202 < (/.f64 #s(literal 1 binary64) n) Initial program 19.3%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified77.7%
Final simplification84.4%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n))))
(if (<= (/ 1.0 n) -4e-14)
(/ t_0 (* n x))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 5e+202)
(- (/ (+ n x) n) t_0)
(+ (exp (/ x n)) -1.0))))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - t_0;
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = x ** (1.0d0 / n)
if ((1.0d0 / n) <= (-4d-14)) then
tmp = t_0 / (n * x)
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else if ((1.0d0 / n) <= 5d+202) then
tmp = ((n + x) / n) - t_0
else
tmp = exp((x / n)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = ((n + x) / n) - t_0;
} else {
tmp = Math.exp((x / n)) + -1.0;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) tmp = 0 if (1.0 / n) <= -4e-14: tmp = t_0 / (n * x) elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n elif (1.0 / n) <= 5e+202: tmp = ((n + x) / n) - t_0 else: tmp = math.exp((x / n)) + -1.0 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(t_0 / Float64(n * x)); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 5e+202) tmp = Float64(Float64(Float64(n + x) / n) - t_0); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = t_0 / (n * x); elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; elseif ((1.0 / n) <= 5e+202) tmp = ((n + x) / n) - t_0; else tmp = exp((x / n)) + -1.0; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(t$95$0 / N[(n * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 5e+202], N[(N[(N[(n + x), $MachinePrecision] / n), $MachinePrecision] - t$95$0), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{t\_0}{n \cdot x}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 5 \cdot 10^{+202}:\\
\;\;\;\;\frac{n + x}{n} - t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 4.9999999999999999e202Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified68.3%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6474.0
Simplified74.0%
Taylor expanded in n around 0
lower-/.f64N/A
+-commutativeN/A
lower-+.f6474.0
Simplified74.0%
if 4.9999999999999999e202 < (/.f64 #s(literal 1 binary64) n) Initial program 19.3%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified77.7%
Final simplification84.4%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n))))
(if (<= (/ 1.0 n) -4e-14)
(/ t_0 (* n x))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 5e+202)
(- (+ 1.0 (/ x n)) t_0)
(+ (exp (/ x n)) -1.0))))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = (1.0 + (x / n)) - t_0;
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = x ** (1.0d0 / n)
if ((1.0d0 / n) <= (-4d-14)) then
tmp = t_0 / (n * x)
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else if ((1.0d0 / n) <= 5d+202) then
tmp = (1.0d0 + (x / n)) - t_0
else
tmp = exp((x / n)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 5e+202) {
tmp = (1.0 + (x / n)) - t_0;
} else {
tmp = Math.exp((x / n)) + -1.0;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) tmp = 0 if (1.0 / n) <= -4e-14: tmp = t_0 / (n * x) elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n elif (1.0 / n) <= 5e+202: tmp = (1.0 + (x / n)) - t_0 else: tmp = math.exp((x / n)) + -1.0 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(t_0 / Float64(n * x)); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 5e+202) tmp = Float64(Float64(1.0 + Float64(x / n)) - t_0); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = t_0 / (n * x); elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; elseif ((1.0 / n) <= 5e+202) tmp = (1.0 + (x / n)) - t_0; else tmp = exp((x / n)) + -1.0; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(t$95$0 / N[(n * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 5e+202], N[(N[(1.0 + N[(x / n), $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{t\_0}{n \cdot x}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 5 \cdot 10^{+202}:\\
\;\;\;\;\left(1 + \frac{x}{n}\right) - t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 4.9999999999999999e202Initial program 77.6%
Taylor expanded in x around 0
*-rgt-identityN/A
associate-*r/N/A
lower-+.f64N/A
associate-*r/N/A
*-rgt-identityN/A
lower-/.f6474.0
Simplified74.0%
if 4.9999999999999999e202 < (/.f64 #s(literal 1 binary64) n) Initial program 19.3%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified77.7%
Final simplification84.4%
(FPCore (x n)
:precision binary64
(let* ((t_0 (pow x (/ 1.0 n))))
(if (<= (/ 1.0 n) -4e-14)
(/ t_0 (* n x))
(if (<= (/ 1.0 n) 1e-16)
(/ (log (/ x (+ 1.0 x))) (- n))
(if (<= (/ 1.0 n) 2e+165) (- 1.0 t_0) (+ (exp (/ x n)) -1.0))))))
double code(double x, double n) {
double t_0 = pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 2e+165) {
tmp = 1.0 - t_0;
} else {
tmp = exp((x / n)) + -1.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = x ** (1.0d0 / n)
if ((1.0d0 / n) <= (-4d-14)) then
tmp = t_0 / (n * x)
else if ((1.0d0 / n) <= 1d-16) then
tmp = log((x / (1.0d0 + x))) / -n
else if ((1.0d0 / n) <= 2d+165) then
tmp = 1.0d0 - t_0
else
tmp = exp((x / n)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x, double n) {
double t_0 = Math.pow(x, (1.0 / n));
double tmp;
if ((1.0 / n) <= -4e-14) {
tmp = t_0 / (n * x);
} else if ((1.0 / n) <= 1e-16) {
tmp = Math.log((x / (1.0 + x))) / -n;
} else if ((1.0 / n) <= 2e+165) {
tmp = 1.0 - t_0;
} else {
tmp = Math.exp((x / n)) + -1.0;
}
return tmp;
}
def code(x, n): t_0 = math.pow(x, (1.0 / n)) tmp = 0 if (1.0 / n) <= -4e-14: tmp = t_0 / (n * x) elif (1.0 / n) <= 1e-16: tmp = math.log((x / (1.0 + x))) / -n elif (1.0 / n) <= 2e+165: tmp = 1.0 - t_0 else: tmp = math.exp((x / n)) + -1.0 return tmp
function code(x, n) t_0 = x ^ Float64(1.0 / n) tmp = 0.0 if (Float64(1.0 / n) <= -4e-14) tmp = Float64(t_0 / Float64(n * x)); elseif (Float64(1.0 / n) <= 1e-16) tmp = Float64(log(Float64(x / Float64(1.0 + x))) / Float64(-n)); elseif (Float64(1.0 / n) <= 2e+165) tmp = Float64(1.0 - t_0); else tmp = Float64(exp(Float64(x / n)) + -1.0); end return tmp end
function tmp_2 = code(x, n) t_0 = x ^ (1.0 / n); tmp = 0.0; if ((1.0 / n) <= -4e-14) tmp = t_0 / (n * x); elseif ((1.0 / n) <= 1e-16) tmp = log((x / (1.0 + x))) / -n; elseif ((1.0 / n) <= 2e+165) tmp = 1.0 - t_0; else tmp = exp((x / n)) + -1.0; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 / n), $MachinePrecision], -4e-14], N[(t$95$0 / N[(n * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 1e-16], N[(N[Log[N[(x / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / (-n)), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 2e+165], N[(1.0 - t$95$0), $MachinePrecision], N[(N[Exp[N[(x / n), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{if}\;\frac{1}{n} \leq -4 \cdot 10^{-14}:\\
\;\;\;\;\frac{t\_0}{n \cdot x}\\
\mathbf{elif}\;\frac{1}{n} \leq 10^{-16}:\\
\;\;\;\;\frac{\log \left(\frac{x}{1 + x}\right)}{-n}\\
\mathbf{elif}\;\frac{1}{n} \leq 2 \cdot 10^{+165}:\\
\;\;\;\;1 - t\_0\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{x}{n}} + -1\\
\end{array}
\end{array}
if (/.f64 #s(literal 1 binary64) n) < -4e-14Initial program 95.5%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6499.9
Simplified99.9%
if -4e-14 < (/.f64 #s(literal 1 binary64) n) < 9.9999999999999998e-17Initial program 31.6%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6477.0
Simplified77.0%
diff-logN/A
clear-numN/A
log-recN/A
lower-neg.f64N/A
lower-log.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f6477.2
Applied egg-rr77.2%
if 9.9999999999999998e-17 < (/.f64 #s(literal 1 binary64) n) < 1.9999999999999998e165Initial program 81.7%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6476.6
Simplified76.6%
if 1.9999999999999998e165 < (/.f64 #s(literal 1 binary64) n) Initial program 27.4%
lift-+.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-exp.f64N/A
lift-/.f64N/A
un-div-invN/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-log1p.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around 0
lower-/.f64100.0
Simplified100.0%
Taylor expanded in n around inf
Simplified71.2%
Final simplification84.3%
(FPCore (x n)
:precision binary64
(if (<= x 5e-235)
(- 1.0 (pow x (/ 1.0 n)))
(if (<= x 3.3e-174)
(- (/ (log x) n))
(if (<= x 8e-61)
(+
(/ 1.0 (* n x))
(/ (+ (/ 0.3333333333333333 (* n x)) (/ -0.5 n)) (* x x)))
(if (<= x 0.88)
(/ 1.0 (/ n (- x (log x))))
(if (<= x 6.4e+81)
(/
(/
(+
1.0
(-
(/ (+ -0.5 (/ 0.3333333333333333 x)) x)
(/ 0.25 (* x (* x x)))))
n)
x)
0.0))))))
double code(double x, double n) {
double tmp;
if (x <= 5e-235) {
tmp = 1.0 - pow(x, (1.0 / n));
} else if (x <= 3.3e-174) {
tmp = -(log(x) / n);
} else if (x <= 8e-61) {
tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = 1.0 / (n / (x - log(x)));
} else if (x <= 6.4e+81) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 5d-235) then
tmp = 1.0d0 - (x ** (1.0d0 / n))
else if (x <= 3.3d-174) then
tmp = -(log(x) / n)
else if (x <= 8d-61) then
tmp = (1.0d0 / (n * x)) + (((0.3333333333333333d0 / (n * x)) + ((-0.5d0) / n)) / (x * x))
else if (x <= 0.88d0) then
tmp = 1.0d0 / (n / (x - log(x)))
else if (x <= 6.4d+81) then
tmp = ((1.0d0 + ((((-0.5d0) + (0.3333333333333333d0 / x)) / x) - (0.25d0 / (x * (x * x))))) / n) / x
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 5e-235) {
tmp = 1.0 - Math.pow(x, (1.0 / n));
} else if (x <= 3.3e-174) {
tmp = -(Math.log(x) / n);
} else if (x <= 8e-61) {
tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = 1.0 / (n / (x - Math.log(x)));
} else if (x <= 6.4e+81) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 5e-235: tmp = 1.0 - math.pow(x, (1.0 / n)) elif x <= 3.3e-174: tmp = -(math.log(x) / n) elif x <= 8e-61: tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)) elif x <= 0.88: tmp = 1.0 / (n / (x - math.log(x))) elif x <= 6.4e+81: tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 5e-235) tmp = Float64(1.0 - (x ^ Float64(1.0 / n))); elseif (x <= 3.3e-174) tmp = Float64(-Float64(log(x) / n)); elseif (x <= 8e-61) tmp = Float64(Float64(1.0 / Float64(n * x)) + Float64(Float64(Float64(0.3333333333333333 / Float64(n * x)) + Float64(-0.5 / n)) / Float64(x * x))); elseif (x <= 0.88) tmp = Float64(1.0 / Float64(n / Float64(x - log(x)))); elseif (x <= 6.4e+81) tmp = Float64(Float64(Float64(1.0 + Float64(Float64(Float64(-0.5 + Float64(0.3333333333333333 / x)) / x) - Float64(0.25 / Float64(x * Float64(x * x))))) / n) / x); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 5e-235) tmp = 1.0 - (x ^ (1.0 / n)); elseif (x <= 3.3e-174) tmp = -(log(x) / n); elseif (x <= 8e-61) tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)); elseif (x <= 0.88) tmp = 1.0 / (n / (x - log(x))); elseif (x <= 6.4e+81) tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 5e-235], N[(1.0 - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.3e-174], (-N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision]), If[LessEqual[x, 8e-61], N[(N[(1.0 / N[(n * x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.3333333333333333 / N[(n * x), $MachinePrecision]), $MachinePrecision] + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 0.88], N[(1.0 / N[(n / N[(x - N[Log[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 6.4e+81], N[(N[(N[(1.0 + N[(N[(N[(-0.5 + N[(0.3333333333333333 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] - N[(0.25 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] / x), $MachinePrecision], 0.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5 \cdot 10^{-235}:\\
\;\;\;\;1 - {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{elif}\;x \leq 3.3 \cdot 10^{-174}:\\
\;\;\;\;-\frac{\log x}{n}\\
\mathbf{elif}\;x \leq 8 \cdot 10^{-61}:\\
\;\;\;\;\frac{1}{n \cdot x} + \frac{\frac{0.3333333333333333}{n \cdot x} + \frac{-0.5}{n}}{x \cdot x}\\
\mathbf{elif}\;x \leq 0.88:\\
\;\;\;\;\frac{1}{\frac{n}{x - \log x}}\\
\mathbf{elif}\;x \leq 6.4 \cdot 10^{+81}:\\
\;\;\;\;\frac{\frac{1 + \left(\frac{-0.5 + \frac{0.3333333333333333}{x}}{x} - \frac{0.25}{x \cdot \left(x \cdot x\right)}\right)}{n}}{x}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 4.9999999999999998e-235Initial program 67.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6467.6
Simplified67.6%
if 4.9999999999999998e-235 < x < 3.3000000000000001e-174Initial program 25.4%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6425.4
Simplified25.4%
Taylor expanded in n around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-log.f6480.1
Simplified80.1%
if 3.3000000000000001e-174 < x < 8.0000000000000003e-61Initial program 54.0%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6436.5
Simplified36.5%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f64N/A
div-subN/A
associate-/r*N/A
lower--.f64N/A
Simplified57.4%
if 8.0000000000000003e-61 < x < 0.880000000000000004Initial program 38.8%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified25.8%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6435.3
Simplified35.3%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log.f6448.1
Simplified48.1%
lift-log.f64N/A
lift--.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6448.1
Applied egg-rr48.1%
if 0.880000000000000004 < x < 6.4e81Initial program 39.0%
Taylor expanded in x around inf
Simplified79.8%
Taylor expanded in n around inf
lower-/.f64N/A
Simplified70.0%
if 6.4e81 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
Final simplification68.0%
(FPCore (x n)
:precision binary64
(if (<= x 4.8e-235)
(- 1.0 (pow x (/ 1.0 n)))
(if (<= x 5.2e-174)
(- (/ (log x) n))
(if (<= x 7.5e-61)
(+
(/ 1.0 (* n x))
(/ (+ (/ 0.3333333333333333 (* n x)) (/ -0.5 n)) (* x x)))
(if (<= x 0.88)
(/ (- x (log x)) n)
(if (<= x 3.7e+82)
(/
(/
(+
1.0
(-
(/ (+ -0.5 (/ 0.3333333333333333 x)) x)
(/ 0.25 (* x (* x x)))))
n)
x)
0.0))))))
double code(double x, double n) {
double tmp;
if (x <= 4.8e-235) {
tmp = 1.0 - pow(x, (1.0 / n));
} else if (x <= 5.2e-174) {
tmp = -(log(x) / n);
} else if (x <= 7.5e-61) {
tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = (x - log(x)) / n;
} else if (x <= 3.7e+82) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 4.8d-235) then
tmp = 1.0d0 - (x ** (1.0d0 / n))
else if (x <= 5.2d-174) then
tmp = -(log(x) / n)
else if (x <= 7.5d-61) then
tmp = (1.0d0 / (n * x)) + (((0.3333333333333333d0 / (n * x)) + ((-0.5d0) / n)) / (x * x))
else if (x <= 0.88d0) then
tmp = (x - log(x)) / n
else if (x <= 3.7d+82) then
tmp = ((1.0d0 + ((((-0.5d0) + (0.3333333333333333d0 / x)) / x) - (0.25d0 / (x * (x * x))))) / n) / x
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 4.8e-235) {
tmp = 1.0 - Math.pow(x, (1.0 / n));
} else if (x <= 5.2e-174) {
tmp = -(Math.log(x) / n);
} else if (x <= 7.5e-61) {
tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = (x - Math.log(x)) / n;
} else if (x <= 3.7e+82) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 4.8e-235: tmp = 1.0 - math.pow(x, (1.0 / n)) elif x <= 5.2e-174: tmp = -(math.log(x) / n) elif x <= 7.5e-61: tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)) elif x <= 0.88: tmp = (x - math.log(x)) / n elif x <= 3.7e+82: tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 4.8e-235) tmp = Float64(1.0 - (x ^ Float64(1.0 / n))); elseif (x <= 5.2e-174) tmp = Float64(-Float64(log(x) / n)); elseif (x <= 7.5e-61) tmp = Float64(Float64(1.0 / Float64(n * x)) + Float64(Float64(Float64(0.3333333333333333 / Float64(n * x)) + Float64(-0.5 / n)) / Float64(x * x))); elseif (x <= 0.88) tmp = Float64(Float64(x - log(x)) / n); elseif (x <= 3.7e+82) tmp = Float64(Float64(Float64(1.0 + Float64(Float64(Float64(-0.5 + Float64(0.3333333333333333 / x)) / x) - Float64(0.25 / Float64(x * Float64(x * x))))) / n) / x); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 4.8e-235) tmp = 1.0 - (x ^ (1.0 / n)); elseif (x <= 5.2e-174) tmp = -(log(x) / n); elseif (x <= 7.5e-61) tmp = (1.0 / (n * x)) + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)); elseif (x <= 0.88) tmp = (x - log(x)) / n; elseif (x <= 3.7e+82) tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 4.8e-235], N[(1.0 - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 5.2e-174], (-N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision]), If[LessEqual[x, 7.5e-61], N[(N[(1.0 / N[(n * x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.3333333333333333 / N[(n * x), $MachinePrecision]), $MachinePrecision] + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 0.88], N[(N[(x - N[Log[x], $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision], If[LessEqual[x, 3.7e+82], N[(N[(N[(1.0 + N[(N[(N[(-0.5 + N[(0.3333333333333333 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] - N[(0.25 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] / x), $MachinePrecision], 0.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.8 \cdot 10^{-235}:\\
\;\;\;\;1 - {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{elif}\;x \leq 5.2 \cdot 10^{-174}:\\
\;\;\;\;-\frac{\log x}{n}\\
\mathbf{elif}\;x \leq 7.5 \cdot 10^{-61}:\\
\;\;\;\;\frac{1}{n \cdot x} + \frac{\frac{0.3333333333333333}{n \cdot x} + \frac{-0.5}{n}}{x \cdot x}\\
\mathbf{elif}\;x \leq 0.88:\\
\;\;\;\;\frac{x - \log x}{n}\\
\mathbf{elif}\;x \leq 3.7 \cdot 10^{+82}:\\
\;\;\;\;\frac{\frac{1 + \left(\frac{-0.5 + \frac{0.3333333333333333}{x}}{x} - \frac{0.25}{x \cdot \left(x \cdot x\right)}\right)}{n}}{x}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 4.80000000000000022e-235Initial program 67.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6467.6
Simplified67.6%
if 4.80000000000000022e-235 < x < 5.2000000000000004e-174Initial program 25.4%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6425.4
Simplified25.4%
Taylor expanded in n around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-log.f6480.1
Simplified80.1%
if 5.2000000000000004e-174 < x < 7.50000000000000047e-61Initial program 54.0%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6436.5
Simplified36.5%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f64N/A
div-subN/A
associate-/r*N/A
lower--.f64N/A
Simplified57.4%
if 7.50000000000000047e-61 < x < 0.880000000000000004Initial program 38.8%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified25.8%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6435.3
Simplified35.3%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log.f6448.1
Simplified48.1%
if 0.880000000000000004 < x < 3.7000000000000002e82Initial program 39.0%
Taylor expanded in x around inf
Simplified79.8%
Taylor expanded in n around inf
lower-/.f64N/A
Simplified70.0%
if 3.7000000000000002e82 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
Final simplification68.0%
(FPCore (x n)
:precision binary64
(let* ((t_0 (/ 1.0 (* n x))))
(if (<= x 8.5e-265)
t_0
(if (<= x 5.2e-174)
(- (/ (log x) n))
(if (<= x 7.5e-61)
(+ t_0 (/ (+ (/ 0.3333333333333333 (* n x)) (/ -0.5 n)) (* x x)))
(if (<= x 0.88)
(/ (- x (log x)) n)
(if (<= x 3.4e+82)
(/
(/
(+
1.0
(-
(/ (+ -0.5 (/ 0.3333333333333333 x)) x)
(/ 0.25 (* x (* x x)))))
n)
x)
0.0)))))))
double code(double x, double n) {
double t_0 = 1.0 / (n * x);
double tmp;
if (x <= 8.5e-265) {
tmp = t_0;
} else if (x <= 5.2e-174) {
tmp = -(log(x) / n);
} else if (x <= 7.5e-61) {
tmp = t_0 + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = (x - log(x)) / n;
} else if (x <= 3.4e+82) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 / (n * x)
if (x <= 8.5d-265) then
tmp = t_0
else if (x <= 5.2d-174) then
tmp = -(log(x) / n)
else if (x <= 7.5d-61) then
tmp = t_0 + (((0.3333333333333333d0 / (n * x)) + ((-0.5d0) / n)) / (x * x))
else if (x <= 0.88d0) then
tmp = (x - log(x)) / n
else if (x <= 3.4d+82) then
tmp = ((1.0d0 + ((((-0.5d0) + (0.3333333333333333d0 / x)) / x) - (0.25d0 / (x * (x * x))))) / n) / x
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double t_0 = 1.0 / (n * x);
double tmp;
if (x <= 8.5e-265) {
tmp = t_0;
} else if (x <= 5.2e-174) {
tmp = -(Math.log(x) / n);
} else if (x <= 7.5e-61) {
tmp = t_0 + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x));
} else if (x <= 0.88) {
tmp = (x - Math.log(x)) / n;
} else if (x <= 3.4e+82) {
tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): t_0 = 1.0 / (n * x) tmp = 0 if x <= 8.5e-265: tmp = t_0 elif x <= 5.2e-174: tmp = -(math.log(x) / n) elif x <= 7.5e-61: tmp = t_0 + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)) elif x <= 0.88: tmp = (x - math.log(x)) / n elif x <= 3.4e+82: tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x else: tmp = 0.0 return tmp
function code(x, n) t_0 = Float64(1.0 / Float64(n * x)) tmp = 0.0 if (x <= 8.5e-265) tmp = t_0; elseif (x <= 5.2e-174) tmp = Float64(-Float64(log(x) / n)); elseif (x <= 7.5e-61) tmp = Float64(t_0 + Float64(Float64(Float64(0.3333333333333333 / Float64(n * x)) + Float64(-0.5 / n)) / Float64(x * x))); elseif (x <= 0.88) tmp = Float64(Float64(x - log(x)) / n); elseif (x <= 3.4e+82) tmp = Float64(Float64(Float64(1.0 + Float64(Float64(Float64(-0.5 + Float64(0.3333333333333333 / x)) / x) - Float64(0.25 / Float64(x * Float64(x * x))))) / n) / x); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) t_0 = 1.0 / (n * x); tmp = 0.0; if (x <= 8.5e-265) tmp = t_0; elseif (x <= 5.2e-174) tmp = -(log(x) / n); elseif (x <= 7.5e-61) tmp = t_0 + (((0.3333333333333333 / (n * x)) + (-0.5 / n)) / (x * x)); elseif (x <= 0.88) tmp = (x - log(x)) / n; elseif (x <= 3.4e+82) tmp = ((1.0 + (((-0.5 + (0.3333333333333333 / x)) / x) - (0.25 / (x * (x * x))))) / n) / x; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := Block[{t$95$0 = N[(1.0 / N[(n * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 8.5e-265], t$95$0, If[LessEqual[x, 5.2e-174], (-N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision]), If[LessEqual[x, 7.5e-61], N[(t$95$0 + N[(N[(N[(0.3333333333333333 / N[(n * x), $MachinePrecision]), $MachinePrecision] + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 0.88], N[(N[(x - N[Log[x], $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision], If[LessEqual[x, 3.4e+82], N[(N[(N[(1.0 + N[(N[(N[(-0.5 + N[(0.3333333333333333 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] - N[(0.25 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] / x), $MachinePrecision], 0.0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{n \cdot x}\\
\mathbf{if}\;x \leq 8.5 \cdot 10^{-265}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 5.2 \cdot 10^{-174}:\\
\;\;\;\;-\frac{\log x}{n}\\
\mathbf{elif}\;x \leq 7.5 \cdot 10^{-61}:\\
\;\;\;\;t\_0 + \frac{\frac{0.3333333333333333}{n \cdot x} + \frac{-0.5}{n}}{x \cdot x}\\
\mathbf{elif}\;x \leq 0.88:\\
\;\;\;\;\frac{x - \log x}{n}\\
\mathbf{elif}\;x \leq 3.4 \cdot 10^{+82}:\\
\;\;\;\;\frac{\frac{1 + \left(\frac{-0.5 + \frac{0.3333333333333333}{x}}{x} - \frac{0.25}{x \cdot \left(x \cdot x\right)}\right)}{n}}{x}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 8.4999999999999997e-265Initial program 74.7%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6439.2
Simplified39.2%
Taylor expanded in n around inf
lower-/.f64N/A
*-commutativeN/A
lower-*.f6447.9
Simplified47.9%
if 8.4999999999999997e-265 < x < 5.2000000000000004e-174Initial program 36.9%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6436.9
Simplified36.9%
Taylor expanded in n around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-log.f6468.4
Simplified68.4%
if 5.2000000000000004e-174 < x < 7.50000000000000047e-61Initial program 54.0%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6436.5
Simplified36.5%
Taylor expanded in x around -inf
mul-1-negN/A
lower-neg.f64N/A
div-subN/A
associate-/r*N/A
lower--.f64N/A
Simplified57.4%
if 7.50000000000000047e-61 < x < 0.880000000000000004Initial program 38.8%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
Simplified25.8%
Taylor expanded in x around 0
+-commutativeN/A
lower-+.f64N/A
lower-/.f6435.3
Simplified35.3%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log.f6448.1
Simplified48.1%
if 0.880000000000000004 < x < 3.39999999999999994e82Initial program 39.0%
Taylor expanded in x around inf
Simplified79.8%
Taylor expanded in n around inf
lower-/.f64N/A
Simplified70.0%
if 3.39999999999999994e82 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
Final simplification65.0%
(FPCore (x n)
:precision binary64
(if (<= x 9.2e-265)
(/ 1.0 (* n x))
(if (<= x 8.5e-175)
(- (/ (log x) n))
(if (<= x 2.3e+81)
(/ (/ (+ 1.0 (/ (+ -0.5 (/ 0.3333333333333333 x)) x)) x) n)
0.0))))
double code(double x, double n) {
double tmp;
if (x <= 9.2e-265) {
tmp = 1.0 / (n * x);
} else if (x <= 8.5e-175) {
tmp = -(log(x) / n);
} else if (x <= 2.3e+81) {
tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 9.2d-265) then
tmp = 1.0d0 / (n * x)
else if (x <= 8.5d-175) then
tmp = -(log(x) / n)
else if (x <= 2.3d+81) then
tmp = ((1.0d0 + (((-0.5d0) + (0.3333333333333333d0 / x)) / x)) / x) / n
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 9.2e-265) {
tmp = 1.0 / (n * x);
} else if (x <= 8.5e-175) {
tmp = -(Math.log(x) / n);
} else if (x <= 2.3e+81) {
tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 9.2e-265: tmp = 1.0 / (n * x) elif x <= 8.5e-175: tmp = -(math.log(x) / n) elif x <= 2.3e+81: tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 9.2e-265) tmp = Float64(1.0 / Float64(n * x)); elseif (x <= 8.5e-175) tmp = Float64(-Float64(log(x) / n)); elseif (x <= 2.3e+81) tmp = Float64(Float64(Float64(1.0 + Float64(Float64(-0.5 + Float64(0.3333333333333333 / x)) / x)) / x) / n); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 9.2e-265) tmp = 1.0 / (n * x); elseif (x <= 8.5e-175) tmp = -(log(x) / n); elseif (x <= 2.3e+81) tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 9.2e-265], N[(1.0 / N[(n * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 8.5e-175], (-N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision]), If[LessEqual[x, 2.3e+81], N[(N[(N[(1.0 + N[(N[(-0.5 + N[(0.3333333333333333 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] / n), $MachinePrecision], 0.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.2 \cdot 10^{-265}:\\
\;\;\;\;\frac{1}{n \cdot x}\\
\mathbf{elif}\;x \leq 8.5 \cdot 10^{-175}:\\
\;\;\;\;-\frac{\log x}{n}\\
\mathbf{elif}\;x \leq 2.3 \cdot 10^{+81}:\\
\;\;\;\;\frac{\frac{1 + \frac{-0.5 + \frac{0.3333333333333333}{x}}{x}}{x}}{n}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 9.1999999999999996e-265Initial program 74.7%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6439.2
Simplified39.2%
Taylor expanded in n around inf
lower-/.f64N/A
*-commutativeN/A
lower-*.f6447.9
Simplified47.9%
if 9.1999999999999996e-265 < x < 8.5000000000000005e-175Initial program 36.9%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6436.9
Simplified36.9%
Taylor expanded in n around inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
lower-log.f6468.4
Simplified68.4%
if 8.5000000000000005e-175 < x < 2.2999999999999999e81Initial program 45.8%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6441.2
Simplified41.2%
Taylor expanded in x around inf
lower-/.f64N/A
Simplified50.6%
if 2.2999999999999999e81 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
Final simplification61.6%
(FPCore (x n) :precision binary64 (if (<= x 2.5e+82) (/ (/ (+ 1.0 (/ (+ -0.5 (/ 0.3333333333333333 x)) x)) x) n) 0.0))
double code(double x, double n) {
double tmp;
if (x <= 2.5e+82) {
tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 2.5d+82) then
tmp = ((1.0d0 + (((-0.5d0) + (0.3333333333333333d0 / x)) / x)) / x) / n
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 2.5e+82) {
tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 2.5e+82: tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 2.5e+82) tmp = Float64(Float64(Float64(1.0 + Float64(Float64(-0.5 + Float64(0.3333333333333333 / x)) / x)) / x) / n); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 2.5e+82) tmp = ((1.0 + ((-0.5 + (0.3333333333333333 / x)) / x)) / x) / n; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 2.5e+82], N[(N[(N[(1.0 + N[(N[(-0.5 + N[(0.3333333333333333 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] / n), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.5 \cdot 10^{+82}:\\
\;\;\;\;\frac{\frac{1 + \frac{-0.5 + \frac{0.3333333333333333}{x}}{x}}{x}}{n}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 2.50000000000000008e82Initial program 48.1%
Taylor expanded in n around inf
lower-/.f64N/A
lower--.f64N/A
lower-log1p.f64N/A
lower-log.f6444.2
Simplified44.2%
Taylor expanded in x around inf
lower-/.f64N/A
Simplified42.8%
if 2.50000000000000008e82 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
(FPCore (x n) :precision binary64 (if (<= x 5.8e+81) (/ (/ 1.0 n) x) 0.0))
double code(double x, double n) {
double tmp;
if (x <= 5.8e+81) {
tmp = (1.0 / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 5.8d+81) then
tmp = (1.0d0 / n) / x
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 5.8e+81) {
tmp = (1.0 / n) / x;
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 5.8e+81: tmp = (1.0 / n) / x else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 5.8e+81) tmp = Float64(Float64(1.0 / n) / x); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 5.8e+81) tmp = (1.0 / n) / x; else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 5.8e+81], N[(N[(1.0 / n), $MachinePrecision] / x), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+81}:\\
\;\;\;\;\frac{\frac{1}{n}}{x}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 5.7999999999999999e81Initial program 48.1%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6442.9
Simplified42.9%
lift-/.f64N/A
lift-pow.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6443.8
Applied egg-rr43.8%
Taylor expanded in n around inf
lower-/.f6433.8
Simplified33.8%
if 5.7999999999999999e81 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
(FPCore (x n) :precision binary64 (if (<= x 7.5e+82) (/ 1.0 (* n x)) 0.0))
double code(double x, double n) {
double tmp;
if (x <= 7.5e+82) {
tmp = 1.0 / (n * x);
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
real(8) :: tmp
if (x <= 7.5d+82) then
tmp = 1.0d0 / (n * x)
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double x, double n) {
double tmp;
if (x <= 7.5e+82) {
tmp = 1.0 / (n * x);
} else {
tmp = 0.0;
}
return tmp;
}
def code(x, n): tmp = 0 if x <= 7.5e+82: tmp = 1.0 / (n * x) else: tmp = 0.0 return tmp
function code(x, n) tmp = 0.0 if (x <= 7.5e+82) tmp = Float64(1.0 / Float64(n * x)); else tmp = 0.0; end return tmp end
function tmp_2 = code(x, n) tmp = 0.0; if (x <= 7.5e+82) tmp = 1.0 / (n * x); else tmp = 0.0; end tmp_2 = tmp; end
code[x_, n_] := If[LessEqual[x, 7.5e+82], N[(1.0 / N[(n * x), $MachinePrecision]), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 7.5 \cdot 10^{+82}:\\
\;\;\;\;\frac{1}{n \cdot x}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if x < 7.4999999999999999e82Initial program 48.1%
Taylor expanded in x around inf
lower-/.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6442.9
Simplified42.9%
Taylor expanded in n around inf
lower-/.f64N/A
*-commutativeN/A
lower-*.f6433.2
Simplified33.2%
if 7.4999999999999999e82 < x Initial program 77.6%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.8
Simplified42.8%
Taylor expanded in n around inf
Simplified77.6%
metadata-eval77.6
Applied egg-rr77.6%
Final simplification47.3%
(FPCore (x n) :precision binary64 0.0)
double code(double x, double n) {
return 0.0;
}
real(8) function code(x, n)
real(8), intent (in) :: x
real(8), intent (in) :: n
code = 0.0d0
end function
public static double code(double x, double n) {
return 0.0;
}
def code(x, n): return 0.0
function code(x, n) return 0.0 end
function tmp = code(x, n) tmp = 0.0; end
code[x_, n_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 57.4%
Taylor expanded in x around 0
remove-double-negN/A
mul-1-negN/A
distribute-neg-fracN/A
mul-1-negN/A
log-recN/A
mul-1-negN/A
lower--.f64N/A
log-recN/A
mul-1-negN/A
associate-*r/N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
associate-/l*N/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f6442.9
Simplified42.9%
Taylor expanded in n around inf
Simplified31.2%
metadata-eval31.2
Applied egg-rr31.2%
herbie shell --seed 2024207
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))