
(FPCore (x) :precision binary64 (/ (exp x) (- (exp x) 1.0)))
double code(double x) {
return exp(x) / (exp(x) - 1.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp(x) / (exp(x) - 1.0d0)
end function
public static double code(double x) {
return Math.exp(x) / (Math.exp(x) - 1.0);
}
def code(x): return math.exp(x) / (math.exp(x) - 1.0)
function code(x) return Float64(exp(x) / Float64(exp(x) - 1.0)) end
function tmp = code(x) tmp = exp(x) / (exp(x) - 1.0); end
code[x_] := N[(N[Exp[x], $MachinePrecision] / N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x}}{e^{x} - 1}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ (exp x) (- (exp x) 1.0)))
double code(double x) {
return exp(x) / (exp(x) - 1.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp(x) / (exp(x) - 1.0d0)
end function
public static double code(double x) {
return Math.exp(x) / (Math.exp(x) - 1.0);
}
def code(x): return math.exp(x) / (math.exp(x) - 1.0)
function code(x) return Float64(exp(x) / Float64(exp(x) - 1.0)) end
function tmp = code(x) tmp = exp(x) / (exp(x) - 1.0); end
code[x_] := N[(N[Exp[x], $MachinePrecision] / N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x}}{e^{x} - 1}
\end{array}
(FPCore (x) :precision binary64 (/ (exp x) (expm1 x)))
double code(double x) {
return exp(x) / expm1(x);
}
public static double code(double x) {
return Math.exp(x) / Math.expm1(x);
}
def code(x): return math.exp(x) / math.expm1(x)
function code(x) return Float64(exp(x) / expm1(x)) end
code[x_] := N[(N[Exp[x], $MachinePrecision] / N[(Exp[x] - 1), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x}}{\mathsf{expm1}\left(x\right)}
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
(FPCore (x) :precision binary64 (if (<= (exp x) 0.95) (/ 1.0 (+ 1.0 (/ -1.0 (exp x)))) (+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333)))))
double code(double x) {
double tmp;
if (exp(x) <= 0.95) {
tmp = 1.0 / (1.0 + (-1.0 / exp(x)));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (exp(x) <= 0.95d0) then
tmp = 1.0d0 / (1.0d0 + ((-1.0d0) / exp(x)))
else
tmp = (1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (Math.exp(x) <= 0.95) {
tmp = 1.0 / (1.0 + (-1.0 / Math.exp(x)));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
def code(x): tmp = 0 if math.exp(x) <= 0.95: tmp = 1.0 / (1.0 + (-1.0 / math.exp(x))) else: tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)) return tmp
function code(x) tmp = 0.0 if (exp(x) <= 0.95) tmp = Float64(1.0 / Float64(1.0 + Float64(-1.0 / exp(x)))); else tmp = Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (exp(x) <= 0.95) tmp = 1.0 / (1.0 + (-1.0 / exp(x))); else tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[N[Exp[x], $MachinePrecision], 0.95], N[(1.0 / N[(1.0 + N[(-1.0 / N[Exp[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{x} \leq 0.95:\\
\;\;\;\;\frac{1}{1 + \frac{-1}{e^{x}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\\
\end{array}
\end{array}
if (exp.f64 x) < 0.94999999999999996Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
clear-numN/A
/-lowering-/.f64N/A
div-subN/A
*-inversesN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6499.9%
Applied egg-rr99.9%
if 0.94999999999999996 < (exp.f64 x) Initial program 7.5%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.1%
Simplified99.1%
Final simplification99.4%
(FPCore (x)
:precision binary64
(if (<= x -3.8)
(/ (exp x) x)
(+
(+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333)))
(* x (* x (* x -0.001388888888888889))))))
double code(double x) {
double tmp;
if (x <= -3.8) {
tmp = exp(x) / x;
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-3.8d0)) then
tmp = exp(x) / x
else
tmp = ((1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))) + (x * (x * (x * (-0.001388888888888889d0))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -3.8) {
tmp = Math.exp(x) / x;
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
def code(x): tmp = 0 if x <= -3.8: tmp = math.exp(x) / x else: tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))) return tmp
function code(x) tmp = 0.0 if (x <= -3.8) tmp = Float64(exp(x) / x); else tmp = Float64(Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))) + Float64(x * Float64(x * Float64(x * -0.001388888888888889)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -3.8) tmp = exp(x) / x; else tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -3.8], N[(N[Exp[x], $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * N[(x * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.8:\\
\;\;\;\;\frac{e^{x}}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\right) + x \cdot \left(x \cdot \left(x \cdot -0.001388888888888889\right)\right)\\
\end{array}
\end{array}
if x < -3.7999999999999998Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified100.0%
if -3.7999999999999998 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified99.0%
distribute-rgt-inN/A
*-commutativeN/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.0%
Applied egg-rr99.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x x))) (t_1 (* (* x x) (* x t_0))))
(if (<= x -5.6e+102)
(/ 1.0 (* x (+ 1.0 (* x (+ -0.5 (* x 0.16666666666666666))))))
(if (<= x -2.35e+51)
(* t_0 (/ (* x x) t_1))
(if (<= x -3.9)
(* t_1 (/ (* x x) (* t_0 t_1)))
(+
(+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333)))
(* x (* x (* x -0.001388888888888889)))))))))
double code(double x) {
double t_0 = x * (x * x);
double t_1 = (x * x) * (x * t_0);
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -2.35e+51) {
tmp = t_0 * ((x * x) / t_1);
} else if (x <= -3.9) {
tmp = t_1 * ((x * x) / (t_0 * t_1));
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = x * (x * x)
t_1 = (x * x) * (x * t_0)
if (x <= (-5.6d+102)) then
tmp = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * 0.16666666666666666d0)))))
else if (x <= (-2.35d+51)) then
tmp = t_0 * ((x * x) / t_1)
else if (x <= (-3.9d0)) then
tmp = t_1 * ((x * x) / (t_0 * t_1))
else
tmp = ((1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))) + (x * (x * (x * (-0.001388888888888889d0))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * x);
double t_1 = (x * x) * (x * t_0);
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -2.35e+51) {
tmp = t_0 * ((x * x) / t_1);
} else if (x <= -3.9) {
tmp = t_1 * ((x * x) / (t_0 * t_1));
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
def code(x): t_0 = x * (x * x) t_1 = (x * x) * (x * t_0) tmp = 0 if x <= -5.6e+102: tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))) elif x <= -2.35e+51: tmp = t_0 * ((x * x) / t_1) elif x <= -3.9: tmp = t_1 * ((x * x) / (t_0 * t_1)) else: tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))) return tmp
function code(x) t_0 = Float64(x * Float64(x * x)) t_1 = Float64(Float64(x * x) * Float64(x * t_0)) tmp = 0.0 if (x <= -5.6e+102) tmp = Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * 0.16666666666666666)))))); elseif (x <= -2.35e+51) tmp = Float64(t_0 * Float64(Float64(x * x) / t_1)); elseif (x <= -3.9) tmp = Float64(t_1 * Float64(Float64(x * x) / Float64(t_0 * t_1))); else tmp = Float64(Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))) + Float64(x * Float64(x * Float64(x * -0.001388888888888889)))); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * x); t_1 = (x * x) * (x * t_0); tmp = 0.0; if (x <= -5.6e+102) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))); elseif (x <= -2.35e+51) tmp = t_0 * ((x * x) / t_1); elseif (x <= -3.9) tmp = t_1 * ((x * x) / (t_0 * t_1)); else tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x * x), $MachinePrecision] * N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -5.6e+102], N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -2.35e+51], N[(t$95$0 * N[(N[(x * x), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.9], N[(t$95$1 * N[(N[(x * x), $MachinePrecision] / N[(t$95$0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * N[(x * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
t_1 := \left(x \cdot x\right) \cdot \left(x \cdot t\_0\right)\\
\mathbf{if}\;x \leq -5.6 \cdot 10^{+102}:\\
\;\;\;\;\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot 0.16666666666666666\right)\right)}\\
\mathbf{elif}\;x \leq -2.35 \cdot 10^{+51}:\\
\;\;\;\;t\_0 \cdot \frac{x \cdot x}{t\_1}\\
\mathbf{elif}\;x \leq -3.9:\\
\;\;\;\;t\_1 \cdot \frac{x \cdot x}{t\_0 \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\right) + x \cdot \left(x \cdot \left(x \cdot -0.001388888888888889\right)\right)\\
\end{array}
\end{array}
if x < -5.60000000000000037e102Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.0%
Simplified2.0%
Taylor expanded in x around 0
Simplified1.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f641.8%
Applied egg-rr1.8%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
if -5.60000000000000037e102 < x < -2.3500000000000001e51Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified1.9%
Taylor expanded in x around -inf
Simplified1.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.8%
Simplified3.8%
Applied egg-rr100.0%
if -2.3500000000000001e51 < x < -3.89999999999999991Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified2.5%
Taylor expanded in x around -inf
Simplified2.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.4%
Simplified3.4%
Applied egg-rr48.8%
if -3.89999999999999991 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified99.0%
distribute-rgt-inN/A
*-commutativeN/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.0%
Applied egg-rr99.0%
Final simplification95.9%
(FPCore (x)
:precision binary64
(if (<= x -5.6e+102)
(/ 1.0 (* x (+ 1.0 (* x (+ -0.5 (* x 0.16666666666666666))))))
(if (<= x -3.9)
(/ (* x (* x x)) (* (* x x) (* x x)))
(+
(+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333)))
(* x (* x (* x -0.001388888888888889)))))))
double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -3.9) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-5.6d+102)) then
tmp = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * 0.16666666666666666d0)))))
else if (x <= (-3.9d0)) then
tmp = (x * (x * x)) / ((x * x) * (x * x))
else
tmp = ((1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))) + (x * (x * (x * (-0.001388888888888889d0))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -3.9) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889)));
}
return tmp;
}
def code(x): tmp = 0 if x <= -5.6e+102: tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))) elif x <= -3.9: tmp = (x * (x * x)) / ((x * x) * (x * x)) else: tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))) return tmp
function code(x) tmp = 0.0 if (x <= -5.6e+102) tmp = Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * 0.16666666666666666)))))); elseif (x <= -3.9) tmp = Float64(Float64(x * Float64(x * x)) / Float64(Float64(x * x) * Float64(x * x))); else tmp = Float64(Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))) + Float64(x * Float64(x * Float64(x * -0.001388888888888889)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -5.6e+102) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))); elseif (x <= -3.9) tmp = (x * (x * x)) / ((x * x) * (x * x)); else tmp = ((1.0 / x) + (0.5 + (x * 0.08333333333333333))) + (x * (x * (x * -0.001388888888888889))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -5.6e+102], N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.9], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * N[(x * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.6 \cdot 10^{+102}:\\
\;\;\;\;\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot 0.16666666666666666\right)\right)}\\
\mathbf{elif}\;x \leq -3.9:\\
\;\;\;\;\frac{x \cdot \left(x \cdot x\right)}{\left(x \cdot x\right) \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\right) + x \cdot \left(x \cdot \left(x \cdot -0.001388888888888889\right)\right)\\
\end{array}
\end{array}
if x < -5.60000000000000037e102Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.0%
Simplified2.0%
Taylor expanded in x around 0
Simplified1.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f641.8%
Applied egg-rr1.8%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
if -5.60000000000000037e102 < x < -3.89999999999999991Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified2.2%
Taylor expanded in x around -inf
Simplified2.2%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.6%
Simplified3.6%
*-commutativeN/A
flip3--N/A
metadata-evalN/A
cube-unmultN/A
frac-timesN/A
neg-mul-1N/A
sub0-negN/A
remove-double-negN/A
metadata-evalN/A
+-lft-identityN/A
distribute-rgt-outN/A
+-commutativeN/A
+-lft-identityN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.9%
Applied egg-rr31.9%
if -3.89999999999999991 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified99.0%
distribute-rgt-inN/A
*-commutativeN/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.0%
Applied egg-rr99.0%
(FPCore (x)
:precision binary64
(if (<= x -5.6e+102)
(/ 1.0 (* x (+ 1.0 (* x (+ -0.5 (* x 0.16666666666666666))))))
(if (<= x -3.9)
(/ (* x (* x x)) (* (* x x) (* x x)))
(+
(+ (/ 1.0 x) 0.5)
(* x (+ 0.08333333333333333 (* x (* x -0.001388888888888889))))))))
double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -3.9) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = ((1.0 / x) + 0.5) + (x * (0.08333333333333333 + (x * (x * -0.001388888888888889))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-5.6d+102)) then
tmp = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * 0.16666666666666666d0)))))
else if (x <= (-3.9d0)) then
tmp = (x * (x * x)) / ((x * x) * (x * x))
else
tmp = ((1.0d0 / x) + 0.5d0) + (x * (0.08333333333333333d0 + (x * (x * (-0.001388888888888889d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -3.9) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = ((1.0 / x) + 0.5) + (x * (0.08333333333333333 + (x * (x * -0.001388888888888889))));
}
return tmp;
}
def code(x): tmp = 0 if x <= -5.6e+102: tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))) elif x <= -3.9: tmp = (x * (x * x)) / ((x * x) * (x * x)) else: tmp = ((1.0 / x) + 0.5) + (x * (0.08333333333333333 + (x * (x * -0.001388888888888889)))) return tmp
function code(x) tmp = 0.0 if (x <= -5.6e+102) tmp = Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * 0.16666666666666666)))))); elseif (x <= -3.9) tmp = Float64(Float64(x * Float64(x * x)) / Float64(Float64(x * x) * Float64(x * x))); else tmp = Float64(Float64(Float64(1.0 / x) + 0.5) + Float64(x * Float64(0.08333333333333333 + Float64(x * Float64(x * -0.001388888888888889))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -5.6e+102) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))); elseif (x <= -3.9) tmp = (x * (x * x)) / ((x * x) * (x * x)); else tmp = ((1.0 / x) + 0.5) + (x * (0.08333333333333333 + (x * (x * -0.001388888888888889)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -5.6e+102], N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.9], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 / x), $MachinePrecision] + 0.5), $MachinePrecision] + N[(x * N[(0.08333333333333333 + N[(x * N[(x * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.6 \cdot 10^{+102}:\\
\;\;\;\;\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot 0.16666666666666666\right)\right)}\\
\mathbf{elif}\;x \leq -3.9:\\
\;\;\;\;\frac{x \cdot \left(x \cdot x\right)}{\left(x \cdot x\right) \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x} + 0.5\right) + x \cdot \left(0.08333333333333333 + x \cdot \left(x \cdot -0.001388888888888889\right)\right)\\
\end{array}
\end{array}
if x < -5.60000000000000037e102Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.0%
Simplified2.0%
Taylor expanded in x around 0
Simplified1.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f641.8%
Applied egg-rr1.8%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
if -5.60000000000000037e102 < x < -3.89999999999999991Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified2.2%
Taylor expanded in x around -inf
Simplified2.2%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.6%
Simplified3.6%
*-commutativeN/A
flip3--N/A
metadata-evalN/A
cube-unmultN/A
frac-timesN/A
neg-mul-1N/A
sub0-negN/A
remove-double-negN/A
metadata-evalN/A
+-lft-identityN/A
distribute-rgt-outN/A
+-commutativeN/A
+-lft-identityN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.9%
Applied egg-rr31.9%
if -3.89999999999999991 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified99.0%
(FPCore (x)
:precision binary64
(if (<= x -5.6e+102)
(/ 1.0 (* x (+ 1.0 (* x (+ -0.5 (* x 0.16666666666666666))))))
(if (<= x -6.0)
(/ (* x (* x x)) (* (* x x) (* x x)))
(+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333))))))
double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -6.0) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-5.6d+102)) then
tmp = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * 0.16666666666666666d0)))))
else if (x <= (-6.0d0)) then
tmp = (x * (x * x)) / ((x * x) * (x * x))
else
tmp = (1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -5.6e+102) {
tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
} else if (x <= -6.0) {
tmp = (x * (x * x)) / ((x * x) * (x * x));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
def code(x): tmp = 0 if x <= -5.6e+102: tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))) elif x <= -6.0: tmp = (x * (x * x)) / ((x * x) * (x * x)) else: tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)) return tmp
function code(x) tmp = 0.0 if (x <= -5.6e+102) tmp = Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * 0.16666666666666666)))))); elseif (x <= -6.0) tmp = Float64(Float64(x * Float64(x * x)) / Float64(Float64(x * x) * Float64(x * x))); else tmp = Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -5.6e+102) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))); elseif (x <= -6.0) tmp = (x * (x * x)) / ((x * x) * (x * x)); else tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -5.6e+102], N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -6.0], N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.6 \cdot 10^{+102}:\\
\;\;\;\;\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot 0.16666666666666666\right)\right)}\\
\mathbf{elif}\;x \leq -6:\\
\;\;\;\;\frac{x \cdot \left(x \cdot x\right)}{\left(x \cdot x\right) \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\\
\end{array}
\end{array}
if x < -5.60000000000000037e102Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.0%
Simplified2.0%
Taylor expanded in x around 0
Simplified1.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f641.8%
Applied egg-rr1.8%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
if -5.60000000000000037e102 < x < -6Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified2.2%
Taylor expanded in x around -inf
Simplified2.2%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.6%
Simplified3.6%
*-commutativeN/A
flip3--N/A
metadata-evalN/A
cube-unmultN/A
frac-timesN/A
neg-mul-1N/A
sub0-negN/A
remove-double-negN/A
metadata-evalN/A
+-lft-identityN/A
distribute-rgt-outN/A
+-commutativeN/A
+-lft-identityN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.9%
Applied egg-rr31.9%
if -6 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.9%
Simplified98.9%
(FPCore (x)
:precision binary64
(if (<= x -1.35e+154)
(/ x (* x x))
(if (<= x -6.0)
(/ (* x x) (* x (* x x)))
(+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333))))))
double code(double x) {
double tmp;
if (x <= -1.35e+154) {
tmp = x / (x * x);
} else if (x <= -6.0) {
tmp = (x * x) / (x * (x * x));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.35d+154)) then
tmp = x / (x * x)
else if (x <= (-6.0d0)) then
tmp = (x * x) / (x * (x * x))
else
tmp = (1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.35e+154) {
tmp = x / (x * x);
} else if (x <= -6.0) {
tmp = (x * x) / (x * (x * x));
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.35e+154: tmp = x / (x * x) elif x <= -6.0: tmp = (x * x) / (x * (x * x)) else: tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)) return tmp
function code(x) tmp = 0.0 if (x <= -1.35e+154) tmp = Float64(x / Float64(x * x)); elseif (x <= -6.0) tmp = Float64(Float64(x * x) / Float64(x * Float64(x * x))); else tmp = Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.35e+154) tmp = x / (x * x); elseif (x <= -6.0) tmp = (x * x) / (x * (x * x)); else tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.35e+154], N[(x / N[(x * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -6.0], N[(N[(x * x), $MachinePrecision] / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.35 \cdot 10^{+154}:\\
\;\;\;\;\frac{x}{x \cdot x}\\
\mathbf{elif}\;x \leq -6:\\
\;\;\;\;\frac{x \cdot x}{x \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\\
\end{array}
\end{array}
if x < -1.35000000000000003e154Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified1.6%
Taylor expanded in x around -inf
Simplified1.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
frac-2negN/A
metadata-evalN/A
un-div-invN/A
neg-sub0N/A
frac-2negN/A
/-lowering-/.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
if -1.35000000000000003e154 < x < -6Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified2.0%
Taylor expanded in x around -inf
Simplified2.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f643.8%
Simplified3.8%
*-commutativeN/A
flip--N/A
+-lft-identityN/A
times-fracN/A
associate-*r*N/A
metadata-evalN/A
neg-sub0N/A
distribute-rgt-neg-inN/A
neg-sub0N/A
associate-*l*N/A
neg-mul-1N/A
neg-sub0N/A
sqr-negN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.7%
Applied egg-rr31.7%
if -6 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.9%
Simplified98.9%
(FPCore (x)
:precision binary64
(/
1.0
(*
x
(+
1.0
(*
x
(+ -0.5 (* x (+ 0.16666666666666666 (* x -0.041666666666666664)))))))))
double code(double x) {
return 1.0 / (x * (1.0 + (x * (-0.5 + (x * (0.16666666666666666 + (x * -0.041666666666666664)))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * (0.16666666666666666d0 + (x * (-0.041666666666666664d0))))))))
end function
public static double code(double x) {
return 1.0 / (x * (1.0 + (x * (-0.5 + (x * (0.16666666666666666 + (x * -0.041666666666666664)))))));
}
def code(x): return 1.0 / (x * (1.0 + (x * (-0.5 + (x * (0.16666666666666666 + (x * -0.041666666666666664)))))))
function code(x) return Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * Float64(0.16666666666666666 + Float64(x * -0.041666666666666664)))))))) end
function tmp = code(x) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * (0.16666666666666666 + (x * -0.041666666666666664))))))); end
code[x_] := N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * N[(0.16666666666666666 + N[(x * -0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot \left(0.16666666666666666 + x \cdot -0.041666666666666664\right)\right)\right)}
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in x around 0
Simplified66.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6490.3%
Simplified90.3%
(FPCore (x) :precision binary64 (if (<= x -6.0) (/ x (* x x)) (+ (/ 1.0 x) (+ 0.5 (* x 0.08333333333333333)))))
double code(double x) {
double tmp;
if (x <= -6.0) {
tmp = x / (x * x);
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-6.0d0)) then
tmp = x / (x * x)
else
tmp = (1.0d0 / x) + (0.5d0 + (x * 0.08333333333333333d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -6.0) {
tmp = x / (x * x);
} else {
tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333));
}
return tmp;
}
def code(x): tmp = 0 if x <= -6.0: tmp = x / (x * x) else: tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)) return tmp
function code(x) tmp = 0.0 if (x <= -6.0) tmp = Float64(x / Float64(x * x)); else tmp = Float64(Float64(1.0 / x) + Float64(0.5 + Float64(x * 0.08333333333333333))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -6.0) tmp = x / (x * x); else tmp = (1.0 / x) + (0.5 + (x * 0.08333333333333333)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -6.0], N[(x / N[(x * x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / x), $MachinePrecision] + N[(0.5 + N[(x * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6:\\
\;\;\;\;\frac{x}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} + \left(0.5 + x \cdot 0.08333333333333333\right)\\
\end{array}
\end{array}
if x < -6Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified1.8%
Taylor expanded in x around -inf
Simplified1.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6446.9%
Simplified46.9%
frac-2negN/A
metadata-evalN/A
un-div-invN/A
neg-sub0N/A
frac-2negN/A
/-lowering-/.f64N/A
*-lowering-*.f6446.9%
Applied egg-rr46.9%
if -6 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.9%
Simplified98.9%
(FPCore (x) :precision binary64 (/ 1.0 (* x (+ 1.0 (* x (+ -0.5 (* x 0.16666666666666666)))))))
double code(double x) {
return 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (x * (1.0d0 + (x * ((-0.5d0) + (x * 0.16666666666666666d0)))))
end function
public static double code(double x) {
return 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))));
}
def code(x): return 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666)))))
function code(x) return Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * Float64(-0.5 + Float64(x * 0.16666666666666666)))))) end
function tmp = code(x) tmp = 1.0 / (x * (1.0 + (x * (-0.5 + (x * 0.16666666666666666))))); end
code[x_] := N[(1.0 / N[(x * N[(1.0 + N[(x * N[(-0.5 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x \cdot \left(1 + x \cdot \left(-0.5 + x \cdot 0.16666666666666666\right)\right)}
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in x around 0
Simplified66.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6486.5%
Simplified86.5%
(FPCore (x) :precision binary64 (if (<= x -1.7) (/ x (* x x)) (+ (/ 1.0 x) 0.5)))
double code(double x) {
double tmp;
if (x <= -1.7) {
tmp = x / (x * x);
} else {
tmp = (1.0 / x) + 0.5;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.7d0)) then
tmp = x / (x * x)
else
tmp = (1.0d0 / x) + 0.5d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.7) {
tmp = x / (x * x);
} else {
tmp = (1.0 / x) + 0.5;
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.7: tmp = x / (x * x) else: tmp = (1.0 / x) + 0.5 return tmp
function code(x) tmp = 0.0 if (x <= -1.7) tmp = Float64(x / Float64(x * x)); else tmp = Float64(Float64(1.0 / x) + 0.5); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.7) tmp = x / (x * x); else tmp = (1.0 / x) + 0.5; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.7], N[(x / N[(x * x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / x), $MachinePrecision] + 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.7:\\
\;\;\;\;\frac{x}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} + 0.5\\
\end{array}
\end{array}
if x < -1.69999999999999996Initial program 100.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
associate-*l/N/A
*-lft-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified1.8%
Taylor expanded in x around -inf
Simplified1.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6446.9%
Simplified46.9%
frac-2negN/A
metadata-evalN/A
un-div-invN/A
neg-sub0N/A
frac-2negN/A
/-lowering-/.f64N/A
*-lowering-*.f6446.9%
Applied egg-rr46.9%
if -1.69999999999999996 < x Initial program 8.0%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-*l/N/A
distribute-rgt-inN/A
associate-*l*N/A
rgt-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f6498.3%
Simplified98.3%
(FPCore (x) :precision binary64 (/ 1.0 (* x (+ 1.0 (* x -0.5)))))
double code(double x) {
return 1.0 / (x * (1.0 + (x * -0.5)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (x * (1.0d0 + (x * (-0.5d0))))
end function
public static double code(double x) {
return 1.0 / (x * (1.0 + (x * -0.5)));
}
def code(x): return 1.0 / (x * (1.0 + (x * -0.5)))
function code(x) return Float64(1.0 / Float64(x * Float64(1.0 + Float64(x * -0.5)))) end
function tmp = code(x) tmp = 1.0 / (x * (1.0 + (x * -0.5))); end
code[x_] := N[(1.0 / N[(x * N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x \cdot \left(1 + x \cdot -0.5\right)}
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-/l*N/A
associate-*l/N/A
distribute-lft-inN/A
*-rgt-identityN/A
associate-*r*N/A
lft-mult-inverseN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.1%
Simplified66.1%
Taylor expanded in x around 0
Simplified66.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6466.0%
Applied egg-rr66.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6481.0%
Simplified81.0%
(FPCore (x) :precision binary64 (+ (/ 1.0 x) 0.5))
double code(double x) {
return (1.0 / x) + 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (1.0d0 / x) + 0.5d0
end function
public static double code(double x) {
return (1.0 / x) + 0.5;
}
def code(x): return (1.0 / x) + 0.5
function code(x) return Float64(Float64(1.0 / x) + 0.5) end
function tmp = code(x) tmp = (1.0 / x) + 0.5; end
code[x_] := N[(N[(1.0 / x), $MachinePrecision] + 0.5), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x} + 0.5
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-*l/N/A
distribute-rgt-inN/A
associate-*l*N/A
rgt-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f6466.0%
Simplified66.0%
(FPCore (x) :precision binary64 (/ 1.0 x))
double code(double x) {
return 1.0 / x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / x
end function
public static double code(double x) {
return 1.0 / x;
}
def code(x): return 1.0 / x
function code(x) return Float64(1.0 / x) end
function tmp = code(x) tmp = 1.0 / x; end
code[x_] := N[(1.0 / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x}
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
/-lowering-/.f6465.4%
Simplified65.4%
(FPCore (x) :precision binary64 1.0)
double code(double x) {
return 1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double x) {
return 1.0;
}
def code(x): return 1.0
function code(x) return 1.0 end
function tmp = code(x) tmp = 1.0; end
code[x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified97.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f6464.8%
Simplified64.8%
Taylor expanded in x around inf
Simplified3.8%
(FPCore (x) :precision binary64 0.5)
double code(double x) {
return 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.5d0
end function
public static double code(double x) {
return 0.5;
}
def code(x): return 0.5
function code(x) return 0.5 end
function tmp = code(x) tmp = 0.5; end
code[x_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 39.3%
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lft-identityN/A
associate-*l/N/A
distribute-rgt-inN/A
associate-*l*N/A
rgt-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f64N/A
/-lowering-/.f6466.0%
Simplified66.0%
Taylor expanded in x around inf
Simplified3.4%
(FPCore (x) :precision binary64 (/ (- 1.0) (expm1 (- x))))
double code(double x) {
return -1.0 / expm1(-x);
}
public static double code(double x) {
return -1.0 / Math.expm1(-x);
}
def code(x): return -1.0 / math.expm1(-x)
function code(x) return Float64(Float64(-1.0) / expm1(Float64(-x))) end
code[x_] := N[((-1.0) / N[(Exp[(-x)] - 1), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-1}{\mathsf{expm1}\left(-x\right)}
\end{array}
herbie shell --seed 2024288
(FPCore (x)
:name "expq2 (section 3.11)"
:precision binary64
:pre (> 710.0 x)
:alt
(! :herbie-platform default (/ (- 1) (expm1 (- x))))
(/ (exp x) (- (exp x) 1.0)))