
(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 19 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 -5e+125) (exp x) (if (<= t_0 2e+67) (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 <= -5e+125) {
tmp = exp(x);
} else if (t_0 <= 2e+67) {
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 <= (-5d+125)) then
tmp = exp(x)
else if (t_0 <= 2d+67) 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 <= -5e+125) {
tmp = Math.exp(x);
} else if (t_0 <= 2e+67) {
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 <= -5e+125: tmp = math.exp(x) elif t_0 <= 2e+67: 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 <= -5e+125) tmp = exp(x); elseif (t_0 <= 2e+67) 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 <= -5e+125) tmp = exp(x); elseif (t_0 <= 2e+67) 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, -5e+125], N[Exp[x], $MachinePrecision], If[LessEqual[t$95$0, 2e+67], 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 -5 \cdot 10^{+125}:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+67}:\\
\;\;\;\;e^{-z}\\
\mathbf{else}:\\
\;\;\;\;{y}^{y}\\
\end{array}
\end{array}
if (+.f64 x (*.f64 y (log.f64 y))) < -4.99999999999999962e125Initial program 100.0%
Taylor expanded in x around inf
Simplified96.2%
if -4.99999999999999962e125 < (+.f64 x (*.f64 y (log.f64 y))) < 1.99999999999999997e67Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6484.4
Simplified84.4%
if 1.99999999999999997e67 < (+.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
*-lowering-*.f64N/A
log-lowering-log.f6476.0
Simplified76.0%
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6476.0
Applied egg-rr76.0%
(FPCore (x y z) :precision binary64 (let* ((t_0 (- (+ x (* y (log y))) z)) (t_1 (* (* x x) 0.5))) (if (<= t_0 -2e+26) t_1 (if (<= t_0 2e+20) (+ x 1.0) t_1))))
double code(double x, double y, double z) {
double t_0 = (x + (y * log(y))) - z;
double t_1 = (x * x) * 0.5;
double tmp;
if (t_0 <= -2e+26) {
tmp = t_1;
} else if (t_0 <= 2e+20) {
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 = (x * x) * 0.5d0
if (t_0 <= (-2d+26)) then
tmp = t_1
else if (t_0 <= 2d+20) 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 = (x * x) * 0.5;
double tmp;
if (t_0 <= -2e+26) {
tmp = t_1;
} else if (t_0 <= 2e+20) {
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 = (x * x) * 0.5 tmp = 0 if t_0 <= -2e+26: tmp = t_1 elif t_0 <= 2e+20: 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(x * x) * 0.5) tmp = 0.0 if (t_0 <= -2e+26) tmp = t_1; elseif (t_0 <= 2e+20) 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 = (x * x) * 0.5; tmp = 0.0; if (t_0 <= -2e+26) tmp = t_1; elseif (t_0 <= 2e+20) 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[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]}, If[LessEqual[t$95$0, -2e+26], t$95$1, If[LessEqual[t$95$0, 2e+20], 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(x \cdot x\right) \cdot 0.5\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+20}:\\
\;\;\;\;x + 1\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) < -2.0000000000000001e26 or 2e20 < (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) Initial program 100.0%
Taylor expanded in x around inf
Simplified45.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6422.5
Simplified22.5%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
if -2.0000000000000001e26 < (-.f64 (+.f64 x (*.f64 y (log.f64 y))) z) < 2e20Initial program 100.0%
Taylor expanded in x around inf
Simplified88.8%
Taylor expanded in x around 0
+-lowering-+.f6475.5
Simplified75.5%
Final simplification34.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (fma (* z (* z z)) -0.004629629629629629 0.125)))
(if (<= z -1e+62)
(fma z (fma z (/ t_0 0.25) -1.0) 1.0)
(if (<= z -2.35e-178)
(exp x)
(if (<= z -1.25e-304)
(fma z (fma z (/ t_0 (* z (* z 0.027777777777777776))) -1.0) 1.0)
(if (<= z 5e+80) (exp x) (* (* x x) (* x 0.16666666666666666))))))))
double code(double x, double y, double z) {
double t_0 = fma((z * (z * z)), -0.004629629629629629, 0.125);
double tmp;
if (z <= -1e+62) {
tmp = fma(z, fma(z, (t_0 / 0.25), -1.0), 1.0);
} else if (z <= -2.35e-178) {
tmp = exp(x);
} else if (z <= -1.25e-304) {
tmp = fma(z, fma(z, (t_0 / (z * (z * 0.027777777777777776))), -1.0), 1.0);
} else if (z <= 5e+80) {
tmp = exp(x);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) t_0 = fma(Float64(z * Float64(z * z)), -0.004629629629629629, 0.125) tmp = 0.0 if (z <= -1e+62) tmp = fma(z, fma(z, Float64(t_0 / 0.25), -1.0), 1.0); elseif (z <= -2.35e-178) tmp = exp(x); elseif (z <= -1.25e-304) tmp = fma(z, fma(z, Float64(t_0 / Float64(z * Float64(z * 0.027777777777777776))), -1.0), 1.0); elseif (z <= 5e+80) tmp = exp(x); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.004629629629629629 + 0.125), $MachinePrecision]}, If[LessEqual[z, -1e+62], N[(z * N[(z * N[(t$95$0 / 0.25), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, -2.35e-178], N[Exp[x], $MachinePrecision], If[LessEqual[z, -1.25e-304], N[(z * N[(z * N[(t$95$0 / N[(z * N[(z * 0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 5e+80], N[Exp[x], $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(z \cdot \left(z \cdot z\right), -0.004629629629629629, 0.125\right)\\
\mathbf{if}\;z \leq -1 \cdot 10^{+62}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{0.25}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq -2.35 \cdot 10^{-178}:\\
\;\;\;\;e^{x}\\
\mathbf{elif}\;z \leq -1.25 \cdot 10^{-304}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{z \cdot \left(z \cdot 0.027777777777777776\right)}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq 5 \cdot 10^{+80}:\\
\;\;\;\;e^{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.00000000000000004e62Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6491.4
Simplified91.4%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6483.2
Simplified83.2%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr14.6%
Taylor expanded in z around 0
Simplified91.4%
if -1.00000000000000004e62 < z < -2.35e-178 or -1.24999999999999991e-304 < z < 4.99999999999999961e80Initial program 100.0%
Taylor expanded in x around inf
Simplified62.9%
if -2.35e-178 < z < -1.24999999999999991e-304Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f648.2
Simplified8.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f648.2
Simplified8.2%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr8.2%
Taylor expanded in z around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6494.9
Simplified94.9%
if 4.99999999999999961e80 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified32.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6411.8
Simplified11.8%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6445.6
Simplified45.6%
(FPCore (x y z) :precision binary64 (if (<= y 140.0) (exp x) (pow y y)))
double code(double x, double y, double z) {
double tmp;
if (y <= 140.0) {
tmp = exp(x);
} 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 <= 140.0d0) then
tmp = exp(x)
else
tmp = y ** y
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= 140.0) {
tmp = Math.exp(x);
} else {
tmp = Math.pow(y, y);
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= 140.0: tmp = math.exp(x) else: tmp = math.pow(y, y) return tmp
function code(x, y, z) tmp = 0.0 if (y <= 140.0) tmp = exp(x); else tmp = y ^ y; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= 140.0) tmp = exp(x); else tmp = y ^ y; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, 140.0], N[Exp[x], $MachinePrecision], N[Power[y, y], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 140:\\
\;\;\;\;e^{x}\\
\mathbf{else}:\\
\;\;\;\;{y}^{y}\\
\end{array}
\end{array}
if y < 140Initial program 100.0%
Taylor expanded in x around inf
Simplified63.9%
if 140 < 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
*-lowering-*.f64N/A
log-lowering-log.f6482.0
Simplified82.0%
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6482.0
Applied egg-rr82.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (fma (* z (* z z)) -0.004629629629629629 0.125)))
(if (<= z -1.32e+23)
(fma z (fma z (/ t_0 0.25) -1.0) 1.0)
(if (<= z -1.75e-164)
(fma x (fma x 0.5 1.0) 1.0)
(if (<= z 1.7e-153)
(fma z (fma z (/ t_0 (* z (* z 0.027777777777777776))) -1.0) 1.0)
(if (<= z 2.3e-68)
(fma
(* z (fma 0.25 (* z z) -1.0))
(/ (- 2.0 (/ (+ 4.0 (/ (+ (/ 16.0 z) -8.0) z)) z)) z)
1.0)
(* (* x x) (* x 0.16666666666666666))))))))
double code(double x, double y, double z) {
double t_0 = fma((z * (z * z)), -0.004629629629629629, 0.125);
double tmp;
if (z <= -1.32e+23) {
tmp = fma(z, fma(z, (t_0 / 0.25), -1.0), 1.0);
} else if (z <= -1.75e-164) {
tmp = fma(x, fma(x, 0.5, 1.0), 1.0);
} else if (z <= 1.7e-153) {
tmp = fma(z, fma(z, (t_0 / (z * (z * 0.027777777777777776))), -1.0), 1.0);
} else if (z <= 2.3e-68) {
tmp = fma((z * fma(0.25, (z * z), -1.0)), ((2.0 - ((4.0 + (((16.0 / z) + -8.0) / z)) / z)) / z), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) t_0 = fma(Float64(z * Float64(z * z)), -0.004629629629629629, 0.125) tmp = 0.0 if (z <= -1.32e+23) tmp = fma(z, fma(z, Float64(t_0 / 0.25), -1.0), 1.0); elseif (z <= -1.75e-164) tmp = fma(x, fma(x, 0.5, 1.0), 1.0); elseif (z <= 1.7e-153) tmp = fma(z, fma(z, Float64(t_0 / Float64(z * Float64(z * 0.027777777777777776))), -1.0), 1.0); elseif (z <= 2.3e-68) tmp = fma(Float64(z * fma(0.25, Float64(z * z), -1.0)), Float64(Float64(2.0 - Float64(Float64(4.0 + Float64(Float64(Float64(16.0 / z) + -8.0) / z)) / z)) / z), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.004629629629629629 + 0.125), $MachinePrecision]}, If[LessEqual[z, -1.32e+23], N[(z * N[(z * N[(t$95$0 / 0.25), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, -1.75e-164], N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 1.7e-153], N[(z * N[(z * N[(t$95$0 / N[(z * N[(z * 0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 2.3e-68], N[(N[(z * N[(0.25 * N[(z * z), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 - N[(N[(4.0 + N[(N[(N[(16.0 / z), $MachinePrecision] + -8.0), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(z \cdot \left(z \cdot z\right), -0.004629629629629629, 0.125\right)\\
\mathbf{if}\;z \leq -1.32 \cdot 10^{+23}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{0.25}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq -1.75 \cdot 10^{-164}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), 1\right)\\
\mathbf{elif}\;z \leq 1.7 \cdot 10^{-153}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{z \cdot \left(z \cdot 0.027777777777777776\right)}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq 2.3 \cdot 10^{-68}:\\
\;\;\;\;\mathsf{fma}\left(z \cdot \mathsf{fma}\left(0.25, z \cdot z, -1\right), \frac{2 - \frac{4 + \frac{\frac{16}{z} + -8}{z}}{z}}{z}, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.3199999999999999e23Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.2
Simplified89.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6474.5
Simplified74.5%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr13.4%
Taylor expanded in z around 0
Simplified81.8%
if -1.3199999999999999e23 < z < -1.75e-164Initial program 100.0%
Taylor expanded in x around inf
Simplified55.5%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6443.4
Simplified43.4%
if -1.75e-164 < z < 1.6999999999999999e-153Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6418.7
Simplified18.7%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6418.7
Simplified18.7%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr18.7%
Taylor expanded in z around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6468.9
Simplified68.9%
if 1.6999999999999999e-153 < z < 2.29999999999999997e-68Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6414.5
Simplified14.5%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f6414.5
Simplified14.5%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
sub-negN/A
swap-sqrN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6414.5
Applied egg-rr14.5%
Taylor expanded in z around -inf
Simplified59.7%
if 2.29999999999999997e-68 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified41.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6418.1
Simplified18.1%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6439.0
Simplified39.0%
Final simplification58.3%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (fma (* z (* z z)) -0.004629629629629629 0.125)))
(if (<= z -3.5e+22)
(fma z (fma z (/ t_0 0.25) -1.0) 1.0)
(if (<= z -4.5e-164)
(fma x (fma x 0.5 1.0) 1.0)
(if (<= z 1.45e-154)
(fma z (fma z (/ t_0 (* z (* z 0.027777777777777776))) -1.0) 1.0)
(if (<= z 660.0)
(fma x (* x (fma x 0.16666666666666666 0.5)) 1.0)
(* (* x x) (* x 0.16666666666666666))))))))
double code(double x, double y, double z) {
double t_0 = fma((z * (z * z)), -0.004629629629629629, 0.125);
double tmp;
if (z <= -3.5e+22) {
tmp = fma(z, fma(z, (t_0 / 0.25), -1.0), 1.0);
} else if (z <= -4.5e-164) {
tmp = fma(x, fma(x, 0.5, 1.0), 1.0);
} else if (z <= 1.45e-154) {
tmp = fma(z, fma(z, (t_0 / (z * (z * 0.027777777777777776))), -1.0), 1.0);
} else if (z <= 660.0) {
tmp = fma(x, (x * fma(x, 0.16666666666666666, 0.5)), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) t_0 = fma(Float64(z * Float64(z * z)), -0.004629629629629629, 0.125) tmp = 0.0 if (z <= -3.5e+22) tmp = fma(z, fma(z, Float64(t_0 / 0.25), -1.0), 1.0); elseif (z <= -4.5e-164) tmp = fma(x, fma(x, 0.5, 1.0), 1.0); elseif (z <= 1.45e-154) tmp = fma(z, fma(z, Float64(t_0 / Float64(z * Float64(z * 0.027777777777777776))), -1.0), 1.0); elseif (z <= 660.0) tmp = fma(x, Float64(x * fma(x, 0.16666666666666666, 0.5)), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.004629629629629629 + 0.125), $MachinePrecision]}, If[LessEqual[z, -3.5e+22], N[(z * N[(z * N[(t$95$0 / 0.25), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, -4.5e-164], N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 1.45e-154], N[(z * N[(z * N[(t$95$0 / N[(z * N[(z * 0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 660.0], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(z \cdot \left(z \cdot z\right), -0.004629629629629629, 0.125\right)\\
\mathbf{if}\;z \leq -3.5 \cdot 10^{+22}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{0.25}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq -4.5 \cdot 10^{-164}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), 1\right)\\
\mathbf{elif}\;z \leq 1.45 \cdot 10^{-154}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{t\_0}{z \cdot \left(z \cdot 0.027777777777777776\right)}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq 660:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -3.5e22Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.2
Simplified89.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6474.5
Simplified74.5%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr13.4%
Taylor expanded in z around 0
Simplified81.8%
if -3.5e22 < z < -4.4999999999999997e-164Initial program 100.0%
Taylor expanded in x around inf
Simplified55.5%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6443.4
Simplified43.4%
if -4.4999999999999997e-164 < z < 1.45e-154Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6417.2
Simplified17.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6417.2
Simplified17.2%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr17.2%
Taylor expanded in z around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6469.8
Simplified69.8%
if 1.45e-154 < z < 660Initial program 100.0%
Taylor expanded in x around inf
Simplified69.6%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6442.1
Simplified42.1%
Taylor expanded in x around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6442.1
Simplified42.1%
if 660 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
(FPCore (x y z)
:precision binary64
(if (<= z -1.32e+23)
(fma
z
(fma z (/ (fma (* z (* z z)) -0.004629629629629629 0.125) 0.25) -1.0)
1.0)
(if (<= z 2.9e-228)
(fma x (fma x (fma x 0.16666666666666666 0.5) 1.0) 1.0)
(if (<= z 1.45e-154)
(fma (* z (fma 0.25 (* z z) -1.0)) (/ (+ 2.0 (/ -4.0 z)) z) 1.0)
(if (<= z 580.0)
(fma x (* x (fma x 0.16666666666666666 0.5)) 1.0)
(* (* x x) (* x 0.16666666666666666)))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.32e+23) {
tmp = fma(z, fma(z, (fma((z * (z * z)), -0.004629629629629629, 0.125) / 0.25), -1.0), 1.0);
} else if (z <= 2.9e-228) {
tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0);
} else if (z <= 1.45e-154) {
tmp = fma((z * fma(0.25, (z * z), -1.0)), ((2.0 + (-4.0 / z)) / z), 1.0);
} else if (z <= 580.0) {
tmp = fma(x, (x * fma(x, 0.16666666666666666, 0.5)), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.32e+23) tmp = fma(z, fma(z, Float64(fma(Float64(z * Float64(z * z)), -0.004629629629629629, 0.125) / 0.25), -1.0), 1.0); elseif (z <= 2.9e-228) tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0); elseif (z <= 1.45e-154) tmp = fma(Float64(z * fma(0.25, Float64(z * z), -1.0)), Float64(Float64(2.0 + Float64(-4.0 / z)) / z), 1.0); elseif (z <= 580.0) tmp = fma(x, Float64(x * fma(x, 0.16666666666666666, 0.5)), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.32e+23], N[(z * N[(z * N[(N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.004629629629629629 + 0.125), $MachinePrecision] / 0.25), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 2.9e-228], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 1.45e-154], N[(N[(z * N[(0.25 * N[(z * z), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 + N[(-4.0 / z), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 580.0], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.32 \cdot 10^{+23}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{\mathsf{fma}\left(z \cdot \left(z \cdot z\right), -0.004629629629629629, 0.125\right)}{0.25}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq 2.9 \cdot 10^{-228}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{elif}\;z \leq 1.45 \cdot 10^{-154}:\\
\;\;\;\;\mathsf{fma}\left(z \cdot \mathsf{fma}\left(0.25, z \cdot z, -1\right), \frac{2 + \frac{-4}{z}}{z}, 1\right)\\
\mathbf{elif}\;z \leq 580:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.3199999999999999e23Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.2
Simplified89.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6474.5
Simplified74.5%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr13.4%
Taylor expanded in z around 0
Simplified81.8%
if -1.3199999999999999e23 < z < 2.9000000000000001e-228Initial program 100.0%
Taylor expanded in x around inf
Simplified59.0%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6442.7
Simplified42.7%
if 2.9000000000000001e-228 < z < 1.45e-154Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6412.4
Simplified12.4%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f6412.4
Simplified12.4%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
sub-negN/A
swap-sqrN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6412.4
Applied egg-rr12.4%
Taylor expanded in z around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval72.0
Simplified72.0%
if 1.45e-154 < z < 580Initial program 100.0%
Taylor expanded in x around inf
Simplified69.6%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6442.1
Simplified42.1%
Taylor expanded in x around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6442.1
Simplified42.1%
if 580 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
Final simplification54.7%
(FPCore (x y z)
:precision binary64
(if (<= z -6.5e+20)
(fma
z
(fma z (/ (fma (* z (* z z)) -0.004629629629629629 0.125) 0.25) -1.0)
1.0)
(if (<= z 1020.0)
(fma x (fma x (fma x 0.16666666666666666 0.5) 1.0) 1.0)
(* (* x x) (* x 0.16666666666666666)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -6.5e+20) {
tmp = fma(z, fma(z, (fma((z * (z * z)), -0.004629629629629629, 0.125) / 0.25), -1.0), 1.0);
} else if (z <= 1020.0) {
tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -6.5e+20) tmp = fma(z, fma(z, Float64(fma(Float64(z * Float64(z * z)), -0.004629629629629629, 0.125) / 0.25), -1.0), 1.0); elseif (z <= 1020.0) tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -6.5e+20], N[(z * N[(z * N[(N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.004629629629629629 + 0.125), $MachinePrecision] / 0.25), $MachinePrecision] + -1.0), $MachinePrecision] + 1.0), $MachinePrecision], If[LessEqual[z, 1020.0], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6.5 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(z, \mathsf{fma}\left(z, \frac{\mathsf{fma}\left(z \cdot \left(z \cdot z\right), -0.004629629629629629, 0.125\right)}{0.25}, -1\right), 1\right)\\
\mathbf{elif}\;z \leq 1020:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -6.5e20Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.2
Simplified89.2%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6474.5
Simplified74.5%
+-commutativeN/A
flip3-+N/A
+-commutativeN/A
/-lowering-/.f64N/A
unpow-prod-downN/A
accelerator-lowering-fma.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-evalN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-commutativeN/A
sub-negN/A
swap-sqrN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
Applied egg-rr13.4%
Taylor expanded in z around 0
Simplified81.8%
if -6.5e20 < z < 1020Initial program 100.0%
Taylor expanded in x around inf
Simplified63.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6440.6
Simplified40.6%
if 1020 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
(FPCore (x y z)
:precision binary64
(if (<= z -1.5e+88)
(* (* z (* z z)) -0.16666666666666666)
(if (<= z 520.0)
(fma x (fma x (fma x 0.16666666666666666 0.5) 1.0) 1.0)
(* (* x x) (* x 0.16666666666666666)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.5e+88) {
tmp = (z * (z * z)) * -0.16666666666666666;
} else if (z <= 520.0) {
tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.5e+88) tmp = Float64(Float64(z * Float64(z * z)) * -0.16666666666666666); elseif (z <= 520.0) tmp = fma(x, fma(x, fma(x, 0.16666666666666666, 0.5), 1.0), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.5e+88], N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision], If[LessEqual[z, 520.0], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+88}:\\
\;\;\;\;\left(z \cdot \left(z \cdot z\right)\right) \cdot -0.16666666666666666\\
\mathbf{elif}\;z \leq 520:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.50000000000000003e88Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6490.9
Simplified90.9%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.5
Simplified87.5%
Taylor expanded in z around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.5
Simplified87.5%
if -1.50000000000000003e88 < z < 520Initial program 100.0%
Taylor expanded in x around inf
Simplified62.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6439.2
Simplified39.2%
if 520 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
Final simplification50.2%
(FPCore (x y z)
:precision binary64
(if (<= z -1.5e+88)
(* (* z (* z z)) -0.16666666666666666)
(if (<= z 580.0)
(fma x (* x (fma x 0.16666666666666666 0.5)) 1.0)
(* (* x x) (* x 0.16666666666666666)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.5e+88) {
tmp = (z * (z * z)) * -0.16666666666666666;
} else if (z <= 580.0) {
tmp = fma(x, (x * fma(x, 0.16666666666666666, 0.5)), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.5e+88) tmp = Float64(Float64(z * Float64(z * z)) * -0.16666666666666666); elseif (z <= 580.0) tmp = fma(x, Float64(x * fma(x, 0.16666666666666666, 0.5)), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.5e+88], N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision], If[LessEqual[z, 580.0], N[(x * N[(x * N[(x * 0.16666666666666666 + 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+88}:\\
\;\;\;\;\left(z \cdot \left(z \cdot z\right)\right) \cdot -0.16666666666666666\\
\mathbf{elif}\;z \leq 580:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot \mathsf{fma}\left(x, 0.16666666666666666, 0.5\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.50000000000000003e88Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6490.9
Simplified90.9%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.5
Simplified87.5%
Taylor expanded in z around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.5
Simplified87.5%
if -1.50000000000000003e88 < z < 580Initial program 100.0%
Taylor expanded in x around inf
Simplified62.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6439.2
Simplified39.2%
Taylor expanded in x around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6439.1
Simplified39.1%
if 580 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
Final simplification50.1%
(FPCore (x y z)
:precision binary64
(if (<= z -1.5e+88)
(* (* z (* z z)) -0.16666666666666666)
(if (<= z 15000000.0)
(fma x (* x (* x 0.16666666666666666)) 1.0)
(* (* x x) (* x 0.16666666666666666)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.5e+88) {
tmp = (z * (z * z)) * -0.16666666666666666;
} else if (z <= 15000000.0) {
tmp = fma(x, (x * (x * 0.16666666666666666)), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.5e+88) tmp = Float64(Float64(z * Float64(z * z)) * -0.16666666666666666); elseif (z <= 15000000.0) tmp = fma(x, Float64(x * Float64(x * 0.16666666666666666)), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.5e+88], N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision], If[LessEqual[z, 15000000.0], N[(x * N[(x * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+88}:\\
\;\;\;\;\left(z \cdot \left(z \cdot z\right)\right) \cdot -0.16666666666666666\\
\mathbf{elif}\;z \leq 15000000:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot \left(x \cdot 0.16666666666666666\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.50000000000000003e88Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6490.9
Simplified90.9%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.5
Simplified87.5%
Taylor expanded in z around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.5
Simplified87.5%
if -1.50000000000000003e88 < z < 1.5e7Initial program 100.0%
Taylor expanded in x around inf
Simplified62.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6439.2
Simplified39.2%
Taylor expanded in x around inf
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6439.1
Simplified39.1%
if 1.5e7 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.7
Simplified12.7%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3
Simplified42.3%
Final simplification50.1%
(FPCore (x y z)
:precision binary64
(if (<= z -1.5e+88)
(* (* z (* z z)) -0.16666666666666666)
(if (<= z 1.7e+39)
(fma x (fma x 0.5 1.0) 1.0)
(* (* x x) (* x 0.16666666666666666)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.5e+88) {
tmp = (z * (z * z)) * -0.16666666666666666;
} else if (z <= 1.7e+39) {
tmp = fma(x, fma(x, 0.5, 1.0), 1.0);
} else {
tmp = (x * x) * (x * 0.16666666666666666);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.5e+88) tmp = Float64(Float64(z * Float64(z * z)) * -0.16666666666666666); elseif (z <= 1.7e+39) tmp = fma(x, fma(x, 0.5, 1.0), 1.0); else tmp = Float64(Float64(x * x) * Float64(x * 0.16666666666666666)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.5e+88], N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision], If[LessEqual[z, 1.7e+39], N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+88}:\\
\;\;\;\;\left(z \cdot \left(z \cdot z\right)\right) \cdot -0.16666666666666666\\
\mathbf{elif}\;z \leq 1.7 \cdot 10^{+39}:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot \left(x \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if z < -1.50000000000000003e88Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6490.9
Simplified90.9%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.5
Simplified87.5%
Taylor expanded in z around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.5
Simplified87.5%
if -1.50000000000000003e88 < z < 1.6999999999999999e39Initial program 100.0%
Taylor expanded in x around inf
Simplified61.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6438.7
Simplified38.7%
if 1.6999999999999999e39 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified35.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6411.8
Simplified11.8%
Taylor expanded in x around inf
cube-multN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6441.8
Simplified41.8%
Final simplification49.8%
(FPCore (x y z) :precision binary64 (if (<= z -1.5e+88) (* (* z (* z z)) -0.16666666666666666) (if (<= z 390.0) (fma x (fma x 0.5 1.0) 1.0) (* (* x x) 0.5))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.5e+88) {
tmp = (z * (z * z)) * -0.16666666666666666;
} else if (z <= 390.0) {
tmp = fma(x, fma(x, 0.5, 1.0), 1.0);
} else {
tmp = (x * x) * 0.5;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -1.5e+88) tmp = Float64(Float64(z * Float64(z * z)) * -0.16666666666666666); elseif (z <= 390.0) tmp = fma(x, fma(x, 0.5, 1.0), 1.0); else tmp = Float64(Float64(x * x) * 0.5); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -1.5e+88], N[(N[(z * N[(z * z), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision], If[LessEqual[z, 390.0], N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+88}:\\
\;\;\;\;\left(z \cdot \left(z \cdot z\right)\right) \cdot -0.16666666666666666\\
\mathbf{elif}\;z \leq 390:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.5\\
\end{array}
\end{array}
if z < -1.50000000000000003e88Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6490.9
Simplified90.9%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6487.5
Simplified87.5%
Taylor expanded in z around inf
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.5
Simplified87.5%
if -1.50000000000000003e88 < z < 390Initial program 100.0%
Taylor expanded in x around inf
Simplified62.2%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6438.4
Simplified38.4%
if 390 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.9
Simplified12.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.0
Simplified35.0%
Final simplification47.8%
(FPCore (x y z) :precision binary64 (if (<= z -2e+137) (* (* z z) 0.5) (if (<= z 700000000.0) (fma x (fma x 0.5 1.0) 1.0) (* (* x x) 0.5))))
double code(double x, double y, double z) {
double tmp;
if (z <= -2e+137) {
tmp = (z * z) * 0.5;
} else if (z <= 700000000.0) {
tmp = fma(x, fma(x, 0.5, 1.0), 1.0);
} else {
tmp = (x * x) * 0.5;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -2e+137) tmp = Float64(Float64(z * z) * 0.5); elseif (z <= 700000000.0) tmp = fma(x, fma(x, 0.5, 1.0), 1.0); else tmp = Float64(Float64(x * x) * 0.5); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -2e+137], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[z, 700000000.0], N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2 \cdot 10^{+137}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;z \leq 700000000:\\
\;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.5\\
\end{array}
\end{array}
if z < -2.0000000000000001e137Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.7
Simplified89.7%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f6482.1
Simplified82.1%
Taylor expanded in z around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.1
Simplified82.1%
if -2.0000000000000001e137 < z < 7e8Initial program 100.0%
Taylor expanded in x around inf
Simplified61.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6437.6
Simplified37.6%
if 7e8 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.9
Simplified12.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.0
Simplified35.0%
Final simplification45.2%
(FPCore (x y z) :precision binary64 (if (<= z -2e+137) (* (* z z) 0.5) (if (<= z 1200.0) (fma x (* x 0.5) 1.0) (* (* x x) 0.5))))
double code(double x, double y, double z) {
double tmp;
if (z <= -2e+137) {
tmp = (z * z) * 0.5;
} else if (z <= 1200.0) {
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 <= -2e+137) tmp = Float64(Float64(z * z) * 0.5); elseif (z <= 1200.0) tmp = fma(x, Float64(x * 0.5), 1.0); else tmp = Float64(Float64(x * x) * 0.5); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -2e+137], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[z, 1200.0], N[(x * N[(x * 0.5), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2 \cdot 10^{+137}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{elif}\;z \leq 1200:\\
\;\;\;\;\mathsf{fma}\left(x, x \cdot 0.5, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.5\\
\end{array}
\end{array}
if z < -2.0000000000000001e137Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.7
Simplified89.7%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f6482.1
Simplified82.1%
Taylor expanded in z around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.1
Simplified82.1%
if -2.0000000000000001e137 < z < 1200Initial program 100.0%
Taylor expanded in x around inf
Simplified61.1%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6437.6
Simplified37.6%
Taylor expanded in x around inf
*-lowering-*.f6437.5
Simplified37.5%
if 1200 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified34.7%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6412.9
Simplified12.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.0
Simplified35.0%
Final simplification45.2%
(FPCore (x y z) :precision binary64 (if (<= z -2e+137) (* (* z z) 0.5) (* (* x x) 0.5)))
double code(double x, double y, double z) {
double tmp;
if (z <= -2e+137) {
tmp = (z * z) * 0.5;
} 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) :: tmp
if (z <= (-2d+137)) then
tmp = (z * z) * 0.5d0
else
tmp = (x * x) * 0.5d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -2e+137) {
tmp = (z * z) * 0.5;
} else {
tmp = (x * x) * 0.5;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -2e+137: tmp = (z * z) * 0.5 else: tmp = (x * x) * 0.5 return tmp
function code(x, y, z) tmp = 0.0 if (z <= -2e+137) tmp = Float64(Float64(z * z) * 0.5); else tmp = Float64(Float64(x * x) * 0.5); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -2e+137) tmp = (z * z) * 0.5; else tmp = (x * x) * 0.5; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -2e+137], N[(N[(z * z), $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2 \cdot 10^{+137}:\\
\;\;\;\;\left(z \cdot z\right) \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.5\\
\end{array}
\end{array}
if z < -2.0000000000000001e137Initial program 100.0%
Taylor expanded in z around inf
mul-1-negN/A
neg-lowering-neg.f6489.7
Simplified89.7%
Taylor expanded in z around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f6482.1
Simplified82.1%
Taylor expanded in z around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.1
Simplified82.1%
if -2.0000000000000001e137 < z Initial program 100.0%
Taylor expanded in x around inf
Simplified52.6%
Taylor expanded in x around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6429.6
Simplified29.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6426.5
Simplified26.5%
Final simplification36.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 x around inf
Simplified50.1%
Taylor expanded in x around 0
+-lowering-+.f6411.6
Simplified11.6%
Final simplification11.6%
(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 x around inf
Simplified50.1%
Taylor expanded in x around 0
Simplified11.3%
(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 2024199
(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)))