
(FPCore (x y z) :precision binary64 (exp (- (+ x (* y (log y))) z)))
double code(double x, double y, double z) {
return exp(((x + (y * log(y))) - z));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = exp(((x + (y * log(y))) - z))
end function
public static double code(double x, double y, double z) {
return Math.exp(((x + (y * Math.log(y))) - z));
}
def code(x, y, z): return math.exp(((x + (y * math.log(y))) - z))
function code(x, y, z) return exp(Float64(Float64(x + Float64(y * log(y))) - z)) end
function tmp = code(x, y, z) tmp = exp(((x + (y * log(y))) - z)); end
code[x_, y_, z_] := N[Exp[N[(N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - z), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x + y \cdot \log y\right) - z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (exp (- (+ x (* y (log y))) z)))
double code(double x, double y, double z) {
return exp(((x + (y * log(y))) - z));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = exp(((x + (y * log(y))) - z))
end function
public static double code(double x, double y, double z) {
return Math.exp(((x + (y * Math.log(y))) - z));
}
def code(x, y, z): return math.exp(((x + (y * math.log(y))) - z))
function code(x, y, z) return exp(Float64(Float64(x + Float64(y * log(y))) - z)) end
function tmp = code(x, y, z) tmp = exp(((x + (y * log(y))) - z)); end
code[x_, y_, z_] := N[Exp[N[(N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - z), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x + y \cdot \log y\right) - z}
\end{array}
(FPCore (x y z) :precision binary64 (exp (- (+ x (* y (log y))) z)))
double code(double x, double y, double z) {
return exp(((x + (y * log(y))) - z));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = exp(((x + (y * log(y))) - z))
end function
public static double code(double x, double y, double z) {
return Math.exp(((x + (y * Math.log(y))) - z));
}
def code(x, y, z): return math.exp(((x + (y * math.log(y))) - z))
function code(x, y, z) return exp(Float64(Float64(x + Float64(y * log(y))) - z)) end
function tmp = code(x, y, z) tmp = exp(((x + (y * log(y))) - z)); end
code[x_, y_, z_] := N[Exp[N[(N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - z), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x + y \cdot \log y\right) - z}
\end{array}
Initial program 100.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (* y (log y)))))
(if (<= t_0 -2e-10)
(exp x)
(if (<= t_0 -1e-122)
(/ 1.0 (fma z (fma z (fma -0.25 (* z z) 0.5) 1.0) 1.0))
(if (<= t_0 5e-39)
(fma
(fma z 0.5 -1.0)
(* z (* (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0)))))
1.0)
(pow y y))))))
double code(double x, double y, double z) {
double t_0 = x + (y * log(y));
double tmp;
if (t_0 <= -2e-10) {
tmp = exp(x);
} else if (t_0 <= -1e-122) {
tmp = 1.0 / fma(z, fma(z, fma(-0.25, (z * z), 0.5), 1.0), 1.0);
} else if (t_0 <= 5e-39) {
tmp = fma(fma(z, 0.5, -1.0), (z * (fma(z, 0.5, -1.0) * (z * (z * fma(z, 0.5, -1.0))))), 1.0);
} else {
tmp = pow(y, y);
}
return tmp;
}
function code(x, y, z) t_0 = Float64(x + Float64(y * log(y))) tmp = 0.0 if (t_0 <= -2e-10) tmp = exp(x); elseif (t_0 <= -1e-122) tmp = Float64(1.0 / fma(z, fma(z, fma(-0.25, Float64(z * z), 0.5), 1.0), 1.0)); elseif (t_0 <= 5e-39) tmp = fma(fma(z, 0.5, -1.0), Float64(z * Float64(fma(z, 0.5, -1.0) * Float64(z * Float64(z * fma(z, 0.5, -1.0))))), 1.0); else tmp = y ^ y; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-10], N[Exp[x], $MachinePrecision], If[LessEqual[t$95$0, -1e-122], N[(1.0 / N[(z * N[(z * N[(-0.25 * N[(z * z), $MachinePrecision] + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 5e-39], N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[Power[y, y], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + y \cdot \log y\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-10}:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;t\_0 \leq -1 \cdot 10^{-122}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(-0.25, z \cdot z, 0.5\right), 1\right), 1\right)}\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{-39}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(\mathsf{fma}\left(z, 0.5, -1\right) \cdot \left(z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right)\right)\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;{y}^{y}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 y (log.f64 y))) < -2.00000000000000007e-10Initial program 100.0%
Taylor expanded in y around 0
lower--.f6497.4
Simplified97.4%
Taylor expanded in z around 0
lower-exp.f6487.8
Simplified87.8%
if -2.00000000000000007e-10 < (+.f64 x (*.f64 y (log.f64 y))) < -1.00000000000000006e-122Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6499.0
Simplified99.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6454.3
Simplified54.3%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6454.3
Applied egg-rr54.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6472.4
Simplified72.4%
if -1.00000000000000006e-122 < (+.f64 x (*.f64 y (log.f64 y))) < 4.9999999999999998e-39Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f64100.0
Simplified100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6480.3
Simplified80.3%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr58.9%
Taylor expanded in z around 0
Simplified86.8%
if 4.9999999999999998e-39 < (+.f64 x (*.f64 y (log.f64 y))) Initial program 100.0%
Taylor expanded in y around inf
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
lower-*.f64N/A
lower-log.f6474.5
Simplified74.5%
lift-log.f64N/A
*-commutativeN/A
lift-log.f64N/A
exp-to-powN/A
lower-pow.f6474.5
Applied egg-rr74.5%
Final simplification78.3%
(FPCore (x y z) :precision binary64 (let* ((t_0 (+ x (* y (log y))))) (if (<= t_0 -4e+23) (exp x) (if (<= t_0 1e+29) (exp (- z)) (pow y y)))))
double code(double x, double y, double z) {
double t_0 = x + (y * log(y));
double tmp;
if (t_0 <= -4e+23) {
tmp = exp(x);
} else if (t_0 <= 1e+29) {
tmp = exp(-z);
} else {
tmp = pow(y, y);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = x + (y * log(y))
if (t_0 <= (-4d+23)) then
tmp = exp(x)
else if (t_0 <= 1d+29) then
tmp = exp(-z)
else
tmp = y ** y
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = x + (y * Math.log(y));
double tmp;
if (t_0 <= -4e+23) {
tmp = Math.exp(x);
} else if (t_0 <= 1e+29) {
tmp = Math.exp(-z);
} else {
tmp = Math.pow(y, y);
}
return tmp;
}
def code(x, y, z): t_0 = x + (y * math.log(y)) tmp = 0 if t_0 <= -4e+23: tmp = math.exp(x) elif t_0 <= 1e+29: tmp = math.exp(-z) else: tmp = math.pow(y, y) return tmp
function code(x, y, z) t_0 = Float64(x + Float64(y * log(y))) tmp = 0.0 if (t_0 <= -4e+23) tmp = exp(x); elseif (t_0 <= 1e+29) tmp = exp(Float64(-z)); else tmp = y ^ y; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x + (y * log(y)); tmp = 0.0; if (t_0 <= -4e+23) tmp = exp(x); elseif (t_0 <= 1e+29) tmp = exp(-z); else tmp = y ^ y; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -4e+23], N[Exp[x], $MachinePrecision], If[LessEqual[t$95$0, 1e+29], N[Exp[(-z)], $MachinePrecision], N[Power[y, y], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + y \cdot \log y\\
\mathbf{if}\;t\_0 \leq -4 \cdot 10^{+23}:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;t\_0 \leq 10^{+29}:\\
\;\;\;\;e^{-z}\\
\mathbf{else}:\\
\;\;\;\;{y}^{y}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 y (log.f64 y))) < -3.9999999999999997e23Initial program 100.0%
Taylor expanded in y around 0
lower--.f64100.0
Simplified100.0%
Taylor expanded in z around 0
lower-exp.f6493.3
Simplified93.3%
if -3.9999999999999997e23 < (+.f64 x (*.f64 y (log.f64 y))) < 9.99999999999999914e28Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6492.2
Simplified92.2%
if 9.99999999999999914e28 < (+.f64 x (*.f64 y (log.f64 y))) Initial program 100.0%
Taylor expanded in y around inf
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
lower-*.f64N/A
lower-log.f6476.6
Simplified76.6%
lift-log.f64N/A
*-commutativeN/A
lift-log.f64N/A
exp-to-powN/A
lower-pow.f6476.6
Applied egg-rr76.6%
(FPCore (x y z) :precision binary64 (if (<= (* y (log y)) 20.0) (exp (- x z)) (exp (fma y (log y) x))))
double code(double x, double y, double z) {
double tmp;
if ((y * log(y)) <= 20.0) {
tmp = exp((x - z));
} else {
tmp = exp(fma(y, log(y), x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (Float64(y * log(y)) <= 20.0) tmp = exp(Float64(x - z)); else tmp = exp(fma(y, log(y), x)); end return tmp end
code[x_, y_, z_] := If[LessEqual[N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision], 20.0], N[Exp[N[(x - z), $MachinePrecision]], $MachinePrecision], N[Exp[N[(y * N[Log[y], $MachinePrecision] + x), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \cdot \log y \leq 20:\\
\;\;\;\;e^{x - z}\\
\mathbf{else}:\\
\;\;\;\;e^{\mathsf{fma}\left(y, \log y, x\right)}\\
\end{array}
\end{array}
if (*.f64 y (log.f64 y)) < 20Initial program 100.0%
Taylor expanded in y around 0
lower--.f6498.7
Simplified98.7%
if 20 < (*.f64 y (log.f64 y)) Initial program 100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower-log.f6491.4
Simplified91.4%
(FPCore (x y z) :precision binary64 (let* ((t_0 (- (+ x (* y (log y))) z)) (t_1 (* (* z z) 0.5))) (if (<= t_0 -2e+127) t_1 (if (<= t_0 2e+85) (+ x 1.0) t_1))))
double code(double x, double y, double z) {
double t_0 = (x + (y * log(y))) - z;
double t_1 = (z * z) * 0.5;
double tmp;
if (t_0 <= -2e+127) {
tmp = t_1;
} else if (t_0 <= 2e+85) {
tmp = x + 1.0;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x + (y * log(y))) - z
t_1 = (z * z) * 0.5d0
if (t_0 <= (-2d+127)) then
tmp = t_1
else if (t_0 <= 2d+85) then
tmp = x + 1.0d0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x + (y * Math.log(y))) - z;
double t_1 = (z * z) * 0.5;
double tmp;
if (t_0 <= -2e+127) {
tmp = t_1;
} else if (t_0 <= 2e+85) {
tmp = x + 1.0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = (x + (y * math.log(y))) - z t_1 = (z * z) * 0.5 tmp = 0 if t_0 <= -2e+127: tmp = t_1 elif t_0 <= 2e+85: tmp = x + 1.0 else: tmp = t_1 return tmp
function code(x, y, z) t_0 = Float64(Float64(x + Float64(y * log(y))) - z) t_1 = Float64(Float64(z * z) * 0.5) tmp = 0.0 if (t_0 <= -2e+127) tmp = t_1; elseif (t_0 <= 2e+85) tmp = Float64(x + 1.0); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + (y * log(y))) - z; t_1 = (z * z) * 0.5; tmp = 0.0; if (t_0 <= -2e+127) tmp = t_1; elseif (t_0 <= 2e+85) tmp = x + 1.0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision]}, If[LessEqual[t$95$0, -2e+127], t$95$1, If[LessEqual[t$95$0, 2e+85], N[(x + 1.0), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + y \cdot \log y\right) - z\\
t_1 := \left(z \cdot z\right) \cdot 0.5\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{+127}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+85}:\\
\;\;\;\;x + 1\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) < -1.99999999999999991e127 or 2e85 < (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6444.8
Simplified44.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6420.1
Simplified20.1%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6425.2
Simplified25.2%
if -1.99999999999999991e127 < (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) < 2e85Initial program 100.0%
Taylor expanded in y around 0
lower--.f6486.5
Simplified86.5%
Taylor expanded in z around 0
lower-exp.f6466.6
Simplified66.6%
Taylor expanded in x around 0
lower-+.f6448.8
Simplified48.8%
Final simplification30.9%
(FPCore (x y z) :precision binary64 (if (<= (* y (log y)) 5e+127) (exp (- x z)) (pow y y)))
double code(double x, double y, double z) {
double tmp;
if ((y * log(y)) <= 5e+127) {
tmp = exp((x - z));
} else {
tmp = pow(y, y);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((y * log(y)) <= 5d+127) then
tmp = exp((x - z))
else
tmp = y ** y
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((y * Math.log(y)) <= 5e+127) {
tmp = Math.exp((x - z));
} else {
tmp = Math.pow(y, y);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (y * math.log(y)) <= 5e+127: tmp = math.exp((x - z)) else: tmp = math.pow(y, y) return tmp
function code(x, y, z) tmp = 0.0 if (Float64(y * log(y)) <= 5e+127) tmp = exp(Float64(x - z)); else tmp = y ^ y; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((y * log(y)) <= 5e+127) tmp = exp((x - z)); else tmp = y ^ y; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision], 5e+127], N[Exp[N[(x - z), $MachinePrecision]], $MachinePrecision], N[Power[y, y], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \cdot \log y \leq 5 \cdot 10^{+127}:\\
\;\;\;\;e^{x - z}\\
\mathbf{else}:\\
\;\;\;\;{y}^{y}\\
\end{array}
\end{array}
if (*.f64 y (log.f64 y)) < 5.0000000000000004e127Initial program 100.0%
Taylor expanded in y around 0
lower--.f6492.2
Simplified92.2%
if 5.0000000000000004e127 < (*.f64 y (log.f64 y)) Initial program 100.0%
Taylor expanded in y around inf
mul-1-negN/A
distribute-rgt-neg-inN/A
log-recN/A
remove-double-negN/A
lower-*.f64N/A
lower-log.f6493.3
Simplified93.3%
lift-log.f64N/A
*-commutativeN/A
lift-log.f64N/A
exp-to-powN/A
lower-pow.f6493.3
Applied egg-rr93.3%
(FPCore (x y z)
:precision binary64
(if (<= (- (+ x (* y (log y))) z) 0.002)
(/ 1.0 (fma z (fma z (fma -0.25 (* z z) 0.5) 1.0) 1.0))
(fma
(fma z 0.5 -1.0)
(* z (* (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0)))))
1.0)))
double code(double x, double y, double z) {
double tmp;
if (((x + (y * log(y))) - z) <= 0.002) {
tmp = 1.0 / fma(z, fma(z, fma(-0.25, (z * z), 0.5), 1.0), 1.0);
} else {
tmp = fma(fma(z, 0.5, -1.0), (z * (fma(z, 0.5, -1.0) * (z * (z * fma(z, 0.5, -1.0))))), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (Float64(Float64(x + Float64(y * log(y))) - z) <= 0.002) tmp = Float64(1.0 / fma(z, fma(z, fma(-0.25, Float64(z * z), 0.5), 1.0), 1.0)); else tmp = fma(fma(z, 0.5, -1.0), Float64(z * Float64(fma(z, 0.5, -1.0) * Float64(z * Float64(z * fma(z, 0.5, -1.0))))), 1.0); end return tmp end
code[x_, y_, z_] := If[LessEqual[N[(N[(x + N[(y * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - z), $MachinePrecision], 0.002], N[(1.0 / N[(z * N[(z * N[(-0.25 * N[(z * z), $MachinePrecision] + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(x + y \cdot \log y\right) - z \leq 0.002:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(-0.25, z \cdot z, 0.5\right), 1\right), 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(\mathsf{fma}\left(z, 0.5, -1\right) \cdot \left(z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right)\right)\right), 1\right)\\
\end{array}
\end{array}
if (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) < 2e-3Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6471.7
Simplified71.7%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6429.4
Simplified29.4%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6429.4
Applied egg-rr29.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6471.5
Simplified71.5%
if 2e-3 < (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6437.0
Simplified37.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6425.0
Simplified25.0%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr6.9%
Taylor expanded in z around 0
Simplified48.2%
Final simplification57.6%
(FPCore (x y z)
:precision binary64
(if (<= z -3.6e+55)
(fma
(fma z 0.5 -1.0)
(* z (* (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0)))))
1.0)
(if (<= z 2150000000000.0)
(exp x)
(if (<= z 6.8e+155) (exp z) (/ 1.0 (fma z (fma 0.5 z 1.0) 1.0))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -3.6e+55) {
tmp = fma(fma(z, 0.5, -1.0), (z * (fma(z, 0.5, -1.0) * (z * (z * fma(z, 0.5, -1.0))))), 1.0);
} else if (z <= 2150000000000.0) {
tmp = exp(x);
} else if (z <= 6.8e+155) {
tmp = exp(z);
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -3.6e+55) tmp = fma(fma(z, 0.5, -1.0), Float64(z * Float64(fma(z, 0.5, -1.0) * Float64(z * Float64(z * fma(z, 0.5, -1.0))))), 1.0); elseif (z <= 2150000000000.0) tmp = exp(x); elseif (z <= 6.8e+155) tmp = exp(z); else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -3.6e+55], N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 2150000000000.0], N[Exp[x], $MachinePrecision], If[LessEqual[z, 6.8e+155], N[Exp[z], $MachinePrecision], N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.6 \cdot 10^{+55}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(\mathsf{fma}\left(z, 0.5, -1\right) \cdot \left(z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right)\right)\right), 1\right)\\
\mathbf{elif}\;z \leq 2150000000000:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;e^{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -3.59999999999999987e55Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6492.1
Simplified92.1%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6456.9
Simplified56.9%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr4.0%
Taylor expanded in z around 0
Simplified92.1%
if -3.59999999999999987e55 < z < 2.15e12Initial program 100.0%
Taylor expanded in y around 0
lower--.f6469.1
Simplified69.1%
Taylor expanded in z around 0
lower-exp.f6464.4
Simplified64.4%
if 2.15e12 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6434.4
Simplified34.4%
neg-sub0N/A
flip3--N/A
metadata-evalN/A
neg-sub0N/A
cube-negN/A
lift-neg.f64N/A
sqr-powN/A
unpow-prod-downN/A
lift-neg.f64N/A
lift-neg.f64N/A
sqr-negN/A
pow-prod-downN/A
sqr-powN/A
metadata-evalN/A
+-lft-identityN/A
distribute-rgt-outN/A
+-commutativeN/A
+-lft-identityN/A
pow2N/A
pow-divN/A
metadata-evalN/A
unpow1N/A
lower-exp.f6467.2
Applied egg-rr67.2%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
Final simplification73.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0
(fma
(fma z 0.5 -1.0)
(* z (* (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0)))))
1.0)))
(if (<= z -3.6e+55)
t_0
(if (<= z 3.3e+51)
(exp x)
(if (<= z 6.8e+155) t_0 (/ 1.0 (fma z (fma 0.5 z 1.0) 1.0)))))))
double code(double x, double y, double z) {
double t_0 = fma(fma(z, 0.5, -1.0), (z * (fma(z, 0.5, -1.0) * (z * (z * fma(z, 0.5, -1.0))))), 1.0);
double tmp;
if (z <= -3.6e+55) {
tmp = t_0;
} else if (z <= 3.3e+51) {
tmp = exp(x);
} else if (z <= 6.8e+155) {
tmp = t_0;
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) t_0 = fma(fma(z, 0.5, -1.0), Float64(z * Float64(fma(z, 0.5, -1.0) * Float64(z * Float64(z * fma(z, 0.5, -1.0))))), 1.0) tmp = 0.0 if (z <= -3.6e+55) tmp = t_0; elseif (z <= 3.3e+51) tmp = exp(x); elseif (z <= 6.8e+155) tmp = t_0; else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[z, -3.6e+55], t$95$0, If[LessEqual[z, 3.3e+51], N[Exp[x], $MachinePrecision], If[LessEqual[z, 6.8e+155], t$95$0, N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(\mathsf{fma}\left(z, 0.5, -1\right) \cdot \left(z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right)\right)\right), 1\right)\\
\mathbf{if}\;z \leq -3.6 \cdot 10^{+55}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{+51}:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -3.59999999999999987e55 or 3.2999999999999997e51 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6472.2
Simplified72.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6439.6
Simplified39.6%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr7.7%
Taylor expanded in z around 0
Simplified82.3%
if -3.59999999999999987e55 < z < 3.2999999999999997e51Initial program 100.0%
Taylor expanded in y around 0
lower--.f6468.8
Simplified68.8%
Taylor expanded in z around 0
lower-exp.f6463.6
Simplified63.6%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
Final simplification72.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0
(fma
(fma z 0.5 -1.0)
(* z (* (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0)))))
1.0)))
(if (<= z -5e-83)
t_0
(if (<= z 2.75e-275)
(/ t_0 (* (* z (* z (* z z))) 0.25))
(if (<= z 9e+45)
(fma (* x (fma (* x x) 0.25 -1.0)) (/ 1.0 (fma x 0.5 -1.0)) 1.0)
(if (<= z 6.8e+155) t_0 (/ 1.0 (fma z (fma 0.5 z 1.0) 1.0))))))))
double code(double x, double y, double z) {
double t_0 = fma(fma(z, 0.5, -1.0), (z * (fma(z, 0.5, -1.0) * (z * (z * fma(z, 0.5, -1.0))))), 1.0);
double tmp;
if (z <= -5e-83) {
tmp = t_0;
} else if (z <= 2.75e-275) {
tmp = t_0 / ((z * (z * (z * z))) * 0.25);
} else if (z <= 9e+45) {
tmp = fma((x * fma((x * x), 0.25, -1.0)), (1.0 / fma(x, 0.5, -1.0)), 1.0);
} else if (z <= 6.8e+155) {
tmp = t_0;
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) t_0 = fma(fma(z, 0.5, -1.0), Float64(z * Float64(fma(z, 0.5, -1.0) * Float64(z * Float64(z * fma(z, 0.5, -1.0))))), 1.0) tmp = 0.0 if (z <= -5e-83) tmp = t_0; elseif (z <= 2.75e-275) tmp = Float64(t_0 / Float64(Float64(z * Float64(z * Float64(z * z))) * 0.25)); elseif (z <= 9e+45) tmp = fma(Float64(x * fma(Float64(x * x), 0.25, -1.0)), Float64(1.0 / fma(x, 0.5, -1.0)), 1.0); elseif (z <= 6.8e+155) tmp = t_0; else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[z, -5e-83], t$95$0, If[LessEqual[z, 2.75e-275], N[(t$95$0 / N[(N[(z * N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 9e+45], N[(N[(x * N[(N[(x * x), $MachinePrecision] * 0.25 + -1.0), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(x * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 6.8e+155], t$95$0, N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(\mathsf{fma}\left(z, 0.5, -1\right) \cdot \left(z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right)\right)\right), 1\right)\\
\mathbf{if}\;z \leq -5 \cdot 10^{-83}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 2.75 \cdot 10^{-275}:\\
\;\;\;\;\frac{t\_0}{\left(z \cdot \left(z \cdot \left(z \cdot z\right)\right)\right) \cdot 0.25}\\
\mathbf{elif}\;z \leq 9 \cdot 10^{+45}:\\
\;\;\;\;\mathsf{fma}\left(x \cdot \mathsf{fma}\left(x \cdot x, 0.25, -1\right), \frac{1}{\mathsf{fma}\left(x, 0.5, -1\right)}, 1\right)\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -5e-83 or 8.9999999999999997e45 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6468.8
Simplified68.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6435.4
Simplified35.4%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr12.4%
Taylor expanded in z around 0
Simplified70.4%
if -5e-83 < z < 2.74999999999999994e-275Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6415.9
Simplified15.9%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.9
Simplified15.9%
lift-fma.f64N/A
flip3-+N/A
lower-/.f64N/A
Applied egg-rr15.9%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-plusN/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.1
Simplified64.1%
if 2.74999999999999994e-275 < z < 8.9999999999999997e45Initial program 100.0%
Taylor expanded in y around 0
lower--.f6464.6
Simplified64.6%
Taylor expanded in z around 0
lower-exp.f6463.3
Simplified63.3%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6445.6
Simplified45.6%
lift-fma.f64N/A
*-commutativeN/A
lift-fma.f64N/A
flip-+N/A
associate-*l/N/A
div-invN/A
lower-fma.f64N/A
Applied egg-rr46.9%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
Final simplification64.9%
(FPCore (x y z)
:precision binary64
(if (<= z -2e+100)
(fma (* z (fma z (* z 0.25) -1.0)) (/ -1.0 (fma z -0.5 -1.0)) 1.0)
(if (<= z -1.2e+56)
(/
(fma (fma z 0.5 -1.0) (* z (* z (fma z 0.5 -1.0))) -1.0)
(fma z (fma z 0.5 -1.0) -1.0))
(if (<= z 6.8e+155)
(fma (* x (fma (* x x) 0.25 -1.0)) (/ 1.0 (fma x 0.5 -1.0)) 1.0)
(/ 1.0 (fma z (fma 0.5 z 1.0) 1.0))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -2e+100) {
tmp = fma((z * fma(z, (z * 0.25), -1.0)), (-1.0 / fma(z, -0.5, -1.0)), 1.0);
} else if (z <= -1.2e+56) {
tmp = fma(fma(z, 0.5, -1.0), (z * (z * fma(z, 0.5, -1.0))), -1.0) / fma(z, fma(z, 0.5, -1.0), -1.0);
} else if (z <= 6.8e+155) {
tmp = fma((x * fma((x * x), 0.25, -1.0)), (1.0 / fma(x, 0.5, -1.0)), 1.0);
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -2e+100) tmp = fma(Float64(z * fma(z, Float64(z * 0.25), -1.0)), Float64(-1.0 / fma(z, -0.5, -1.0)), 1.0); elseif (z <= -1.2e+56) tmp = Float64(fma(fma(z, 0.5, -1.0), Float64(z * Float64(z * fma(z, 0.5, -1.0))), -1.0) / fma(z, fma(z, 0.5, -1.0), -1.0)); elseif (z <= 6.8e+155) tmp = fma(Float64(x * fma(Float64(x * x), 0.25, -1.0)), Float64(1.0 / fma(x, 0.5, -1.0)), 1.0); else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -2e+100], N[(N[(z * N[(z * N[(z * 0.25), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision] * N[(-1.0 / N[(z * -0.5 + -1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, -1.2e+56], N[(N[(N[(z * 0.5 + -1.0), $MachinePrecision] * N[(z * N[(z * N[(z * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] / N[(z * N[(z * 0.5 + -1.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 6.8e+155], N[(N[(x * N[(N[(x * x), $MachinePrecision] * 0.25 + -1.0), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(x * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2 \cdot 10^{+100}:\\
\;\;\;\;\mathsf{fma}\left(z \cdot \mathsf{fma}\left(z, z \cdot 0.25, -1\right), \frac{-1}{\mathsf{fma}\left(z, -0.5, -1\right)}, 1\right)\\
\mathbf{elif}\;z \leq -1.2 \cdot 10^{+56}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(z, 0.5, -1\right), z \cdot \left(z \cdot \mathsf{fma}\left(z, 0.5, -1\right)\right), -1\right)}{\mathsf{fma}\left(z, \mathsf{fma}\left(z, 0.5, -1\right), -1\right)}\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;\mathsf{fma}\left(x \cdot \mathsf{fma}\left(x \cdot x, 0.25, -1\right), \frac{1}{\mathsf{fma}\left(x, 0.5, -1\right)}, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -2.00000000000000003e100Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6490.4
Simplified90.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6468.3
Simplified68.3%
lift-fma.f64N/A
*-commutativeN/A
lift-fma.f64N/A
flip-+N/A
associate-*l/N/A
div-invN/A
lower-fma.f64N/A
Applied egg-rr90.4%
if -2.00000000000000003e100 < z < -1.20000000000000007e56Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f64100.0
Simplified100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f645.0
Simplified5.0%
lift-fma.f64N/A
flip-+N/A
lower-/.f64N/A
metadata-evalN/A
sub-negN/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f6478.7
Applied egg-rr78.7%
if -1.20000000000000007e56 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in y around 0
lower--.f6468.1
Simplified68.1%
Taylor expanded in z around 0
lower-exp.f6460.4
Simplified60.4%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6437.3
Simplified37.3%
lift-fma.f64N/A
*-commutativeN/A
lift-fma.f64N/A
flip-+N/A
associate-*l/N/A
div-invN/A
lower-fma.f64N/A
Applied egg-rr41.0%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
Final simplification56.6%
(FPCore (x y z)
:precision binary64
(if (<= z -1.2e+56)
(fma z (fma z (fma z -0.16666666666666666 0.5) -1.0) 1.0)
(if (<= z 6.8e+155)
(fma (* x (fma (* x x) 0.25 -1.0)) (/ 1.0 (fma x 0.5 -1.0)) 1.0)
(/ 1.0 (fma z (fma 0.5 z 1.0) 1.0)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.2e+56) {
tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0);
} else if (z <= 6.8e+155) {
tmp = fma((x * fma((x * x), 0.25, -1.0)), (1.0 / fma(x, 0.5, -1.0)), 1.0);
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.2e+56) tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0); elseif (z <= 6.8e+155) tmp = fma(Float64(x * fma(Float64(x * x), 0.25, -1.0)), Float64(1.0 / fma(x, 0.5, -1.0)), 1.0); else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.2e+56], N[(z * N[(z * N[(z * -0.16666666666666666 + 0.5), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 6.8e+155], N[(N[(x * N[(N[(x * x), $MachinePrecision] * 0.25 + -1.0), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(x * 0.5 + -1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.2 \cdot 10^{+56}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, -0.16666666666666666, 0.5\right), -1\right), 1\right)\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;\mathsf{fma}\left(x \cdot \mathsf{fma}\left(x \cdot x, 0.25, -1\right), \frac{1}{\mathsf{fma}\left(x, 0.5, -1\right)}, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -1.20000000000000007e56Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6492.1
Simplified92.1%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6475.5
Simplified75.5%
if -1.20000000000000007e56 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in y around 0
lower--.f6468.1
Simplified68.1%
Taylor expanded in z around 0
lower-exp.f6460.4
Simplified60.4%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6437.3
Simplified37.3%
lift-fma.f64N/A
*-commutativeN/A
lift-fma.f64N/A
flip-+N/A
associate-*l/N/A
div-invN/A
lower-fma.f64N/A
Applied egg-rr41.0%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
Final simplification54.1%
(FPCore (x y z)
:precision binary64
(if (<= z -1.2e+56)
(fma z (fma z (fma z -0.16666666666666666 0.5) -1.0) 1.0)
(if (<= z 2.55e+81)
(fma x (fma x (fma 0.16666666666666666 x 0.5) 1.0) 1.0)
(/ 1.0 (fma z (fma z (fma -0.25 (* z z) 0.5) 1.0) 1.0)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.2e+56) {
tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0);
} else if (z <= 2.55e+81) {
tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0);
} else {
tmp = 1.0 / fma(z, fma(z, fma(-0.25, (z * z), 0.5), 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.2e+56) tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0); elseif (z <= 2.55e+81) tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0); else tmp = Float64(1.0 / fma(z, fma(z, fma(-0.25, Float64(z * z), 0.5), 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.2e+56], N[(z * N[(z * N[(z * -0.16666666666666666 + 0.5), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 2.55e+81], N[(x * N[(x * N[(0.16666666666666666 * x + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(1.0 / N[(z * N[(z * N[(-0.25 * N[(z * z), $MachinePrecision] + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.2 \cdot 10^{+56}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, -0.16666666666666666, 0.5\right), -1\right), 1\right)\\
\mathbf{elif}\;z \leq 2.55 \cdot 10^{+81}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(0.16666666666666666, x, 0.5\right), 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(-0.25, z \cdot z, 0.5\right), 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -1.20000000000000007e56Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6492.1
Simplified92.1%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6475.5
Simplified75.5%
if -1.20000000000000007e56 < z < 2.5500000000000001e81Initial program 100.0%
Taylor expanded in y around 0
lower--.f6468.1
Simplified68.1%
Taylor expanded in z around 0
lower-exp.f6463.1
Simplified63.1%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6438.7
Simplified38.7%
if 2.5500000000000001e81 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6469.5
Simplified69.5%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6413.2
Simplified13.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6413.2
Applied egg-rr13.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6469.5
Simplified69.5%
(FPCore (x y z)
:precision binary64
(if (<= z -1.2e+56)
(fma z (fma z (fma z -0.16666666666666666 0.5) -1.0) 1.0)
(if (<= z 6.8e+155)
(fma x (fma x (fma 0.16666666666666666 x 0.5) 1.0) 1.0)
(/ 1.0 (fma z (fma 0.5 z 1.0) 1.0)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.2e+56) {
tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0);
} else if (z <= 6.8e+155) {
tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0);
} else {
tmp = 1.0 / fma(z, fma(0.5, z, 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.2e+56) tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0); elseif (z <= 6.8e+155) tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0); else tmp = Float64(1.0 / fma(z, fma(0.5, z, 1.0), 1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.2e+56], N[(z * N[(z * N[(z * -0.16666666666666666 + 0.5), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 6.8e+155], N[(x * N[(x * N[(0.16666666666666666 * x + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(1.0 / N[(z * N[(0.5 * z + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.2 \cdot 10^{+56}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, -0.16666666666666666, 0.5\right), -1\right), 1\right)\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+155}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(0.16666666666666666, x, 0.5\right), 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(z, \mathsf{fma}\left(0.5, z, 1\right), 1\right)}\\
\end{array}
\end{array}
if z < -1.20000000000000007e56Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6492.1
Simplified92.1%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6475.5
Simplified75.5%
if -1.20000000000000007e56 < z < 6.8000000000000002e155Initial program 100.0%
Taylor expanded in y around 0
lower--.f6468.1
Simplified68.1%
Taylor expanded in z around 0
lower-exp.f6460.4
Simplified60.4%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6437.4
Simplified37.4%
if 6.8000000000000002e155 < z Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6486.3
Simplified86.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6415.2
Simplified15.2%
lift-fma.f64N/A
flip-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
lift-fma.f64N/A
lower-/.f6415.2
Applied egg-rr15.2%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6486.3
Simplified86.3%
(FPCore (x y z)
:precision binary64
(if (<= x -2.3e+19)
(* (* z z) 0.5)
(if (<= x 8.5e+40)
(fma z (fma z (fma z -0.16666666666666666 0.5) -1.0) 1.0)
(fma x (fma x (fma 0.16666666666666666 x 0.5) 1.0) 1.0))))
double code(double x, double y, double z) {
double tmp;
if (x <= -2.3e+19) {
tmp = (z * z) * 0.5;
} else if (x <= 8.5e+40) {
tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0);
} else {
tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -2.3e+19) tmp = Float64(Float64(z * z) * 0.5); elseif (x <= 8.5e+40) tmp = fma(z, fma(z, fma(z, -0.16666666666666666, 0.5), -1.0), 1.0); else tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -2.3e+19], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 8.5e+40], N[(z * N[(z * N[(z * -0.16666666666666666 + 0.5), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(x * N[(x * N[(0.16666666666666666 * x + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.3 \cdot 10^{+19}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 8.5 \cdot 10^{+40}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, -0.16666666666666666, 0.5\right), -1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(0.16666666666666666, x, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if x < -2.3e19Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6411.5
Simplified11.5%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6428.8
Simplified28.8%
if -2.3e19 < x < 8.49999999999999996e40Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6466.5
Simplified66.5%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6441.8
Simplified41.8%
if 8.49999999999999996e40 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f6496.5
Simplified96.5%
Taylor expanded in z around 0
lower-exp.f6486.2
Simplified86.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6468.1
Simplified68.1%
Final simplification44.8%
(FPCore (x y z)
:precision binary64
(if (<= x -3.15e+19)
(* (* z z) 0.5)
(if (<= x 4.5e+96)
(- (fma z (* z 0.5) 1.0) z)
(fma x (fma x (fma 0.16666666666666666 x 0.5) 1.0) 1.0))))
double code(double x, double y, double z) {
double tmp;
if (x <= -3.15e+19) {
tmp = (z * z) * 0.5;
} else if (x <= 4.5e+96) {
tmp = fma(z, (z * 0.5), 1.0) - z;
} else {
tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -3.15e+19) tmp = Float64(Float64(z * z) * 0.5); elseif (x <= 4.5e+96) tmp = Float64(fma(z, Float64(z * 0.5), 1.0) - z); else tmp = fma(x, fma(x, fma(0.16666666666666666, x, 0.5), 1.0), 1.0); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -3.15e+19], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 4.5e+96], N[(N[(z * N[(z * 0.5), $MachinePrecision] + 1.0), $MachinePrecision] - z), $MachinePrecision], N[(x * N[(x * N[(0.16666666666666666 * x + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.15 \cdot 10^{+19}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 4.5 \cdot 10^{+96}:\\
\;\;\;\;\mathsf{fma}\left(z, z \cdot 0.5, 1\right) - z\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(0.16666666666666666, x, 0.5\right), 1\right), 1\right)\\
\end{array}
\end{array}
if x < -3.15e19Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6411.5
Simplified11.5%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6428.8
Simplified28.8%
if -3.15e19 < x < 4.49999999999999957e96Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6463.1
Simplified63.1%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6434.9
Simplified34.9%
lift-fma.f64N/A
+-commutativeN/A
lift-fma.f64N/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6434.9
Applied egg-rr34.9%
if 4.49999999999999957e96 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f6497.9
Simplified97.9%
Taylor expanded in z around 0
lower-exp.f6485.0
Simplified85.0%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.1
Simplified83.1%
Final simplification42.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (* z z) 0.5)))
(if (<= x -7.2e-147)
t_0
(if (<= x 1.1e-139)
(- 1.0 z)
(if (<= x 1.85e+152) t_0 (* x (* x 0.5)))))))
double code(double x, double y, double z) {
double t_0 = (z * z) * 0.5;
double tmp;
if (x <= -7.2e-147) {
tmp = t_0;
} else if (x <= 1.1e-139) {
tmp = 1.0 - z;
} else if (x <= 1.85e+152) {
tmp = t_0;
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (z * z) * 0.5d0
if (x <= (-7.2d-147)) then
tmp = t_0
else if (x <= 1.1d-139) then
tmp = 1.0d0 - z
else if (x <= 1.85d+152) then
tmp = t_0
else
tmp = x * (x * 0.5d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (z * z) * 0.5;
double tmp;
if (x <= -7.2e-147) {
tmp = t_0;
} else if (x <= 1.1e-139) {
tmp = 1.0 - z;
} else if (x <= 1.85e+152) {
tmp = t_0;
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
def code(x, y, z): t_0 = (z * z) * 0.5 tmp = 0 if x <= -7.2e-147: tmp = t_0 elif x <= 1.1e-139: tmp = 1.0 - z elif x <= 1.85e+152: tmp = t_0 else: tmp = x * (x * 0.5) return tmp
function code(x, y, z) t_0 = Float64(Float64(z * z) * 0.5) tmp = 0.0 if (x <= -7.2e-147) tmp = t_0; elseif (x <= 1.1e-139) tmp = Float64(1.0 - z); elseif (x <= 1.85e+152) tmp = t_0; else tmp = Float64(x * Float64(x * 0.5)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (z * z) * 0.5; tmp = 0.0; if (x <= -7.2e-147) tmp = t_0; elseif (x <= 1.1e-139) tmp = 1.0 - z; elseif (x <= 1.85e+152) tmp = t_0; else tmp = x * (x * 0.5); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision]}, If[LessEqual[x, -7.2e-147], t$95$0, If[LessEqual[x, 1.1e-139], N[(1.0 - z), $MachinePrecision], If[LessEqual[x, 1.85e+152], t$95$0, N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(z \cdot z\right) \cdot 0.5\\
\mathbf{if}\;x \leq -7.2 \cdot 10^{-147}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 1.1 \cdot 10^{-139}:\\
\;\;\;\;1 - z\\
\mathbf{elif}\;x \leq 1.85 \cdot 10^{+152}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.5\right)\\
\end{array}
\end{array}
if x < -7.20000000000000023e-147 or 1.10000000000000005e-139 < x < 1.84999999999999998e152Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6449.3
Simplified49.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6421.3
Simplified21.3%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6423.5
Simplified23.5%
if -7.20000000000000023e-147 < x < 1.10000000000000005e-139Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6469.2
Simplified69.2%
Taylor expanded in z around 0
neg-mul-1N/A
unsub-negN/A
lower--.f6433.2
Simplified33.2%
if 1.84999999999999998e152 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f64100.0
Simplified100.0%
Taylor expanded in z around 0
lower-exp.f6489.9
Simplified89.9%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6487.7
Simplified87.7%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6487.7
Simplified87.7%
Final simplification35.9%
(FPCore (x y z) :precision binary64 (if (<= x -3.15e+19) (* (* z z) 0.5) (if (<= x 1.85e+152) (- (fma z (* z 0.5) 1.0) z) (* x (* x 0.5)))))
double code(double x, double y, double z) {
double tmp;
if (x <= -3.15e+19) {
tmp = (z * z) * 0.5;
} else if (x <= 1.85e+152) {
tmp = fma(z, (z * 0.5), 1.0) - z;
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -3.15e+19) tmp = Float64(Float64(z * z) * 0.5); elseif (x <= 1.85e+152) tmp = Float64(fma(z, Float64(z * 0.5), 1.0) - z); else tmp = Float64(x * Float64(x * 0.5)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -3.15e+19], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 1.85e+152], N[(N[(z * N[(z * 0.5), $MachinePrecision] + 1.0), $MachinePrecision] - z), $MachinePrecision], N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.15 \cdot 10^{+19}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 1.85 \cdot 10^{+152}:\\
\;\;\;\;\mathsf{fma}\left(z, z \cdot 0.5, 1\right) - z\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.5\right)\\
\end{array}
\end{array}
if x < -3.15e19Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6411.5
Simplified11.5%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6428.8
Simplified28.8%
if -3.15e19 < x < 1.84999999999999998e152Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6462.3
Simplified62.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6434.1
Simplified34.1%
lift-fma.f64N/A
+-commutativeN/A
lift-fma.f64N/A
distribute-lft-inN/A
associate-+r+N/A
*-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6434.1
Applied egg-rr34.1%
if 1.84999999999999998e152 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f64100.0
Simplified100.0%
Taylor expanded in z around 0
lower-exp.f6489.9
Simplified89.9%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6487.7
Simplified87.7%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6487.7
Simplified87.7%
Final simplification41.1%
(FPCore (x y z) :precision binary64 (if (<= x -3.15e+19) (* (* z z) 0.5) (if (<= x 1.85e+152) (fma z (fma z 0.5 -1.0) 1.0) (* x (* x 0.5)))))
double code(double x, double y, double z) {
double tmp;
if (x <= -3.15e+19) {
tmp = (z * z) * 0.5;
} else if (x <= 1.85e+152) {
tmp = fma(z, fma(z, 0.5, -1.0), 1.0);
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -3.15e+19) tmp = Float64(Float64(z * z) * 0.5); elseif (x <= 1.85e+152) tmp = fma(z, fma(z, 0.5, -1.0), 1.0); else tmp = Float64(x * Float64(x * 0.5)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -3.15e+19], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 1.85e+152], N[(z * N[(z * 0.5 + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.15 \cdot 10^{+19}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 1.85 \cdot 10^{+152}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, 0.5, -1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.5\right)\\
\end{array}
\end{array}
if x < -3.15e19Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6411.5
Simplified11.5%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6428.8
Simplified28.8%
if -3.15e19 < x < 1.84999999999999998e152Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6462.3
Simplified62.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6434.1
Simplified34.1%
if 1.84999999999999998e152 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f64100.0
Simplified100.0%
Taylor expanded in z around 0
lower-exp.f6489.9
Simplified89.9%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6487.7
Simplified87.7%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6487.7
Simplified87.7%
Final simplification41.1%
(FPCore (x y z) :precision binary64 (if (<= x -3.15e+19) (* (* z z) 0.5) (if (<= x 1.85e+152) (fma z (* z 0.5) 1.0) (* x (* x 0.5)))))
double code(double x, double y, double z) {
double tmp;
if (x <= -3.15e+19) {
tmp = (z * z) * 0.5;
} else if (x <= 1.85e+152) {
tmp = fma(z, (z * 0.5), 1.0);
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -3.15e+19) tmp = Float64(Float64(z * z) * 0.5); elseif (x <= 1.85e+152) tmp = fma(z, Float64(z * 0.5), 1.0); else tmp = Float64(x * Float64(x * 0.5)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -3.15e+19], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[x, 1.85e+152], N[(z * N[(z * 0.5), $MachinePrecision] + 1.0), $MachinePrecision], N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.15 \cdot 10^{+19}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;x \leq 1.85 \cdot 10^{+152}:\\
\;\;\;\;\mathsf{fma}\left(z, z \cdot 0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.5\right)\\
\end{array}
\end{array}
if x < -3.15e19Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6411.5
Simplified11.5%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6428.8
Simplified28.8%
if -3.15e19 < x < 1.84999999999999998e152Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6462.3
Simplified62.3%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6434.1
Simplified34.1%
Taylor expanded in z around inf
lower-*.f6433.9
Simplified33.9%
if 1.84999999999999998e152 < x Initial program 100.0%
Taylor expanded in y around 0
lower--.f64100.0
Simplified100.0%
Taylor expanded in z around 0
lower-exp.f6489.9
Simplified89.9%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6487.7
Simplified87.7%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6487.7
Simplified87.7%
Final simplification41.0%
(FPCore (x y z) :precision binary64 (if (<= z -1.7e+152) (* (* z z) 0.5) (if (<= z 1.25e+15) (fma x (* x 0.5) 1.0) (* x (* x 0.5)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.7e+152) {
tmp = (z * z) * 0.5;
} else if (z <= 1.25e+15) {
tmp = fma(x, (x * 0.5), 1.0);
} else {
tmp = x * (x * 0.5);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.7e+152) tmp = Float64(Float64(z * z) * 0.5); elseif (z <= 1.25e+15) tmp = fma(x, Float64(x * 0.5), 1.0); else tmp = Float64(x * Float64(x * 0.5)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.7e+152], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[z, 1.25e+15], N[(x * N[(x * 0.5), $MachinePrecision] + 1.0), $MachinePrecision], N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.7 \cdot 10^{+152}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;z \leq 1.25 \cdot 10^{+15}:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot 0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.5\right)\\
\end{array}
\end{array}
if z < -1.7000000000000001e152Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6496.7
Simplified96.7%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6491.1
Simplified91.1%
Taylor expanded in z around inf
lower-*.f64N/A
unpow2N/A
lower-*.f6493.8
Simplified93.8%
if -1.7000000000000001e152 < z < 1.25e15Initial program 100.0%
Taylor expanded in y around 0
lower--.f6472.6
Simplified72.6%
Taylor expanded in z around 0
lower-exp.f6460.5
Simplified60.5%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6434.8
Simplified34.8%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f6433.9
Simplified33.9%
if 1.25e15 < z Initial program 100.0%
Taylor expanded in y around 0
lower--.f6479.7
Simplified79.7%
Taylor expanded in z around 0
lower-exp.f6441.2
Simplified41.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6418.9
Simplified18.9%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6429.9
Simplified29.9%
Final simplification39.9%
(FPCore (x y z) :precision binary64 (+ x 1.0))
double code(double x, double y, double z) {
return x + 1.0;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x + 1.0d0
end function
public static double code(double x, double y, double z) {
return x + 1.0;
}
def code(x, y, z): return x + 1.0
function code(x, y, z) return Float64(x + 1.0) end
function tmp = code(x, y, z) tmp = x + 1.0; end
code[x_, y_, z_] := N[(x + 1.0), $MachinePrecision]
\begin{array}{l}
\\
x + 1
\end{array}
Initial program 100.0%
Taylor expanded in y around 0
lower--.f6477.7
Simplified77.7%
Taylor expanded in z around 0
lower-exp.f6450.5
Simplified50.5%
Taylor expanded in x around 0
lower-+.f6414.4
Simplified14.4%
Final simplification14.4%
(FPCore (x y z) :precision binary64 1.0)
double code(double x, double y, double z) {
return 1.0;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 1.0d0
end function
public static double code(double x, double y, double z) {
return 1.0;
}
def code(x, y, z): return 1.0
function code(x, y, z) return 1.0 end
function tmp = code(x, y, z) tmp = 1.0; end
code[x_, y_, z_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
lower-neg.f6451.0
Simplified51.0%
Taylor expanded in z around 0
Simplified13.8%
(FPCore (x y z) :precision binary64 (exp (+ (- x z) (* (log y) y))))
double code(double x, double y, double z) {
return exp(((x - z) + (log(y) * y)));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = exp(((x - z) + (log(y) * y)))
end function
public static double code(double x, double y, double z) {
return Math.exp(((x - z) + (Math.log(y) * y)));
}
def code(x, y, z): return math.exp(((x - z) + (math.log(y) * y)))
function code(x, y, z) return exp(Float64(Float64(x - z) + Float64(log(y) * y))) end
function tmp = code(x, y, z) tmp = exp(((x - z) + (log(y) * y))); end
code[x_, y_, z_] := N[Exp[N[(N[(x - z), $MachinePrecision] + N[(N[Log[y], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x - z\right) + \log y \cdot y}
\end{array}
herbie shell --seed 2024219
(FPCore (x y z)
:name "Statistics.Distribution.Poisson.Internal:probability from math-functions-0.1.5.2"
:precision binary64
:alt
(! :herbie-platform default (exp (+ (- x z) (* (log y) y))))
(exp (- (+ x (* y (log y))) z)))