
(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 25 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 1e-78)
(+
t_0
(+
(*
(/ 1.0 x)
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651))))))
(- 0.91893853320467 x)))
(+
(+ t_0 (+ 0.91893853320467 (/ 0.083333333333333 x)))
(-
(*
z
(+
(/ (+ -0.0027777777777778 (* z 0.0007936500793651)) x)
(* z (/ y x))))
x)))))
double code(double x, double y, double z) {
double t_0 = (x + -0.5) * log(x);
double tmp;
if (x <= 1e-78) {
tmp = t_0 + (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + (0.91893853320467 - x));
} else {
tmp = (t_0 + (0.91893853320467 + (0.083333333333333 / x))) + ((z * (((-0.0027777777777778 + (z * 0.0007936500793651)) / x) + (z * (y / 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) :: t_0
real(8) :: tmp
t_0 = (x + (-0.5d0)) * log(x)
if (x <= 1d-78) then
tmp = t_0 + (((1.0d0 / x) * (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0)))))) + (0.91893853320467d0 - x))
else
tmp = (t_0 + (0.91893853320467d0 + (0.083333333333333d0 / x))) + ((z * ((((-0.0027777777777778d0) + (z * 0.0007936500793651d0)) / x) + (z * (y / x)))) - x)
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 <= 1e-78) {
tmp = t_0 + (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + (0.91893853320467 - x));
} else {
tmp = (t_0 + (0.91893853320467 + (0.083333333333333 / x))) + ((z * (((-0.0027777777777778 + (z * 0.0007936500793651)) / x) + (z * (y / x)))) - x);
}
return tmp;
}
def code(x, y, z): t_0 = (x + -0.5) * math.log(x) tmp = 0 if x <= 1e-78: tmp = t_0 + (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + (0.91893853320467 - x)) else: tmp = (t_0 + (0.91893853320467 + (0.083333333333333 / x))) + ((z * (((-0.0027777777777778 + (z * 0.0007936500793651)) / x) + (z * (y / x)))) - x) return tmp
function code(x, y, z) t_0 = Float64(Float64(x + -0.5) * log(x)) tmp = 0.0 if (x <= 1e-78) tmp = Float64(t_0 + Float64(Float64(Float64(1.0 / x) * Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651)))))) + Float64(0.91893853320467 - x))); else tmp = Float64(Float64(t_0 + Float64(0.91893853320467 + Float64(0.083333333333333 / x))) + Float64(Float64(z * Float64(Float64(Float64(-0.0027777777777778 + Float64(z * 0.0007936500793651)) / x) + Float64(z * Float64(y / x)))) - x)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + -0.5) * log(x); tmp = 0.0; if (x <= 1e-78) tmp = t_0 + (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + (0.91893853320467 - x)); else tmp = (t_0 + (0.91893853320467 + (0.083333333333333 / x))) + ((z * (((-0.0027777777777778 + (z * 0.0007936500793651)) / x) + (z * (y / x)))) - x); 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, 1e-78], N[(t$95$0 + N[(N[(N[(1.0 / x), $MachinePrecision] * N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.91893853320467 - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$0 + N[(0.91893853320467 + N[(0.083333333333333 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(N[(-0.0027777777777778 + N[(z * 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + N[(z * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + -0.5\right) \cdot \log x\\
\mathbf{if}\;x \leq 10^{-78}:\\
\;\;\;\;t\_0 + \left(\frac{1}{x} \cdot \left(0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)\right) + \left(0.91893853320467 - x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(t\_0 + \left(0.91893853320467 + \frac{0.083333333333333}{x}\right)\right) + \left(z \cdot \left(\frac{-0.0027777777777778 + z \cdot 0.0007936500793651}{x} + z \cdot \frac{y}{x}\right) - x\right)\\
\end{array}
\end{array}
if x < 9.99999999999999999e-79Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6499.8%
Applied egg-rr99.8%
if 9.99999999999999999e-79 < x Initial program 91.2%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified91.2%
Taylor expanded in z around 0
Simplified99.6%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (* (+ x -0.5) (log x)) x)))
(if (<= x 5000000000000.0)
(+
(+
(*
(/ 1.0 x)
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651))))))
0.91893853320467)
t_0)
(+ (* z (* (/ z x) (+ y 0.0007936500793651))) t_0))))
double code(double x, double y, double z) {
double t_0 = ((x + -0.5) * log(x)) - x;
double tmp;
if (x <= 5000000000000.0) {
tmp = (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + 0.91893853320467) + t_0;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + 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 = ((x + (-0.5d0)) * log(x)) - x
if (x <= 5000000000000.0d0) then
tmp = (((1.0d0 / x) * (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0)))))) + 0.91893853320467d0) + t_0
else
tmp = (z * ((z / x) * (y + 0.0007936500793651d0))) + t_0
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)) - x;
double tmp;
if (x <= 5000000000000.0) {
tmp = (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + 0.91893853320467) + t_0;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0;
}
return tmp;
}
def code(x, y, z): t_0 = ((x + -0.5) * math.log(x)) - x tmp = 0 if x <= 5000000000000.0: tmp = (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + 0.91893853320467) + t_0 else: tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(Float64(x + -0.5) * log(x)) - x) tmp = 0.0 if (x <= 5000000000000.0) tmp = Float64(Float64(Float64(Float64(1.0 / x) * Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651)))))) + 0.91893853320467) + t_0); else tmp = Float64(Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) + t_0); end return tmp end
function tmp_2 = code(x, y, z) t_0 = ((x + -0.5) * log(x)) - x; tmp = 0.0; if (x <= 5000000000000.0) tmp = (((1.0 / x) * (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))))) + 0.91893853320467) + t_0; else tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]}, If[LessEqual[x, 5000000000000.0], N[(N[(N[(N[(1.0 / x), $MachinePrecision] * N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.91893853320467), $MachinePrecision] + t$95$0), $MachinePrecision], N[(N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + -0.5\right) \cdot \log x - x\\
\mathbf{if}\;x \leq 5000000000000:\\
\;\;\;\;\left(\frac{1}{x} \cdot \left(0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)\right) + 0.91893853320467\right) + t\_0\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right) + t\_0\\
\end{array}
\end{array}
if x < 5e12Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6499.7%
Applied egg-rr99.7%
if 5e12 < x Initial program 88.9%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.9%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr88.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6488.9%
Applied egg-rr88.9%
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%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (* (+ x -0.5) (log x)) x)))
(if (<= x 8600000000000.0)
(+
(+
0.91893853320467
(/
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651)))))
x))
t_0)
(+ (* z (* (/ z x) (+ y 0.0007936500793651))) t_0))))
double code(double x, double y, double z) {
double t_0 = ((x + -0.5) * log(x)) - x;
double tmp;
if (x <= 8600000000000.0) {
tmp = (0.91893853320467 + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)) + t_0;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + 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 = ((x + (-0.5d0)) * log(x)) - x
if (x <= 8600000000000.0d0) then
tmp = (0.91893853320467d0 + ((0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0))))) / x)) + t_0
else
tmp = (z * ((z / x) * (y + 0.0007936500793651d0))) + t_0
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)) - x;
double tmp;
if (x <= 8600000000000.0) {
tmp = (0.91893853320467 + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)) + t_0;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0;
}
return tmp;
}
def code(x, y, z): t_0 = ((x + -0.5) * math.log(x)) - x tmp = 0 if x <= 8600000000000.0: tmp = (0.91893853320467 + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)) + t_0 else: tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(Float64(x + -0.5) * log(x)) - x) tmp = 0.0 if (x <= 8600000000000.0) tmp = Float64(Float64(0.91893853320467 + Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651))))) / x)) + t_0); else tmp = Float64(Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) + t_0); end return tmp end
function tmp_2 = code(x, y, z) t_0 = ((x + -0.5) * log(x)) - x; tmp = 0.0; if (x <= 8600000000000.0) tmp = (0.91893853320467 + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)) + t_0; else tmp = (z * ((z / x) * (y + 0.0007936500793651))) + t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]}, If[LessEqual[x, 8600000000000.0], N[(N[(0.91893853320467 + N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + t$95$0), $MachinePrecision], N[(N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + -0.5\right) \cdot \log x - x\\
\mathbf{if}\;x \leq 8600000000000:\\
\;\;\;\;\left(0.91893853320467 + \frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)}{x}\right) + t\_0\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right) + t\_0\\
\end{array}
\end{array}
if x < 8.6e12Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.7%
if 8.6e12 < x Initial program 88.9%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.9%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr88.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6488.9%
Applied egg-rr88.9%
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%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (+ x -0.5) (log x))))
(if (<= x 9000000000000.0)
(+
t_0
(+
(- 0.91893853320467 x)
(/
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651)))))
x)))
(+ (* z (* (/ z x) (+ y 0.0007936500793651))) (- t_0 x)))))
double code(double x, double y, double z) {
double t_0 = (x + -0.5) * log(x);
double tmp;
if (x <= 9000000000000.0) {
tmp = t_0 + ((0.91893853320467 - x) + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x));
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (t_0 - 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) :: t_0
real(8) :: tmp
t_0 = (x + (-0.5d0)) * log(x)
if (x <= 9000000000000.0d0) then
tmp = t_0 + ((0.91893853320467d0 - x) + ((0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0))))) / x))
else
tmp = (z * ((z / x) * (y + 0.0007936500793651d0))) + (t_0 - x)
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 <= 9000000000000.0) {
tmp = t_0 + ((0.91893853320467 - x) + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x));
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (t_0 - x);
}
return tmp;
}
def code(x, y, z): t_0 = (x + -0.5) * math.log(x) tmp = 0 if x <= 9000000000000.0: tmp = t_0 + ((0.91893853320467 - x) + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)) else: tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (t_0 - x) return tmp
function code(x, y, z) t_0 = Float64(Float64(x + -0.5) * log(x)) tmp = 0.0 if (x <= 9000000000000.0) tmp = Float64(t_0 + Float64(Float64(0.91893853320467 - x) + Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651))))) / x))); else tmp = Float64(Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) + Float64(t_0 - x)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + -0.5) * log(x); tmp = 0.0; if (x <= 9000000000000.0) tmp = t_0 + ((0.91893853320467 - x) + ((0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x)); else tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (t_0 - x); 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, 9000000000000.0], N[(t$95$0 + N[(N[(0.91893853320467 - x), $MachinePrecision] + N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$0 - x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x + -0.5\right) \cdot \log x\\
\mathbf{if}\;x \leq 9000000000000:\\
\;\;\;\;t\_0 + \left(\left(0.91893853320467 - x\right) + \frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right) + \left(t\_0 - x\right)\\
\end{array}
\end{array}
if x < 9e12Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
if 9e12 < x Initial program 88.9%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.9%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr88.9%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6488.9%
Applied egg-rr88.9%
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%
Final simplification99.6%
(FPCore (x y z)
:precision binary64
(if (<= x 0.0125)
(/
(+
0.083333333333333
(+
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651))))
(* x (+ 0.91893853320467 (* -0.5 (log x))))))
x)
(+
(* z (* (/ z x) (+ y 0.0007936500793651)))
(- (* (+ x -0.5) (log x)) x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 0.0125) {
tmp = (0.083333333333333 + ((z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))) + (x * (0.91893853320467 + (-0.5 * log(x)))))) / x;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * 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 <= 0.0125d0) then
tmp = (0.083333333333333d0 + ((z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0)))) + (x * (0.91893853320467d0 + ((-0.5d0) * log(x)))))) / x
else
tmp = (z * ((z / x) * (y + 0.0007936500793651d0))) + (((x + (-0.5d0)) * log(x)) - x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 0.0125) {
tmp = (0.083333333333333 + ((z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))) + (x * (0.91893853320467 + (-0.5 * Math.log(x)))))) / x;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * Math.log(x)) - x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 0.0125: tmp = (0.083333333333333 + ((z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))) + (x * (0.91893853320467 + (-0.5 * math.log(x)))))) / x else: tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * math.log(x)) - x) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 0.0125) tmp = Float64(Float64(0.083333333333333 + Float64(Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651)))) + Float64(x * Float64(0.91893853320467 + Float64(-0.5 * log(x)))))) / x); else tmp = Float64(Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) + Float64(Float64(Float64(x + -0.5) * log(x)) - x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 0.0125) tmp = (0.083333333333333 + ((z * (-0.0027777777777778 + (z * (y + 0.0007936500793651)))) + (x * (0.91893853320467 + (-0.5 * log(x)))))) / x; else tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * log(x)) - x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 0.0125], N[(N[(0.083333333333333 + N[(N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(0.91893853320467 + N[(-0.5 * N[Log[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.0125:\\
\;\;\;\;\frac{0.083333333333333 + \left(z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right) + x \cdot \left(0.91893853320467 + -0.5 \cdot \log x\right)\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right) + \left(\left(x + -0.5\right) \cdot \log x - x\right)\\
\end{array}
\end{array}
if x < 0.012500000000000001Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
Simplified99.2%
if 0.012500000000000001 < x Initial program 89.5%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified89.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr89.5%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6489.5%
Applied egg-rr89.5%
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 simplification99.0%
(FPCore (x y z)
:precision binary64
(if (<= z -8e-16)
(*
z
(*
z
(+
(/ y x)
(* (/ 1.0 x) (+ 0.0007936500793651 (/ -0.0027777777777778 z))))))
(if (<= z 4.5e+16)
(+
(+ 0.91893853320467 (/ 0.083333333333333 x))
(- (* (+ x -0.5) (log x)) x))
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -8e-16) {
tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z)))));
} else if (z <= 4.5e+16) {
tmp = (0.91893853320467 + (0.083333333333333 / x)) + (((x + -0.5) * log(x)) - 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 <= (-8d-16)) then
tmp = z * (z * ((y / x) + ((1.0d0 / x) * (0.0007936500793651d0 + ((-0.0027777777777778d0) / z)))))
else if (z <= 4.5d+16) then
tmp = (0.91893853320467d0 + (0.083333333333333d0 / x)) + (((x + (-0.5d0)) * log(x)) - 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 <= -8e-16) {
tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z)))));
} else if (z <= 4.5e+16) {
tmp = (0.91893853320467 + (0.083333333333333 / x)) + (((x + -0.5) * Math.log(x)) - x);
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -8e-16: tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z))))) elif z <= 4.5e+16: tmp = (0.91893853320467 + (0.083333333333333 / x)) + (((x + -0.5) * math.log(x)) - x) else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -8e-16) tmp = Float64(z * Float64(z * Float64(Float64(y / x) + Float64(Float64(1.0 / x) * Float64(0.0007936500793651 + Float64(-0.0027777777777778 / z)))))); elseif (z <= 4.5e+16) tmp = Float64(Float64(0.91893853320467 + Float64(0.083333333333333 / x)) + Float64(Float64(Float64(x + -0.5) * log(x)) - 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 <= -8e-16) tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z))))); elseif (z <= 4.5e+16) tmp = (0.91893853320467 + (0.083333333333333 / x)) + (((x + -0.5) * log(x)) - x); else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -8e-16], N[(z * N[(z * N[(N[(y / x), $MachinePrecision] + N[(N[(1.0 / x), $MachinePrecision] * N[(0.0007936500793651 + N[(-0.0027777777777778 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.5e+16], N[(N[(0.91893853320467 + N[(0.083333333333333 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $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 -8 \cdot 10^{-16}:\\
\;\;\;\;z \cdot \left(z \cdot \left(\frac{y}{x} + \frac{1}{x} \cdot \left(0.0007936500793651 + \frac{-0.0027777777777778}{z}\right)\right)\right)\\
\mathbf{elif}\;z \leq 4.5 \cdot 10^{+16}:\\
\;\;\;\;\left(0.91893853320467 + \frac{0.083333333333333}{x}\right) + \left(\left(x + -0.5\right) \cdot \log x - x\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -7.9999999999999998e-16Initial program 93.1%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.1%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
Simplified79.8%
if -7.9999999999999998e-16 < z < 4.5e16Initial program 99.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.5%
Taylor expanded in z around 0
/-lowering-/.f6494.7%
Simplified94.7%
if 4.5e16 < z Initial program 84.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified84.4%
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-/.f6490.6%
Simplified90.6%
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-+.f6490.7%
Simplified90.7%
Final simplification89.9%
(FPCore (x y z)
:precision binary64
(if (<= z -8e-16)
(*
z
(*
z
(+
(/ y x)
(* (/ 1.0 x) (+ 0.0007936500793651 (/ -0.0027777777777778 z))))))
(if (<= z 5.6e+16)
(+
(* (+ x -0.5) (log x))
(- (+ 0.91893853320467 (/ 0.083333333333333 x)) x))
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -8e-16) {
tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z)))));
} else if (z <= 5.6e+16) {
tmp = ((x + -0.5) * log(x)) + ((0.91893853320467 + (0.083333333333333 / x)) - 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 <= (-8d-16)) then
tmp = z * (z * ((y / x) + ((1.0d0 / x) * (0.0007936500793651d0 + ((-0.0027777777777778d0) / z)))))
else if (z <= 5.6d+16) then
tmp = ((x + (-0.5d0)) * log(x)) + ((0.91893853320467d0 + (0.083333333333333d0 / x)) - 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 <= -8e-16) {
tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z)))));
} else if (z <= 5.6e+16) {
tmp = ((x + -0.5) * Math.log(x)) + ((0.91893853320467 + (0.083333333333333 / x)) - x);
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -8e-16: tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z))))) elif z <= 5.6e+16: tmp = ((x + -0.5) * math.log(x)) + ((0.91893853320467 + (0.083333333333333 / x)) - x) else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -8e-16) tmp = Float64(z * Float64(z * Float64(Float64(y / x) + Float64(Float64(1.0 / x) * Float64(0.0007936500793651 + Float64(-0.0027777777777778 / z)))))); elseif (z <= 5.6e+16) tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(Float64(0.91893853320467 + Float64(0.083333333333333 / x)) - 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 <= -8e-16) tmp = z * (z * ((y / x) + ((1.0 / x) * (0.0007936500793651 + (-0.0027777777777778 / z))))); elseif (z <= 5.6e+16) tmp = ((x + -0.5) * log(x)) + ((0.91893853320467 + (0.083333333333333 / x)) - x); else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -8e-16], N[(z * N[(z * N[(N[(y / x), $MachinePrecision] + N[(N[(1.0 / x), $MachinePrecision] * N[(0.0007936500793651 + N[(-0.0027777777777778 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 5.6e+16], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(0.91893853320467 + N[(0.083333333333333 / x), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $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 -8 \cdot 10^{-16}:\\
\;\;\;\;z \cdot \left(z \cdot \left(\frac{y}{x} + \frac{1}{x} \cdot \left(0.0007936500793651 + \frac{-0.0027777777777778}{z}\right)\right)\right)\\
\mathbf{elif}\;z \leq 5.6 \cdot 10^{+16}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(\left(0.91893853320467 + \frac{0.083333333333333}{x}\right) - x\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -7.9999999999999998e-16Initial program 93.1%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.1%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
Simplified79.8%
if -7.9999999999999998e-16 < z < 5.6e16Initial program 99.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
Taylor expanded in z around 0
associate-+r+N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
remove-double-negN/A
log-recN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
log-recN/A
remove-double-negN/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6494.7%
Simplified94.7%
if 5.6e16 < z Initial program 84.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified84.4%
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-/.f6490.6%
Simplified90.6%
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-+.f6490.7%
Simplified90.7%
Final simplification89.9%
(FPCore (x y z)
:precision binary64
(if (<= x 0.0125)
(/
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651)))))
x)
(+
(* z (* (/ z x) (+ y 0.0007936500793651)))
(- (* (+ x -0.5) (log x)) x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 0.0125) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * 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 <= 0.0125d0) then
tmp = (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0))))) / x
else
tmp = (z * ((z / x) * (y + 0.0007936500793651d0))) + (((x + (-0.5d0)) * log(x)) - x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 0.0125) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * Math.log(x)) - x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 0.0125: tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x else: tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * math.log(x)) - x) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 0.0125) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651))))) / x); else tmp = Float64(Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) + Float64(Float64(Float64(x + -0.5) * log(x)) - x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 0.0125) tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x; else tmp = (z * ((z / x) * (y + 0.0007936500793651))) + (((x + -0.5) * log(x)) - x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 0.0125], N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.0125:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right) + \left(\left(x + -0.5\right) \cdot \log x - x\right)\\
\end{array}
\end{array}
if x < 0.012500000000000001Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.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-+.f6498.1%
Simplified98.1%
if 0.012500000000000001 < x Initial program 89.5%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified89.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr89.5%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6489.5%
Applied egg-rr89.5%
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.4%
(FPCore (x y z)
:precision binary64
(if (<= x 16.5)
(/
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651)))))
x)
(+ (* (+ x -0.5) (log x)) (+ (- 0.91893853320467 x) (* (/ z x) (* z y))))))
double code(double x, double y, double z) {
double tmp;
if (x <= 16.5) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = ((x + -0.5) * log(x)) + ((0.91893853320467 - x) + ((z / x) * (z * y)));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= 16.5d0) then
tmp = (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0))))) / x
else
tmp = ((x + (-0.5d0)) * log(x)) + ((0.91893853320467d0 - x) + ((z / x) * (z * y)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 16.5) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = ((x + -0.5) * Math.log(x)) + ((0.91893853320467 - x) + ((z / x) * (z * y)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 16.5: tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x else: tmp = ((x + -0.5) * math.log(x)) + ((0.91893853320467 - x) + ((z / x) * (z * y))) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 16.5) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651))))) / x); else tmp = Float64(Float64(Float64(x + -0.5) * log(x)) + Float64(Float64(0.91893853320467 - x) + Float64(Float64(z / x) * Float64(z * y)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 16.5) tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x; else tmp = ((x + -0.5) * log(x)) + ((0.91893853320467 - x) + ((z / x) * (z * y))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 16.5], N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] + N[(N[(0.91893853320467 - x), $MachinePrecision] + N[(N[(z / x), $MachinePrecision] * N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 16.5:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -0.5\right) \cdot \log x + \left(\left(0.91893853320467 - x\right) + \frac{z}{x} \cdot \left(z \cdot y\right)\right)\\
\end{array}
\end{array}
if x < 16.5Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.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-+.f6498.1%
Simplified98.1%
if 16.5 < x Initial program 89.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified89.4%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.3%
Simplified82.3%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6489.6%
Applied egg-rr89.6%
Final simplification94.2%
(FPCore (x y z)
:precision binary64
(if (<= x 920.0)
(/
(+
0.083333333333333
(* z (+ -0.0027777777777778 (* z (+ y 0.0007936500793651)))))
x)
(+ (/ (* z (* z y)) x) (- (* (+ x -0.5) (log x)) x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 920.0) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = ((z * (z * y)) / x) + (((x + -0.5) * 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 <= 920.0d0) then
tmp = (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + (z * (y + 0.0007936500793651d0))))) / x
else
tmp = ((z * (z * y)) / x) + (((x + (-0.5d0)) * log(x)) - x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 920.0) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x;
} else {
tmp = ((z * (z * y)) / x) + (((x + -0.5) * Math.log(x)) - x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 920.0: tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x else: tmp = ((z * (z * y)) / x) + (((x + -0.5) * math.log(x)) - x) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 920.0) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + Float64(z * Float64(y + 0.0007936500793651))))) / x); else tmp = Float64(Float64(Float64(z * Float64(z * y)) / x) + Float64(Float64(Float64(x + -0.5) * log(x)) - x)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 920.0) tmp = (0.083333333333333 + (z * (-0.0027777777777778 + (z * (y + 0.0007936500793651))))) / x; else tmp = ((z * (z * y)) / x) + (((x + -0.5) * log(x)) - x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 920.0], N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(N[(z * N[(z * y), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 920:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + z \cdot \left(y + 0.0007936500793651\right)\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{z \cdot \left(z \cdot y\right)}{x} + \left(\left(x + -0.5\right) \cdot \log x - x\right)\\
\end{array}
\end{array}
if x < 920Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.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-+.f6498.1%
Simplified98.1%
if 920 < x Initial program 89.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified89.4%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr89.4%
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-*.f6485.1%
Simplified85.1%
Final simplification92.1%
(FPCore (x y z)
:precision binary64
(if (<= x 1.75e+32)
(/
(+
0.083333333333333
(* y (* z (+ z (/ (+ -0.0027777777777778 (* z 0.0007936500793651)) y)))))
x)
(if (<= x 7.8e+107)
(* z (+ (* z (/ y x)) (/ z (/ x 0.0007936500793651))))
(* x (+ (log x) -1.0)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 1.75e+32) {
tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x;
} else if (x <= 7.8e+107) {
tmp = z * ((z * (y / x)) + (z / (x / 0.0007936500793651)));
} 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 <= 1.75d+32) then
tmp = (0.083333333333333d0 + (y * (z * (z + (((-0.0027777777777778d0) + (z * 0.0007936500793651d0)) / y))))) / x
else if (x <= 7.8d+107) then
tmp = z * ((z * (y / x)) + (z / (x / 0.0007936500793651d0)))
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 <= 1.75e+32) {
tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x;
} else if (x <= 7.8e+107) {
tmp = z * ((z * (y / x)) + (z / (x / 0.0007936500793651)));
} else {
tmp = x * (Math.log(x) + -1.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 1.75e+32: tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x elif x <= 7.8e+107: tmp = z * ((z * (y / x)) + (z / (x / 0.0007936500793651))) else: tmp = x * (math.log(x) + -1.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 1.75e+32) tmp = Float64(Float64(0.083333333333333 + Float64(y * Float64(z * Float64(z + Float64(Float64(-0.0027777777777778 + Float64(z * 0.0007936500793651)) / y))))) / x); elseif (x <= 7.8e+107) tmp = Float64(z * Float64(Float64(z * Float64(y / x)) + Float64(z / Float64(x / 0.0007936500793651)))); 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 <= 1.75e+32) tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x; elseif (x <= 7.8e+107) tmp = z * ((z * (y / x)) + (z / (x / 0.0007936500793651))); else tmp = x * (log(x) + -1.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 1.75e+32], N[(N[(0.083333333333333 + N[(y * N[(z * N[(z + N[(N[(-0.0027777777777778 + N[(z * 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, 7.8e+107], N[(z * N[(N[(z * N[(y / x), $MachinePrecision]), $MachinePrecision] + N[(z / N[(x / 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.75 \cdot 10^{+32}:\\
\;\;\;\;\frac{0.083333333333333 + y \cdot \left(z \cdot \left(z + \frac{-0.0027777777777778 + z \cdot 0.0007936500793651}{y}\right)\right)}{x}\\
\mathbf{elif}\;x \leq 7.8 \cdot 10^{+107}:\\
\;\;\;\;z \cdot \left(z \cdot \frac{y}{x} + \frac{z}{\frac{x}{0.0007936500793651}}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 1.75e32Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
Taylor expanded in y around -inf
Simplified69.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
Simplified95.1%
if 1.75e32 < x < 7.7999999999999997e107Initial program 85.8%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified85.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-/.f6464.3%
Simplified64.3%
distribute-lft-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6464.4%
Applied egg-rr64.4%
if 7.7999999999999997e107 < x Initial program 89.6%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified89.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.f6488.1%
Simplified88.1%
Final simplification89.1%
(FPCore (x y z)
:precision binary64
(if (<= x 5.2e+32)
(/
(+
0.083333333333333
(* y (* z (+ z (/ (+ -0.0027777777777778 (* z 0.0007936500793651)) y)))))
x)
(* z (/ (* z (+ y 0.0007936500793651)) x))))
double code(double x, double y, double z) {
double tmp;
if (x <= 5.2e+32) {
tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / 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 (x <= 5.2d+32) then
tmp = (0.083333333333333d0 + (y * (z * (z + (((-0.0027777777777778d0) + (z * 0.0007936500793651d0)) / y))))) / 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 (x <= 5.2e+32) {
tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 5.2e+32: tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 5.2e+32) tmp = Float64(Float64(0.083333333333333 + Float64(y * Float64(z * Float64(z + Float64(Float64(-0.0027777777777778 + Float64(z * 0.0007936500793651)) / y))))) / 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 (x <= 5.2e+32) tmp = (0.083333333333333 + (y * (z * (z + ((-0.0027777777777778 + (z * 0.0007936500793651)) / y))))) / x; else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 5.2e+32], N[(N[(0.083333333333333 + N[(y * N[(z * N[(z + N[(N[(-0.0027777777777778 + N[(z * 0.0007936500793651), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $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}\;x \leq 5.2 \cdot 10^{+32}:\\
\;\;\;\;\frac{0.083333333333333 + y \cdot \left(z \cdot \left(z + \frac{-0.0027777777777778 + z \cdot 0.0007936500793651}{y}\right)\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if x < 5.2000000000000004e32Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
Taylor expanded in y around -inf
Simplified69.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
Simplified95.1%
if 5.2000000000000004e32 < x Initial program 88.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.4%
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-/.f6429.4%
Simplified29.4%
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-+.f6429.5%
Simplified29.5%
Final simplification67.7%
(FPCore (x y z)
:precision binary64
(if (<= z -110.0)
(* z (* (/ z x) (+ y 0.0007936500793651)))
(if (<= z 2.7e-5)
(/ (+ 0.083333333333333 (* z (* z y))) x)
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -110.0) {
tmp = z * ((z / x) * (y + 0.0007936500793651));
} else if (z <= 2.7e-5) {
tmp = (0.083333333333333 + (z * (z * y))) / 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 <= (-110.0d0)) then
tmp = z * ((z / x) * (y + 0.0007936500793651d0))
else if (z <= 2.7d-5) then
tmp = (0.083333333333333d0 + (z * (z * y))) / 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 <= -110.0) {
tmp = z * ((z / x) * (y + 0.0007936500793651));
} else if (z <= 2.7e-5) {
tmp = (0.083333333333333 + (z * (z * y))) / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -110.0: tmp = z * ((z / x) * (y + 0.0007936500793651)) elif z <= 2.7e-5: tmp = (0.083333333333333 + (z * (z * y))) / x else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -110.0) tmp = Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))); elseif (z <= 2.7e-5) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(z * y))) / 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 <= -110.0) tmp = z * ((z / x) * (y + 0.0007936500793651)); elseif (z <= 2.7e-5) tmp = (0.083333333333333 + (z * (z * y))) / x; else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -110.0], N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.7e-5], N[(N[(0.083333333333333 + N[(z * N[(z * y), $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 -110:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right)\\
\mathbf{elif}\;z \leq 2.7 \cdot 10^{-5}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot y\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -110Initial program 92.8%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified92.8%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr92.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6492.7%
Applied egg-rr92.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-+.f6480.2%
Simplified80.2%
if -110 < z < 2.6999999999999999e-5Initial program 99.5%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
Taylor expanded in x around 0
Simplified75.7%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6455.5%
Simplified55.5%
if 2.6999999999999999e-5 < z Initial program 86.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified86.4%
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-/.f6482.8%
Simplified82.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-+.f6482.8%
Simplified82.8%
Final simplification67.6%
(FPCore (x y z)
:precision binary64
(if (<= z -1.9e-19)
(* z (* (/ z x) (+ y 0.0007936500793651)))
(if (<= z 3.1e-13)
(/ (+ 0.083333333333333 (* z -0.0027777777777778)) x)
(* z (/ (* z (+ y 0.0007936500793651)) x)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.9e-19) {
tmp = z * ((z / x) * (y + 0.0007936500793651));
} else if (z <= 3.1e-13) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / 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 <= (-1.9d-19)) then
tmp = z * ((z / x) * (y + 0.0007936500793651d0))
else if (z <= 3.1d-13) then
tmp = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / 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 <= -1.9e-19) {
tmp = z * ((z / x) * (y + 0.0007936500793651));
} else if (z <= 3.1e-13) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = z * ((z * (y + 0.0007936500793651)) / x);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -1.9e-19: tmp = z * ((z / x) * (y + 0.0007936500793651)) elif z <= 3.1e-13: tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x else: tmp = z * ((z * (y + 0.0007936500793651)) / x) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -1.9e-19) tmp = Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))); elseif (z <= 3.1e-13) tmp = Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / 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 <= -1.9e-19) tmp = z * ((z / x) * (y + 0.0007936500793651)); elseif (z <= 3.1e-13) tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; else tmp = z * ((z * (y + 0.0007936500793651)) / x); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -1.9e-19], N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.1e-13], N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $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 -1.9 \cdot 10^{-19}:\\
\;\;\;\;z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right)\\
\mathbf{elif}\;z \leq 3.1 \cdot 10^{-13}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right)}{x}\\
\end{array}
\end{array}
if z < -1.9e-19Initial program 93.1%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.1%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr93.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6493.0%
Applied egg-rr93.0%
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-+.f6479.7%
Simplified79.7%
if -1.9e-19 < z < 3.0999999999999999e-13Initial program 99.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.5%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6495.9%
Simplified95.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6451.6%
Simplified51.6%
if 3.0999999999999999e-13 < z Initial program 86.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified86.4%
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-/.f6482.8%
Simplified82.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-+.f6482.8%
Simplified82.8%
Final simplification65.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (* (/ z x) (+ y 0.0007936500793651)))))
(if (<= z -3.5e-17)
t_0
(if (<= z 2.8e-10)
(/ (+ 0.083333333333333 (* z -0.0027777777777778)) x)
t_0))))
double code(double x, double y, double z) {
double t_0 = z * ((z / x) * (y + 0.0007936500793651));
double tmp;
if (z <= -3.5e-17) {
tmp = t_0;
} else if (z <= 2.8e-10) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / 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 * ((z / x) * (y + 0.0007936500793651d0))
if (z <= (-3.5d-17)) then
tmp = t_0
else if (z <= 2.8d-10) then
tmp = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / 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 * ((z / x) * (y + 0.0007936500793651));
double tmp;
if (z <= -3.5e-17) {
tmp = t_0;
} else if (z <= 2.8e-10) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = z * ((z / x) * (y + 0.0007936500793651)) tmp = 0 if z <= -3.5e-17: tmp = t_0 elif z <= 2.8e-10: tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(z * Float64(Float64(z / x) * Float64(y + 0.0007936500793651))) tmp = 0.0 if (z <= -3.5e-17) tmp = t_0; elseif (z <= 2.8e-10) tmp = Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / x); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * ((z / x) * (y + 0.0007936500793651)); tmp = 0.0; if (z <= -3.5e-17) tmp = t_0; elseif (z <= 2.8e-10) tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(N[(z / x), $MachinePrecision] * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.5e-17], t$95$0, If[LessEqual[z, 2.8e-10], N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(\frac{z}{x} \cdot \left(y + 0.0007936500793651\right)\right)\\
\mathbf{if}\;z \leq -3.5 \cdot 10^{-17}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 2.8 \cdot 10^{-10}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -3.5000000000000002e-17 or 2.80000000000000015e-10 < z Initial program 90.1%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified90.1%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr90.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6490.1%
Applied egg-rr90.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-+.f6481.0%
Simplified81.0%
if -3.5000000000000002e-17 < z < 2.80000000000000015e-10Initial program 99.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.5%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6495.9%
Simplified95.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6451.6%
Simplified51.6%
Final simplification65.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (+ y 0.0007936500793651))))
(if (<= x 1.16e+33)
(/ (+ 0.083333333333333 (* z (+ -0.0027777777777778 t_0))) x)
(* z (/ t_0 x)))))
double code(double x, double y, double z) {
double t_0 = z * (y + 0.0007936500793651);
double tmp;
if (x <= 1.16e+33) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + t_0))) / x;
} else {
tmp = z * (t_0 / 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) :: t_0
real(8) :: tmp
t_0 = z * (y + 0.0007936500793651d0)
if (x <= 1.16d+33) then
tmp = (0.083333333333333d0 + (z * ((-0.0027777777777778d0) + t_0))) / x
else
tmp = z * (t_0 / x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = z * (y + 0.0007936500793651);
double tmp;
if (x <= 1.16e+33) {
tmp = (0.083333333333333 + (z * (-0.0027777777777778 + t_0))) / x;
} else {
tmp = z * (t_0 / x);
}
return tmp;
}
def code(x, y, z): t_0 = z * (y + 0.0007936500793651) tmp = 0 if x <= 1.16e+33: tmp = (0.083333333333333 + (z * (-0.0027777777777778 + t_0))) / x else: tmp = z * (t_0 / x) return tmp
function code(x, y, z) t_0 = Float64(z * Float64(y + 0.0007936500793651)) tmp = 0.0 if (x <= 1.16e+33) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(-0.0027777777777778 + t_0))) / x); else tmp = Float64(z * Float64(t_0 / x)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * (y + 0.0007936500793651); tmp = 0.0; if (x <= 1.16e+33) tmp = (0.083333333333333 + (z * (-0.0027777777777778 + t_0))) / x; else tmp = z * (t_0 / x); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.16e+33], N[(N[(0.083333333333333 + N[(z * N[(-0.0027777777777778 + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(t$95$0 / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(y + 0.0007936500793651\right)\\
\mathbf{if}\;x \leq 1.16 \cdot 10^{+33}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(-0.0027777777777778 + t\_0\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{t\_0}{x}\\
\end{array}
\end{array}
if x < 1.16000000000000001e33Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.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-+.f6495.1%
Simplified95.1%
if 1.16000000000000001e33 < x Initial program 88.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.4%
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-/.f6429.4%
Simplified29.4%
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-+.f6429.5%
Simplified29.5%
Final simplification67.7%
(FPCore (x y z)
:precision binary64
(if (<= z -3.85e-16)
(* (/ z x) (* z y))
(if (<= z 3.2e-8)
(/ (+ 0.083333333333333 (* z -0.0027777777777778)) x)
(* y (/ z (/ x z))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -3.85e-16) {
tmp = (z / x) * (z * y);
} else if (z <= 3.2e-8) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = y * (z / (x / z));
}
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 <= (-3.85d-16)) then
tmp = (z / x) * (z * y)
else if (z <= 3.2d-8) then
tmp = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / x
else
tmp = y * (z / (x / z))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -3.85e-16) {
tmp = (z / x) * (z * y);
} else if (z <= 3.2e-8) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = y * (z / (x / z));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -3.85e-16: tmp = (z / x) * (z * y) elif z <= 3.2e-8: tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x else: tmp = y * (z / (x / z)) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -3.85e-16) tmp = Float64(Float64(z / x) * Float64(z * y)); elseif (z <= 3.2e-8) tmp = Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / x); else tmp = Float64(y * Float64(z / Float64(x / z))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -3.85e-16) tmp = (z / x) * (z * y); elseif (z <= 3.2e-8) tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; else tmp = y * (z / (x / z)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -3.85e-16], N[(N[(z / x), $MachinePrecision] * N[(z * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.2e-8], N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(y * N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.85 \cdot 10^{-16}:\\
\;\;\;\;\frac{z}{x} \cdot \left(z \cdot y\right)\\
\mathbf{elif}\;z \leq 3.2 \cdot 10^{-8}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z}{\frac{x}{z}}\\
\end{array}
\end{array}
if z < -3.84999999999999994e-16Initial program 93.1%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.1%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.8%
Simplified54.8%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6457.5%
Applied egg-rr57.5%
if -3.84999999999999994e-16 < z < 3.2000000000000002e-8Initial program 99.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.5%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr99.5%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6495.9%
Simplified95.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6451.6%
Simplified51.6%
if 3.2000000000000002e-8 < z Initial program 86.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified86.4%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.3%
Simplified55.3%
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6457.3%
Applied egg-rr57.3%
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6460.8%
Applied egg-rr60.8%
Final simplification55.1%
(FPCore (x y z)
:precision binary64
(if (<= y -2.4e+16)
(/ y (/ x (* z z)))
(if (<= y 0.0008)
(* z (* z (/ 0.0007936500793651 x)))
(* (/ z x) (* z y)))))
double code(double x, double y, double z) {
double tmp;
if (y <= -2.4e+16) {
tmp = y / (x / (z * z));
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = (z / x) * (z * y);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (y <= (-2.4d+16)) then
tmp = y / (x / (z * z))
else if (y <= 0.0008d0) then
tmp = z * (z * (0.0007936500793651d0 / x))
else
tmp = (z / x) * (z * y)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -2.4e+16) {
tmp = y / (x / (z * z));
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = (z / x) * (z * y);
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -2.4e+16: tmp = y / (x / (z * z)) elif y <= 0.0008: tmp = z * (z * (0.0007936500793651 / x)) else: tmp = (z / x) * (z * y) return tmp
function code(x, y, z) tmp = 0.0 if (y <= -2.4e+16) tmp = Float64(y / Float64(x / Float64(z * z))); elseif (y <= 0.0008) tmp = Float64(z * Float64(z * Float64(0.0007936500793651 / x))); else tmp = Float64(Float64(z / x) * Float64(z * y)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -2.4e+16) tmp = y / (x / (z * z)); elseif (y <= 0.0008) tmp = z * (z * (0.0007936500793651 / x)); else tmp = (z / x) * (z * y); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -2.4e+16], N[(y / N[(x / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.0008], N[(z * N[(z * N[(0.0007936500793651 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z / x), $MachinePrecision] * N[(z * y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.4 \cdot 10^{+16}:\\
\;\;\;\;\frac{y}{\frac{x}{z \cdot z}}\\
\mathbf{elif}\;y \leq 0.0008:\\
\;\;\;\;z \cdot \left(z \cdot \frac{0.0007936500793651}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{z}{x} \cdot \left(z \cdot y\right)\\
\end{array}
\end{array}
if y < -2.4e16Initial program 98.0%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified98.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.8%
Simplified50.8%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6452.4%
Applied egg-rr52.4%
if -2.4e16 < y < 8.00000000000000038e-4Initial program 93.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.7%
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-/.f6432.5%
Simplified32.5%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6432.5%
Simplified32.5%
if 8.00000000000000038e-4 < y Initial program 95.2%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified95.2%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.6%
Simplified53.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6455.8%
Applied egg-rr55.8%
Final simplification42.3%
(FPCore (x y z)
:precision binary64
(if (<= y -2.4e+16)
(* y (/ z (/ x z)))
(if (<= y 0.0008)
(* z (* z (/ 0.0007936500793651 x)))
(* (/ z x) (* z y)))))
double code(double x, double y, double z) {
double tmp;
if (y <= -2.4e+16) {
tmp = y * (z / (x / z));
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = (z / x) * (z * y);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (y <= (-2.4d+16)) then
tmp = y * (z / (x / z))
else if (y <= 0.0008d0) then
tmp = z * (z * (0.0007936500793651d0 / x))
else
tmp = (z / x) * (z * y)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -2.4e+16) {
tmp = y * (z / (x / z));
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = (z / x) * (z * y);
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -2.4e+16: tmp = y * (z / (x / z)) elif y <= 0.0008: tmp = z * (z * (0.0007936500793651 / x)) else: tmp = (z / x) * (z * y) return tmp
function code(x, y, z) tmp = 0.0 if (y <= -2.4e+16) tmp = Float64(y * Float64(z / Float64(x / z))); elseif (y <= 0.0008) tmp = Float64(z * Float64(z * Float64(0.0007936500793651 / x))); else tmp = Float64(Float64(z / x) * Float64(z * y)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -2.4e+16) tmp = y * (z / (x / z)); elseif (y <= 0.0008) tmp = z * (z * (0.0007936500793651 / x)); else tmp = (z / x) * (z * y); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -2.4e+16], N[(y * N[(z / N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.0008], N[(z * N[(z * N[(0.0007936500793651 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z / x), $MachinePrecision] * N[(z * y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.4 \cdot 10^{+16}:\\
\;\;\;\;y \cdot \frac{z}{\frac{x}{z}}\\
\mathbf{elif}\;y \leq 0.0008:\\
\;\;\;\;z \cdot \left(z \cdot \frac{0.0007936500793651}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{z}{x} \cdot \left(z \cdot y\right)\\
\end{array}
\end{array}
if y < -2.4e16Initial program 98.0%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified98.0%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.8%
Simplified50.8%
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6452.3%
Applied egg-rr52.3%
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6452.4%
Applied egg-rr52.4%
if -2.4e16 < y < 8.00000000000000038e-4Initial program 93.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.7%
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-/.f6432.5%
Simplified32.5%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6432.5%
Simplified32.5%
if 8.00000000000000038e-4 < y Initial program 95.2%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified95.2%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.6%
Simplified53.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6455.8%
Applied egg-rr55.8%
Final simplification42.3%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (/ z x) (* z y))))
(if (<= y -2.4e+16)
t_0
(if (<= y 0.0008) (* z (* z (/ 0.0007936500793651 x))) t_0))))
double code(double x, double y, double z) {
double t_0 = (z / x) * (z * y);
double tmp;
if (y <= -2.4e+16) {
tmp = t_0;
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / 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) * (z * y)
if (y <= (-2.4d+16)) then
tmp = t_0
else if (y <= 0.0008d0) then
tmp = z * (z * (0.0007936500793651d0 / 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) * (z * y);
double tmp;
if (y <= -2.4e+16) {
tmp = t_0;
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (z / x) * (z * y) tmp = 0 if y <= -2.4e+16: tmp = t_0 elif y <= 0.0008: tmp = z * (z * (0.0007936500793651 / x)) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(z / x) * Float64(z * y)) tmp = 0.0 if (y <= -2.4e+16) tmp = t_0; elseif (y <= 0.0008) tmp = Float64(z * Float64(z * Float64(0.0007936500793651 / x))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (z / x) * (z * y); tmp = 0.0; if (y <= -2.4e+16) tmp = t_0; elseif (y <= 0.0008) tmp = z * (z * (0.0007936500793651 / x)); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z / x), $MachinePrecision] * N[(z * y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.4e+16], t$95$0, If[LessEqual[y, 0.0008], N[(z * N[(z * N[(0.0007936500793651 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{z}{x} \cdot \left(z \cdot y\right)\\
\mathbf{if}\;y \leq -2.4 \cdot 10^{+16}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 0.0008:\\
\;\;\;\;z \cdot \left(z \cdot \frac{0.0007936500793651}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -2.4e16 or 8.00000000000000038e-4 < y Initial program 96.5%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified96.5%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.3%
Simplified52.3%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6454.2%
Applied egg-rr54.2%
if -2.4e16 < y < 8.00000000000000038e-4Initial program 93.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.7%
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-/.f6432.5%
Simplified32.5%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6432.5%
Simplified32.5%
Final simplification42.3%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (* y (/ z x)))))
(if (<= y -2.4e+16)
t_0
(if (<= y 0.0008) (* z (* z (/ 0.0007936500793651 x))) t_0))))
double code(double x, double y, double z) {
double t_0 = z * (y * (z / x));
double tmp;
if (y <= -2.4e+16) {
tmp = t_0;
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / 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 * (z / x))
if (y <= (-2.4d+16)) then
tmp = t_0
else if (y <= 0.0008d0) then
tmp = z * (z * (0.0007936500793651d0 / 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 * (z / x));
double tmp;
if (y <= -2.4e+16) {
tmp = t_0;
} else if (y <= 0.0008) {
tmp = z * (z * (0.0007936500793651 / x));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = z * (y * (z / x)) tmp = 0 if y <= -2.4e+16: tmp = t_0 elif y <= 0.0008: tmp = z * (z * (0.0007936500793651 / x)) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(z * Float64(y * Float64(z / x))) tmp = 0.0 if (y <= -2.4e+16) tmp = t_0; elseif (y <= 0.0008) tmp = Float64(z * Float64(z * Float64(0.0007936500793651 / x))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * (y * (z / x)); tmp = 0.0; if (y <= -2.4e+16) tmp = t_0; elseif (y <= 0.0008) tmp = z * (z * (0.0007936500793651 / x)); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(y * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.4e+16], t$95$0, If[LessEqual[y, 0.0008], N[(z * N[(z * N[(0.0007936500793651 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(y \cdot \frac{z}{x}\right)\\
\mathbf{if}\;y \leq -2.4 \cdot 10^{+16}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 0.0008:\\
\;\;\;\;z \cdot \left(z \cdot \frac{0.0007936500793651}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -2.4e16 or 8.00000000000000038e-4 < y Initial program 96.5%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified96.5%
Taylor expanded in y around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.3%
Simplified52.3%
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6453.3%
Applied egg-rr53.3%
if -2.4e16 < y < 8.00000000000000038e-4Initial program 93.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified93.7%
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-/.f6432.5%
Simplified32.5%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6432.5%
Simplified32.5%
Final simplification42.0%
(FPCore (x y z) :precision binary64 (let* ((t_0 (* z (+ y 0.0007936500793651)))) (if (<= x 9.5e+32) (/ (+ 0.083333333333333 (* z t_0)) x) (* z (/ t_0 x)))))
double code(double x, double y, double z) {
double t_0 = z * (y + 0.0007936500793651);
double tmp;
if (x <= 9.5e+32) {
tmp = (0.083333333333333 + (z * t_0)) / x;
} else {
tmp = z * (t_0 / 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) :: t_0
real(8) :: tmp
t_0 = z * (y + 0.0007936500793651d0)
if (x <= 9.5d+32) then
tmp = (0.083333333333333d0 + (z * t_0)) / x
else
tmp = z * (t_0 / x)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = z * (y + 0.0007936500793651);
double tmp;
if (x <= 9.5e+32) {
tmp = (0.083333333333333 + (z * t_0)) / x;
} else {
tmp = z * (t_0 / x);
}
return tmp;
}
def code(x, y, z): t_0 = z * (y + 0.0007936500793651) tmp = 0 if x <= 9.5e+32: tmp = (0.083333333333333 + (z * t_0)) / x else: tmp = z * (t_0 / x) return tmp
function code(x, y, z) t_0 = Float64(z * Float64(y + 0.0007936500793651)) tmp = 0.0 if (x <= 9.5e+32) tmp = Float64(Float64(0.083333333333333 + Float64(z * t_0)) / x); else tmp = Float64(z * Float64(t_0 / x)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * (y + 0.0007936500793651); tmp = 0.0; if (x <= 9.5e+32) tmp = (0.083333333333333 + (z * t_0)) / x; else tmp = z * (t_0 / x); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 9.5e+32], N[(N[(0.083333333333333 + N[(z * t$95$0), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(t$95$0 / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(y + 0.0007936500793651\right)\\
\mathbf{if}\;x \leq 9.5 \cdot 10^{+32}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot t\_0}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{t\_0}{x}\\
\end{array}
\end{array}
if x < 9.50000000000000006e32Initial program 99.7%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified99.7%
Taylor expanded in x around 0
Simplified99.7%
Taylor expanded in z around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6494.7%
Simplified94.7%
if 9.50000000000000006e32 < x Initial program 88.4%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified88.4%
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-/.f6429.4%
Simplified29.4%
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-+.f6429.5%
Simplified29.5%
Final simplification67.4%
(FPCore (x y z) :precision binary64 (* z (* z (/ 0.0007936500793651 x))))
double code(double x, double y, double z) {
return z * (z * (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 = z * (z * (0.0007936500793651d0 / x))
end function
public static double code(double x, double y, double z) {
return z * (z * (0.0007936500793651 / x));
}
def code(x, y, z): return z * (z * (0.0007936500793651 / x))
function code(x, y, z) return Float64(z * Float64(z * Float64(0.0007936500793651 / x))) end
function tmp = code(x, y, z) tmp = z * (z * (0.0007936500793651 / x)); end
code[x_, y_, z_] := N[(z * N[(z * N[(0.0007936500793651 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
z \cdot \left(z \cdot \frac{0.0007936500793651}{x}\right)
\end{array}
Initial program 95.0%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified95.0%
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-/.f6440.9%
Simplified40.9%
Taylor expanded in y around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6426.9%
Simplified26.9%
(FPCore (x y z) :precision binary64 (* z (/ -0.0027777777777778 x)))
double code(double x, double y, double z) {
return z * (-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 = z * ((-0.0027777777777778d0) / x)
end function
public static double code(double x, double y, double z) {
return z * (-0.0027777777777778 / x);
}
def code(x, y, z): return z * (-0.0027777777777778 / x)
function code(x, y, z) return Float64(z * Float64(-0.0027777777777778 / x)) end
function tmp = code(x, y, z) tmp = z * (-0.0027777777777778 / x); end
code[x_, y_, z_] := N[(z * N[(-0.0027777777777778 / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
z \cdot \frac{-0.0027777777777778}{x}
\end{array}
Initial program 95.0%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified95.0%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr95.0%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6467.3%
Simplified67.3%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6410.1%
Simplified10.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6410.1%
Applied egg-rr10.1%
Final simplification10.1%
(FPCore (x y z) :precision binary64 (* -0.0027777777777778 (/ z x)))
double code(double x, double y, double z) {
return -0.0027777777777778 * (z / 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.0027777777777778d0) * (z / x)
end function
public static double code(double x, double y, double z) {
return -0.0027777777777778 * (z / x);
}
def code(x, y, z): return -0.0027777777777778 * (z / x)
function code(x, y, z) return Float64(-0.0027777777777778 * Float64(z / x)) end
function tmp = code(x, y, z) tmp = -0.0027777777777778 * (z / x); end
code[x_, y_, z_] := N[(-0.0027777777777778 * N[(z / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.0027777777777778 \cdot \frac{z}{x}
\end{array}
Initial program 95.0%
associate-+l+N/A
associate-+l-N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
log-lowering-log.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
Simplified95.0%
+-commutativeN/A
associate-+r-N/A
associate-+l-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
Applied egg-rr95.0%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6467.3%
Simplified67.3%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6410.1%
Simplified10.1%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6410.1%
Applied egg-rr10.1%
(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 2024138
(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)))