
(FPCore (x eps) :precision binary64 (- (pow (+ x eps) 5.0) (pow x 5.0)))
double code(double x, double eps) {
return pow((x + eps), 5.0) - pow(x, 5.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
end function
public static double code(double x, double eps) {
return Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
}
def code(x, eps): return math.pow((x + eps), 5.0) - math.pow(x, 5.0)
function code(x, eps) return Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) end
function tmp = code(x, eps) tmp = ((x + eps) ^ 5.0) - (x ^ 5.0); end
code[x_, eps_] := N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
{\left(x + \varepsilon\right)}^{5} - {x}^{5}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x eps) :precision binary64 (- (pow (+ x eps) 5.0) (pow x 5.0)))
double code(double x, double eps) {
return pow((x + eps), 5.0) - pow(x, 5.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
end function
public static double code(double x, double eps) {
return Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
}
def code(x, eps): return math.pow((x + eps), 5.0) - math.pow(x, 5.0)
function code(x, eps) return Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) end
function tmp = code(x, eps) tmp = ((x + eps) ^ 5.0) - (x ^ 5.0); end
code[x_, eps_] := N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
{\left(x + \varepsilon\right)}^{5} - {x}^{5}
\end{array}
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (+ x eps) 5.0))
(t_1 (- t_0 (pow x 5.0)))
(t_2 (- t_0 (* (* x x) (* x (* x x))))))
(if (<= t_1 -2e-304)
t_2
(if (<= t_1 0.0) (* (pow x 4.0) (* eps 5.0)) t_2))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0);
double t_1 = t_0 - pow(x, 5.0);
double t_2 = t_0 - ((x * x) * (x * (x * x)));
double tmp;
if (t_1 <= -2e-304) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = pow(x, 4.0) * (eps * 5.0);
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = (x + eps) ** 5.0d0
t_1 = t_0 - (x ** 5.0d0)
t_2 = t_0 - ((x * x) * (x * (x * x)))
if (t_1 <= (-2d-304)) then
tmp = t_2
else if (t_1 <= 0.0d0) then
tmp = (x ** 4.0d0) * (eps * 5.0d0)
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double eps) {
double t_0 = Math.pow((x + eps), 5.0);
double t_1 = t_0 - Math.pow(x, 5.0);
double t_2 = t_0 - ((x * x) * (x * (x * x)));
double tmp;
if (t_1 <= -2e-304) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = Math.pow(x, 4.0) * (eps * 5.0);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, eps): t_0 = math.pow((x + eps), 5.0) t_1 = t_0 - math.pow(x, 5.0) t_2 = t_0 - ((x * x) * (x * (x * x))) tmp = 0 if t_1 <= -2e-304: tmp = t_2 elif t_1 <= 0.0: tmp = math.pow(x, 4.0) * (eps * 5.0) else: tmp = t_2 return tmp
function code(x, eps) t_0 = Float64(x + eps) ^ 5.0 t_1 = Float64(t_0 - (x ^ 5.0)) t_2 = Float64(t_0 - Float64(Float64(x * x) * Float64(x * Float64(x * x)))) tmp = 0.0 if (t_1 <= -2e-304) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64((x ^ 4.0) * Float64(eps * 5.0)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, eps) t_0 = (x + eps) ^ 5.0; t_1 = t_0 - (x ^ 5.0); t_2 = t_0 - ((x * x) * (x * (x * x))); tmp = 0.0; if (t_1 <= -2e-304) tmp = t_2; elseif (t_1 <= 0.0) tmp = (x ^ 4.0) * (eps * 5.0); else tmp = t_2; end tmp_2 = tmp; end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 - N[(N[(x * x), $MachinePrecision] * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e-304], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[Power[x, 4.0], $MachinePrecision] * N[(eps * 5.0), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5}\\
t_1 := t\_0 - {x}^{5}\\
t_2 := t\_0 - \left(x \cdot x\right) \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;{x}^{4} \cdot \left(\varepsilon \cdot 5\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in eps around 0
*-lowering-*.f64100.0
Simplified100.0%
Final simplification99.3%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(/
1.0
(/
1.0
(fma
(* (* eps eps) (* eps eps))
(fma 5.0 x eps)
(* (* eps eps) (* (* x 10.0) (* x (+ x eps)))))))
(if (<= t_0 0.0) (* (pow x 4.0) (* eps 5.0)) (pow eps 5.0)))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = 1.0 / (1.0 / fma(((eps * eps) * (eps * eps)), fma(5.0, x, eps), ((eps * eps) * ((x * 10.0) * (x * (x + eps))))));
} else if (t_0 <= 0.0) {
tmp = pow(x, 4.0) * (eps * 5.0);
} else {
tmp = pow(eps, 5.0);
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(1.0 / Float64(1.0 / fma(Float64(Float64(eps * eps) * Float64(eps * eps)), fma(5.0, x, eps), Float64(Float64(eps * eps) * Float64(Float64(x * 10.0) * Float64(x * Float64(x + eps))))))); elseif (t_0 <= 0.0) tmp = Float64((x ^ 4.0) * Float64(eps * 5.0)); else tmp = eps ^ 5.0; end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(1.0 / N[(1.0 / N[(N[(N[(eps * eps), $MachinePrecision] * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(5.0 * x + eps), $MachinePrecision] + N[(N[(eps * eps), $MachinePrecision] * N[(N[(x * 10.0), $MachinePrecision] * N[(x * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[Power[x, 4.0], $MachinePrecision] * N[(eps * 5.0), $MachinePrecision]), $MachinePrecision], N[Power[eps, 5.0], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\frac{1}{\frac{1}{\mathsf{fma}\left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right), \mathsf{fma}\left(5, x, \varepsilon\right), \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(x \cdot 10\right) \cdot \left(x \cdot \left(x + \varepsilon\right)\right)\right)\right)}}\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;{x}^{4} \cdot \left(\varepsilon \cdot 5\right)\\
\mathbf{else}:\\
\;\;\;\;{\varepsilon}^{5}\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
flip3-+N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr95.5%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in eps around 0
*-lowering-*.f64100.0
Simplified100.0%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
Taylor expanded in x around 0
pow-lowering-pow.f6491.6
Simplified91.6%
Final simplification98.8%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(/
1.0
(/
1.0
(fma
(* (* eps eps) (* eps eps))
(fma 5.0 x eps)
(* (* eps eps) (* (* x 10.0) (* x (+ x eps)))))))
(if (<= t_0 0.0)
(* (fma eps 5.0 (/ (* (* eps eps) 10.0) x)) (* x (* x (* x x))))
(pow eps 5.0)))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = 1.0 / (1.0 / fma(((eps * eps) * (eps * eps)), fma(5.0, x, eps), ((eps * eps) * ((x * 10.0) * (x * (x + eps))))));
} else if (t_0 <= 0.0) {
tmp = fma(eps, 5.0, (((eps * eps) * 10.0) / x)) * (x * (x * (x * x)));
} else {
tmp = pow(eps, 5.0);
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(1.0 / Float64(1.0 / fma(Float64(Float64(eps * eps) * Float64(eps * eps)), fma(5.0, x, eps), Float64(Float64(eps * eps) * Float64(Float64(x * 10.0) * Float64(x * Float64(x + eps))))))); elseif (t_0 <= 0.0) tmp = Float64(fma(eps, 5.0, Float64(Float64(Float64(eps * eps) * 10.0) / x)) * Float64(x * Float64(x * Float64(x * x)))); else tmp = eps ^ 5.0; end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(1.0 / N[(1.0 / N[(N[(N[(eps * eps), $MachinePrecision] * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(5.0 * x + eps), $MachinePrecision] + N[(N[(eps * eps), $MachinePrecision] * N[(N[(x * 10.0), $MachinePrecision] * N[(x * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(eps * 5.0 + N[(N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[eps, 5.0], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\frac{1}{\frac{1}{\mathsf{fma}\left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right), \mathsf{fma}\left(5, x, \varepsilon\right), \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(x \cdot 10\right) \cdot \left(x \cdot \left(x + \varepsilon\right)\right)\right)\right)}}\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(\varepsilon, 5, \frac{\left(\varepsilon \cdot \varepsilon\right) \cdot 10}{x}\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\varepsilon}^{5}\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
flip3-+N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr95.5%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-plusN/A
cube-unmultN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
Taylor expanded in x around 0
pow-lowering-pow.f6491.6
Simplified91.6%
Final simplification98.8%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(/
1.0
(/
1.0
(fma
(* (* eps eps) (* eps eps))
(fma 5.0 x eps)
(* (* eps eps) (* (* x 10.0) (* x (+ x eps)))))))
(if (<= t_0 0.0)
(* (fma eps 5.0 (/ (* (* eps eps) 10.0) x)) (* x (* x (* x x))))
(* (* eps eps) (* eps (fma eps (fma x 5.0 eps) (* (* x x) 10.0))))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = 1.0 / (1.0 / fma(((eps * eps) * (eps * eps)), fma(5.0, x, eps), ((eps * eps) * ((x * 10.0) * (x * (x + eps))))));
} else if (t_0 <= 0.0) {
tmp = fma(eps, 5.0, (((eps * eps) * 10.0) / x)) * (x * (x * (x * x)));
} else {
tmp = (eps * eps) * (eps * fma(eps, fma(x, 5.0, eps), ((x * x) * 10.0)));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(1.0 / Float64(1.0 / fma(Float64(Float64(eps * eps) * Float64(eps * eps)), fma(5.0, x, eps), Float64(Float64(eps * eps) * Float64(Float64(x * 10.0) * Float64(x * Float64(x + eps))))))); elseif (t_0 <= 0.0) tmp = Float64(fma(eps, 5.0, Float64(Float64(Float64(eps * eps) * 10.0) / x)) * Float64(x * Float64(x * Float64(x * x)))); else tmp = Float64(Float64(eps * eps) * Float64(eps * fma(eps, fma(x, 5.0, eps), Float64(Float64(x * x) * 10.0)))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(1.0 / N[(1.0 / N[(N[(N[(eps * eps), $MachinePrecision] * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(5.0 * x + eps), $MachinePrecision] + N[(N[(eps * eps), $MachinePrecision] * N[(N[(x * 10.0), $MachinePrecision] * N[(x * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(eps * 5.0 + N[(N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * N[(x * 5.0 + eps), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\frac{1}{\frac{1}{\mathsf{fma}\left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right), \mathsf{fma}\left(5, x, \varepsilon\right), \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(x \cdot 10\right) \cdot \left(x \cdot \left(x + \varepsilon\right)\right)\right)\right)}}\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(\varepsilon, 5, \frac{\left(\varepsilon \cdot \varepsilon\right) \cdot 10}{x}\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon, \mathsf{fma}\left(x, 5, \varepsilon\right), \left(x \cdot x\right) \cdot 10\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
flip3-+N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr95.5%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-plusN/A
cube-unmultN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-+r+N/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-inN/A
+-commutativeN/A
Simplified90.7%
Final simplification98.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(fma
(fma 5.0 x eps)
(* (* eps eps) (* eps eps))
(* (* eps eps) (* (* x 10.0) (* x (+ x eps)))))
(if (<= t_0 0.0)
(* (fma eps 5.0 (/ (* (* eps eps) 10.0) x)) (* x (* x (* x x))))
(* (* eps eps) (* eps (fma eps (fma x 5.0 eps) (* (* x x) 10.0))))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = fma(fma(5.0, x, eps), ((eps * eps) * (eps * eps)), ((eps * eps) * ((x * 10.0) * (x * (x + eps)))));
} else if (t_0 <= 0.0) {
tmp = fma(eps, 5.0, (((eps * eps) * 10.0) / x)) * (x * (x * (x * x)));
} else {
tmp = (eps * eps) * (eps * fma(eps, fma(x, 5.0, eps), ((x * x) * 10.0)));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = fma(fma(5.0, x, eps), Float64(Float64(eps * eps) * Float64(eps * eps)), Float64(Float64(eps * eps) * Float64(Float64(x * 10.0) * Float64(x * Float64(x + eps))))); elseif (t_0 <= 0.0) tmp = Float64(fma(eps, 5.0, Float64(Float64(Float64(eps * eps) * 10.0) / x)) * Float64(x * Float64(x * Float64(x * x)))); else tmp = Float64(Float64(eps * eps) * Float64(eps * fma(eps, fma(x, 5.0, eps), Float64(Float64(x * x) * 10.0)))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(N[(5.0 * x + eps), $MachinePrecision] * N[(N[(eps * eps), $MachinePrecision] * N[(eps * eps), $MachinePrecision]), $MachinePrecision] + N[(N[(eps * eps), $MachinePrecision] * N[(N[(x * 10.0), $MachinePrecision] * N[(x * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(eps * 5.0 + N[(N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * N[(x * 5.0 + eps), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(5, x, \varepsilon\right), \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right), \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(x \cdot 10\right) \cdot \left(x \cdot \left(x + \varepsilon\right)\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(\varepsilon, 5, \frac{\left(\varepsilon \cdot \varepsilon\right) \cdot 10}{x}\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon, \mathsf{fma}\left(x, 5, \varepsilon\right), \left(x \cdot x\right) \cdot 10\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
*-commutativeN/A
accelerator-lowering-fma.f64N/A
accelerator-lowering-fma.f64N/A
metadata-evalN/A
pow-plusN/A
unpow3N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6495.5
Applied egg-rr95.5%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-plusN/A
cube-unmultN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-+r+N/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-inN/A
+-commutativeN/A
Simplified90.7%
Final simplification98.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(*
eps
(* eps (fma (* eps eps) (fma 5.0 x eps) (* x (* x (* 10.0 (+ x eps)))))))
(if (<= t_0 0.0)
(* (fma eps 5.0 (/ (* (* eps eps) 10.0) x)) (* x (* x (* x x))))
(* (* eps eps) (* eps (fma eps (fma x 5.0 eps) (* (* x x) 10.0))))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = eps * (eps * fma((eps * eps), fma(5.0, x, eps), (x * (x * (10.0 * (x + eps))))));
} else if (t_0 <= 0.0) {
tmp = fma(eps, 5.0, (((eps * eps) * 10.0) / x)) * (x * (x * (x * x)));
} else {
tmp = (eps * eps) * (eps * fma(eps, fma(x, 5.0, eps), ((x * x) * 10.0)));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(eps * Float64(eps * fma(Float64(eps * eps), fma(5.0, x, eps), Float64(x * Float64(x * Float64(10.0 * Float64(x + eps))))))); elseif (t_0 <= 0.0) tmp = Float64(fma(eps, 5.0, Float64(Float64(Float64(eps * eps) * 10.0) / x)) * Float64(x * Float64(x * Float64(x * x)))); else tmp = Float64(Float64(eps * eps) * Float64(eps * fma(eps, fma(x, 5.0, eps), Float64(Float64(x * x) * 10.0)))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(eps * N[(eps * N[(N[(eps * eps), $MachinePrecision] * N[(5.0 * x + eps), $MachinePrecision] + N[(x * N[(x * N[(10.0 * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(eps * 5.0 + N[(N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * N[(x * 5.0 + eps), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\varepsilon \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon \cdot \varepsilon, \mathsf{fma}\left(5, x, \varepsilon\right), x \cdot \left(x \cdot \left(10 \cdot \left(x + \varepsilon\right)\right)\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(\varepsilon, 5, \frac{\left(\varepsilon \cdot \varepsilon\right) \cdot 10}{x}\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon, \mathsf{fma}\left(x, 5, \varepsilon\right), \left(x \cdot x\right) \cdot 10\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
metadata-evalN/A
pow-plusN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6495.4
Applied egg-rr95.4%
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-lft-inN/A
Applied egg-rr95.4%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-plusN/A
cube-unmultN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-+r+N/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-inN/A
+-commutativeN/A
Simplified90.7%
Final simplification98.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(*
eps
(* eps (fma (* eps eps) (fma 5.0 x eps) (* x (* x (* 10.0 (+ x eps)))))))
(if (<= t_0 0.0)
(* (* x (* x x)) (fma (* x eps) 5.0 (* (* eps eps) 10.0)))
(* (* eps eps) (* eps (fma eps (fma x 5.0 eps) (* (* x x) 10.0))))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = eps * (eps * fma((eps * eps), fma(5.0, x, eps), (x * (x * (10.0 * (x + eps))))));
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * fma((x * eps), 5.0, ((eps * eps) * 10.0));
} else {
tmp = (eps * eps) * (eps * fma(eps, fma(x, 5.0, eps), ((x * x) * 10.0)));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(eps * Float64(eps * fma(Float64(eps * eps), fma(5.0, x, eps), Float64(x * Float64(x * Float64(10.0 * Float64(x + eps))))))); elseif (t_0 <= 0.0) tmp = Float64(Float64(x * Float64(x * x)) * fma(Float64(x * eps), 5.0, Float64(Float64(eps * eps) * 10.0))); else tmp = Float64(Float64(eps * eps) * Float64(eps * fma(eps, fma(x, 5.0, eps), Float64(Float64(x * x) * 10.0)))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(eps * N[(eps * N[(N[(eps * eps), $MachinePrecision] * N[(5.0 * x + eps), $MachinePrecision] + N[(x * N[(x * N[(10.0 * N[(x + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(N[(x * eps), $MachinePrecision] * 5.0 + N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * N[(x * 5.0 + eps), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\varepsilon \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon \cdot \varepsilon, \mathsf{fma}\left(5, x, \varepsilon\right), x \cdot \left(x \cdot \left(10 \cdot \left(x + \varepsilon\right)\right)\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot \mathsf{fma}\left(x \cdot \varepsilon, 5, \left(\varepsilon \cdot \varepsilon\right) \cdot 10\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon, \mathsf{fma}\left(x, 5, \varepsilon\right), \left(x \cdot x\right) \cdot 10\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
metadata-evalN/A
pow-plusN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6495.4
Applied egg-rr95.4%
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
distribute-lft-inN/A
Applied egg-rr95.4%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
distribute-lft-inN/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-+r+N/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-inN/A
+-commutativeN/A
Simplified90.7%
Final simplification98.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0)))
(t_1
(* (* eps eps) (* eps (fma eps (fma x 5.0 eps) (* (* x x) 10.0))))))
(if (<= t_0 -2e-304)
t_1
(if (<= t_0 0.0)
(* (* x (* x x)) (fma (* x eps) 5.0 (* (* eps eps) 10.0)))
t_1))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double t_1 = (eps * eps) * (eps * fma(eps, fma(x, 5.0, eps), ((x * x) * 10.0)));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * fma((x * eps), 5.0, ((eps * eps) * 10.0));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) t_1 = Float64(Float64(eps * eps) * Float64(eps * fma(eps, fma(x, 5.0, eps), Float64(Float64(x * x) * 10.0)))) tmp = 0.0 if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(Float64(x * Float64(x * x)) * fma(Float64(x * eps), 5.0, Float64(Float64(eps * eps) * 10.0))); else tmp = t_1; end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * N[(x * 5.0 + eps), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], t$95$1, If[LessEqual[t$95$0, 0.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(N[(x * eps), $MachinePrecision] * 5.0 + N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
t_1 := \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(\varepsilon, \mathsf{fma}\left(x, 5, \varepsilon\right), \left(x \cdot x\right) \cdot 10\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot \mathsf{fma}\left(x \cdot \varepsilon, 5, \left(\varepsilon \cdot \varepsilon\right) \cdot 10\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
Taylor expanded in x around 0
Simplified93.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified92.7%
Taylor expanded in x around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-+r+N/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-inN/A
+-commutativeN/A
Simplified93.0%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
distribute-lft-inN/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
Final simplification98.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(* (* eps eps) (* (* eps eps) (fma x 5.0 eps)))
(if (<= t_0 0.0)
(* (* x (* x x)) (fma (* x eps) 5.0 (* (* eps eps) 10.0)))
(* (* eps eps) (* eps (* eps eps)))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = (eps * eps) * ((eps * eps) * fma(x, 5.0, eps));
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * fma((x * eps), 5.0, ((eps * eps) * 10.0));
} else {
tmp = (eps * eps) * (eps * (eps * eps));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(Float64(eps * eps) * Float64(Float64(eps * eps) * fma(x, 5.0, eps))); elseif (t_0 <= 0.0) tmp = Float64(Float64(x * Float64(x * x)) * fma(Float64(x * eps), 5.0, Float64(Float64(eps * eps) * 10.0))); else tmp = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(N[(eps * eps), $MachinePrecision] * N[(N[(eps * eps), $MachinePrecision] * N[(x * 5.0 + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(N[(x * eps), $MachinePrecision] * 5.0 + N[(N[(eps * eps), $MachinePrecision] * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \mathsf{fma}\left(x, 5, \varepsilon\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot \mathsf{fma}\left(x \cdot \varepsilon, 5, \left(\varepsilon \cdot \varepsilon\right) \cdot 10\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified95.1%
Taylor expanded in x around 0
*-commutativeN/A
associate-*r*N/A
cube-multN/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6493.9
Simplified93.9%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
distribute-lft-inN/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.6
Simplified90.6%
Final simplification98.6%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(* (* eps eps) (* (* eps eps) (fma x 5.0 eps)))
(if (<= t_0 0.0)
(* (* x (* x x)) (* eps (fma 5.0 x (* eps 10.0))))
(* (* eps eps) (* eps (* eps eps)))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = (eps * eps) * ((eps * eps) * fma(x, 5.0, eps));
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * (eps * fma(5.0, x, (eps * 10.0)));
} else {
tmp = (eps * eps) * (eps * (eps * eps));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(Float64(eps * eps) * Float64(Float64(eps * eps) * fma(x, 5.0, eps))); elseif (t_0 <= 0.0) tmp = Float64(Float64(x * Float64(x * x)) * Float64(eps * fma(5.0, x, Float64(eps * 10.0)))); else tmp = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(N[(eps * eps), $MachinePrecision] * N[(N[(eps * eps), $MachinePrecision] * N[(x * 5.0 + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(eps * N[(5.0 * x + N[(eps * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \mathsf{fma}\left(x, 5, \varepsilon\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(\varepsilon \cdot \mathsf{fma}\left(5, x, \varepsilon \cdot 10\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified95.1%
Taylor expanded in x around 0
*-commutativeN/A
associate-*r*N/A
cube-multN/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6493.9
Simplified93.9%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.6
Simplified90.6%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0))))
(if (<= t_0 -2e-304)
(* (* eps eps) (* (* eps eps) (fma x 5.0 eps)))
(if (<= t_0 0.0)
(* eps (* (* x (* x x)) (* x 5.0)))
(* (* eps eps) (* eps (* eps eps)))))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double tmp;
if (t_0 <= -2e-304) {
tmp = (eps * eps) * ((eps * eps) * fma(x, 5.0, eps));
} else if (t_0 <= 0.0) {
tmp = eps * ((x * (x * x)) * (x * 5.0));
} else {
tmp = (eps * eps) * (eps * (eps * eps));
}
return tmp;
}
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) tmp = 0.0 if (t_0 <= -2e-304) tmp = Float64(Float64(eps * eps) * Float64(Float64(eps * eps) * fma(x, 5.0, eps))); elseif (t_0 <= 0.0) tmp = Float64(eps * Float64(Float64(x * Float64(x * x)) * Float64(x * 5.0))); else tmp = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))); end return tmp end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], N[(N[(eps * eps), $MachinePrecision] * N[(N[(eps * eps), $MachinePrecision] * N[(x * 5.0 + eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[(eps * N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \mathsf{fma}\left(x, 5, \varepsilon\right)\right)\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\varepsilon \cdot \left(\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot 5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304Initial program 97.3%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.3
Applied egg-rr97.3%
Taylor expanded in x around 0
Simplified95.8%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified95.1%
Taylor expanded in x around 0
*-commutativeN/A
associate-*r*N/A
cube-multN/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6493.9
Simplified93.9%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
if 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 94.7%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.1
Applied egg-rr95.1%
Taylor expanded in x around 0
Simplified90.7%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified90.1%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.6
Simplified90.6%
Final simplification98.6%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0)))
(t_1 (* (* eps eps) (* eps (* eps eps)))))
(if (<= t_0 -2e-304)
t_1
(if (<= t_0 0.0) (* eps (* (* x (* x x)) (* x 5.0))) t_1))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = eps * ((x * (x * x)) * (x * 5.0));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
t_1 = (eps * eps) * (eps * (eps * eps))
if (t_0 <= (-2d-304)) then
tmp = t_1
else if (t_0 <= 0.0d0) then
tmp = eps * ((x * (x * x)) * (x * 5.0d0))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double eps) {
double t_0 = Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = eps * ((x * (x * x)) * (x * 5.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, eps): t_0 = math.pow((x + eps), 5.0) - math.pow(x, 5.0) t_1 = (eps * eps) * (eps * (eps * eps)) tmp = 0 if t_0 <= -2e-304: tmp = t_1 elif t_0 <= 0.0: tmp = eps * ((x * (x * x)) * (x * 5.0)) else: tmp = t_1 return tmp
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) t_1 = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))) tmp = 0.0 if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(eps * Float64(Float64(x * Float64(x * x)) * Float64(x * 5.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, eps) t_0 = ((x + eps) ^ 5.0) - (x ^ 5.0); t_1 = (eps * eps) * (eps * (eps * eps)); tmp = 0.0; if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = eps * ((x * (x * x)) * (x * 5.0)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], t$95$1, If[LessEqual[t$95$0, 0.0], N[(eps * N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
t_1 := \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\varepsilon \cdot \left(\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot 5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
Taylor expanded in x around 0
Simplified93.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified92.7%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.7
Simplified91.7%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
Final simplification98.4%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0)))
(t_1 (* (* eps eps) (* eps (* eps eps)))))
(if (<= t_0 -2e-304)
t_1
(if (<= t_0 0.0) (* (* x (* x x)) (* eps (* x 5.0))) t_1))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * (eps * (x * 5.0));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
t_1 = (eps * eps) * (eps * (eps * eps))
if (t_0 <= (-2d-304)) then
tmp = t_1
else if (t_0 <= 0.0d0) then
tmp = (x * (x * x)) * (eps * (x * 5.0d0))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double eps) {
double t_0 = Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (x * (x * x)) * (eps * (x * 5.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, eps): t_0 = math.pow((x + eps), 5.0) - math.pow(x, 5.0) t_1 = (eps * eps) * (eps * (eps * eps)) tmp = 0 if t_0 <= -2e-304: tmp = t_1 elif t_0 <= 0.0: tmp = (x * (x * x)) * (eps * (x * 5.0)) else: tmp = t_1 return tmp
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) t_1 = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))) tmp = 0.0 if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(Float64(x * Float64(x * x)) * Float64(eps * Float64(x * 5.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, eps) t_0 = ((x + eps) ^ 5.0) - (x ^ 5.0); t_1 = (eps * eps) * (eps * (eps * eps)); tmp = 0.0; if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = (x * (x * x)) * (eps * (x * 5.0)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], t$95$1, If[LessEqual[t$95$0, 0.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(eps * N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
t_1 := \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(\varepsilon \cdot \left(x \cdot 5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
Taylor expanded in x around 0
Simplified93.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified92.7%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.7
Simplified91.7%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
Taylor expanded in eps around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0)))
(t_1 (* (* eps eps) (* eps (* eps eps)))))
(if (<= t_0 -2e-304)
t_1
(if (<= t_0 0.0) (* (* x x) (* x (* eps (* x 5.0)))) t_1))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (x * x) * (x * (eps * (x * 5.0)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
t_1 = (eps * eps) * (eps * (eps * eps))
if (t_0 <= (-2d-304)) then
tmp = t_1
else if (t_0 <= 0.0d0) then
tmp = (x * x) * (x * (eps * (x * 5.0d0)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double eps) {
double t_0 = Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = (x * x) * (x * (eps * (x * 5.0)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, eps): t_0 = math.pow((x + eps), 5.0) - math.pow(x, 5.0) t_1 = (eps * eps) * (eps * (eps * eps)) tmp = 0 if t_0 <= -2e-304: tmp = t_1 elif t_0 <= 0.0: tmp = (x * x) * (x * (eps * (x * 5.0))) else: tmp = t_1 return tmp
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) t_1 = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))) tmp = 0.0 if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(Float64(x * x) * Float64(x * Float64(eps * Float64(x * 5.0)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, eps) t_0 = ((x + eps) ^ 5.0) - (x ^ 5.0); t_1 = (eps * eps) * (eps * (eps * eps)); tmp = 0.0; if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = (x * x) * (x * (eps * (x * 5.0))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], t$95$1, If[LessEqual[t$95$0, 0.0], N[(N[(x * x), $MachinePrecision] * N[(x * N[(eps * N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
t_1 := \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot \left(\varepsilon \cdot \left(x \cdot 5\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
Taylor expanded in x around 0
Simplified93.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified92.7%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.7
Simplified91.7%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
pow2N/A
pow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
Taylor expanded in eps around 0
unpow2N/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
Final simplification98.4%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (- (pow (+ x eps) 5.0) (pow x 5.0)))
(t_1 (* (* eps eps) (* eps (* eps eps)))))
(if (<= t_0 -2e-304)
t_1
(if (<= t_0 0.0) (* eps (* x (* 5.0 (* x (* x x))))) t_1))))
double code(double x, double eps) {
double t_0 = pow((x + eps), 5.0) - pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = eps * (x * (5.0 * (x * (x * x))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = ((x + eps) ** 5.0d0) - (x ** 5.0d0)
t_1 = (eps * eps) * (eps * (eps * eps))
if (t_0 <= (-2d-304)) then
tmp = t_1
else if (t_0 <= 0.0d0) then
tmp = eps * (x * (5.0d0 * (x * (x * x))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double eps) {
double t_0 = Math.pow((x + eps), 5.0) - Math.pow(x, 5.0);
double t_1 = (eps * eps) * (eps * (eps * eps));
double tmp;
if (t_0 <= -2e-304) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = eps * (x * (5.0 * (x * (x * x))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, eps): t_0 = math.pow((x + eps), 5.0) - math.pow(x, 5.0) t_1 = (eps * eps) * (eps * (eps * eps)) tmp = 0 if t_0 <= -2e-304: tmp = t_1 elif t_0 <= 0.0: tmp = eps * (x * (5.0 * (x * (x * x)))) else: tmp = t_1 return tmp
function code(x, eps) t_0 = Float64((Float64(x + eps) ^ 5.0) - (x ^ 5.0)) t_1 = Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))) tmp = 0.0 if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(eps * Float64(x * Float64(5.0 * Float64(x * Float64(x * x))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, eps) t_0 = ((x + eps) ^ 5.0) - (x ^ 5.0); t_1 = (eps * eps) * (eps * (eps * eps)); tmp = 0.0; if (t_0 <= -2e-304) tmp = t_1; elseif (t_0 <= 0.0) tmp = eps * (x * (5.0 * (x * (x * x)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[(x + eps), $MachinePrecision], 5.0], $MachinePrecision] - N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-304], t$95$1, If[LessEqual[t$95$0, 0.0], N[(eps * N[(x * N[(5.0 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(x + \varepsilon\right)}^{5} - {x}^{5}\\
t_1 := \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-304}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\varepsilon \cdot \left(x \cdot \left(5 \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < -1.99999999999999994e-304 or 0.0 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) Initial program 96.0%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.2
Applied egg-rr96.2%
Taylor expanded in x around 0
Simplified93.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified92.7%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.7
Simplified91.7%
if -1.99999999999999994e-304 < (-.f64 (pow.f64 (+.f64 x eps) #s(literal 5 binary64)) (pow.f64 x #s(literal 5 binary64))) < 0.0Initial program 87.7%
Taylor expanded in x around -inf
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
+-commutativeN/A
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Taylor expanded in x around 0
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9
Simplified99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f6499.9
Applied egg-rr99.9%
Taylor expanded in x around inf
*-commutativeN/A
metadata-evalN/A
pow-plusN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9
Simplified99.9%
Final simplification98.4%
(FPCore (x eps) :precision binary64 (* (* eps eps) (* eps (* eps eps))))
double code(double x, double eps) {
return (eps * eps) * (eps * (eps * eps));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (eps * eps) * (eps * (eps * eps))
end function
public static double code(double x, double eps) {
return (eps * eps) * (eps * (eps * eps));
}
def code(x, eps): return (eps * eps) * (eps * (eps * eps))
function code(x, eps) return Float64(Float64(eps * eps) * Float64(eps * Float64(eps * eps))) end
function tmp = code(x, eps) tmp = (eps * eps) * (eps * (eps * eps)); end
code[x_, eps_] := N[(N[(eps * eps), $MachinePrecision] * N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right)
\end{array}
Initial program 89.2%
--lowering--.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6486.8
Applied egg-rr86.8%
Taylor expanded in x around 0
Simplified88.8%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-outN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified88.6%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.4
Simplified88.4%
herbie shell --seed 2024199
(FPCore (x eps)
:name "ENA, Section 1.4, Exercise 4b, n=5"
:precision binary64
:pre (and (and (<= -1000000000.0 x) (<= x 1000000000.0)) (and (<= -1.0 eps) (<= eps 1.0)))
(- (pow (+ x eps) 5.0) (pow x 5.0)))