
(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 21 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
(let* ((t_0 (* (+ x -0.5) (log x))))
(if (<= x 1500000.0)
(+
(- (+ 0.91893853320467 t_0) x)
(/
(+
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778))
0.083333333333333)
x))
(+
t_0
(-
(* (/ z (/ x z)) (+ y 0.0007936500793651))
(+ x -0.91893853320467))))))
double code(double x, double y, double z) {
double t_0 = (x + -0.5) * log(x);
double tmp;
if (x <= 1500000.0) {
tmp = ((0.91893853320467 + t_0) - x) + (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x);
} else {
tmp = t_0 + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
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 + (-0.5d0)) * log(x)
if (x <= 1500000.0d0) then
tmp = ((0.91893853320467d0 + t_0) - x) + (((z * ((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0)) + 0.083333333333333d0) / x)
else
tmp = t_0 + (((z / (x / z)) * (y + 0.0007936500793651d0)) - (x + (-0.91893853320467d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x + -0.5) * Math.log(x);
double tmp;
if (x <= 1500000.0) {
tmp = ((0.91893853320467 + t_0) - x) + (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x);
} else {
tmp = t_0 + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
return tmp;
}
def code(x, y, z): t_0 = (x + -0.5) * math.log(x) tmp = 0 if x <= 1500000.0: tmp = ((0.91893853320467 + t_0) - x) + (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) else: tmp = t_0 + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)) return tmp
function code(x, y, z) t_0 = Float64(Float64(x + -0.5) * log(x)) tmp = 0.0 if (x <= 1500000.0) tmp = Float64(Float64(Float64(0.91893853320467 + t_0) - x) + Float64(Float64(Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x)); else tmp = Float64(t_0 + Float64(Float64(Float64(z / Float64(x / z)) * Float64(y + 0.0007936500793651)) - Float64(x + -0.91893853320467))); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + -0.5) * log(x); tmp = 0.0; if (x <= 1500000.0) tmp = ((0.91893853320467 + t_0) - x) + (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x); else tmp = t_0 + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1500000.0], N[(N[(N[(0.91893853320467 + t$95$0), $MachinePrecision] - x), $MachinePrecision] + N[(N[(N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], N[(t$95$0 + N[(N[(N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - N[(x + -0.91893853320467), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + -0.5\right) \cdot \log x\\
\mathbf{if}\;x \leq 1500000:\\
\;\;\;\;\left(\left(0.91893853320467 + t\_0\right) - x\right) + \frac{z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right) + 0.083333333333333}{x}\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(\frac{z}{\frac{x}{z}} \cdot \left(y + 0.0007936500793651\right) - \left(x + -0.91893853320467\right)\right)\\
\end{array}
\end{array}
if x < 1.5e6Initial program 99.7%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f6499.8%
Applied egg-rr99.8%
if 1.5e6 < x Initial program 84.2%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr84.2%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.6%
Simplified99.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6499.6%
Simplified99.6%
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6499.6%
Applied egg-rr99.6%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(if (<= x 1500000.0)
(+
(/
(+
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778))
0.083333333333333)
x)
(+ 0.91893853320467 (- (* (log x) (- x 0.5)) x)))
(+
(* (+ x -0.5) (log x))
(- (* (/ z (/ x z)) (+ y 0.0007936500793651)) (+ x -0.91893853320467)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 1500000.0) {
tmp = (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) + (0.91893853320467 + ((log(x) * (x - 0.5)) - x));
} else {
tmp = ((x + -0.5) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
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 (x <= 1500000.0d0) then
tmp = (((z * ((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0)) + 0.083333333333333d0) / x) + (0.91893853320467d0 + ((log(x) * (x - 0.5d0)) - x))
else
tmp = ((x + (-0.5d0)) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651d0)) - (x + (-0.91893853320467d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 1500000.0) {
tmp = (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) + (0.91893853320467 + ((Math.log(x) * (x - 0.5)) - x));
} else {
tmp = ((x + -0.5) * Math.log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 1500000.0: tmp = (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) + (0.91893853320467 + ((math.log(x) * (x - 0.5)) - x)) else: tmp = ((x + -0.5) * math.log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 1500000.0) tmp = Float64(Float64(Float64(Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) + Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x))); else tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(Float64(Float64(z / Float64(x / z)) * Float64(y + 0.0007936500793651)) - Float64(x + -0.91893853320467))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 1500000.0) tmp = (((z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778)) + 0.083333333333333) / x) + (0.91893853320467 + ((log(x) * (x - 0.5)) - x)); else tmp = ((x + -0.5) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 1500000.0], N[(N[(N[(N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision] + N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - N[(x + -0.91893853320467), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1500000:\\
\;\;\;\;\frac{z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right) + 0.083333333333333}{x} + \left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(\frac{z}{\frac{x}{z}} \cdot \left(y + 0.0007936500793651\right) - \left(x + -0.91893853320467\right)\right)\\
\end{array}
\end{array}
if x < 1.5e6Initial program 99.7%
if 1.5e6 < x Initial program 84.2%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr84.2%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.6%
Simplified99.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6499.6%
Simplified99.6%
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6499.6%
Applied egg-rr99.6%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(if (<= x 2550000.0)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(if (<= x 3.8e+274)
(+ (+ 0.91893853320467 (- (* (log x) (- x 0.5)) x)) (/ (* z (* y z)) x))
(* x (+ (log x) -1.0)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 2550000.0) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else if (x <= 3.8e+274) {
tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + ((z * (y * z)) / x);
} else {
tmp = x * (log(x) + -1.0);
}
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 (x <= 2550000.0d0) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else if (x <= 3.8d+274) then
tmp = (0.91893853320467d0 + ((log(x) * (x - 0.5d0)) - x)) + ((z * (y * z)) / x)
else
tmp = x * (log(x) + (-1.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 2550000.0) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else if (x <= 3.8e+274) {
tmp = (0.91893853320467 + ((Math.log(x) * (x - 0.5)) - x)) + ((z * (y * z)) / x);
} else {
tmp = x * (Math.log(x) + -1.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 2550000.0: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x elif x <= 3.8e+274: tmp = (0.91893853320467 + ((math.log(x) * (x - 0.5)) - x)) + ((z * (y * z)) / x) else: tmp = x * (math.log(x) + -1.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 2550000.0) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); elseif (x <= 3.8e+274) tmp = Float64(Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x)) + Float64(Float64(z * Float64(y * z)) / x)); else tmp = Float64(x * Float64(log(x) + -1.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 2550000.0) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; elseif (x <= 3.8e+274) tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + ((z * (y * z)) / x); else tmp = x * (log(x) + -1.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 2550000.0], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 3.8e+274], N[(N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(y * z), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2550000:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{elif}\;x \leq 3.8 \cdot 10^{+274}:\\
\;\;\;\;\left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right) + \frac{z \cdot \left(y \cdot z\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 2.55e6Initial program 99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6497.5%
Simplified97.5%
if 2.55e6 < x < 3.7999999999999998e274Initial program 88.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6483.0%
Simplified83.0%
if 3.7999999999999998e274 < x Initial program 56.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6497.6%
Simplified97.6%
Final simplification90.7%
(FPCore (x y z)
:precision binary64
(if (<= x 28000000.0)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(if (<= x 1.45e+274)
(- (* x (log x)) (- (+ x -0.91893853320467) (/ (* z (* y z)) x)))
(* x (+ (log x) -1.0)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 28000000.0) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else if (x <= 1.45e+274) {
tmp = (x * log(x)) - ((x + -0.91893853320467) - ((z * (y * z)) / x));
} else {
tmp = x * (log(x) + -1.0);
}
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 (x <= 28000000.0d0) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else if (x <= 1.45d+274) then
tmp = (x * log(x)) - ((x + (-0.91893853320467d0)) - ((z * (y * z)) / x))
else
tmp = x * (log(x) + (-1.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 28000000.0) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else if (x <= 1.45e+274) {
tmp = (x * Math.log(x)) - ((x + -0.91893853320467) - ((z * (y * z)) / x));
} else {
tmp = x * (Math.log(x) + -1.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 28000000.0: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x elif x <= 1.45e+274: tmp = (x * math.log(x)) - ((x + -0.91893853320467) - ((z * (y * z)) / x)) else: tmp = x * (math.log(x) + -1.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 28000000.0) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); elseif (x <= 1.45e+274) tmp = Float64(Float64(x * log(x)) - Float64(Float64(x + -0.91893853320467) - Float64(Float64(z * Float64(y * z)) / x))); else tmp = Float64(x * Float64(log(x) + -1.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 28000000.0) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; elseif (x <= 1.45e+274) tmp = (x * log(x)) - ((x + -0.91893853320467) - ((z * (y * z)) / x)); else tmp = x * (log(x) + -1.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 28000000.0], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 1.45e+274], N[(N[(x * N[Log[x], $MachinePrecision]), $MachinePrecision] - N[(N[(x + -0.91893853320467), $MachinePrecision] - N[(N[(z * N[(y * z), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 28000000:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{elif}\;x \leq 1.45 \cdot 10^{+274}:\\
\;\;\;\;x \cdot \log x - \left(\left(x + -0.91893853320467\right) - \frac{z \cdot \left(y \cdot z\right)}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 2.8e7Initial program 99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6497.5%
Simplified97.5%
if 2.8e7 < x < 1.45e274Initial program 88.0%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr88.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6483.0%
Simplified83.0%
Taylor expanded in x 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.5%
Simplified82.5%
if 1.45e274 < x Initial program 56.6%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6497.6%
Simplified97.6%
Final simplification90.4%
(FPCore (x y z)
:precision binary64
(if (<= x 0.00032)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(+
(* (+ x -0.5) (log x))
(- (* (/ z (/ x z)) (+ y 0.0007936500793651)) (+ x -0.91893853320467)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
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 (x <= 0.00032d0) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = ((x + (-0.5d0)) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651d0)) - (x + (-0.91893853320467d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * Math.log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 0.00032: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = ((x + -0.5) * math.log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 0.00032) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(Float64(Float64(z / Float64(x / z)) * Float64(y + 0.0007936500793651)) - Float64(x + -0.91893853320467))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 0.00032) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = ((x + -0.5) * log(x)) + (((z / (x / z)) * (y + 0.0007936500793651)) - (x + -0.91893853320467)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 0.00032], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - N[(x + -0.91893853320467), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.00032:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(\frac{z}{\frac{x}{z}} \cdot \left(y + 0.0007936500793651\right) - \left(x + -0.91893853320467\right)\right)\\
\end{array}
\end{array}
if x < 3.20000000000000026e-4Initial program 99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6499.1%
Simplified99.1%
if 3.20000000000000026e-4 < x Initial program 85.0%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.6%
Simplified99.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.8%
Simplified98.8%
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.8%
Applied egg-rr98.8%
Final simplification98.9%
(FPCore (x y z)
:precision binary64
(if (<= x 0.00032)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(+
(* (+ x -0.5) (log x))
(- (* z (* (+ y 0.0007936500793651) (/ z x))) (+ x -0.91893853320467)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467));
}
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 (x <= 0.00032d0) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = ((x + (-0.5d0)) * log(x)) + ((z * ((y + 0.0007936500793651d0) * (z / x))) - (x + (-0.91893853320467d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * Math.log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 0.00032: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = ((x + -0.5) * math.log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 0.00032) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))) - Float64(x + -0.91893853320467))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 0.00032) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = ((x + -0.5) * log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 0.00032], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x + -0.91893853320467), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.00032:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right) - \left(x + -0.91893853320467\right)\right)\\
\end{array}
\end{array}
if x < 3.20000000000000026e-4Initial program 99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6499.1%
Simplified99.1%
if 3.20000000000000026e-4 < x Initial program 85.0%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.6%
Simplified99.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.8%
Simplified98.8%
Final simplification98.9%
(FPCore (x y z)
:precision binary64
(if (<= x 0.00032)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(+
(* x (log x))
(- (* z (* (+ y 0.0007936500793651) (/ z x))) (+ x -0.91893853320467)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = (x * log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467));
}
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 (x <= 0.00032d0) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = (x * log(x)) + ((z * ((y + 0.0007936500793651d0) * (z / x))) - (x + (-0.91893853320467d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 0.00032) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = (x * Math.log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 0.00032: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = (x * math.log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 0.00032) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(Float64(x * log(x)) + Float64(Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))) - Float64(x + -0.91893853320467))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 0.00032) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = (x * log(x)) + ((z * ((y + 0.0007936500793651) * (z / x))) - (x + -0.91893853320467)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 0.00032], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x + -0.91893853320467), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.00032:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \log x + \left(z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right) - \left(x + -0.91893853320467\right)\right)\\
\end{array}
\end{array}
if x < 3.20000000000000026e-4Initial program 99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6499.1%
Simplified99.1%
if 3.20000000000000026e-4 < x Initial program 85.0%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.6%
Simplified99.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.8%
Simplified98.8%
Taylor expanded in x around inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6498.2%
Simplified98.2%
Final simplification98.6%
(FPCore (x y z)
:precision binary64
(if (<= x 1.05e+68)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(+ (* (+ x -0.5) (log x)) (- 0.91893853320467 x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 1.05e+68) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * log(x)) + (0.91893853320467 - x);
}
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 (x <= 1.05d+68) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = ((x + (-0.5d0)) * log(x)) + (0.91893853320467d0 - x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 1.05e+68) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = ((x + -0.5) * Math.log(x)) + (0.91893853320467 - x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 1.05e+68: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = ((x + -0.5) * math.log(x)) + (0.91893853320467 - x) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 1.05e+68) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(0.91893853320467 - x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 1.05e+68) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = ((x + -0.5) * log(x)) + (0.91893853320467 - x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 1.05e+68], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(0.91893853320467 - x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.05 \cdot 10^{+68}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(0.91893853320467 - x\right)\\
\end{array}
\end{array}
if x < 1.05e68Initial program 99.1%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6490.3%
Simplified90.3%
if 1.05e68 < x Initial program 81.1%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr81.1%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6482.5%
Simplified82.5%
Taylor expanded in z around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f6475.6%
Simplified75.6%
Final simplification83.9%
(FPCore (x y z)
:precision binary64
(if (<= x 6.2e+67)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(- (* x (log x)) x)))
double code(double x, double y, double z) {
double tmp;
if (x <= 6.2e+67) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = (x * log(x)) - x;
}
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 (x <= 6.2d+67) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = (x * log(x)) - x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 6.2e+67) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = (x * Math.log(x)) - x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 6.2e+67: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = (x * math.log(x)) - x return tmp
function code(x, y, z) tmp = 0.0 if (x <= 6.2e+67) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(Float64(x * log(x)) - x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 6.2e+67) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = (x * log(x)) - x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 6.2e+67], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6.2 \cdot 10^{+67}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \log x - x\\
\end{array}
\end{array}
if x < 6.19999999999999992e67Initial program 99.1%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6490.3%
Simplified90.3%
if 6.19999999999999992e67 < x Initial program 81.1%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6475.6%
Simplified75.6%
distribute-rgt-inN/A
neg-mul-1N/A
+-lowering-+.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6475.6%
Applied egg-rr75.6%
Final simplification83.9%
(FPCore (x y z)
:precision binary64
(if (<= x 4e+67)
(/
(+
0.083333333333333
(* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778)))
x)
(* x (+ (log x) -1.0))))
double code(double x, double y, double z) {
double tmp;
if (x <= 4e+67) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = x * (log(x) + -1.0);
}
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 (x <= 4d+67) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
else
tmp = x * (log(x) + (-1.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 4e+67) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
} else {
tmp = x * (Math.log(x) + -1.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 4e+67: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x else: tmp = x * (math.log(x) + -1.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 4e+67) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x); else tmp = Float64(x * Float64(log(x) + -1.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 4e+67) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; else tmp = x * (log(x) + -1.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 4e+67], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4 \cdot 10^{+67}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 3.99999999999999993e67Initial program 99.1%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6490.3%
Simplified90.3%
if 3.99999999999999993e67 < x Initial program 81.1%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f6475.6%
Simplified75.6%
Final simplification83.9%
(FPCore (x y z)
:precision binary64
(if (<= z -2050000.0)
(* z (* (+ y 0.0007936500793651) (/ z x)))
(if (<= z 11.0)
(/ (+ 0.083333333333333 (* z (* y z))) x)
(* z (* z (+ (/ 0.0007936500793651 x) (/ y x)))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -2050000.0) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 11.0) {
tmp = (0.083333333333333 + (z * (y * z))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 / x) + (y / x)));
}
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 <= (-2050000.0d0)) then
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
else if (z <= 11.0d0) then
tmp = (0.083333333333333d0 + (z * (y * z))) / x
else
tmp = z * (z * ((0.0007936500793651d0 / x) + (y / x)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -2050000.0) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 11.0) {
tmp = (0.083333333333333 + (z * (y * z))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 / x) + (y / x)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -2050000.0: tmp = z * ((y + 0.0007936500793651) * (z / x)) elif z <= 11.0: tmp = (0.083333333333333 + (z * (y * z))) / x else: tmp = z * (z * ((0.0007936500793651 / x) + (y / x))) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -2050000.0) tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); elseif (z <= 11.0) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(y * z))) / x); else tmp = Float64(z * Float64(z * Float64(Float64(0.0007936500793651 / x) + Float64(y / x)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -2050000.0) tmp = z * ((y + 0.0007936500793651) * (z / x)); elseif (z <= 11.0) tmp = (0.083333333333333 + (z * (y * z))) / x; else tmp = z * (z * ((0.0007936500793651 / x) + (y / x))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -2050000.0], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 11.0], N[(N[(0.083333333333333 + N[(z * N[(y * z), $MachinePrecision]), $MachinePrecision]), $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}\;z \leq -2050000:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\mathbf{elif}\;z \leq 11:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(y \cdot z\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(z \cdot \left(\frac{0.0007936500793651}{x} + \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if z < -2.05e6Initial program 81.1%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr81.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.8%
Simplified99.8%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6473.4%
Simplified73.4%
if -2.05e6 < z < 11Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6450.6%
Simplified50.6%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
if 11 < z Initial program 85.8%
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-/.f6472.2%
Simplified72.2%
Final simplification61.8%
(FPCore (x y z)
:precision binary64
(if (<= z -1850000.0)
(* z (* (+ y 0.0007936500793651) (/ z x)))
(if (<= z 9.5)
(/ (+ 0.083333333333333 (* z (* y z))) x)
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1850000.0) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 9.5) {
tmp = (0.083333333333333 + (z * (y * z))) / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
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 <= (-1850000.0d0)) then
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
else if (z <= 9.5d0) then
tmp = (0.083333333333333d0 + (z * (y * z))) / x
else
tmp = z * ((z * (y + 0.0007936500793651d0)) / x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -1850000.0) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 9.5) {
tmp = (0.083333333333333 + (z * (y * z))) / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -1850000.0: tmp = z * ((y + 0.0007936500793651) * (z / x)) elif z <= 9.5: tmp = (0.083333333333333 + (z * (y * z))) / x else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -1850000.0) tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); elseif (z <= 9.5) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(y * z))) / x); else tmp = Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) / x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -1850000.0) tmp = z * ((y + 0.0007936500793651) * (z / x)); elseif (z <= 9.5) tmp = (0.083333333333333 + (z * (y * z))) / x; else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -1850000.0], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 9.5], N[(N[(0.083333333333333 + N[(z * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1850000:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\mathbf{elif}\;z \leq 9.5:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(y \cdot z\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -1.85e6Initial program 81.1%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr81.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.8%
Simplified99.8%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6473.4%
Simplified73.4%
if -1.85e6 < z < 9.5Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6450.6%
Simplified50.6%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
if 9.5 < z Initial program 85.8%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.8%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.7%
Simplified99.7%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in x around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6472.2%
Simplified72.2%
Final simplification61.8%
(FPCore (x y z)
:precision binary64
(if (<= z -5.5e-80)
(* z (* (+ y 0.0007936500793651) (/ z x)))
(if (<= z 8.2e-44)
(/ 0.083333333333333 x)
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -5.5e-80) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 8.2e-44) {
tmp = 0.083333333333333 / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
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 <= (-5.5d-80)) then
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
else if (z <= 8.2d-44) then
tmp = 0.083333333333333d0 / x
else
tmp = z * ((z * (y + 0.0007936500793651d0)) / x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -5.5e-80) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else if (z <= 8.2e-44) {
tmp = 0.083333333333333 / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -5.5e-80: tmp = z * ((y + 0.0007936500793651) * (z / x)) elif z <= 8.2e-44: tmp = 0.083333333333333 / x else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -5.5e-80) tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); elseif (z <= 8.2e-44) tmp = Float64(0.083333333333333 / x); else tmp = Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) / x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -5.5e-80) tmp = z * ((y + 0.0007936500793651) * (z / x)); elseif (z <= 8.2e-44) tmp = 0.083333333333333 / x; else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -5.5e-80], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 8.2e-44], N[(0.083333333333333 / x), $MachinePrecision], N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.5 \cdot 10^{-80}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\mathbf{elif}\;z \leq 8.2 \cdot 10^{-44}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -5.4999999999999997e-80Initial program 85.1%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.1%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6498.6%
Simplified98.6%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6466.3%
Simplified66.3%
if -5.4999999999999997e-80 < z < 8.19999999999999984e-44Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6449.0%
Simplified49.0%
Taylor expanded in z around 0
/-lowering-/.f6446.5%
Simplified46.5%
if 8.19999999999999984e-44 < z Initial program 86.8%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr86.8%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.7%
Simplified99.7%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6496.8%
Simplified96.8%
Taylor expanded in x around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6470.1%
Simplified70.1%
Final simplification59.5%
(FPCore (x y z) :precision binary64 (let* ((t_0 (* z (* (+ y 0.0007936500793651) (/ z x))))) (if (<= z -4.8e-80) t_0 (if (<= z 3.8e-42) (/ 0.083333333333333 x) t_0))))
double code(double x, double y, double z) {
double t_0 = z * ((y + 0.0007936500793651) * (z / x));
double tmp;
if (z <= -4.8e-80) {
tmp = t_0;
} else if (z <= 3.8e-42) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = z * ((y + 0.0007936500793651d0) * (z / x))
if (z <= (-4.8d-80)) then
tmp = t_0
else if (z <= 3.8d-42) then
tmp = 0.083333333333333d0 / x
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = z * ((y + 0.0007936500793651) * (z / x));
double tmp;
if (z <= -4.8e-80) {
tmp = t_0;
} else if (z <= 3.8e-42) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = z * ((y + 0.0007936500793651) * (z / x)) tmp = 0 if z <= -4.8e-80: tmp = t_0 elif z <= 3.8e-42: tmp = 0.083333333333333 / x else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))) tmp = 0.0 if (z <= -4.8e-80) tmp = t_0; elseif (z <= 3.8e-42) tmp = Float64(0.083333333333333 / x); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * ((y + 0.0007936500793651) * (z / x)); tmp = 0.0; if (z <= -4.8e-80) tmp = t_0; elseif (z <= 3.8e-42) tmp = 0.083333333333333 / x; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -4.8e-80], t$95$0, If[LessEqual[z, 3.8e-42], N[(0.083333333333333 / x), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\mathbf{if}\;z \leq -4.8 \cdot 10^{-80}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{-42}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -4.7999999999999998e-80 or 3.80000000000000017e-42 < z Initial program 85.9%
associate-+l-N/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr85.9%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6499.1%
Simplified99.1%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6468.1%
Simplified68.1%
if -4.7999999999999998e-80 < z < 3.80000000000000017e-42Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6449.0%
Simplified49.0%
Taylor expanded in z around 0
/-lowering-/.f6446.5%
Simplified46.5%
Final simplification59.5%
(FPCore (x y z) :precision binary64 (let* ((t_0 (* y (/ z (/ x z))))) (if (<= z -3.7e-84) t_0 (if (<= z 3.6e-44) (/ 0.083333333333333 x) t_0))))
double code(double x, double y, double z) {
double t_0 = y * (z / (x / z));
double tmp;
if (z <= -3.7e-84) {
tmp = t_0;
} else if (z <= 3.6e-44) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
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 = y * (z / (x / z))
if (z <= (-3.7d-84)) then
tmp = t_0
else if (z <= 3.6d-44) then
tmp = 0.083333333333333d0 / x
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = y * (z / (x / z));
double tmp;
if (z <= -3.7e-84) {
tmp = t_0;
} else if (z <= 3.6e-44) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = y * (z / (x / z)) tmp = 0 if z <= -3.7e-84: tmp = t_0 elif z <= 3.6e-44: tmp = 0.083333333333333 / x else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(y * Float64(z / Float64(x / z))) tmp = 0.0 if (z <= -3.7e-84) tmp = t_0; elseif (z <= 3.6e-44) tmp = Float64(0.083333333333333 / x); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = y * (z / (x / z)); tmp = 0.0; if (z <= -3.7e-84) tmp = t_0; elseif (z <= 3.6e-44) tmp = 0.083333333333333 / x; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(y * N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.7e-84], t$95$0, If[LessEqual[z, 3.6e-44], N[(0.083333333333333 / x), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \frac{z}{\frac{x}{z}}\\
\mathbf{if}\;z \leq -3.7 \cdot 10^{-84}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{-44}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -3.6999999999999999e-84 or 3.5999999999999999e-44 < z Initial program 86.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6445.6%
Simplified45.6%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6447.3%
Applied egg-rr47.3%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6450.4%
Applied egg-rr50.4%
if -3.6999999999999999e-84 < z < 3.5999999999999999e-44Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6448.5%
Simplified48.5%
Taylor expanded in z around 0
/-lowering-/.f6446.9%
Simplified46.9%
Final simplification49.0%
(FPCore (x y z) :precision binary64 (let* ((t_0 (* (/ z x) (* y z)))) (if (<= z -3.7e-84) t_0 (if (<= z 8.5e-38) (/ 0.083333333333333 x) t_0))))
double code(double x, double y, double z) {
double t_0 = (z / x) * (y * z);
double tmp;
if (z <= -3.7e-84) {
tmp = t_0;
} else if (z <= 8.5e-38) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (z / x) * (y * z)
if (z <= (-3.7d-84)) then
tmp = t_0
else if (z <= 8.5d-38) then
tmp = 0.083333333333333d0 / x
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (z / x) * (y * z);
double tmp;
if (z <= -3.7e-84) {
tmp = t_0;
} else if (z <= 8.5e-38) {
tmp = 0.083333333333333 / x;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (z / x) * (y * z) tmp = 0 if z <= -3.7e-84: tmp = t_0 elif z <= 8.5e-38: tmp = 0.083333333333333 / x else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(z / x) * Float64(y * z)) tmp = 0.0 if (z <= -3.7e-84) tmp = t_0; elseif (z <= 8.5e-38) tmp = Float64(0.083333333333333 / x); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (z / x) * (y * z); tmp = 0.0; if (z <= -3.7e-84) tmp = t_0; elseif (z <= 8.5e-38) tmp = 0.083333333333333 / x; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z / x), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.7e-84], t$95$0, If[LessEqual[z, 8.5e-38], N[(0.083333333333333 / x), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{z}{x} \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -3.7 \cdot 10^{-84}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 8.5 \cdot 10^{-38}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -3.6999999999999999e-84 or 8.50000000000000046e-38 < z Initial program 86.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6445.6%
Simplified45.6%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6447.3%
Applied egg-rr47.3%
if -3.6999999999999999e-84 < z < 8.50000000000000046e-38Initial program 99.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6448.5%
Simplified48.5%
Taylor expanded in z around 0
/-lowering-/.f6446.9%
Simplified46.9%
Final simplification47.1%
(FPCore (x y z) :precision binary64 (+ (/ 0.083333333333333 x) (* z (+ (* (+ y 0.0007936500793651) (/ z x)) (/ -0.0027777777777778 x)))))
double code(double x, double y, double z) {
return (0.083333333333333 / x) + (z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / 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) + (z * (((y + 0.0007936500793651d0) * (z / x)) + ((-0.0027777777777778d0) / x)))
end function
public static double code(double x, double y, double z) {
return (0.083333333333333 / x) + (z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x)));
}
def code(x, y, z): return (0.083333333333333 / x) + (z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x)))
function code(x, y, z) return Float64(Float64(0.083333333333333 / x) + Float64(z * Float64(Float64(Float64(y + 0.0007936500793651) * Float64(z / x)) + Float64(-0.0027777777777778 / x)))) end
function tmp = code(x, y, z) tmp = (0.083333333333333 / x) + (z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x))); end
code[x_, y_, z_] := N[(N[(0.083333333333333 / x), $MachinePrecision] + N[(z * N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision] + N[(-0.0027777777777778 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333}{x} + z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x} + \frac{-0.0027777777777778}{x}\right)
\end{array}
Initial program 91.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6460.8%
Simplified60.8%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt-inN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-*r/N/A
associate-*l/N/A
associate-/l*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
Simplified62.0%
Final simplification62.0%
(FPCore (x y z) :precision binary64 (/ (+ 0.083333333333333 (* z (+ (* z (+ y 0.0007936500793651)) -0.0027777777777778))) x))
double code(double x, double y, double z) {
return (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / 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 + (z * ((z * (y + 0.0007936500793651d0)) + (-0.0027777777777778d0)))) / x
end function
public static double code(double x, double y, double z) {
return (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x;
}
def code(x, y, z): return (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x
function code(x, y, z) return Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) + -0.0027777777777778))) / x) end
function tmp = code(x, y, z) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) + -0.0027777777777778))) / x; end
code[x_, y_, z_] := N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] + -0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) + -0.0027777777777778\right)}{x}
\end{array}
Initial program 91.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6460.8%
Simplified60.8%
Final simplification60.8%
(FPCore (x y z) :precision binary64 (/ (+ 0.083333333333333 (* z (* z (+ y 0.0007936500793651)))) x))
double code(double x, double y, double z) {
return (0.083333333333333 + (z * (z * (y + 0.0007936500793651)))) / 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 + (z * (z * (y + 0.0007936500793651d0)))) / x
end function
public static double code(double x, double y, double z) {
return (0.083333333333333 + (z * (z * (y + 0.0007936500793651)))) / x;
}
def code(x, y, z): return (0.083333333333333 + (z * (z * (y + 0.0007936500793651)))) / x
function code(x, y, z) return Float64(Float64(0.083333333333333 + Float64(z * Float64(z * Float64(y + 0.0007936500793651)))) / x) end
function tmp = code(x, y, z) tmp = (0.083333333333333 + (z * (z * (y + 0.0007936500793651)))) / x; end
code[x_, y_, z_] := N[(N[(0.083333333333333 + N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right)\right)}{x}
\end{array}
Initial program 91.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6460.8%
Simplified60.8%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6460.1%
Simplified60.1%
Final simplification60.1%
(FPCore (x y z) :precision binary64 (if (<= z -1.35e+19) (* z (/ -0.0027777777777778 x)) (/ 0.083333333333333 x)))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.35e+19) {
tmp = z * (-0.0027777777777778 / x);
} else {
tmp = 0.083333333333333 / x;
}
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 <= (-1.35d+19)) then
tmp = z * ((-0.0027777777777778d0) / x)
else
tmp = 0.083333333333333d0 / x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -1.35e+19) {
tmp = z * (-0.0027777777777778 / x);
} else {
tmp = 0.083333333333333 / x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -1.35e+19: tmp = z * (-0.0027777777777778 / x) else: tmp = 0.083333333333333 / x return tmp
function code(x, y, z) tmp = 0.0 if (z <= -1.35e+19) tmp = Float64(z * Float64(-0.0027777777777778 / x)); else tmp = Float64(0.083333333333333 / x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -1.35e+19) tmp = z * (-0.0027777777777778 / x); else tmp = 0.083333333333333 / x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -1.35e+19], N[(z * N[(-0.0027777777777778 / x), $MachinePrecision]), $MachinePrecision], N[(0.083333333333333 / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.35 \cdot 10^{+19}:\\
\;\;\;\;z \cdot \frac{-0.0027777777777778}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\end{array}
\end{array}
if z < -1.35e19Initial program 80.5%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6468.3%
Simplified68.3%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6468.3%
Simplified68.3%
Taylor expanded in z around 0
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6415.9%
Simplified15.9%
if -1.35e19 < z Initial program 94.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6458.4%
Simplified58.4%
Taylor expanded in z around 0
/-lowering-/.f6427.8%
Simplified27.8%
(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 91.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6460.8%
Simplified60.8%
Taylor expanded in z around 0
/-lowering-/.f6421.8%
Simplified21.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 2024158
(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)))