
(FPCore (x) :precision binary64 (/ (- (exp x) 1.0) x))
double code(double x) {
return (exp(x) - 1.0) / x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (exp(x) - 1.0d0) / x
end function
public static double code(double x) {
return (Math.exp(x) - 1.0) / x;
}
def code(x): return (math.exp(x) - 1.0) / x
function code(x) return Float64(Float64(exp(x) - 1.0) / x) end
function tmp = code(x) tmp = (exp(x) - 1.0) / x; end
code[x_] := N[(N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x} - 1}{x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ (- (exp x) 1.0) x))
double code(double x) {
return (exp(x) - 1.0) / x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (exp(x) - 1.0d0) / x
end function
public static double code(double x) {
return (Math.exp(x) - 1.0) / x;
}
def code(x): return (math.exp(x) - 1.0) / x
function code(x) return Float64(Float64(exp(x) - 1.0) / x) end
function tmp = code(x) tmp = (exp(x) - 1.0) / x; end
code[x_] := N[(N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{x} - 1}{x}
\end{array}
(FPCore (x) :precision binary64 (/ (expm1 x) x))
double code(double x) {
return expm1(x) / x;
}
public static double code(double x) {
return Math.expm1(x) / x;
}
def code(x): return math.expm1(x) / x
function code(x) return Float64(expm1(x) / x) end
code[x_] := N[(N[(Exp[x] - 1), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{\mathsf{expm1}\left(x\right)}{x}
\end{array}
Initial program 56.6%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ 0.16666666666666666 (* x 0.041666666666666664)))
(t_1 (* x t_0))
(t_2 (* t_0 t_0))
(t_3 (* (* (* x x) t_0) (* x t_2))))
(if (<= x -1.55)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 5e+51)
(+
1.0
(/
(/
(* x (- 0.015625 (* (* x (* x x)) (* t_3 (* t_0 t_2)))))
(- 0.125 t_3))
(+ 0.25 (* t_1 (+ -0.5 t_1)))))
(*
(+ 0.004629629629629629 (* x (* (* x x) 7.233796296296296e-5)))
(* x (* x (+ 36.0 (* x 9.0)))))))))
double code(double x) {
double t_0 = 0.16666666666666666 + (x * 0.041666666666666664);
double t_1 = x * t_0;
double t_2 = t_0 * t_0;
double t_3 = ((x * x) * t_0) * (x * t_2);
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5e+51) {
tmp = 1.0 + (((x * (0.015625 - ((x * (x * x)) * (t_3 * (t_0 * t_2))))) / (0.125 - t_3)) / (0.25 + (t_1 * (-0.5 + t_1))));
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_0 = 0.16666666666666666d0 + (x * 0.041666666666666664d0)
t_1 = x * t_0
t_2 = t_0 * t_0
t_3 = ((x * x) * t_0) * (x * t_2)
if (x <= (-1.55d0)) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 5d+51) then
tmp = 1.0d0 + (((x * (0.015625d0 - ((x * (x * x)) * (t_3 * (t_0 * t_2))))) / (0.125d0 - t_3)) / (0.25d0 + (t_1 * ((-0.5d0) + t_1))))
else
tmp = (0.004629629629629629d0 + (x * ((x * x) * 7.233796296296296d-5))) * (x * (x * (36.0d0 + (x * 9.0d0))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 0.16666666666666666 + (x * 0.041666666666666664);
double t_1 = x * t_0;
double t_2 = t_0 * t_0;
double t_3 = ((x * x) * t_0) * (x * t_2);
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5e+51) {
tmp = 1.0 + (((x * (0.015625 - ((x * (x * x)) * (t_3 * (t_0 * t_2))))) / (0.125 - t_3)) / (0.25 + (t_1 * (-0.5 + t_1))));
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
def code(x): t_0 = 0.16666666666666666 + (x * 0.041666666666666664) t_1 = x * t_0 t_2 = t_0 * t_0 t_3 = ((x * x) * t_0) * (x * t_2) tmp = 0 if x <= -1.55: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 5e+51: tmp = 1.0 + (((x * (0.015625 - ((x * (x * x)) * (t_3 * (t_0 * t_2))))) / (0.125 - t_3)) / (0.25 + (t_1 * (-0.5 + t_1)))) else: tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))) return tmp
function code(x) t_0 = Float64(0.16666666666666666 + Float64(x * 0.041666666666666664)) t_1 = Float64(x * t_0) t_2 = Float64(t_0 * t_0) t_3 = Float64(Float64(Float64(x * x) * t_0) * Float64(x * t_2)) tmp = 0.0 if (x <= -1.55) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 5e+51) tmp = Float64(1.0 + Float64(Float64(Float64(x * Float64(0.015625 - Float64(Float64(x * Float64(x * x)) * Float64(t_3 * Float64(t_0 * t_2))))) / Float64(0.125 - t_3)) / Float64(0.25 + Float64(t_1 * Float64(-0.5 + t_1))))); else tmp = Float64(Float64(0.004629629629629629 + Float64(x * Float64(Float64(x * x) * 7.233796296296296e-5))) * Float64(x * Float64(x * Float64(36.0 + Float64(x * 9.0))))); end return tmp end
function tmp_2 = code(x) t_0 = 0.16666666666666666 + (x * 0.041666666666666664); t_1 = x * t_0; t_2 = t_0 * t_0; t_3 = ((x * x) * t_0) * (x * t_2); tmp = 0.0; if (x <= -1.55) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 5e+51) tmp = 1.0 + (((x * (0.015625 - ((x * (x * x)) * (t_3 * (t_0 * t_2))))) / (0.125 - t_3)) / (0.25 + (t_1 * (-0.5 + t_1)))); else tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 * t$95$0), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(x * x), $MachinePrecision] * t$95$0), $MachinePrecision] * N[(x * t$95$2), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.55], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 5e+51], N[(1.0 + N[(N[(N[(x * N[(0.015625 - N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(t$95$3 * N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.125 - t$95$3), $MachinePrecision]), $MachinePrecision] / N[(0.25 + N[(t$95$1 * N[(-0.5 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.004629629629629629 + N[(x * N[(N[(x * x), $MachinePrecision] * 7.233796296296296e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(36.0 + N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.16666666666666666 + x \cdot 0.041666666666666664\\
t_1 := x \cdot t\_0\\
t_2 := t\_0 \cdot t\_0\\
t_3 := \left(\left(x \cdot x\right) \cdot t\_0\right) \cdot \left(x \cdot t\_2\right)\\
\mathbf{if}\;x \leq -1.55:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 5 \cdot 10^{+51}:\\
\;\;\;\;1 + \frac{\frac{x \cdot \left(0.015625 - \left(x \cdot \left(x \cdot x\right)\right) \cdot \left(t\_3 \cdot \left(t\_0 \cdot t\_2\right)\right)\right)}{0.125 - t\_3}}{0.25 + t\_1 \cdot \left(-0.5 + t\_1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(0.004629629629629629 + x \cdot \left(\left(x \cdot x\right) \cdot 7.233796296296296 \cdot 10^{-5}\right)\right) \cdot \left(x \cdot \left(x \cdot \left(36 + x \cdot 9\right)\right)\right)\\
\end{array}
\end{array}
if x < -1.55000000000000004Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f641.1%
Simplified1.1%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr1.1%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6418.8%
Simplified18.8%
if -1.55000000000000004 < x < 5e51Initial program 16.5%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.2%
Simplified88.2%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr89.6%
Applied egg-rr94.0%
if 5e51 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6487.4%
Simplified87.4%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.4%
Simplified87.4%
Applied egg-rr20.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification76.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (+ (* x (+ 0.16666666666666666 (* x 0.041666666666666664))) 0.5))
(t_1 (* t_0 (* (* x x) t_0))))
(if (<= x -1.55)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 5e+34)
(/ (- 1.0 (* t_1 t_1)) (* (+ 1.0 t_1) (- 1.0 (* x t_0))))
(*
(+ 0.004629629629629629 (* x (* (* x x) 7.233796296296296e-5)))
(* x (* x (+ 36.0 (* x 9.0)))))))))
double code(double x) {
double t_0 = (x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5;
double t_1 = t_0 * ((x * x) * t_0);
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5e+34) {
tmp = (1.0 - (t_1 * t_1)) / ((1.0 + t_1) * (1.0 - (x * t_0)));
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x * (0.16666666666666666d0 + (x * 0.041666666666666664d0))) + 0.5d0
t_1 = t_0 * ((x * x) * t_0)
if (x <= (-1.55d0)) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 5d+34) then
tmp = (1.0d0 - (t_1 * t_1)) / ((1.0d0 + t_1) * (1.0d0 - (x * t_0)))
else
tmp = (0.004629629629629629d0 + (x * ((x * x) * 7.233796296296296d-5))) * (x * (x * (36.0d0 + (x * 9.0d0))))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = (x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5;
double t_1 = t_0 * ((x * x) * t_0);
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5e+34) {
tmp = (1.0 - (t_1 * t_1)) / ((1.0 + t_1) * (1.0 - (x * t_0)));
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
def code(x): t_0 = (x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5 t_1 = t_0 * ((x * x) * t_0) tmp = 0 if x <= -1.55: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 5e+34: tmp = (1.0 - (t_1 * t_1)) / ((1.0 + t_1) * (1.0 - (x * t_0))) else: tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))) return tmp
function code(x) t_0 = Float64(Float64(x * Float64(0.16666666666666666 + Float64(x * 0.041666666666666664))) + 0.5) t_1 = Float64(t_0 * Float64(Float64(x * x) * t_0)) tmp = 0.0 if (x <= -1.55) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 5e+34) tmp = Float64(Float64(1.0 - Float64(t_1 * t_1)) / Float64(Float64(1.0 + t_1) * Float64(1.0 - Float64(x * t_0)))); else tmp = Float64(Float64(0.004629629629629629 + Float64(x * Float64(Float64(x * x) * 7.233796296296296e-5))) * Float64(x * Float64(x * Float64(36.0 + Float64(x * 9.0))))); end return tmp end
function tmp_2 = code(x) t_0 = (x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5; t_1 = t_0 * ((x * x) * t_0); tmp = 0.0; if (x <= -1.55) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 5e+34) tmp = (1.0 - (t_1 * t_1)) / ((1.0 + t_1) * (1.0 - (x * t_0))); else tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[(x * N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[(N[(x * x), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.55], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 5e+34], N[(N[(1.0 - N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + t$95$1), $MachinePrecision] * N[(1.0 - N[(x * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.004629629629629629 + N[(x * N[(N[(x * x), $MachinePrecision] * 7.233796296296296e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(36.0 + N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(0.16666666666666666 + x \cdot 0.041666666666666664\right) + 0.5\\
t_1 := t\_0 \cdot \left(\left(x \cdot x\right) \cdot t\_0\right)\\
\mathbf{if}\;x \leq -1.55:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 5 \cdot 10^{+34}:\\
\;\;\;\;\frac{1 - t\_1 \cdot t\_1}{\left(1 + t\_1\right) \cdot \left(1 - x \cdot t\_0\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(0.004629629629629629 + x \cdot \left(\left(x \cdot x\right) \cdot 7.233796296296296 \cdot 10^{-5}\right)\right) \cdot \left(x \cdot \left(x \cdot \left(36 + x \cdot 9\right)\right)\right)\\
\end{array}
\end{array}
if x < -1.55000000000000004Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f641.1%
Simplified1.1%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr1.1%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6418.8%
Simplified18.8%
if -1.55000000000000004 < x < 4.9999999999999998e34Initial program 14.6%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6490.1%
Simplified90.1%
clear-numN/A
associate-/r*N/A
*-inversesN/A
flip-+N/A
clear-numN/A
clear-numN/A
div-invN/A
Applied egg-rr93.1%
if 4.9999999999999998e34 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6483.5%
Simplified83.5%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.5%
Simplified83.5%
Applied egg-rr19.2%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6495.7%
Simplified95.7%
Final simplification75.4%
(FPCore (x)
:precision binary64
(if (<= x 1.75)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(*
(+ 0.004629629629629629 (* x (* (* x x) 7.233796296296296e-5)))
(* x (* x (+ 36.0 (* x 9.0)))))))
double code(double x) {
double tmp;
if (x <= 1.75) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.75d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = (0.004629629629629629d0 + (x * ((x * x) * 7.233796296296296d-5))) * (x * (x * (36.0d0 + (x * 9.0d0))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.75) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.75: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))) return tmp
function code(x) tmp = 0.0 if (x <= 1.75) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(0.004629629629629629 + Float64(x * Float64(Float64(x * x) * 7.233796296296296e-5))) * Float64(x * Float64(x * Float64(36.0 + Float64(x * 9.0))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.75) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x * (36.0 + (x * 9.0)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.75], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.004629629629629629 + N[(x * N[(N[(x * x), $MachinePrecision] * 7.233796296296296e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x * N[(36.0 + N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.75:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(0.004629629629629629 + x \cdot \left(\left(x \cdot x\right) \cdot 7.233796296296296 \cdot 10^{-5}\right)\right) \cdot \left(x \cdot \left(x \cdot \left(36 + x \cdot 9\right)\right)\right)\\
\end{array}
\end{array}
if x < 1.75Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.75 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
Applied egg-rr16.6%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6480.1%
Simplified80.1%
Final simplification73.9%
(FPCore (x)
:precision binary64
(if (<= x 1.8)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(*
(+ 0.004629629629629629 (* x (* (* x x) 7.233796296296296e-5)))
(* x (/ x 0.027777777777777776)))))
double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x / 0.027777777777777776));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.8d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = (0.004629629629629629d0 + (x * ((x * x) * 7.233796296296296d-5))) * (x * (x / 0.027777777777777776d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x / 0.027777777777777776));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.8: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x / 0.027777777777777776)) return tmp
function code(x) tmp = 0.0 if (x <= 1.8) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(0.004629629629629629 + Float64(x * Float64(Float64(x * x) * 7.233796296296296e-5))) * Float64(x * Float64(x / 0.027777777777777776))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.8) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = (0.004629629629629629 + (x * ((x * x) * 7.233796296296296e-5))) * (x * (x / 0.027777777777777776)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.8], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.004629629629629629 + N[(x * N[(N[(x * x), $MachinePrecision] * 7.233796296296296e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x * N[(x / 0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.8:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(0.004629629629629629 + x \cdot \left(\left(x \cdot x\right) \cdot 7.233796296296296 \cdot 10^{-5}\right)\right) \cdot \left(x \cdot \frac{x}{0.027777777777777776}\right)\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.80000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
Applied egg-rr16.6%
Taylor expanded in x around 0
Simplified78.7%
Final simplification73.5%
(FPCore (x) :precision binary64 (if (<= x 1.8) (/ (/ x (+ 1.0 (* x -0.5))) x) (/ (* (* x (* x x)) (+ 0.16666666666666666 (* x 0.041666666666666664))) x)))
double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = ((x * (x * x)) * (0.16666666666666666 + (x * 0.041666666666666664))) / x;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.8d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = ((x * (x * x)) * (0.16666666666666666d0 + (x * 0.041666666666666664d0))) / x
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = ((x * (x * x)) * (0.16666666666666666 + (x * 0.041666666666666664))) / x;
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.8: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = ((x * (x * x)) * (0.16666666666666666 + (x * 0.041666666666666664))) / x return tmp
function code(x) tmp = 0.0 if (x <= 1.8) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(Float64(x * Float64(x * x)) * Float64(0.16666666666666666 + Float64(x * 0.041666666666666664))) / x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.8) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = ((x * (x * x)) * (0.16666666666666666 + (x * 0.041666666666666664))) / x; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.8], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.8:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(0.16666666666666666 + x \cdot 0.041666666666666664\right)}{x}\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.80000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
Final simplification70.8%
(FPCore (x) :precision binary64 (if (<= x 1.6) (/ (/ x (+ 1.0 (* x -0.5))) x) (/ (* x (+ 1.0 (* x (* x (* x 0.041666666666666664))))) x)))
double code(double x) {
double tmp;
if (x <= 1.6) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (x * (1.0 + (x * (x * (x * 0.041666666666666664))))) / x;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.6d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = (x * (1.0d0 + (x * (x * (x * 0.041666666666666664d0))))) / x
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.6) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (x * (1.0 + (x * (x * (x * 0.041666666666666664))))) / x;
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.6: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = (x * (1.0 + (x * (x * (x * 0.041666666666666664))))) / x return tmp
function code(x) tmp = 0.0 if (x <= 1.6) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(x * 0.041666666666666664))))) / x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.6) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = (x * (1.0 + (x * (x * (x * 0.041666666666666664))))) / x; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.6], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x * N[(1.0 + N[(x * N[(x * N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.6:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \left(1 + x \cdot \left(x \cdot \left(x \cdot 0.041666666666666664\right)\right)\right)}{x}\\
\end{array}
\end{array}
if x < 1.6000000000000001Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.6000000000000001 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
(FPCore (x) :precision binary64 (if (<= x 1.95) (/ (/ x (+ 1.0 (* x -0.5))) x) (/ (* 0.041666666666666664 (* x (* x (* x x)))) x)))
double code(double x) {
double tmp;
if (x <= 1.95) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.041666666666666664 * (x * (x * (x * x)))) / x;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.95d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = (0.041666666666666664d0 * (x * (x * (x * x)))) / x
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.95) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (0.041666666666666664 * (x * (x * (x * x)))) / x;
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.95: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = (0.041666666666666664 * (x * (x * (x * x)))) / x return tmp
function code(x) tmp = 0.0 if (x <= 1.95) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(0.041666666666666664 * Float64(x * Float64(x * Float64(x * x)))) / x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.95) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = (0.041666666666666664 * (x * (x * (x * x)))) / x; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.95], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.041666666666666664 * N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.95:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.041666666666666664 \cdot \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right)}{x}\\
\end{array}
\end{array}
if x < 1.94999999999999996Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.94999999999999996 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
(FPCore (x) :precision binary64 (if (<= x 4.0) (/ (/ x (+ 1.0 (* x -0.5))) x) (* x (/ (* x x) (+ 24.0 (/ -96.0 x))))))
double code(double x) {
double tmp;
if (x <= 4.0) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = x * ((x * x) / (24.0 + (-96.0 / x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 4.0d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = x * ((x * x) / (24.0d0 + ((-96.0d0) / x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 4.0) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = x * ((x * x) / (24.0 + (-96.0 / x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= 4.0: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = x * ((x * x) / (24.0 + (-96.0 / x))) return tmp
function code(x) tmp = 0.0 if (x <= 4.0) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(x * Float64(Float64(x * x) / Float64(24.0 + Float64(-96.0 / x)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 4.0) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = x * ((x * x) / (24.0 + (-96.0 / x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 4.0], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(N[(x * x), $MachinePrecision] / N[(24.0 + N[(-96.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{x \cdot x}{24 + \frac{-96}{x}}\\
\end{array}
\end{array}
if x < 4Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 4 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr69.8%
Taylor expanded in x around inf
/-lowering-/.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
div-invN/A
associate-/r/N/A
associate-*l*N/A
cube-unmultN/A
inv-powN/A
pow-prod-upN/A
metadata-evalN/A
pow2N/A
clear-numN/A
associate-/r/N/A
associate-/r*N/A
clear-numN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval62.6%
Applied egg-rr62.6%
(FPCore (x) :precision binary64 (if (<= x 1.55) (/ (/ x (+ 1.0 (* x -0.5))) x) (* x (+ (* x (+ 0.16666666666666666 (* x 0.041666666666666664))) 0.5))))
double code(double x) {
double tmp;
if (x <= 1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = x * ((x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.55d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = x * ((x * (0.16666666666666666d0 + (x * 0.041666666666666664d0))) + 0.5d0)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = x * ((x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5);
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.55: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = x * ((x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5) return tmp
function code(x) tmp = 0.0 if (x <= 1.55) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(x * Float64(Float64(x * Float64(0.16666666666666666 + Float64(x * 0.041666666666666664))) + 0.5)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.55) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = x * ((x * (0.16666666666666666 + (x * 0.041666666666666664))) + 0.5); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.55], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(N[(x * N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.55:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(0.16666666666666666 + x \cdot 0.041666666666666664\right) + 0.5\right)\\
\end{array}
\end{array}
if x < 1.55000000000000004Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.55000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
Simplified62.6%
Final simplification68.7%
(FPCore (x) :precision binary64 (if (<= x 1.8) (/ (/ x (+ 1.0 (* x -0.5))) x) (* (* x x) (+ 0.16666666666666666 (* x 0.041666666666666664)))))
double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.8d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = (x * x) * (0.16666666666666666d0 + (x * 0.041666666666666664d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.8: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664)) return tmp
function code(x) tmp = 0.0 if (x <= 1.8) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * 0.041666666666666664))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.8) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.8], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.8:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot 0.041666666666666664\right)\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.2%
Simplified65.2%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr65.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if 1.80000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.6%
Simplified62.6%
Final simplification68.7%
(FPCore (x) :precision binary64 (if (<= x 2.0) 1.0 (* (* x x) (+ 0.16666666666666666 (* x 0.041666666666666664)))))
double code(double x) {
double tmp;
if (x <= 2.0) {
tmp = 1.0;
} else {
tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 2.0d0) then
tmp = 1.0d0
else
tmp = (x * x) * (0.16666666666666666d0 + (x * 0.041666666666666664d0))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 2.0) {
tmp = 1.0;
} else {
tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664));
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.0: tmp = 1.0 else: tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664)) return tmp
function code(x) tmp = 0.0 if (x <= 2.0) tmp = 1.0; else tmp = Float64(Float64(x * x) * Float64(0.16666666666666666 + Float64(x * 0.041666666666666664))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.0) tmp = 1.0; else tmp = (x * x) * (0.16666666666666666 + (x * 0.041666666666666664)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.0], 1.0, N[(N[(x * x), $MachinePrecision] * N[(0.16666666666666666 + N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(0.16666666666666666 + x \cdot 0.041666666666666664\right)\\
\end{array}
\end{array}
if x < 2Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified66.2%
if 2 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.6%
Simplified62.6%
Final simplification65.1%
(FPCore (x) :precision binary64 (if (<= x 2.9) 1.0 (* (* x (* x x)) 0.041666666666666664)))
double code(double x) {
double tmp;
if (x <= 2.9) {
tmp = 1.0;
} else {
tmp = (x * (x * x)) * 0.041666666666666664;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 2.9d0) then
tmp = 1.0d0
else
tmp = (x * (x * x)) * 0.041666666666666664d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 2.9) {
tmp = 1.0;
} else {
tmp = (x * (x * x)) * 0.041666666666666664;
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.9: tmp = 1.0 else: tmp = (x * (x * x)) * 0.041666666666666664 return tmp
function code(x) tmp = 0.0 if (x <= 2.9) tmp = 1.0; else tmp = Float64(Float64(x * Float64(x * x)) * 0.041666666666666664); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.9) tmp = 1.0; else tmp = (x * (x * x)) * 0.041666666666666664; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.9], 1.0, N[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * 0.041666666666666664), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.9:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot \left(x \cdot x\right)\right) \cdot 0.041666666666666664\\
\end{array}
\end{array}
if x < 2.89999999999999991Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified66.2%
if 2.89999999999999991 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8%
Simplified69.8%
Taylor expanded in x around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.6%
Simplified62.6%
Final simplification65.1%
(FPCore (x) :precision binary64 (if (<= x 2.4) 1.0 (* (* x x) 0.16666666666666666)))
double code(double x) {
double tmp;
if (x <= 2.4) {
tmp = 1.0;
} else {
tmp = (x * x) * 0.16666666666666666;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 2.4d0) then
tmp = 1.0d0
else
tmp = (x * x) * 0.16666666666666666d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 2.4) {
tmp = 1.0;
} else {
tmp = (x * x) * 0.16666666666666666;
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.4: tmp = 1.0 else: tmp = (x * x) * 0.16666666666666666 return tmp
function code(x) tmp = 0.0 if (x <= 2.4) tmp = 1.0; else tmp = Float64(Float64(x * x) * 0.16666666666666666); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.4) tmp = 1.0; else tmp = (x * x) * 0.16666666666666666; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.4], 1.0, N[(N[(x * x), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.4:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.16666666666666666\\
\end{array}
\end{array}
if x < 2.39999999999999991Initial program 38.3%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified66.2%
if 2.39999999999999991 < x Initial program 100.0%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
remove-double-negN/A
distribute-lft-neg-outN/A
unsub-negN/A
--lowering--.f64N/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-eval48.9%
Simplified48.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.9%
Simplified48.9%
Final simplification61.0%
(FPCore (x) :precision binary64 (+ 1.0 (* x (* x (* x 0.041666666666666664)))))
double code(double x) {
return 1.0 + (x * (x * (x * 0.041666666666666664)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 + (x * (x * (x * 0.041666666666666664d0)))
end function
public static double code(double x) {
return 1.0 + (x * (x * (x * 0.041666666666666664)));
}
def code(x): return 1.0 + (x * (x * (x * 0.041666666666666664)))
function code(x) return Float64(1.0 + Float64(x * Float64(x * Float64(x * 0.041666666666666664)))) end
function tmp = code(x) tmp = 1.0 + (x * (x * (x * 0.041666666666666664))); end
code[x_] := N[(1.0 + N[(x * N[(x * N[(x * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + x \cdot \left(x \cdot \left(x \cdot 0.041666666666666664\right)\right)
\end{array}
Initial program 56.6%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6464.5%
Simplified64.5%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
(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 56.6%
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified47.5%
(FPCore (x) :precision binary64 (let* ((t_0 (- (exp x) 1.0))) (if (and (< x 1.0) (> x -1.0)) (/ t_0 (log (exp x))) (/ t_0 x))))
double code(double x) {
double t_0 = exp(x) - 1.0;
double tmp;
if ((x < 1.0) && (x > -1.0)) {
tmp = t_0 / log(exp(x));
} else {
tmp = t_0 / x;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = exp(x) - 1.0d0
if ((x < 1.0d0) .and. (x > (-1.0d0))) then
tmp = t_0 / log(exp(x))
else
tmp = t_0 / x
end if
code = tmp
end function
public static double code(double x) {
double t_0 = Math.exp(x) - 1.0;
double tmp;
if ((x < 1.0) && (x > -1.0)) {
tmp = t_0 / Math.log(Math.exp(x));
} else {
tmp = t_0 / x;
}
return tmp;
}
def code(x): t_0 = math.exp(x) - 1.0 tmp = 0 if (x < 1.0) and (x > -1.0): tmp = t_0 / math.log(math.exp(x)) else: tmp = t_0 / x return tmp
function code(x) t_0 = Float64(exp(x) - 1.0) tmp = 0.0 if ((x < 1.0) && (x > -1.0)) tmp = Float64(t_0 / log(exp(x))); else tmp = Float64(t_0 / x); end return tmp end
function tmp_2 = code(x) t_0 = exp(x) - 1.0; tmp = 0.0; if ((x < 1.0) && (x > -1.0)) tmp = t_0 / log(exp(x)); else tmp = t_0 / x; end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision]}, If[And[Less[x, 1.0], Greater[x, -1.0]], N[(t$95$0 / N[Log[N[Exp[x], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$0 / x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{x} - 1\\
\mathbf{if}\;x < 1 \land x > -1:\\
\;\;\;\;\frac{t\_0}{\log \left(e^{x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_0}{x}\\
\end{array}
\end{array}
herbie shell --seed 2024149
(FPCore (x)
:name "Kahan's exp quotient"
:precision binary64
:alt
(! :herbie-platform default (if (and (< x 1) (> x -1)) (/ (- (exp x) 1) (log (exp x))) (/ (- (exp x) 1) x)))
(/ (- (exp x) 1.0) x))