
(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 14 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 53.2%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(let* ((t_0
(+ -0.5 (* x (+ (* x -0.041666666666666664) -0.16666666666666666))))
(t_1 (* x t_0)))
(if (<= x -1.55)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 1.65e+103)
(/ (/ (* x (- 1.0 (* x (* t_0 t_1)))) (+ 1.0 t_1)) x)
(* x (* 0.041666666666666664 (* x x)))))))
double code(double x) {
double t_0 = -0.5 + (x * ((x * -0.041666666666666664) + -0.16666666666666666));
double t_1 = x * t_0;
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 1.65e+103) {
tmp = ((x * (1.0 - (x * (t_0 * t_1)))) / (1.0 + t_1)) / x;
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (-0.5d0) + (x * ((x * (-0.041666666666666664d0)) + (-0.16666666666666666d0)))
t_1 = x * t_0
if (x <= (-1.55d0)) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 1.65d+103) then
tmp = ((x * (1.0d0 - (x * (t_0 * t_1)))) / (1.0d0 + t_1)) / x
else
tmp = x * (0.041666666666666664d0 * (x * x))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = -0.5 + (x * ((x * -0.041666666666666664) + -0.16666666666666666));
double t_1 = x * t_0;
double tmp;
if (x <= -1.55) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 1.65e+103) {
tmp = ((x * (1.0 - (x * (t_0 * t_1)))) / (1.0 + t_1)) / x;
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
return tmp;
}
def code(x): t_0 = -0.5 + (x * ((x * -0.041666666666666664) + -0.16666666666666666)) t_1 = x * t_0 tmp = 0 if x <= -1.55: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 1.65e+103: tmp = ((x * (1.0 - (x * (t_0 * t_1)))) / (1.0 + t_1)) / x else: tmp = x * (0.041666666666666664 * (x * x)) return tmp
function code(x) t_0 = Float64(-0.5 + Float64(x * Float64(Float64(x * -0.041666666666666664) + -0.16666666666666666))) t_1 = Float64(x * t_0) tmp = 0.0 if (x <= -1.55) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 1.65e+103) tmp = Float64(Float64(Float64(x * Float64(1.0 - Float64(x * Float64(t_0 * t_1)))) / Float64(1.0 + t_1)) / x); else tmp = Float64(x * Float64(0.041666666666666664 * Float64(x * x))); end return tmp end
function tmp_2 = code(x) t_0 = -0.5 + (x * ((x * -0.041666666666666664) + -0.16666666666666666)); t_1 = x * t_0; tmp = 0.0; if (x <= -1.55) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 1.65e+103) tmp = ((x * (1.0 - (x * (t_0 * t_1)))) / (1.0 + t_1)) / x; else tmp = x * (0.041666666666666664 * (x * x)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(-0.5 + N[(x * N[(N[(x * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x * t$95$0), $MachinePrecision]}, If[LessEqual[x, -1.55], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 1.65e+103], N[(N[(N[(x * N[(1.0 - N[(x * N[(t$95$0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(0.041666666666666664 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.5 + x \cdot \left(x \cdot -0.041666666666666664 + -0.16666666666666666\right)\\
t_1 := x \cdot t\_0\\
\mathbf{if}\;x \leq -1.55:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 1.65 \cdot 10^{+103}:\\
\;\;\;\;\frac{\frac{x \cdot \left(1 - x \cdot \left(t\_0 \cdot t\_1\right)\right)}{1 + t\_1}}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0.041666666666666664 \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < -1.55000000000000004Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified1.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr1.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6418.8%
Simplified18.8%
if -1.55000000000000004 < x < 1.65000000000000004e103Initial program 22.3%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified85.8%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr95.6%
if 1.65000000000000004e103 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification76.7%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x x))) (t_1 (/ t_0 -24.0)))
(if (<= x 1.65)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 7e+51)
(/
(+ 1.0 (* 7.233796296296296e-5 (* t_0 (* t_0 t_0))))
(+ 1.0 (* t_1 (+ 1.0 t_1))))
(if (<= x 2e+154)
(/
(/
(* t_0 (+ (* 7.233796296296296e-5 t_0) 0.004629629629629629))
(+
0.027777777777777776
(*
(* x 0.041666666666666664)
(- (* x 0.041666666666666664) 0.16666666666666666))))
x)
(* x (* x 0.16666666666666666)))))))
double code(double x) {
double t_0 = x * (x * x);
double t_1 = t_0 / -24.0;
double tmp;
if (x <= 1.65) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 7e+51) {
tmp = (1.0 + (7.233796296296296e-5 * (t_0 * (t_0 * t_0)))) / (1.0 + (t_1 * (1.0 + t_1)));
} else if (x <= 2e+154) {
tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x;
} else {
tmp = x * (x * 0.16666666666666666);
}
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 = t_0 / (-24.0d0)
if (x <= 1.65d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 7d+51) then
tmp = (1.0d0 + (7.233796296296296d-5 * (t_0 * (t_0 * t_0)))) / (1.0d0 + (t_1 * (1.0d0 + t_1)))
else if (x <= 2d+154) then
tmp = ((t_0 * ((7.233796296296296d-5 * t_0) + 0.004629629629629629d0)) / (0.027777777777777776d0 + ((x * 0.041666666666666664d0) * ((x * 0.041666666666666664d0) - 0.16666666666666666d0)))) / x
else
tmp = x * (x * 0.16666666666666666d0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * x);
double t_1 = t_0 / -24.0;
double tmp;
if (x <= 1.65) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 7e+51) {
tmp = (1.0 + (7.233796296296296e-5 * (t_0 * (t_0 * t_0)))) / (1.0 + (t_1 * (1.0 + t_1)));
} else if (x <= 2e+154) {
tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x;
} else {
tmp = x * (x * 0.16666666666666666);
}
return tmp;
}
def code(x): t_0 = x * (x * x) t_1 = t_0 / -24.0 tmp = 0 if x <= 1.65: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 7e+51: tmp = (1.0 + (7.233796296296296e-5 * (t_0 * (t_0 * t_0)))) / (1.0 + (t_1 * (1.0 + t_1))) elif x <= 2e+154: tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x else: tmp = x * (x * 0.16666666666666666) return tmp
function code(x) t_0 = Float64(x * Float64(x * x)) t_1 = Float64(t_0 / -24.0) tmp = 0.0 if (x <= 1.65) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 7e+51) tmp = Float64(Float64(1.0 + Float64(7.233796296296296e-5 * Float64(t_0 * Float64(t_0 * t_0)))) / Float64(1.0 + Float64(t_1 * Float64(1.0 + t_1)))); elseif (x <= 2e+154) tmp = Float64(Float64(Float64(t_0 * Float64(Float64(7.233796296296296e-5 * t_0) + 0.004629629629629629)) / Float64(0.027777777777777776 + Float64(Float64(x * 0.041666666666666664) * Float64(Float64(x * 0.041666666666666664) - 0.16666666666666666)))) / x); else tmp = Float64(x * Float64(x * 0.16666666666666666)); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * x); t_1 = t_0 / -24.0; tmp = 0.0; if (x <= 1.65) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 7e+51) tmp = (1.0 + (7.233796296296296e-5 * (t_0 * (t_0 * t_0)))) / (1.0 + (t_1 * (1.0 + t_1))); elseif (x <= 2e+154) tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x; else tmp = x * (x * 0.16666666666666666); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 / -24.0), $MachinePrecision]}, If[LessEqual[x, 1.65], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 7e+51], N[(N[(1.0 + N[(7.233796296296296e-5 * N[(t$95$0 * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$1 * N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2e+154], N[(N[(N[(t$95$0 * N[(N[(7.233796296296296e-5 * t$95$0), $MachinePrecision] + 0.004629629629629629), $MachinePrecision]), $MachinePrecision] / N[(0.027777777777777776 + N[(N[(x * 0.041666666666666664), $MachinePrecision] * N[(N[(x * 0.041666666666666664), $MachinePrecision] - 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
t_1 := \frac{t\_0}{-24}\\
\mathbf{if}\;x \leq 1.65:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 7 \cdot 10^{+51}:\\
\;\;\;\;\frac{1 + 7.233796296296296 \cdot 10^{-5} \cdot \left(t\_0 \cdot \left(t\_0 \cdot t\_0\right)\right)}{1 + t\_1 \cdot \left(1 + t\_1\right)}\\
\mathbf{elif}\;x \leq 2 \cdot 10^{+154}:\\
\;\;\;\;\frac{\frac{t\_0 \cdot \left(7.233796296296296 \cdot 10^{-5} \cdot t\_0 + 0.004629629629629629\right)}{0.027777777777777776 + \left(x \cdot 0.041666666666666664\right) \cdot \left(x \cdot 0.041666666666666664 - 0.16666666666666666\right)}}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if x < 1.6499999999999999Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.6499999999999999 < x < 7e51Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified4.1%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr4.1%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f644.1%
Simplified4.1%
flip3--N/A
/-lowering-/.f64N/A
Applied egg-rr43.8%
if 7e51 < x < 2.00000000000000007e154Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified56.6%
Taylor expanded in x around inf
*-commutativeN/A
metadata-evalN/A
pow-plusN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6456.6%
Simplified56.6%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr100.0%
if 2.00000000000000007e154 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
distribute-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
metadata-evalN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification76.7%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x x))))
(if (<= x 1.8)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 2e+154)
(/
(/
(* t_0 (+ (* 7.233796296296296e-5 t_0) 0.004629629629629629))
(+
0.027777777777777776
(*
(* x 0.041666666666666664)
(- (* x 0.041666666666666664) 0.16666666666666666))))
x)
(* x (* x 0.16666666666666666))))))
double code(double x) {
double t_0 = x * (x * x);
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 2e+154) {
tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x;
} else {
tmp = x * (x * 0.16666666666666666);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x * (x * x)
if (x <= 1.8d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 2d+154) then
tmp = ((t_0 * ((7.233796296296296d-5 * t_0) + 0.004629629629629629d0)) / (0.027777777777777776d0 + ((x * 0.041666666666666664d0) * ((x * 0.041666666666666664d0) - 0.16666666666666666d0)))) / x
else
tmp = x * (x * 0.16666666666666666d0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * x);
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 2e+154) {
tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x;
} else {
tmp = x * (x * 0.16666666666666666);
}
return tmp;
}
def code(x): t_0 = x * (x * x) tmp = 0 if x <= 1.8: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 2e+154: tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x else: tmp = x * (x * 0.16666666666666666) return tmp
function code(x) t_0 = Float64(x * Float64(x * x)) tmp = 0.0 if (x <= 1.8) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 2e+154) tmp = Float64(Float64(Float64(t_0 * Float64(Float64(7.233796296296296e-5 * t_0) + 0.004629629629629629)) / Float64(0.027777777777777776 + Float64(Float64(x * 0.041666666666666664) * Float64(Float64(x * 0.041666666666666664) - 0.16666666666666666)))) / x); else tmp = Float64(x * Float64(x * 0.16666666666666666)); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * x); tmp = 0.0; if (x <= 1.8) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 2e+154) tmp = ((t_0 * ((7.233796296296296e-5 * t_0) + 0.004629629629629629)) / (0.027777777777777776 + ((x * 0.041666666666666664) * ((x * 0.041666666666666664) - 0.16666666666666666)))) / x; else tmp = x * (x * 0.16666666666666666); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.8], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 2e+154], N[(N[(N[(t$95$0 * N[(N[(7.233796296296296e-5 * t$95$0), $MachinePrecision] + 0.004629629629629629), $MachinePrecision]), $MachinePrecision] / N[(0.027777777777777776 + N[(N[(x * 0.041666666666666664), $MachinePrecision] * N[(N[(x * 0.041666666666666664), $MachinePrecision] - 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 1.8:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 2 \cdot 10^{+154}:\\
\;\;\;\;\frac{\frac{t\_0 \cdot \left(7.233796296296296 \cdot 10^{-5} \cdot t\_0 + 0.004629629629629629\right)}{0.027777777777777776 + \left(x \cdot 0.041666666666666664\right) \cdot \left(x \cdot 0.041666666666666664 - 0.16666666666666666\right)}}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.80000000000000004 < x < 2.00000000000000007e154Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified39.1%
Taylor expanded in x around inf
*-commutativeN/A
metadata-evalN/A
pow-plusN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6439.1%
Simplified39.1%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr68.0%
if 2.00000000000000007e154 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
distribute-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
metadata-evalN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification74.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x x))))
(if (<= x 1.65)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(if (<= x 5.5e+102)
(/
(- 1.0 (* (* t_0 t_0) 0.001736111111111111))
(- 1.0 (* 0.041666666666666664 t_0)))
(* x (* 0.041666666666666664 (* x x)))))))
double code(double x) {
double t_0 = x * (x * x);
double tmp;
if (x <= 1.65) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5.5e+102) {
tmp = (1.0 - ((t_0 * t_0) * 0.001736111111111111)) / (1.0 - (0.041666666666666664 * t_0));
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x * (x * x)
if (x <= 1.65d0) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else if (x <= 5.5d+102) then
tmp = (1.0d0 - ((t_0 * t_0) * 0.001736111111111111d0)) / (1.0d0 - (0.041666666666666664d0 * t_0))
else
tmp = x * (0.041666666666666664d0 * (x * x))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * x);
double tmp;
if (x <= 1.65) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else if (x <= 5.5e+102) {
tmp = (1.0 - ((t_0 * t_0) * 0.001736111111111111)) / (1.0 - (0.041666666666666664 * t_0));
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
return tmp;
}
def code(x): t_0 = x * (x * x) tmp = 0 if x <= 1.65: tmp = (x / (1.0 + (x * -0.5))) / x elif x <= 5.5e+102: tmp = (1.0 - ((t_0 * t_0) * 0.001736111111111111)) / (1.0 - (0.041666666666666664 * t_0)) else: tmp = x * (0.041666666666666664 * (x * x)) return tmp
function code(x) t_0 = Float64(x * Float64(x * x)) tmp = 0.0 if (x <= 1.65) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); elseif (x <= 5.5e+102) tmp = Float64(Float64(1.0 - Float64(Float64(t_0 * t_0) * 0.001736111111111111)) / Float64(1.0 - Float64(0.041666666666666664 * t_0))); else tmp = Float64(x * Float64(0.041666666666666664 * Float64(x * x))); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * x); tmp = 0.0; if (x <= 1.65) tmp = (x / (1.0 + (x * -0.5))) / x; elseif (x <= 5.5e+102) tmp = (1.0 - ((t_0 * t_0) * 0.001736111111111111)) / (1.0 - (0.041666666666666664 * t_0)); else tmp = x * (0.041666666666666664 * (x * x)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.65], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 5.5e+102], N[(N[(1.0 - N[(N[(t$95$0 * t$95$0), $MachinePrecision] * 0.001736111111111111), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(0.041666666666666664 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.041666666666666664 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot x\right)\\
\mathbf{if}\;x \leq 1.65:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{elif}\;x \leq 5.5 \cdot 10^{+102}:\\
\;\;\;\;\frac{1 - \left(t\_0 \cdot t\_0\right) \cdot 0.001736111111111111}{1 - 0.041666666666666664 \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0.041666666666666664 \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < 1.6499999999999999Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.6499999999999999 < x < 5.49999999999999981e102Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified5.2%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr43.9%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f645.2%
Simplified5.2%
sub-negN/A
flip-+N/A
sqr-negN/A
/-lowering-/.f64N/A
Applied egg-rr60.3%
if 5.49999999999999981e102 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(if (<= x -2.2)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(+
1.0
(/
-1.0
(/
(+ (* x -0.16666666666666666) 0.5)
(* x (+ (* x (* x (* (* x x) 0.001736111111111111))) -0.25)))))))
double code(double x) {
double tmp;
if (x <= -2.2) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = 1.0 + (-1.0 / (((x * -0.16666666666666666) + 0.5) / (x * ((x * (x * ((x * x) * 0.001736111111111111))) + -0.25))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-2.2d0)) then
tmp = (x / (1.0d0 + (x * (-0.5d0)))) / x
else
tmp = 1.0d0 + ((-1.0d0) / (((x * (-0.16666666666666666d0)) + 0.5d0) / (x * ((x * (x * ((x * x) * 0.001736111111111111d0))) + (-0.25d0)))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -2.2) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = 1.0 + (-1.0 / (((x * -0.16666666666666666) + 0.5) / (x * ((x * (x * ((x * x) * 0.001736111111111111))) + -0.25))));
}
return tmp;
}
def code(x): tmp = 0 if x <= -2.2: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = 1.0 + (-1.0 / (((x * -0.16666666666666666) + 0.5) / (x * ((x * (x * ((x * x) * 0.001736111111111111))) + -0.25)))) return tmp
function code(x) tmp = 0.0 if (x <= -2.2) tmp = Float64(Float64(x / Float64(1.0 + Float64(x * -0.5))) / x); else tmp = Float64(1.0 + Float64(-1.0 / Float64(Float64(Float64(x * -0.16666666666666666) + 0.5) / Float64(x * Float64(Float64(x * Float64(x * Float64(Float64(x * x) * 0.001736111111111111))) + -0.25))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -2.2) tmp = (x / (1.0 + (x * -0.5))) / x; else tmp = 1.0 + (-1.0 / (((x * -0.16666666666666666) + 0.5) / (x * ((x * (x * ((x * x) * 0.001736111111111111))) + -0.25)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -2.2], N[(N[(x / N[(1.0 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(N[(N[(x * -0.16666666666666666), $MachinePrecision] + 0.5), $MachinePrecision] / N[(x * N[(N[(x * N[(x * N[(N[(x * x), $MachinePrecision] * 0.001736111111111111), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.2:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{\frac{x \cdot -0.16666666666666666 + 0.5}{x \cdot \left(x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot 0.001736111111111111\right)\right) + -0.25\right)}}\\
\end{array}
\end{array}
if x < -2.2000000000000002Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified1.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr1.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6418.8%
Simplified18.8%
if -2.2000000000000002 < x Initial program 37.3%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified85.6%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr75.8%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6491.5%
Simplified91.5%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.5%
Simplified91.5%
Final simplification73.1%
(FPCore (x)
:precision binary64
(if (<= x 1.8)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(/
(/
(*
(* x (* x x))
(- 0.027777777777777776 (* (* x x) 0.001736111111111111)))
(+ (* x -0.041666666666666664) 0.16666666666666666))
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.027777777777777776 - ((x * x) * 0.001736111111111111))) / ((x * -0.041666666666666664) + 0.16666666666666666)) / 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.027777777777777776d0 - ((x * x) * 0.001736111111111111d0))) / ((x * (-0.041666666666666664d0)) + 0.16666666666666666d0)) / 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.027777777777777776 - ((x * x) * 0.001736111111111111))) / ((x * -0.041666666666666664) + 0.16666666666666666)) / 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.027777777777777776 - ((x * x) * 0.001736111111111111))) / ((x * -0.041666666666666664) + 0.16666666666666666)) / 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(Float64(x * Float64(x * x)) * Float64(0.027777777777777776 - Float64(Float64(x * x) * 0.001736111111111111))) / Float64(Float64(x * -0.041666666666666664) + 0.16666666666666666)) / 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.027777777777777776 - ((x * x) * 0.001736111111111111))) / ((x * -0.041666666666666664) + 0.16666666666666666)) / 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[(N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(0.027777777777777776 - N[(N[(x * x), $MachinePrecision] * 0.001736111111111111), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(x * -0.041666666666666664), $MachinePrecision] + 0.16666666666666666), $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{\frac{\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(0.027777777777777776 - \left(x \cdot x\right) \cdot 0.001736111111111111\right)}{x \cdot -0.041666666666666664 + 0.16666666666666666}}{x}\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.80000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.8%
Taylor expanded in x around inf
*-commutativeN/A
metadata-evalN/A
pow-plusN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.8%
Simplified66.8%
*-commutativeN/A
flip-+N/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
*-commutativeN/A
metadata-evalN/A
cancel-sign-sub-invN/A
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr75.4%
Final simplification73.0%
(FPCore (x)
:precision binary64
(if (<= x 1.8)
(/ (/ x (+ 1.0 (* x -0.5))) x)
(*
(* (* x x) (* x x))
(- 0.010416666666666666 (/ -0.11458333333333333 x)))))
double code(double x) {
double tmp;
if (x <= 1.8) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = ((x * x) * (x * x)) * (0.010416666666666666 - (-0.11458333333333333 / 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 * x)) * (0.010416666666666666d0 - ((-0.11458333333333333d0) / 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 * x)) * (0.010416666666666666 - (-0.11458333333333333 / 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 * x)) * (0.010416666666666666 - (-0.11458333333333333 / 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 * x) * Float64(x * x)) * Float64(0.010416666666666666 - Float64(-0.11458333333333333 / 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 * x)) * (0.010416666666666666 - (-0.11458333333333333 / 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 * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(0.010416666666666666 - N[(-0.11458333333333333 / x), $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(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(0.010416666666666666 - \frac{-0.11458333333333333}{x}\right)\\
\end{array}
\end{array}
if x < 1.80000000000000004Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.80000000000000004 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified58.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr29.9%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6475.5%
Simplified75.5%
Taylor expanded in x around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
remove-double-negN/A
sub-negN/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval67.2%
Simplified67.2%
(FPCore (x) :precision binary64 (if (<= x 1.95) (/ (/ x (+ 1.0 (* x -0.5))) x) (* (* (* x x) (* x x)) 0.010416666666666666)))
double code(double x) {
double tmp;
if (x <= 1.95) {
tmp = (x / (1.0 + (x * -0.5))) / x;
} else {
tmp = ((x * x) * (x * x)) * 0.010416666666666666;
}
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 = ((x * x) * (x * x)) * 0.010416666666666666d0
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 = ((x * x) * (x * x)) * 0.010416666666666666;
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.95: tmp = (x / (1.0 + (x * -0.5))) / x else: tmp = ((x * x) * (x * x)) * 0.010416666666666666 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(Float64(x * x) * Float64(x * x)) * 0.010416666666666666); 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 = ((x * x) * (x * x)) * 0.010416666666666666; 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[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * 0.010416666666666666), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.95:\\
\;\;\;\;\frac{\frac{x}{1 + x \cdot -0.5}}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.010416666666666666\\
\end{array}
\end{array}
if x < 1.94999999999999996Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.2%
flip3--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr66.2%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6472.2%
Simplified72.2%
if 1.94999999999999996 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified58.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr29.9%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6475.5%
Simplified75.5%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.2%
Simplified67.2%
(FPCore (x) :precision binary64 (if (<= x 3.1) 1.0 (* (* (* x x) (* x x)) 0.010416666666666666)))
double code(double x) {
double tmp;
if (x <= 3.1) {
tmp = 1.0;
} else {
tmp = ((x * x) * (x * x)) * 0.010416666666666666;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 3.1d0) then
tmp = 1.0d0
else
tmp = ((x * x) * (x * x)) * 0.010416666666666666d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 3.1) {
tmp = 1.0;
} else {
tmp = ((x * x) * (x * x)) * 0.010416666666666666;
}
return tmp;
}
def code(x): tmp = 0 if x <= 3.1: tmp = 1.0 else: tmp = ((x * x) * (x * x)) * 0.010416666666666666 return tmp
function code(x) tmp = 0.0 if (x <= 3.1) tmp = 1.0; else tmp = Float64(Float64(Float64(x * x) * Float64(x * x)) * 0.010416666666666666); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 3.1) tmp = 1.0; else tmp = ((x * x) * (x * x)) * 0.010416666666666666; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 3.1], 1.0, N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * 0.010416666666666666), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.1:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.010416666666666666\\
\end{array}
\end{array}
if x < 3.10000000000000009Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified67.2%
if 3.10000000000000009 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified58.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr29.9%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6475.5%
Simplified75.5%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.2%
Simplified67.2%
(FPCore (x) :precision binary64 (if (<= x 2.0) 1.0 (* (* x x) (+ (* x 0.041666666666666664) 0.16666666666666666))))
double code(double x) {
double tmp;
if (x <= 2.0) {
tmp = 1.0;
} else {
tmp = (x * x) * ((x * 0.041666666666666664) + 0.16666666666666666);
}
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) * ((x * 0.041666666666666664d0) + 0.16666666666666666d0)
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) * ((x * 0.041666666666666664) + 0.16666666666666666);
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.0: tmp = 1.0 else: tmp = (x * x) * ((x * 0.041666666666666664) + 0.16666666666666666) return tmp
function code(x) tmp = 0.0 if (x <= 2.0) tmp = 1.0; else tmp = Float64(Float64(x * x) * Float64(Float64(x * 0.041666666666666664) + 0.16666666666666666)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.0) tmp = 1.0; else tmp = (x * x) * ((x * 0.041666666666666664) + 0.16666666666666666); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.0], 1.0, N[(N[(x * x), $MachinePrecision] * N[(N[(x * 0.041666666666666664), $MachinePrecision] + 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.041666666666666664 + 0.16666666666666666\right)\\
\end{array}
\end{array}
if x < 2Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified67.2%
if 2 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.8%
Taylor expanded in x around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6458.4%
Simplified58.4%
Final simplification64.9%
(FPCore (x) :precision binary64 (if (<= x 2.9) 1.0 (* x (* 0.041666666666666664 (* x x)))))
double code(double x) {
double tmp;
if (x <= 2.9) {
tmp = 1.0;
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
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 * (0.041666666666666664d0 * (x * x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 2.9) {
tmp = 1.0;
} else {
tmp = x * (0.041666666666666664 * (x * x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 2.9: tmp = 1.0 else: tmp = x * (0.041666666666666664 * (x * x)) return tmp
function code(x) tmp = 0.0 if (x <= 2.9) tmp = 1.0; else tmp = Float64(x * Float64(0.041666666666666664 * Float64(x * x))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 2.9) tmp = 1.0; else tmp = x * (0.041666666666666664 * (x * x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 2.9], 1.0, N[(x * N[(0.041666666666666664 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.9:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0.041666666666666664 \cdot \left(x \cdot x\right)\right)\\
\end{array}
\end{array}
if x < 2.89999999999999991Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified67.2%
if 2.89999999999999991 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-outN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified66.8%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.4%
Simplified58.4%
(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(x * Float64(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[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.4:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if x < 2.39999999999999991Initial program 37.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified67.2%
if 2.39999999999999991 < x Initial program 100.0%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
cancel-sign-subN/A
--lowering--.f64N/A
*-rgt-identityN/A
distribute-rgt-neg-inN/A
associate-*r*N/A
*-lowering-*.f64N/A
metadata-evalN/A
mul-1-negN/A
distribute-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
metadata-evalN/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f6448.1%
Simplified48.1%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6448.1%
Simplified48.1%
(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 53.2%
/-lowering-/.f64N/A
accelerator-lowering-expm1.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
Simplified50.7%
(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 2024191
(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))