
(FPCore (x eps) :precision binary64 (/ (- (* (+ 1.0 (/ 1.0 eps)) (exp (- (* (- 1.0 eps) x)))) (* (- (/ 1.0 eps) 1.0) (exp (- (* (+ 1.0 eps) x))))) 2.0))
double code(double x, double eps) {
return (((1.0 + (1.0 / eps)) * exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * exp(-((1.0 + eps) * x)))) / 2.0;
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, eps)
use fmin_fmax_functions
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (((1.0d0 + (1.0d0 / eps)) * exp(-((1.0d0 - eps) * x))) - (((1.0d0 / eps) - 1.0d0) * exp(-((1.0d0 + eps) * x)))) / 2.0d0
end function
public static double code(double x, double eps) {
return (((1.0 + (1.0 / eps)) * Math.exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * Math.exp(-((1.0 + eps) * x)))) / 2.0;
}
def code(x, eps): return (((1.0 + (1.0 / eps)) * math.exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * math.exp(-((1.0 + eps) * x)))) / 2.0
function code(x, eps) return Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps)) * exp(Float64(-Float64(Float64(1.0 - eps) * x)))) - Float64(Float64(Float64(1.0 / eps) - 1.0) * exp(Float64(-Float64(Float64(1.0 + eps) * x))))) / 2.0) end
function tmp = code(x, eps) tmp = (((1.0 + (1.0 / eps)) * exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * exp(-((1.0 + eps) * x)))) / 2.0; end
code[x_, eps_] := N[(N[(N[(N[(1.0 + N[(1.0 / eps), $MachinePrecision]), $MachinePrecision] * N[Exp[(-N[(N[(1.0 - eps), $MachinePrecision] * x), $MachinePrecision])], $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps), $MachinePrecision] - 1.0), $MachinePrecision] * N[Exp[(-N[(N[(1.0 + eps), $MachinePrecision] * x), $MachinePrecision])], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{-\left(1 + \varepsilon\right) \cdot x}}{2}
\end{array}
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x eps) :precision binary64 (/ (- (* (+ 1.0 (/ 1.0 eps)) (exp (- (* (- 1.0 eps) x)))) (* (- (/ 1.0 eps) 1.0) (exp (- (* (+ 1.0 eps) x))))) 2.0))
double code(double x, double eps) {
return (((1.0 + (1.0 / eps)) * exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * exp(-((1.0 + eps) * x)))) / 2.0;
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, eps)
use fmin_fmax_functions
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (((1.0d0 + (1.0d0 / eps)) * exp(-((1.0d0 - eps) * x))) - (((1.0d0 / eps) - 1.0d0) * exp(-((1.0d0 + eps) * x)))) / 2.0d0
end function
public static double code(double x, double eps) {
return (((1.0 + (1.0 / eps)) * Math.exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * Math.exp(-((1.0 + eps) * x)))) / 2.0;
}
def code(x, eps): return (((1.0 + (1.0 / eps)) * math.exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * math.exp(-((1.0 + eps) * x)))) / 2.0
function code(x, eps) return Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps)) * exp(Float64(-Float64(Float64(1.0 - eps) * x)))) - Float64(Float64(Float64(1.0 / eps) - 1.0) * exp(Float64(-Float64(Float64(1.0 + eps) * x))))) / 2.0) end
function tmp = code(x, eps) tmp = (((1.0 + (1.0 / eps)) * exp(-((1.0 - eps) * x))) - (((1.0 / eps) - 1.0) * exp(-((1.0 + eps) * x)))) / 2.0; end
code[x_, eps_] := N[(N[(N[(N[(1.0 + N[(1.0 / eps), $MachinePrecision]), $MachinePrecision] * N[Exp[(-N[(N[(1.0 - eps), $MachinePrecision] * x), $MachinePrecision])], $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps), $MachinePrecision] - 1.0), $MachinePrecision] * N[Exp[(-N[(N[(1.0 + eps), $MachinePrecision] * x), $MachinePrecision])], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{-\left(1 + \varepsilon\right) \cdot x}}{2}
\end{array}
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 (* (- (pow (exp -1.0) (* x (- 1.0 eps_m))) (- (exp (- (fma x eps_m x))))) 0.5))
eps_m = fabs(eps);
double code(double x, double eps_m) {
return (pow(exp(-1.0), (x * (1.0 - eps_m))) - -exp(-fma(x, eps_m, x))) * 0.5;
}
eps_m = abs(eps) function code(x, eps_m) return Float64(Float64((exp(-1.0) ^ Float64(x * Float64(1.0 - eps_m))) - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5) end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := N[(N[(N[Power[N[Exp[-1.0], $MachinePrecision], N[(x * N[(1.0 - eps$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\left({\left(e^{-1}\right)}^{\left(x \cdot \left(1 - eps\_m\right)\right)} - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5
\end{array}
Initial program 73.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.1%
lift-exp.f64N/A
lift-neg.f64N/A
lift--.f64N/A
lift-*.f64N/A
distribute-lft-neg-inN/A
mul-1-negN/A
exp-prodN/A
*-commutativeN/A
lower-pow.f64N/A
lower-exp.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f6499.1
Applied rewrites99.1%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (* (- eps_m 1.0) x)))
(if (<= x 13000000000000.0)
(* (- (exp (* x eps_m)) (- (exp (- (fma x eps_m x))))) 0.5)
(if (<= x 4.7e+116)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (/ (- (* t_0 t_0) 1.0) (- t_0 1.0)))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 (fma x eps_m x) 1.0)))
2.0)
(* (- (exp (* (- x) (- 1.0 eps_m))) -1.0) 0.5)))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (eps_m - 1.0) * x;
double tmp;
if (x <= 13000000000000.0) {
tmp = (exp((x * eps_m)) - -exp(-fma(x, eps_m, x))) * 0.5;
} else if (x <= 4.7e+116) {
tmp = (((1.0 + (1.0 / eps_m)) * (((t_0 * t_0) - 1.0) / (t_0 - 1.0))) - (((1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0;
} else {
tmp = (exp((-x * (1.0 - eps_m))) - -1.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(eps_m - 1.0) * x) tmp = 0.0 if (x <= 13000000000000.0) tmp = Float64(Float64(exp(Float64(x * eps_m)) - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5); elseif (x <= 4.7e+116) tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * Float64(Float64(Float64(t_0 * t_0) - 1.0) / Float64(t_0 - 1.0))) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0); else tmp = Float64(Float64(exp(Float64(Float64(-x) * Float64(1.0 - eps_m))) - -1.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(eps$95$m - 1.0), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[x, 13000000000000.0], N[(N[(N[Exp[N[(x * eps$95$m), $MachinePrecision]], $MachinePrecision] - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 4.7e+116], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t$95$0 * t$95$0), $MachinePrecision] - 1.0), $MachinePrecision] / N[(t$95$0 - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * N[(x * eps$95$m + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(N[Exp[N[((-x) * N[(1.0 - eps$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := \left(eps\_m - 1\right) \cdot x\\
\mathbf{if}\;x \leq 13000000000000:\\
\;\;\;\;\left(e^{x \cdot eps\_m} - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 4.7 \cdot 10^{+116}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \frac{t\_0 \cdot t\_0 - 1}{t\_0 - 1} - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, \mathsf{fma}\left(x, eps\_m, x\right), 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\left(e^{\left(-x\right) \cdot \left(1 - eps\_m\right)} - -1\right) \cdot 0.5\\
\end{array}
\end{array}
if x < 1.3e13Initial program 63.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f6498.5
Applied rewrites98.5%
if 1.3e13 < x < 4.7000000000000003e116Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6439.8
Applied rewrites39.8%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6438.0
Applied rewrites38.0%
lift--.f64N/A
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
*-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lower--.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f6460.6
Applied rewrites60.6%
if 4.7000000000000003e116 < x Initial program 99.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
Applied rewrites52.3%
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 (* (- (exp (* (- x) (- 1.0 eps_m))) (- (exp (- (fma x eps_m x))))) 0.5))
eps_m = fabs(eps);
double code(double x, double eps_m) {
return (exp((-x * (1.0 - eps_m))) - -exp(-fma(x, eps_m, x))) * 0.5;
}
eps_m = abs(eps) function code(x, eps_m) return Float64(Float64(exp(Float64(Float64(-x) * Float64(1.0 - eps_m))) - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5) end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := N[(N[(N[Exp[N[((-x) * N[(1.0 - eps$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\left(e^{\left(-x\right) \cdot \left(1 - eps\_m\right)} - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5
\end{array}
Initial program 73.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.1%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(if (<= x -3e-265)
(* (- 1.0 (- (exp (- (fma x eps_m x))))) 0.5)
(if (<= x 65000.0)
(* (- (exp (* (- x) (- 1.0 eps_m))) (- (fma x eps_m x) 1.0)) 0.5)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (fma (- eps_m 1.0) x 1.0))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 x 1.0)))
2.0))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -3e-265) {
tmp = (1.0 - -exp(-fma(x, eps_m, x))) * 0.5;
} else if (x <= 65000.0) {
tmp = (exp((-x * (1.0 - eps_m))) - (fma(x, eps_m, x) - 1.0)) * 0.5;
} else {
tmp = (((1.0 + (1.0 / eps_m)) * fma((eps_m - 1.0), x, 1.0)) - (((1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -3e-265) tmp = Float64(Float64(1.0 - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5); elseif (x <= 65000.0) tmp = Float64(Float64(exp(Float64(Float64(-x) * Float64(1.0 - eps_m))) - Float64(fma(x, eps_m, x) - 1.0)) * 0.5); else tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * fma(Float64(eps_m - 1.0), x, 1.0)) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -3e-265], N[(N[(1.0 - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 65000.0], N[(N[(N[Exp[N[((-x) * N[(1.0 - eps$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(N[(x * eps$95$m + x), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(eps$95$m - 1.0), $MachinePrecision] * x + 1.0), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3 \cdot 10^{-265}:\\
\;\;\;\;\left(1 - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 65000:\\
\;\;\;\;\left(e^{\left(-x\right) \cdot \left(1 - eps\_m\right)} - \left(\mathsf{fma}\left(x, eps\_m, x\right) - 1\right)\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \mathsf{fma}\left(eps\_m - 1, x, 1\right) - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, x, 1\right)}{2}\\
\end{array}
\end{array}
if x < -2.9999999999999998e-265Initial program 70.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.7%
Taylor expanded in x around 0
*-commutative98.1
distribute-rgt-neg-in98.1
sinh---cosh-rev98.1
Applied rewrites98.1%
if -2.9999999999999998e-265 < x < 65000Initial program 52.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.9%
Taylor expanded in x around 0
lower--.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt1-inN/A
*-commutativeN/A
+-commutativeN/A
lift-fma.f6497.9
Applied rewrites97.9%
if 65000 < x Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6425.5
Applied rewrites25.5%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6423.8
Applied rewrites23.8%
Taylor expanded in eps around 0
Applied rewrites49.7%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (* (- eps_m 1.0) x)) (t_1 (* (- (exp (* x eps_m)) -1.0) 0.5)))
(if (<= x -1.15e-13)
(* (- (exp (- x)) -1.0) 0.5)
(if (<= x -6e-223)
(*
(fma (fma -1.0 (/ (- (* eps_m eps_m) 1.0) (- eps_m 1.0)) -1.0) x 2.0)
0.5)
(if (<= x 13000000000000.0)
t_1
(if (<= x 4.7e+116)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (/ (- (* t_0 t_0) 1.0) (- t_0 1.0)))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 (fma x eps_m x) 1.0)))
2.0)
t_1))))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (eps_m - 1.0) * x;
double t_1 = (exp((x * eps_m)) - -1.0) * 0.5;
double tmp;
if (x <= -1.15e-13) {
tmp = (exp(-x) - -1.0) * 0.5;
} else if (x <= -6e-223) {
tmp = fma(fma(-1.0, (((eps_m * eps_m) - 1.0) / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else if (x <= 13000000000000.0) {
tmp = t_1;
} else if (x <= 4.7e+116) {
tmp = (((1.0 + (1.0 / eps_m)) * (((t_0 * t_0) - 1.0) / (t_0 - 1.0))) - (((1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0;
} else {
tmp = t_1;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(eps_m - 1.0) * x) t_1 = Float64(Float64(exp(Float64(x * eps_m)) - -1.0) * 0.5) tmp = 0.0 if (x <= -1.15e-13) tmp = Float64(Float64(exp(Float64(-x)) - -1.0) * 0.5); elseif (x <= -6e-223) tmp = Float64(fma(fma(-1.0, Float64(Float64(Float64(eps_m * eps_m) - 1.0) / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); elseif (x <= 13000000000000.0) tmp = t_1; elseif (x <= 4.7e+116) tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * Float64(Float64(Float64(t_0 * t_0) - 1.0) / Float64(t_0 - 1.0))) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0); else tmp = t_1; end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(eps$95$m - 1.0), $MachinePrecision] * x), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[Exp[N[(x * eps$95$m), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision]}, If[LessEqual[x, -1.15e-13], N[(N[(N[Exp[(-x)], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, -6e-223], N[(N[(N[(-1.0 * N[(N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 13000000000000.0], t$95$1, If[LessEqual[x, 4.7e+116], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t$95$0 * t$95$0), $MachinePrecision] - 1.0), $MachinePrecision] / N[(t$95$0 - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * N[(x * eps$95$m + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := \left(eps\_m - 1\right) \cdot x\\
t_1 := \left(e^{x \cdot eps\_m} - -1\right) \cdot 0.5\\
\mathbf{if}\;x \leq -1.15 \cdot 10^{-13}:\\
\;\;\;\;\left(e^{-x} - -1\right) \cdot 0.5\\
\mathbf{elif}\;x \leq -6 \cdot 10^{-223}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{eps\_m \cdot eps\_m - 1}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 13000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 4.7 \cdot 10^{+116}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \frac{t\_0 \cdot t\_0 - 1}{t\_0 - 1} - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, \mathsf{fma}\left(x, eps\_m, x\right), 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -1.1499999999999999e-13Initial program 94.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites96.9%
Taylor expanded in x around 0
Applied rewrites5.4%
Taylor expanded in eps around 0
mul-1-negN/A
lift-neg.f6492.7
Applied rewrites92.7%
if -1.1499999999999999e-13 < x < -5.99999999999999983e-223Initial program 53.2%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6470.1
Applied rewrites70.1%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6490.9
Applied rewrites90.9%
Taylor expanded in eps around 0
Applied rewrites90.4%
if -5.99999999999999983e-223 < x < 1.3e13 or 4.7000000000000003e116 < x Initial program 69.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.4%
Taylor expanded in x around 0
Applied rewrites81.7%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f6481.8
Applied rewrites81.8%
if 1.3e13 < x < 4.7000000000000003e116Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6439.8
Applied rewrites39.8%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6438.0
Applied rewrites38.0%
lift--.f64N/A
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
*-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lower--.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f6460.6
Applied rewrites60.6%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (* (- eps_m 1.0) x)))
(if (<= x -2e-266)
(* (- 1.0 (- (exp (- (fma x eps_m x))))) 0.5)
(if (<= x 13000000000000.0)
(* (- (exp (* x eps_m)) -1.0) 0.5)
(if (<= x 4.7e+116)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (/ (- (* t_0 t_0) 1.0) (- t_0 1.0)))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 (fma x eps_m x) 1.0)))
2.0)
(* (- (exp (* (- x) (- 1.0 eps_m))) -1.0) 0.5))))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (eps_m - 1.0) * x;
double tmp;
if (x <= -2e-266) {
tmp = (1.0 - -exp(-fma(x, eps_m, x))) * 0.5;
} else if (x <= 13000000000000.0) {
tmp = (exp((x * eps_m)) - -1.0) * 0.5;
} else if (x <= 4.7e+116) {
tmp = (((1.0 + (1.0 / eps_m)) * (((t_0 * t_0) - 1.0) / (t_0 - 1.0))) - (((1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0;
} else {
tmp = (exp((-x * (1.0 - eps_m))) - -1.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(eps_m - 1.0) * x) tmp = 0.0 if (x <= -2e-266) tmp = Float64(Float64(1.0 - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5); elseif (x <= 13000000000000.0) tmp = Float64(Float64(exp(Float64(x * eps_m)) - -1.0) * 0.5); elseif (x <= 4.7e+116) tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * Float64(Float64(Float64(t_0 * t_0) - 1.0) / Float64(t_0 - 1.0))) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0); else tmp = Float64(Float64(exp(Float64(Float64(-x) * Float64(1.0 - eps_m))) - -1.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(eps$95$m - 1.0), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[x, -2e-266], N[(N[(1.0 - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 13000000000000.0], N[(N[(N[Exp[N[(x * eps$95$m), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 4.7e+116], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t$95$0 * t$95$0), $MachinePrecision] - 1.0), $MachinePrecision] / N[(t$95$0 - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * N[(x * eps$95$m + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(N[Exp[N[((-x) * N[(1.0 - eps$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := \left(eps\_m - 1\right) \cdot x\\
\mathbf{if}\;x \leq -2 \cdot 10^{-266}:\\
\;\;\;\;\left(1 - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 13000000000000:\\
\;\;\;\;\left(e^{x \cdot eps\_m} - -1\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 4.7 \cdot 10^{+116}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \frac{t\_0 \cdot t\_0 - 1}{t\_0 - 1} - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, \mathsf{fma}\left(x, eps\_m, x\right), 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\left(e^{\left(-x\right) \cdot \left(1 - eps\_m\right)} - -1\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -2e-266Initial program 70.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.7%
Taylor expanded in x around 0
*-commutative98.1
distribute-rgt-neg-in98.1
sinh---cosh-rev98.1
Applied rewrites98.1%
if -2e-266 < x < 1.3e13Initial program 53.6%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.9%
Taylor expanded in x around 0
Applied rewrites96.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f6497.3
Applied rewrites97.3%
if 1.3e13 < x < 4.7000000000000003e116Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6439.8
Applied rewrites39.8%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6438.0
Applied rewrites38.0%
lift--.f64N/A
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
*-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lower--.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f6460.6
Applied rewrites60.6%
if 4.7000000000000003e116 < x Initial program 99.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
Applied rewrites52.3%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (* (- (exp (* x eps_m)) -1.0) 0.5)) (t_1 (* (- eps_m 1.0) x)))
(if (<= x -2e-266)
(* (- 1.0 (- (exp (- (fma x eps_m x))))) 0.5)
(if (<= x 13000000000000.0)
t_0
(if (<= x 4.7e+116)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (/ (- (* t_1 t_1) 1.0) (- t_1 1.0)))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 (fma x eps_m x) 1.0)))
2.0)
t_0)))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (exp((x * eps_m)) - -1.0) * 0.5;
double t_1 = (eps_m - 1.0) * x;
double tmp;
if (x <= -2e-266) {
tmp = (1.0 - -exp(-fma(x, eps_m, x))) * 0.5;
} else if (x <= 13000000000000.0) {
tmp = t_0;
} else if (x <= 4.7e+116) {
tmp = (((1.0 + (1.0 / eps_m)) * (((t_1 * t_1) - 1.0) / (t_1 - 1.0))) - (((1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0;
} else {
tmp = t_0;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(exp(Float64(x * eps_m)) - -1.0) * 0.5) t_1 = Float64(Float64(eps_m - 1.0) * x) tmp = 0.0 if (x <= -2e-266) tmp = Float64(Float64(1.0 - Float64(-exp(Float64(-fma(x, eps_m, x))))) * 0.5); elseif (x <= 13000000000000.0) tmp = t_0; elseif (x <= 4.7e+116) tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * Float64(Float64(Float64(t_1 * t_1) - 1.0) / Float64(t_1 - 1.0))) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, fma(x, eps_m, x), 1.0))) / 2.0); else tmp = t_0; end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(N[Exp[N[(x * eps$95$m), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision]}, Block[{t$95$1 = N[(N[(eps$95$m - 1.0), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[x, -2e-266], N[(N[(1.0 - (-N[Exp[(-N[(x * eps$95$m + x), $MachinePrecision])], $MachinePrecision])), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 13000000000000.0], t$95$0, If[LessEqual[x, 4.7e+116], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t$95$1 * t$95$1), $MachinePrecision] - 1.0), $MachinePrecision] / N[(t$95$1 - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * N[(x * eps$95$m + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := \left(e^{x \cdot eps\_m} - -1\right) \cdot 0.5\\
t_1 := \left(eps\_m - 1\right) \cdot x\\
\mathbf{if}\;x \leq -2 \cdot 10^{-266}:\\
\;\;\;\;\left(1 - \left(-e^{-\mathsf{fma}\left(x, eps\_m, x\right)}\right)\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 13000000000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 4.7 \cdot 10^{+116}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \frac{t\_1 \cdot t\_1 - 1}{t\_1 - 1} - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, \mathsf{fma}\left(x, eps\_m, x\right), 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x < -2e-266Initial program 70.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.7%
Taylor expanded in x around 0
*-commutative98.1
distribute-rgt-neg-in98.1
sinh---cosh-rev98.1
Applied rewrites98.1%
if -2e-266 < x < 1.3e13 or 4.7000000000000003e116 < x Initial program 69.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.3%
Taylor expanded in x around 0
Applied rewrites81.3%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f6481.4
Applied rewrites81.4%
if 1.3e13 < x < 4.7000000000000003e116Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6439.8
Applied rewrites39.8%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6438.0
Applied rewrites38.0%
lift--.f64N/A
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
*-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lower--.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f64N/A
*-commutativeN/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift--.f6460.6
Applied rewrites60.6%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(if (<= x -1.15e-13)
(* (- (exp (- x)) -1.0) 0.5)
(if (<= x -5.6e-213)
(*
(fma (fma -1.0 (/ (- (* eps_m eps_m) 1.0) (- eps_m 1.0)) -1.0) x 2.0)
0.5)
(if (<= x 7.2e-268)
1.0
(if (<= x 60000.0)
(*
(fma
(fma
-1.0
(+ eps_m 1.0)
(- (/ (- 1.0 (* eps_m eps_m)) (+ eps_m 1.0))))
x
2.0)
0.5)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (fma (- eps_m 1.0) x 1.0))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 x 1.0)))
2.0))))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -1.15e-13) {
tmp = (exp(-x) - -1.0) * 0.5;
} else if (x <= -5.6e-213) {
tmp = fma(fma(-1.0, (((eps_m * eps_m) - 1.0) / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else if (x <= 7.2e-268) {
tmp = 1.0;
} else if (x <= 60000.0) {
tmp = fma(fma(-1.0, (eps_m + 1.0), -((1.0 - (eps_m * eps_m)) / (eps_m + 1.0))), x, 2.0) * 0.5;
} else {
tmp = (((1.0 + (1.0 / eps_m)) * fma((eps_m - 1.0), x, 1.0)) - (((1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -1.15e-13) tmp = Float64(Float64(exp(Float64(-x)) - -1.0) * 0.5); elseif (x <= -5.6e-213) tmp = Float64(fma(fma(-1.0, Float64(Float64(Float64(eps_m * eps_m) - 1.0) / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); elseif (x <= 7.2e-268) tmp = 1.0; elseif (x <= 60000.0) tmp = Float64(fma(fma(-1.0, Float64(eps_m + 1.0), Float64(-Float64(Float64(1.0 - Float64(eps_m * eps_m)) / Float64(eps_m + 1.0)))), x, 2.0) * 0.5); else tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * fma(Float64(eps_m - 1.0), x, 1.0)) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -1.15e-13], N[(N[(N[Exp[(-x)], $MachinePrecision] - -1.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, -5.6e-213], N[(N[(N[(-1.0 * N[(N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 7.2e-268], 1.0, If[LessEqual[x, 60000.0], N[(N[(N[(-1.0 * N[(eps$95$m + 1.0), $MachinePrecision] + (-N[(N[(1.0 - N[(eps$95$m * eps$95$m), $MachinePrecision]), $MachinePrecision] / N[(eps$95$m + 1.0), $MachinePrecision]), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(eps$95$m - 1.0), $MachinePrecision] * x + 1.0), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.15 \cdot 10^{-13}:\\
\;\;\;\;\left(e^{-x} - -1\right) \cdot 0.5\\
\mathbf{elif}\;x \leq -5.6 \cdot 10^{-213}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{eps\_m \cdot eps\_m - 1}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{-268}:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 60000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, eps\_m + 1, -\frac{1 - eps\_m \cdot eps\_m}{eps\_m + 1}\right), x, 2\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \mathsf{fma}\left(eps\_m - 1, x, 1\right) - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, x, 1\right)}{2}\\
\end{array}
\end{array}
if x < -1.1499999999999999e-13Initial program 94.8%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites96.9%
Taylor expanded in x around 0
Applied rewrites5.4%
Taylor expanded in eps around 0
mul-1-negN/A
lift-neg.f6492.7
Applied rewrites92.7%
if -1.1499999999999999e-13 < x < -5.6e-213Initial program 53.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6469.8
Applied rewrites69.8%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6490.7
Applied rewrites90.7%
Taylor expanded in eps around 0
Applied rewrites90.2%
if -5.6e-213 < x < 7.2000000000000002e-268Initial program 54.7%
Taylor expanded in x around 0
Applied rewrites93.1%
if 7.2000000000000002e-268 < x < 6e4Initial program 52.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.6%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6470.7
Applied rewrites70.7%
lift--.f64N/A
flip--N/A
lower-/.f64N/A
metadata-evalN/A
unpow2N/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lift-+.f6488.0
Applied rewrites88.0%
if 6e4 < x Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6425.5
Applied rewrites25.5%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6423.9
Applied rewrites23.9%
Taylor expanded in eps around 0
Applied rewrites49.7%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(if (<= x -3.8e+184)
(* (/ (- (* (* -2.0 x) (* -2.0 x)) 4.0) (- (* -2.0 x) 2.0)) 0.5)
(if (<= x -5.6e-213)
(*
(fma (fma -1.0 (/ (- (* eps_m eps_m) 1.0) (- eps_m 1.0)) -1.0) x 2.0)
0.5)
(if (<= x 7.2e-268)
1.0
(if (<= x 60000.0)
(*
(fma
(fma
-1.0
(+ eps_m 1.0)
(- (/ (- 1.0 (* eps_m eps_m)) (+ eps_m 1.0))))
x
2.0)
0.5)
(/
(-
(* (+ 1.0 (/ 1.0 eps_m)) (fma (- eps_m 1.0) x 1.0))
(* (- (/ 1.0 eps_m) 1.0) (fma -1.0 x 1.0)))
2.0))))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -3.8e+184) {
tmp = ((((-2.0 * x) * (-2.0 * x)) - 4.0) / ((-2.0 * x) - 2.0)) * 0.5;
} else if (x <= -5.6e-213) {
tmp = fma(fma(-1.0, (((eps_m * eps_m) - 1.0) / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else if (x <= 7.2e-268) {
tmp = 1.0;
} else if (x <= 60000.0) {
tmp = fma(fma(-1.0, (eps_m + 1.0), -((1.0 - (eps_m * eps_m)) / (eps_m + 1.0))), x, 2.0) * 0.5;
} else {
tmp = (((1.0 + (1.0 / eps_m)) * fma((eps_m - 1.0), x, 1.0)) - (((1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -3.8e+184) tmp = Float64(Float64(Float64(Float64(Float64(-2.0 * x) * Float64(-2.0 * x)) - 4.0) / Float64(Float64(-2.0 * x) - 2.0)) * 0.5); elseif (x <= -5.6e-213) tmp = Float64(fma(fma(-1.0, Float64(Float64(Float64(eps_m * eps_m) - 1.0) / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); elseif (x <= 7.2e-268) tmp = 1.0; elseif (x <= 60000.0) tmp = Float64(fma(fma(-1.0, Float64(eps_m + 1.0), Float64(-Float64(Float64(1.0 - Float64(eps_m * eps_m)) / Float64(eps_m + 1.0)))), x, 2.0) * 0.5); else tmp = Float64(Float64(Float64(Float64(1.0 + Float64(1.0 / eps_m)) * fma(Float64(eps_m - 1.0), x, 1.0)) - Float64(Float64(Float64(1.0 / eps_m) - 1.0) * fma(-1.0, x, 1.0))) / 2.0); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -3.8e+184], N[(N[(N[(N[(N[(-2.0 * x), $MachinePrecision] * N[(-2.0 * x), $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(-2.0 * x), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, -5.6e-213], N[(N[(N[(-1.0 * N[(N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 7.2e-268], 1.0, If[LessEqual[x, 60000.0], N[(N[(N[(-1.0 * N[(eps$95$m + 1.0), $MachinePrecision] + (-N[(N[(1.0 - N[(eps$95$m * eps$95$m), $MachinePrecision]), $MachinePrecision] / N[(eps$95$m + 1.0), $MachinePrecision]), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(N[(1.0 + N[(1.0 / eps$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(eps$95$m - 1.0), $MachinePrecision] * x + 1.0), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.0 / eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] * N[(-1.0 * x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.8 \cdot 10^{+184}:\\
\;\;\;\;\frac{\left(-2 \cdot x\right) \cdot \left(-2 \cdot x\right) - 4}{-2 \cdot x - 2} \cdot 0.5\\
\mathbf{elif}\;x \leq -5.6 \cdot 10^{-213}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{eps\_m \cdot eps\_m - 1}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{-268}:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 60000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, eps\_m + 1, -\frac{1 - eps\_m \cdot eps\_m}{eps\_m + 1}\right), x, 2\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(1 + \frac{1}{eps\_m}\right) \cdot \mathsf{fma}\left(eps\_m - 1, x, 1\right) - \left(\frac{1}{eps\_m} - 1\right) \cdot \mathsf{fma}\left(-1, x, 1\right)}{2}\\
\end{array}
\end{array}
if x < -3.8000000000000001e184Initial program 100.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f643.3
Applied rewrites3.3%
Taylor expanded in eps around 0
Applied rewrites7.6%
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
Applied rewrites99.8%
if -3.8000000000000001e184 < x < -5.6e-213Initial program 67.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6447.1
Applied rewrites47.1%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6476.4
Applied rewrites76.4%
Taylor expanded in eps around 0
Applied rewrites76.5%
if -5.6e-213 < x < 7.2000000000000002e-268Initial program 54.7%
Taylor expanded in x around 0
Applied rewrites93.1%
if 7.2000000000000002e-268 < x < 6e4Initial program 52.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.6%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6470.7
Applied rewrites70.7%
lift--.f64N/A
flip--N/A
lower-/.f64N/A
metadata-evalN/A
unpow2N/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lift-+.f6488.0
Applied rewrites88.0%
if 6e4 < x Initial program 100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6425.5
Applied rewrites25.5%
Taylor expanded in x around 0
sinh---cosh-revN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6423.9
Applied rewrites23.9%
Taylor expanded in eps around 0
Applied rewrites49.7%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (- (* eps_m eps_m) 1.0)))
(if (<= x -3.8e+184)
(* (/ (- (* (* -2.0 x) (* -2.0 x)) 4.0) (- (* -2.0 x) 2.0)) 0.5)
(if (<= x -5.6e-213)
(* (fma (fma -1.0 (/ t_0 (- eps_m 1.0)) -1.0) x 2.0) 0.5)
(if (<= x 7.2e-268)
1.0
(* (fma (fma -1.0 (/ t_0 -1.0) (- (- 1.0 eps_m))) x 2.0) 0.5))))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (eps_m * eps_m) - 1.0;
double tmp;
if (x <= -3.8e+184) {
tmp = ((((-2.0 * x) * (-2.0 * x)) - 4.0) / ((-2.0 * x) - 2.0)) * 0.5;
} else if (x <= -5.6e-213) {
tmp = fma(fma(-1.0, (t_0 / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else if (x <= 7.2e-268) {
tmp = 1.0;
} else {
tmp = fma(fma(-1.0, (t_0 / -1.0), -(1.0 - eps_m)), x, 2.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(eps_m * eps_m) - 1.0) tmp = 0.0 if (x <= -3.8e+184) tmp = Float64(Float64(Float64(Float64(Float64(-2.0 * x) * Float64(-2.0 * x)) - 4.0) / Float64(Float64(-2.0 * x) - 2.0)) * 0.5); elseif (x <= -5.6e-213) tmp = Float64(fma(fma(-1.0, Float64(t_0 / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); elseif (x <= 7.2e-268) tmp = 1.0; else tmp = Float64(fma(fma(-1.0, Float64(t_0 / -1.0), Float64(-Float64(1.0 - eps_m))), x, 2.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision]}, If[LessEqual[x, -3.8e+184], N[(N[(N[(N[(N[(-2.0 * x), $MachinePrecision] * N[(-2.0 * x), $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] / N[(N[(-2.0 * x), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, -5.6e-213], N[(N[(N[(-1.0 * N[(t$95$0 / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 7.2e-268], 1.0, N[(N[(N[(-1.0 * N[(t$95$0 / -1.0), $MachinePrecision] + (-N[(1.0 - eps$95$m), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := eps\_m \cdot eps\_m - 1\\
\mathbf{if}\;x \leq -3.8 \cdot 10^{+184}:\\
\;\;\;\;\frac{\left(-2 \cdot x\right) \cdot \left(-2 \cdot x\right) - 4}{-2 \cdot x - 2} \cdot 0.5\\
\mathbf{elif}\;x \leq -5.6 \cdot 10^{-213}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{t\_0}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{-268}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{t\_0}{-1}, -\left(1 - eps\_m\right)\right), x, 2\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -3.8000000000000001e184Initial program 100.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f643.3
Applied rewrites3.3%
Taylor expanded in eps around 0
Applied rewrites7.6%
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
Applied rewrites99.8%
if -3.8000000000000001e184 < x < -5.6e-213Initial program 67.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6447.1
Applied rewrites47.1%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6476.4
Applied rewrites76.4%
Taylor expanded in eps around 0
Applied rewrites76.5%
if -5.6e-213 < x < 7.2000000000000002e-268Initial program 54.7%
Taylor expanded in x around 0
Applied rewrites93.1%
if 7.2000000000000002e-268 < x Initial program 77.7%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6434.8
Applied rewrites34.8%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6432.6
Applied rewrites32.6%
Taylor expanded in eps around 0
Applied rewrites60.1%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(let* ((t_0 (- (* eps_m eps_m) 1.0)))
(if (<= x -5.6e-213)
(* (fma (fma -1.0 (/ t_0 (- eps_m 1.0)) -1.0) x 2.0) 0.5)
(if (<= x 7.2e-268)
1.0
(* (fma (fma -1.0 (/ t_0 -1.0) (- (- 1.0 eps_m))) x 2.0) 0.5)))))eps_m = fabs(eps);
double code(double x, double eps_m) {
double t_0 = (eps_m * eps_m) - 1.0;
double tmp;
if (x <= -5.6e-213) {
tmp = fma(fma(-1.0, (t_0 / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else if (x <= 7.2e-268) {
tmp = 1.0;
} else {
tmp = fma(fma(-1.0, (t_0 / -1.0), -(1.0 - eps_m)), x, 2.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) t_0 = Float64(Float64(eps_m * eps_m) - 1.0) tmp = 0.0 if (x <= -5.6e-213) tmp = Float64(fma(fma(-1.0, Float64(t_0 / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); elseif (x <= 7.2e-268) tmp = 1.0; else tmp = Float64(fma(fma(-1.0, Float64(t_0 / -1.0), Float64(-Float64(1.0 - eps_m))), x, 2.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision]
code[x_, eps$95$m_] := Block[{t$95$0 = N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision]}, If[LessEqual[x, -5.6e-213], N[(N[(N[(-1.0 * N[(t$95$0 / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 7.2e-268], 1.0, N[(N[(N[(-1.0 * N[(t$95$0 / -1.0), $MachinePrecision] + (-N[(1.0 - eps$95$m), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]]]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
t_0 := eps\_m \cdot eps\_m - 1\\
\mathbf{if}\;x \leq -5.6 \cdot 10^{-213}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{t\_0}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{-268}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{t\_0}{-1}, -\left(1 - eps\_m\right)\right), x, 2\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -5.6e-213Initial program 72.4%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.5%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6439.9
Applied rewrites39.9%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6472.9
Applied rewrites72.9%
Taylor expanded in eps around 0
Applied rewrites77.7%
if -5.6e-213 < x < 7.2000000000000002e-268Initial program 54.7%
Taylor expanded in x around 0
Applied rewrites93.1%
if 7.2000000000000002e-268 < x Initial program 77.7%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6434.8
Applied rewrites34.8%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6432.6
Applied rewrites32.6%
Taylor expanded in eps around 0
Applied rewrites60.1%
eps_m = (fabs.f64 eps)
(FPCore (x eps_m)
:precision binary64
(if (<= x -6e-223)
(*
(fma (fma -1.0 (/ (- (* eps_m eps_m) 1.0) (- eps_m 1.0)) -1.0) x 2.0)
0.5)
(* (fma (fma -1.0 (/ -1.0 (- eps_m 1.0)) (- (- 1.0 eps_m))) x 2.0) 0.5)))eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -6e-223) {
tmp = fma(fma(-1.0, (((eps_m * eps_m) - 1.0) / (eps_m - 1.0)), -1.0), x, 2.0) * 0.5;
} else {
tmp = fma(fma(-1.0, (-1.0 / (eps_m - 1.0)), -(1.0 - eps_m)), x, 2.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -6e-223) tmp = Float64(fma(fma(-1.0, Float64(Float64(Float64(eps_m * eps_m) - 1.0) / Float64(eps_m - 1.0)), -1.0), x, 2.0) * 0.5); else tmp = Float64(fma(fma(-1.0, Float64(-1.0 / Float64(eps_m - 1.0)), Float64(-Float64(1.0 - eps_m))), x, 2.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -6e-223], N[(N[(N[(-1.0 * N[(N[(N[(eps$95$m * eps$95$m), $MachinePrecision] - 1.0), $MachinePrecision] / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(-1.0 * N[(-1.0 / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + (-N[(1.0 - eps$95$m), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6 \cdot 10^{-223}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{eps\_m \cdot eps\_m - 1}{eps\_m - 1}, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{-1}{eps\_m - 1}, -\left(1 - eps\_m\right)\right), x, 2\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -5.99999999999999983e-223Initial program 72.1%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.6%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6440.7
Applied rewrites40.7%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6473.3
Applied rewrites73.3%
Taylor expanded in eps around 0
Applied rewrites78.0%
if -5.99999999999999983e-223 < x Initial program 73.5%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.4%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6445.5
Applied rewrites45.5%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6441.4
Applied rewrites41.4%
Taylor expanded in eps around 0
Applied rewrites56.7%
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 (if (<= x -2e-265) (* (fma (fma -1.0 (+ eps_m 1.0) -1.0) x 2.0) 0.5) (* (fma (fma -1.0 (/ -1.0 (- eps_m 1.0)) (- (- 1.0 eps_m))) x 2.0) 0.5)))
eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -2e-265) {
tmp = fma(fma(-1.0, (eps_m + 1.0), -1.0), x, 2.0) * 0.5;
} else {
tmp = fma(fma(-1.0, (-1.0 / (eps_m - 1.0)), -(1.0 - eps_m)), x, 2.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -2e-265) tmp = Float64(fma(fma(-1.0, Float64(eps_m + 1.0), -1.0), x, 2.0) * 0.5); else tmp = Float64(fma(fma(-1.0, Float64(-1.0 / Float64(eps_m - 1.0)), Float64(-Float64(1.0 - eps_m))), x, 2.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -2e-265], N[(N[(N[(-1.0 * N[(eps$95$m + 1.0), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(-1.0 * N[(-1.0 / N[(eps$95$m - 1.0), $MachinePrecision]), $MachinePrecision] + (-N[(1.0 - eps$95$m), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2 \cdot 10^{-265}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, eps\_m + 1, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, \frac{-1}{eps\_m - 1}, -\left(1 - eps\_m\right)\right), x, 2\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -1.99999999999999997e-265Initial program 70.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.7%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6445.1
Applied rewrites45.1%
Taylor expanded in eps around 0
Applied rewrites64.0%
if -1.99999999999999997e-265 < x Initial program 74.4%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.4%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6442.9
Applied rewrites42.9%
lift-+.f64N/A
flip-+N/A
lower-/.f64N/A
unpow2N/A
metadata-evalN/A
lower--.f64N/A
unpow2N/A
lower-*.f64N/A
lift--.f6438.7
Applied rewrites38.7%
Taylor expanded in eps around 0
Applied rewrites54.9%
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 (if (<= x -2e-265) (* (fma (fma -1.0 (+ eps_m 1.0) -1.0) x 2.0) 0.5) (* (fma (fma -1.0 1.0 (- (- 1.0 eps_m))) x 2.0) 0.5)))
eps_m = fabs(eps);
double code(double x, double eps_m) {
double tmp;
if (x <= -2e-265) {
tmp = fma(fma(-1.0, (eps_m + 1.0), -1.0), x, 2.0) * 0.5;
} else {
tmp = fma(fma(-1.0, 1.0, -(1.0 - eps_m)), x, 2.0) * 0.5;
}
return tmp;
}
eps_m = abs(eps) function code(x, eps_m) tmp = 0.0 if (x <= -2e-265) tmp = Float64(fma(fma(-1.0, Float64(eps_m + 1.0), -1.0), x, 2.0) * 0.5); else tmp = Float64(fma(fma(-1.0, 1.0, Float64(-Float64(1.0 - eps_m))), x, 2.0) * 0.5); end return tmp end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := If[LessEqual[x, -2e-265], N[(N[(N[(-1.0 * N[(eps$95$m + 1.0), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(N[(-1.0 * 1.0 + (-N[(1.0 - eps$95$m), $MachinePrecision])), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2 \cdot 10^{-265}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, eps\_m + 1, -1\right), x, 2\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-1, 1, -\left(1 - eps\_m\right)\right), x, 2\right) \cdot 0.5\\
\end{array}
\end{array}
if x < -1.99999999999999997e-265Initial program 70.9%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.7%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6445.1
Applied rewrites45.1%
Taylor expanded in eps around 0
Applied rewrites64.0%
if -1.99999999999999997e-265 < x Initial program 74.4%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.4%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6442.9
Applied rewrites42.9%
Taylor expanded in eps around 0
Applied rewrites54.9%
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 (* (fma (fma -1.0 (+ eps_m 1.0) -1.0) x 2.0) 0.5))
eps_m = fabs(eps);
double code(double x, double eps_m) {
return fma(fma(-1.0, (eps_m + 1.0), -1.0), x, 2.0) * 0.5;
}
eps_m = abs(eps) function code(x, eps_m) return Float64(fma(fma(-1.0, Float64(eps_m + 1.0), -1.0), x, 2.0) * 0.5) end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := N[(N[(N[(-1.0 * N[(eps$95$m + 1.0), $MachinePrecision] + -1.0), $MachinePrecision] * x + 2.0), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
\mathsf{fma}\left(\mathsf{fma}\left(-1, eps\_m + 1, -1\right), x, 2\right) \cdot 0.5
\end{array}
Initial program 73.0%
Taylor expanded in eps around inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.1%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift--.f6443.8
Applied rewrites43.8%
Taylor expanded in eps around 0
Applied rewrites50.5%
eps_m = (fabs.f64 eps) (FPCore (x eps_m) :precision binary64 1.0)
eps_m = fabs(eps);
double code(double x, double eps_m) {
return 1.0;
}
eps_m = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, eps_m)
use fmin_fmax_functions
real(8), intent (in) :: x
real(8), intent (in) :: eps_m
code = 1.0d0
end function
eps_m = Math.abs(eps);
public static double code(double x, double eps_m) {
return 1.0;
}
eps_m = math.fabs(eps) def code(x, eps_m): return 1.0
eps_m = abs(eps) function code(x, eps_m) return 1.0 end
eps_m = abs(eps); function tmp = code(x, eps_m) tmp = 1.0; end
eps_m = N[Abs[eps], $MachinePrecision] code[x_, eps$95$m_] := 1.0
\begin{array}{l}
eps_m = \left|\varepsilon\right|
\\
1
\end{array}
Initial program 73.0%
Taylor expanded in x around 0
Applied rewrites44.2%
herbie shell --seed 2025089
(FPCore (x eps)
:name "NMSE Section 6.1 mentioned, A"
:precision binary64
(/ (- (* (+ 1.0 (/ 1.0 eps)) (exp (- (* (- 1.0 eps) x)))) (* (- (/ 1.0 eps) 1.0) (exp (- (* (+ 1.0 eps) x))))) 2.0))