
(FPCore (x y z)
:precision binary64
(+
(+ (- (* (- x 0.5) (log x)) x) 0.91893853320467)
(/
(+
(* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z)
0.083333333333333)
x)))
double code(double x, double y, double z) {
return ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x - 0.5d0) * log(x)) - x) + 0.91893853320467d0) + ((((((y + 0.0007936500793651d0) * z) - 0.0027777777777778d0) * z) + 0.083333333333333d0) / x)
end function
public static double code(double x, double y, double z) {
return ((((x - 0.5) * Math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
def code(x, y, z): return ((((x - 0.5) * math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x - 0.5) * log(x)) - x) + 0.91893853320467) + Float64(Float64(Float64(Float64(Float64(Float64(y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)) end
function tmp = code(x, y, z) tmp = ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x); end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x - 0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision] + 0.91893853320467), $MachinePrecision] + N[(N[(N[(N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * z), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] * z), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.91893853320467\right) + \frac{\left(\left(y + 0.0007936500793651\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.083333333333333}{x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z)
:precision binary64
(+
(+ (- (* (- x 0.5) (log x)) x) 0.91893853320467)
(/
(+
(* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z)
0.083333333333333)
x)))
double code(double x, double y, double z) {
return ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x - 0.5d0) * log(x)) - x) + 0.91893853320467d0) + ((((((y + 0.0007936500793651d0) * z) - 0.0027777777777778d0) * z) + 0.083333333333333d0) / x)
end function
public static double code(double x, double y, double z) {
return ((((x - 0.5) * Math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
def code(x, y, z): return ((((x - 0.5) * math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x - 0.5) * log(x)) - x) + 0.91893853320467) + Float64(Float64(Float64(Float64(Float64(Float64(y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)) end
function tmp = code(x, y, z) tmp = ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x); end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x - 0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision] + 0.91893853320467), $MachinePrecision] + N[(N[(N[(N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * z), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] * z), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.91893853320467\right) + \frac{\left(\left(y + 0.0007936500793651\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.083333333333333}{x}
\end{array}
(FPCore (x y z)
:precision binary64
(if (<= x 5e+36)
(fma
(fma
z
(fma (+ y 0.0007936500793651) z -0.0027777777777778)
0.083333333333333)
(/ 1.0 x)
(fma (+ x -0.5) (log x) (- 0.91893853320467 x)))
(+
(+ 0.91893853320467 (- (* (log x) (- x 0.5)) x))
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 5e+36) {
tmp = fma(fma(z, fma((y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), (1.0 / x), fma((x + -0.5), log(x), (0.91893853320467 - x)));
} else {
tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + (z * ((z * (y + 0.0007936500793651)) / x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 5e+36) tmp = fma(fma(z, fma(Float64(y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), Float64(1.0 / x), fma(Float64(x + -0.5), log(x), Float64(0.91893853320467 - x))); else tmp = Float64(Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x)) + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) / x))); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 5e+36], N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * z + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] * N[(1.0 / x), $MachinePrecision] + N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision] + N[(0.91893853320467 - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision] + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5 \cdot 10^{+36}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, \mathsf{fma}\left(y + 0.0007936500793651, z, -0.0027777777777778\right), 0.083333333333333\right), \frac{1}{x}, \mathsf{fma}\left(x + -0.5, \log x, 0.91893853320467 - x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right) + z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if x < 4.99999999999999977e36Initial program 99.6%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6499.7
Applied egg-rr99.7%
if 4.99999999999999977e36 < x Initial program 86.8%
Taylor expanded in z around inf
associate-/l*N/A
unpow2N/A
associate-*l*N/A
associate-*r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6497.2
Simplified97.2%
Final simplification98.6%
(FPCore (x y z)
:precision binary64
(if (<=
(+
(+ 0.91893853320467 (- (* (log x) (- x 0.5)) x))
(/
(+
0.083333333333333
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))
x))
-5e+55)
(* (* z y) (/ z x))
(/
(fma z (fma z 0.0007936500793651 -0.0027777777777778) 0.083333333333333)
x)))
double code(double x, double y, double z) {
double tmp;
if (((0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + ((0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x)) <= -5e+55) {
tmp = (z * y) * (z / x);
} else {
tmp = fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (Float64(Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x)) + Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) / x)) <= -5e+55) tmp = Float64(Float64(z * y) * Float64(z / x)); else tmp = Float64(fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x); end return tmp end
code[x_, y_, z_] := If[LessEqual[N[(N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision] + N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], -5e+55], N[(N[(z * y), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision], N[(N[(z * N[(z * 0.0007936500793651 + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right) + \frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)}{x} \leq -5 \cdot 10^{+55}:\\
\;\;\;\;\left(z \cdot y\right) \cdot \frac{z}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\end{array}
\end{array}
if (+.f64 (+.f64 (-.f64 (*.f64 (-.f64 x #s(literal 1/2 binary64)) (log.f64 x)) x) #s(literal 91893853320467/100000000000000 binary64)) (/.f64 (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) x)) < -5.00000000000000046e55Initial program 83.9%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.0
Simplified84.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6489.3
Applied egg-rr89.3%
if -5.00000000000000046e55 < (+.f64 (+.f64 (-.f64 (*.f64 (-.f64 x #s(literal 1/2 binary64)) (log.f64 x)) x) #s(literal 91893853320467/100000000000000 binary64)) (/.f64 (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) x)) Initial program 95.9%
Taylor expanded in y around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified93.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6454.2
Simplified54.2%
Final simplification59.1%
(FPCore (x y z)
:precision binary64
(if (<= x 1.8e-19)
(fma
(fma
z
(fma (+ y 0.0007936500793651) z -0.0027777777777778)
0.083333333333333)
(/ 1.0 x)
(* x (+ (log x) -1.0)))
(+
(+ 0.91893853320467 (- (* (log x) (- x 0.5)) x))
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 1.8e-19) {
tmp = fma(fma(z, fma((y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), (1.0 / x), (x * (log(x) + -1.0)));
} else {
tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + (z * ((z * (y + 0.0007936500793651)) / x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 1.8e-19) tmp = fma(fma(z, fma(Float64(y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), Float64(1.0 / x), Float64(x * Float64(log(x) + -1.0))); else tmp = Float64(Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x)) + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) / x))); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 1.8e-19], N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * z + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] * N[(1.0 / x), $MachinePrecision] + N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision] + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.8 \cdot 10^{-19}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, \mathsf{fma}\left(y + 0.0007936500793651, z, -0.0027777777777778\right), 0.083333333333333\right), \frac{1}{x}, x \cdot \left(\log x + -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right) + z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if x < 1.8000000000000001e-19Initial program 99.6%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6499.7
Simplified99.7%
if 1.8000000000000001e-19 < x Initial program 88.4%
Taylor expanded in z around inf
associate-/l*N/A
unpow2N/A
associate-*l*N/A
associate-*r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6497.4
Simplified97.4%
Final simplification98.6%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ (log x) -1.0)))
(if (<= x 100000000000.0)
(fma
(fma
z
(fma (+ y 0.0007936500793651) z -0.0027777777777778)
0.083333333333333)
(/ 1.0 x)
(* x t_0))
(fma t_0 x (* (+ y 0.0007936500793651) (/ (* z z) x))))))
double code(double x, double y, double z) {
double t_0 = log(x) + -1.0;
double tmp;
if (x <= 100000000000.0) {
tmp = fma(fma(z, fma((y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), (1.0 / x), (x * t_0));
} else {
tmp = fma(t_0, x, ((y + 0.0007936500793651) * ((z * z) / x)));
}
return tmp;
}
function code(x, y, z) t_0 = Float64(log(x) + -1.0) tmp = 0.0 if (x <= 100000000000.0) tmp = fma(fma(z, fma(Float64(y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333), Float64(1.0 / x), Float64(x * t_0)); else tmp = fma(t_0, x, Float64(Float64(y + 0.0007936500793651) * Float64(Float64(z * z) / x))); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]}, If[LessEqual[x, 100000000000.0], N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * z + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] * N[(1.0 / x), $MachinePrecision] + N[(x * t$95$0), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * x + N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(N[(z * z), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log x + -1\\
\mathbf{if}\;x \leq 100000000000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, \mathsf{fma}\left(y + 0.0007936500793651, z, -0.0027777777777778\right), 0.083333333333333\right), \frac{1}{x}, x \cdot t\_0\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_0, x, \left(y + 0.0007936500793651\right) \cdot \frac{z \cdot z}{x}\right)\\
\end{array}
\end{array}
if x < 1e11Initial program 99.6%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6499.5
Simplified99.5%
if 1e11 < x Initial program 87.7%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6487.8
Applied egg-rr87.8%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6487.8
Simplified87.8%
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6487.8
Applied egg-rr87.8%
Taylor expanded in z around inf
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6493.4
Simplified93.4%
Final simplification96.7%
(FPCore (x y z)
:precision binary64
(if (<= x 3e-6)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(fma (+ (log x) -1.0) x (* (+ y 0.0007936500793651) (/ (* z z) x)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 3e-6) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = fma((log(x) + -1.0), x, ((y + 0.0007936500793651) * ((z * z) / x)));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 3e-6) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = fma(Float64(log(x) + -1.0), x, Float64(Float64(y + 0.0007936500793651) * Float64(Float64(z * z) / x))); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 3e-6], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision] * x + N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(N[(z * z), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3 \cdot 10^{-6}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\log x + -1, x, \left(y + 0.0007936500793651\right) \cdot \frac{z \cdot z}{x}\right)\\
\end{array}
\end{array}
if x < 3.0000000000000001e-6Initial program 99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6499.7
Simplified99.7%
if 3.0000000000000001e-6 < x Initial program 88.2%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6488.3
Applied egg-rr88.3%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6488.1
Simplified88.1%
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6488.0
Applied egg-rr88.0%
Taylor expanded in z around inf
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6493.5
Simplified93.5%
Final simplification96.7%
(FPCore (x y z)
:precision binary64
(if (<= x 114.0)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(fma
(+ (log x) -1.0)
x
(/
(fma z (fma z 0.0007936500793651 -0.0027777777777778) 0.083333333333333)
x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 114.0) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = fma((log(x) + -1.0), x, (fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 114.0) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = fma(Float64(log(x) + -1.0), x, Float64(fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 114.0], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision] * x + N[(N[(z * N[(z * 0.0007936500793651 + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 114:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\log x + -1, x, \frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\right)\\
\end{array}
\end{array}
if x < 114Initial program 99.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6499.6
Simplified99.6%
if 114 < x Initial program 88.0%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6488.1
Applied egg-rr88.1%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6487.9
Simplified87.9%
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6487.9
Applied egg-rr87.9%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6486.9
Simplified86.9%
Final simplification93.7%
(FPCore (x y z)
:precision binary64
(if (<= x 4.7e+27)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(fma (+ (log x) -1.0) x (/ (* y (* z z)) x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 4.7e+27) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = fma((log(x) + -1.0), x, ((y * (z * z)) / x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 4.7e+27) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = fma(Float64(log(x) + -1.0), x, Float64(Float64(y * Float64(z * z)) / x)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 4.7e+27], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision] * x + N[(N[(y * N[(z * z), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.7 \cdot 10^{+27}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\log x + -1, x, \frac{y \cdot \left(z \cdot z\right)}{x}\right)\\
\end{array}
\end{array}
if x < 4.69999999999999976e27Initial program 99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6498.4
Simplified98.4%
if 4.69999999999999976e27 < x Initial program 87.3%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6487.3
Applied egg-rr87.3%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6487.3
Simplified87.3%
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6487.3
Applied egg-rr87.3%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.7
Simplified84.7%
Final simplification92.3%
(FPCore (x y z)
:precision binary64
(if (<= x 1.7e+35)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(* x (+ (log x) -1.0))))
double code(double x, double y, double z) {
double tmp;
if (x <= 1.7e+35) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = x * (log(x) + -1.0);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 1.7e+35) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = Float64(x * Float64(log(x) + -1.0)); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 1.7e+35], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.7 \cdot 10^{+35}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 1.7000000000000001e35Initial program 99.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6497.7
Simplified97.7%
if 1.7000000000000001e35 < x Initial program 87.0%
Taylor expanded in x around inf
mul-1-negN/A
log-recN/A
remove-double-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
log-lowering-log.f6480.9
Simplified80.9%
Final simplification90.4%
(FPCore (x y z)
:precision binary64
(let* ((t_0
(+
0.083333333333333
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))))
(if (<= t_0 -5e+36)
(* (* z y) (/ z x))
(if (<= t_0 2e+82)
(/
(fma
z
(fma z 0.0007936500793651 -0.0027777777777778)
0.083333333333333)
x)
(* (+ y 0.0007936500793651) (/ (* z z) x))))))
double code(double x, double y, double z) {
double t_0 = 0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778));
double tmp;
if (t_0 <= -5e+36) {
tmp = (z * y) * (z / x);
} else if (t_0 <= 2e+82) {
tmp = fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = (y + 0.0007936500793651) * ((z * z) / x);
}
return tmp;
}
function code(x, y, z) t_0 = Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) tmp = 0.0 if (t_0 <= -5e+36) tmp = Float64(Float64(z * y) * Float64(z / x)); elseif (t_0 <= 2e+82) tmp = Float64(fma(z, fma(z, 0.0007936500793651, -0.0027777777777778), 0.083333333333333) / x); else tmp = Float64(Float64(y + 0.0007936500793651) * Float64(Float64(z * z) / x)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -5e+36], N[(N[(z * y), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+82], N[(N[(z * N[(z * 0.0007936500793651 + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(N[(z * z), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{+36}:\\
\;\;\;\;\left(z \cdot y\right) \cdot \frac{z}{x}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+82}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(y + 0.0007936500793651\right) \cdot \frac{z \cdot z}{x}\\
\end{array}
\end{array}
if (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) < -4.99999999999999977e36Initial program 86.8%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.0
Simplified69.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6473.2
Applied egg-rr73.2%
if -4.99999999999999977e36 < (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) < 1.9999999999999999e82Initial program 99.4%
Taylor expanded in y around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified98.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
accelerator-lowering-fma.f6451.2
Simplified51.2%
if 1.9999999999999999e82 < (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) Initial program 90.3%
+-commutativeN/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
accelerator-lowering-fma.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f6490.3
Applied egg-rr90.3%
Taylor expanded in x around inf
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6490.3
Simplified90.3%
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6490.3
Applied egg-rr90.3%
Taylor expanded in z around inf
distribute-rgt-inN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f6473.6
Simplified73.6%
Final simplification62.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0
(+
0.083333333333333
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))))
(if (<= t_0 -5e+36)
(* (* z y) (/ z x))
(if (<= t_0 5e+23)
(fma (/ 1.0 x) 0.083333333333333 0.91893853320467)
(* z (* 0.0007936500793651 (/ z x)))))))
double code(double x, double y, double z) {
double t_0 = 0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778));
double tmp;
if (t_0 <= -5e+36) {
tmp = (z * y) * (z / x);
} else if (t_0 <= 5e+23) {
tmp = fma((1.0 / x), 0.083333333333333, 0.91893853320467);
} else {
tmp = z * (0.0007936500793651 * (z / x));
}
return tmp;
}
function code(x, y, z) t_0 = Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) tmp = 0.0 if (t_0 <= -5e+36) tmp = Float64(Float64(z * y) * Float64(z / x)); elseif (t_0 <= 5e+23) tmp = fma(Float64(1.0 / x), 0.083333333333333, 0.91893853320467); else tmp = Float64(z * Float64(0.0007936500793651 * Float64(z / x))); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -5e+36], N[(N[(z * y), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 5e+23], N[(N[(1.0 / x), $MachinePrecision] * 0.083333333333333 + 0.91893853320467), $MachinePrecision], N[(z * N[(0.0007936500793651 * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{+36}:\\
\;\;\;\;\left(z \cdot y\right) \cdot \frac{z}{x}\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{+23}:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{x}, 0.083333333333333, 0.91893853320467\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(0.0007936500793651 \cdot \frac{z}{x}\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) < -4.99999999999999977e36Initial program 86.8%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.0
Simplified69.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6473.2
Applied egg-rr73.2%
if -4.99999999999999977e36 < (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) < 4.9999999999999999e23Initial program 99.4%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6497.9
Simplified97.9%
Taylor expanded in x around 0
/-lowering-/.f6450.1
Simplified50.1%
+-commutativeN/A
clear-numN/A
associate-/r/N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6450.1
Applied egg-rr50.1%
if 4.9999999999999999e23 < (+.f64 (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) #s(literal 83333333333333/1000000000000000 binary64)) Initial program 90.8%
Taylor expanded in y around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.3%
Taylor expanded in z around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Simplified81.4%
Taylor expanded in z around inf
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6464.8
Simplified64.8%
Final simplification59.3%
(FPCore (x y z) :precision binary64 (if (<= (* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)) 5e+23) (fma (/ 1.0 x) 0.083333333333333 0.91893853320467) (* z (* 0.0007936500793651 (/ z x)))))
double code(double x, double y, double z) {
double tmp;
if ((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) <= 5e+23) {
tmp = fma((1.0 / x), 0.083333333333333, 0.91893853320467);
} else {
tmp = z * (0.0007936500793651 * (z / x));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778)) <= 5e+23) tmp = fma(Float64(1.0 / x), 0.083333333333333, 0.91893853320467); else tmp = Float64(z * Float64(0.0007936500793651 * Float64(z / x))); end return tmp end
code[x_, y_, z_] := If[LessEqual[N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision], 5e+23], N[(N[(1.0 / x), $MachinePrecision] * 0.083333333333333 + 0.91893853320467), $MachinePrecision], N[(z * N[(0.0007936500793651 * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right) \leq 5 \cdot 10^{+23}:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{x}, 0.083333333333333, 0.91893853320467\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(0.0007936500793651 \cdot \frac{z}{x}\right)\\
\end{array}
\end{array}
if (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) < 4.9999999999999999e23Initial program 96.0%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6478.5
Simplified78.5%
Taylor expanded in x around 0
/-lowering-/.f6437.1
Simplified37.1%
+-commutativeN/A
clear-numN/A
associate-/r/N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6437.1
Applied egg-rr37.1%
if 4.9999999999999999e23 < (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) Initial program 90.8%
Taylor expanded in y around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.3%
Taylor expanded in z around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Simplified81.4%
Taylor expanded in z around inf
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6464.8
Simplified64.8%
Final simplification47.0%
(FPCore (x y z)
:precision binary64
(if (<= x 7e+103)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(* z (* z (+ (/ 0.0007936500793651 x) (/ y x))))))
double code(double x, double y, double z) {
double tmp;
if (x <= 7e+103) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = z * (z * ((0.0007936500793651 / x) + (y / x)));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 7e+103) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = Float64(z * Float64(z * Float64(Float64(0.0007936500793651 / x) + Float64(y / x)))); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 7e+103], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(z * N[(N[(0.0007936500793651 / x), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 7 \cdot 10^{+103}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(z \cdot \left(\frac{0.0007936500793651}{x} + \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if x < 7e103Initial program 99.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6489.9
Simplified89.9%
if 7e103 < x Initial program 84.1%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f6417.8
Simplified17.8%
Final simplification64.5%
(FPCore (x y z) :precision binary64 (/ (fma z (fma z (+ y 0.0007936500793651) -0.0027777777777778) 0.083333333333333) x))
double code(double x, double y, double z) {
return fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
}
function code(x, y, z) return Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x) end
code[x_, y_, z_] := N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}
\end{array}
Initial program 94.2%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
+-lowering-+.f6461.8
Simplified61.8%
Final simplification61.8%
(FPCore (x y z) :precision binary64 (fma (/ 1.0 x) 0.083333333333333 0.91893853320467))
double code(double x, double y, double z) {
return fma((1.0 / x), 0.083333333333333, 0.91893853320467);
}
function code(x, y, z) return fma(Float64(1.0 / x), 0.083333333333333, 0.91893853320467) end
code[x_, y_, z_] := N[(N[(1.0 / x), $MachinePrecision] * 0.083333333333333 + 0.91893853320467), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\frac{1}{x}, 0.083333333333333, 0.91893853320467\right)
\end{array}
Initial program 94.2%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6460.0
Simplified60.0%
Taylor expanded in x around 0
/-lowering-/.f6425.6
Simplified25.6%
+-commutativeN/A
clear-numN/A
associate-/r/N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6425.6
Applied egg-rr25.6%
(FPCore (x y z) :precision binary64 (+ 0.91893853320467 (/ 0.083333333333333 x)))
double code(double x, double y, double z) {
return 0.91893853320467 + (0.083333333333333 / x);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.91893853320467d0 + (0.083333333333333d0 / x)
end function
public static double code(double x, double y, double z) {
return 0.91893853320467 + (0.083333333333333 / x);
}
def code(x, y, z): return 0.91893853320467 + (0.083333333333333 / x)
function code(x, y, z) return Float64(0.91893853320467 + Float64(0.083333333333333 / x)) end
function tmp = code(x, y, z) tmp = 0.91893853320467 + (0.083333333333333 / x); end
code[x_, y_, z_] := N[(0.91893853320467 + N[(0.083333333333333 / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.91893853320467 + \frac{0.083333333333333}{x}
\end{array}
Initial program 94.2%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6460.0
Simplified60.0%
Taylor expanded in x around 0
/-lowering-/.f6425.6
Simplified25.6%
(FPCore (x y z) :precision binary64 (/ 0.083333333333333 x))
double code(double x, double y, double z) {
return 0.083333333333333 / x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.083333333333333d0 / x
end function
public static double code(double x, double y, double z) {
return 0.083333333333333 / x;
}
def code(x, y, z): return 0.083333333333333 / x
function code(x, y, z) return Float64(0.083333333333333 / x) end
function tmp = code(x, y, z) tmp = 0.083333333333333 / x; end
code[x_, y_, z_] := N[(0.083333333333333 / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333}{x}
\end{array}
Initial program 94.2%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6460.0
Simplified60.0%
Taylor expanded in x around 0
/-lowering-/.f6425.0
Simplified25.0%
(FPCore (x y z) :precision binary64 0.91893853320467)
double code(double x, double y, double z) {
return 0.91893853320467;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.91893853320467d0
end function
public static double code(double x, double y, double z) {
return 0.91893853320467;
}
def code(x, y, z): return 0.91893853320467
function code(x, y, z) return 0.91893853320467 end
function tmp = code(x, y, z) tmp = 0.91893853320467; end
code[x_, y_, z_] := 0.91893853320467
\begin{array}{l}
\\
0.91893853320467
\end{array}
Initial program 94.2%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6460.0
Simplified60.0%
Taylor expanded in x around 0
/-lowering-/.f6425.6
Simplified25.6%
Taylor expanded in x around inf
Simplified3.8%
(FPCore (x y z) :precision binary64 (+ (+ (+ (* (- x 0.5) (log x)) (- 0.91893853320467 x)) (/ 0.083333333333333 x)) (* (/ z x) (- (* z (+ y 0.0007936500793651)) 0.0027777777777778))))
double code(double x, double y, double z) {
return ((((x - 0.5) * log(x)) + (0.91893853320467 - x)) + (0.083333333333333 / x)) + ((z / x) * ((z * (y + 0.0007936500793651)) - 0.0027777777777778));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x - 0.5d0) * log(x)) + (0.91893853320467d0 - x)) + (0.083333333333333d0 / x)) + ((z / x) * ((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0))
end function
public static double code(double x, double y, double z) {
return ((((x - 0.5) * Math.log(x)) + (0.91893853320467 - x)) + (0.083333333333333 / x)) + ((z / x) * ((z * (y + 0.0007936500793651)) - 0.0027777777777778));
}
def code(x, y, z): return ((((x - 0.5) * math.log(x)) + (0.91893853320467 - x)) + (0.083333333333333 / x)) + ((z / x) * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x - 0.5) * log(x)) + Float64(0.91893853320467 - x)) + Float64(0.083333333333333 / x)) + Float64(Float64(z / x) * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) end
function tmp = code(x, y, z) tmp = ((((x - 0.5) * log(x)) + (0.91893853320467 - x)) + (0.083333333333333 / x)) + ((z / x) * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)); end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x - 0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(0.91893853320467 - x), $MachinePrecision]), $MachinePrecision] + N[(0.083333333333333 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(z / x), $MachinePrecision] * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x - 0.5\right) \cdot \log x + \left(0.91893853320467 - x\right)\right) + \frac{0.083333333333333}{x}\right) + \frac{z}{x} \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)
\end{array}
herbie shell --seed 2024195
(FPCore (x y z)
:name "Numeric.SpecFunctions:$slogFactorial from math-functions-0.1.5.2, B"
:precision binary64
:alt
(! :herbie-platform default (+ (+ (+ (* (- x 1/2) (log x)) (- 91893853320467/100000000000000 x)) (/ 83333333333333/1000000000000000 x)) (* (/ z x) (- (* z (+ y 7936500793651/10000000000000000)) 13888888888889/5000000000000000))))
(+ (+ (- (* (- x 0.5) (log x)) x) 0.91893853320467) (/ (+ (* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z) 0.083333333333333) x)))