
(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 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z)
:precision binary64
(+
(+ (- (* (- x 0.5) (log x)) x) 0.91893853320467)
(/
(+
(* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z)
0.083333333333333)
x)))
double code(double x, double y, double z) {
return ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x - 0.5d0) * log(x)) - x) + 0.91893853320467d0) + ((((((y + 0.0007936500793651d0) * z) - 0.0027777777777778d0) * z) + 0.083333333333333d0) / x)
end function
public static double code(double x, double y, double z) {
return ((((x - 0.5) * Math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x);
}
def code(x, y, z): return ((((x - 0.5) * math.log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x - 0.5) * log(x)) - x) + 0.91893853320467) + Float64(Float64(Float64(Float64(Float64(Float64(y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x)) end
function tmp = code(x, y, z) tmp = ((((x - 0.5) * log(x)) - x) + 0.91893853320467) + ((((((y + 0.0007936500793651) * z) - 0.0027777777777778) * z) + 0.083333333333333) / x); end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x - 0.5), $MachinePrecision] * N[Log[x], $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision] + 0.91893853320467), $MachinePrecision] + N[(N[(N[(N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * z), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] * z), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.91893853320467\right) + \frac{\left(\left(y + 0.0007936500793651\right) \cdot z - 0.0027777777777778\right) \cdot z + 0.083333333333333}{x}
\end{array}
(FPCore (x y z)
:precision binary64
(if (<= x 2e+30)
(+
(fma (+ x -0.5) (log x) (- 0.91893853320467 x))
(/
(fma
z
(fma (+ y 0.0007936500793651) z -0.0027777777777778)
0.083333333333333)
x))
(-
(+
0.91893853320467
(+
(* 0.083333333333333 (/ 1.0 x))
(+
(*
z
(+
(* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x)))
(* 0.0027777777777778 (/ -1.0 x))))
(* (log x) (- x 0.5)))))
x)))
double code(double x, double y, double z) {
double tmp;
if (x <= 2e+30) {
tmp = fma((x + -0.5), log(x), (0.91893853320467 - x)) + (fma(z, fma((y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333) / x);
} else {
tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (log(x) * (x - 0.5))))) - x;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 2e+30) tmp = Float64(fma(Float64(x + -0.5), log(x), Float64(0.91893853320467 - x)) + Float64(fma(z, fma(Float64(y + 0.0007936500793651), z, -0.0027777777777778), 0.083333333333333) / x)); else tmp = Float64(Float64(0.91893853320467 + Float64(Float64(0.083333333333333 * Float64(1.0 / x)) + Float64(Float64(z * Float64(Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x))) + Float64(0.0027777777777778 * Float64(-1.0 / x)))) + Float64(log(x) * Float64(x - 0.5))))) - x); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 2e+30], N[(N[(N[(x + -0.5), $MachinePrecision] * N[Log[x], $MachinePrecision] + N[(0.91893853320467 - x), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * z + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], N[(N[(0.91893853320467 + N[(N[(0.083333333333333 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0027777777777778 * N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2 \cdot 10^{+30}:\\
\;\;\;\;\mathsf{fma}\left(x + -0.5, \log x, 0.91893853320467 - x\right) + \frac{\mathsf{fma}\left(z, \mathsf{fma}\left(y + 0.0007936500793651, z, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(0.91893853320467 + \left(0.083333333333333 \cdot \frac{1}{x} + \left(z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right) + 0.0027777777777778 \cdot \frac{-1}{x}\right) + \log x \cdot \left(x - 0.5\right)\right)\right)\right) - x\\
\end{array}
\end{array}
if x < 2e30Initial program 99.8%
remove-double-neg99.8%
distribute-frac-neg299.8%
sub-neg99.8%
associate-+l+99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
if 2e30 < x Initial program 88.7%
remove-double-neg88.7%
distribute-frac-neg288.7%
sub-neg88.7%
associate-+l+88.7%
fma-define88.9%
sub-neg88.9%
metadata-eval88.9%
+-commutative88.9%
unsub-neg88.9%
distribute-frac-neg288.9%
remove-double-neg88.9%
Simplified88.9%
Taylor expanded in z around 0 99.5%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(if (<= x 6.5e-58)
(/
(fma
z
(fma z (+ y 0.0007936500793651) -0.0027777777777778)
0.083333333333333)
x)
(-
(+
0.91893853320467
(+
(* 0.083333333333333 (/ 1.0 x))
(+
(*
z
(+
(* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x)))
(* 0.0027777777777778 (/ -1.0 x))))
(* (log x) (- x 0.5)))))
x)))
double code(double x, double y, double z) {
double tmp;
if (x <= 6.5e-58) {
tmp = fma(z, fma(z, (y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x;
} else {
tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (log(x) * (x - 0.5))))) - x;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= 6.5e-58) tmp = Float64(fma(z, fma(z, Float64(y + 0.0007936500793651), -0.0027777777777778), 0.083333333333333) / x); else tmp = Float64(Float64(0.91893853320467 + Float64(Float64(0.083333333333333 * Float64(1.0 / x)) + Float64(Float64(z * Float64(Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x))) + Float64(0.0027777777777778 * Float64(-1.0 / x)))) + Float64(log(x) * Float64(x - 0.5))))) - x); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, 6.5e-58], N[(N[(z * N[(z * N[(y + 0.0007936500793651), $MachinePrecision] + -0.0027777777777778), $MachinePrecision] + 0.083333333333333), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.91893853320467 + N[(N[(0.083333333333333 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0027777777777778 * N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6.5 \cdot 10^{-58}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(z, y + 0.0007936500793651, -0.0027777777777778\right), 0.083333333333333\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(0.91893853320467 + \left(0.083333333333333 \cdot \frac{1}{x} + \left(z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right) + 0.0027777777777778 \cdot \frac{-1}{x}\right) + \log x \cdot \left(x - 0.5\right)\right)\right)\right) - x\\
\end{array}
\end{array}
if x < 6.49999999999999964e-58Initial program 99.8%
remove-double-neg99.8%
distribute-frac-neg299.8%
sub-neg99.8%
associate-+l+99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 99.8%
+-commutative99.8%
*-commutative99.8%
fmm-def99.8%
metadata-eval99.8%
+-commutative99.8%
fma-define99.8%
fma-define99.8%
+-commutative99.8%
*-commutative99.8%
fma-define99.8%
Simplified99.8%
if 6.49999999999999964e-58 < x Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.9%
Taylor expanded in z around 0 99.6%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(if (<= x 7e-58)
(/
(+
0.083333333333333
(* z (- (* z (* y (+ 1.0 (/ 0.0007936500793651 y)))) 0.0027777777777778)))
x)
(-
(+
0.91893853320467
(+
(* 0.083333333333333 (/ 1.0 x))
(+
(*
z
(+
(* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x)))
(* 0.0027777777777778 (/ -1.0 x))))
(* (log x) (- x 0.5)))))
x)))
double code(double x, double y, double z) {
double tmp;
if (x <= 7e-58) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (log(x) * (x - 0.5))))) - 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 <= 7d-58) then
tmp = (0.083333333333333d0 + (z * ((z * (y * (1.0d0 + (0.0007936500793651d0 / y)))) - 0.0027777777777778d0))) / x
else
tmp = (0.91893853320467d0 + ((0.083333333333333d0 * (1.0d0 / x)) + ((z * ((z * ((0.0007936500793651d0 * (1.0d0 / x)) + (y / x))) + (0.0027777777777778d0 * ((-1.0d0) / x)))) + (log(x) * (x - 0.5d0))))) - x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 7e-58) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (Math.log(x) * (x - 0.5))))) - x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 7e-58: tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x else: tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (math.log(x) * (x - 0.5))))) - x return tmp
function code(x, y, z) tmp = 0.0 if (x <= 7e-58) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y * Float64(1.0 + Float64(0.0007936500793651 / y)))) - 0.0027777777777778))) / x); else tmp = Float64(Float64(0.91893853320467 + Float64(Float64(0.083333333333333 * Float64(1.0 / x)) + Float64(Float64(z * Float64(Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x))) + Float64(0.0027777777777778 * Float64(-1.0 / x)))) + Float64(log(x) * Float64(x - 0.5))))) - x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 7e-58) tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x; else tmp = (0.91893853320467 + ((0.083333333333333 * (1.0 / x)) + ((z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) + (log(x) * (x - 0.5))))) - x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 7e-58], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y * N[(1.0 + N[(0.0007936500793651 / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.91893853320467 + N[(N[(0.083333333333333 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(z * N[(N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0027777777777778 * N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 7 \cdot 10^{-58}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y \cdot \left(1 + \frac{0.0007936500793651}{y}\right)\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(0.91893853320467 + \left(0.083333333333333 \cdot \frac{1}{x} + \left(z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right) + 0.0027777777777778 \cdot \frac{-1}{x}\right) + \log x \cdot \left(x - 0.5\right)\right)\right)\right) - x\\
\end{array}
\end{array}
if x < 6.9999999999999998e-58Initial program 99.8%
remove-double-neg99.8%
distribute-frac-neg299.8%
sub-neg99.8%
associate-+l+99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 99.8%
Taylor expanded in y around inf 99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
if 6.9999999999999998e-58 < x Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.9%
Taylor expanded in z around 0 99.6%
Final simplification99.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (* z (+ y 0.0007936500793651)) 0.0027777777777778))
(t_1 (* z t_0)))
(if (<= t_1 1e+307)
(+
(+ 0.91893853320467 (- (* (log x) (- x 0.5)) x))
(/ (+ 0.083333333333333 t_1) x))
(* z (/ t_0 x)))))
double code(double x, double y, double z) {
double t_0 = (z * (y + 0.0007936500793651)) - 0.0027777777777778;
double t_1 = z * t_0;
double tmp;
if (t_1 <= 1e+307) {
tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + ((0.083333333333333 + t_1) / 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) :: t_1
real(8) :: tmp
t_0 = (z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0
t_1 = z * t_0
if (t_1 <= 1d+307) then
tmp = (0.91893853320467d0 + ((log(x) * (x - 0.5d0)) - x)) + ((0.083333333333333d0 + t_1) / 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)) - 0.0027777777777778;
double t_1 = z * t_0;
double tmp;
if (t_1 <= 1e+307) {
tmp = (0.91893853320467 + ((Math.log(x) * (x - 0.5)) - x)) + ((0.083333333333333 + t_1) / x);
} else {
tmp = z * (t_0 / x);
}
return tmp;
}
def code(x, y, z): t_0 = (z * (y + 0.0007936500793651)) - 0.0027777777777778 t_1 = z * t_0 tmp = 0 if t_1 <= 1e+307: tmp = (0.91893853320467 + ((math.log(x) * (x - 0.5)) - x)) + ((0.083333333333333 + t_1) / x) else: tmp = z * (t_0 / x) return tmp
function code(x, y, z) t_0 = Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) t_1 = Float64(z * t_0) tmp = 0.0 if (t_1 <= 1e+307) tmp = Float64(Float64(0.91893853320467 + Float64(Float64(log(x) * Float64(x - 0.5)) - x)) + Float64(Float64(0.083333333333333 + t_1) / 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)) - 0.0027777777777778; t_1 = z * t_0; tmp = 0.0; if (t_1 <= 1e+307) tmp = (0.91893853320467 + ((log(x) * (x - 0.5)) - x)) + ((0.083333333333333 + t_1) / x); else tmp = z * (t_0 / x); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]}, Block[{t$95$1 = N[(z * t$95$0), $MachinePrecision]}, If[LessEqual[t$95$1, 1e+307], N[(N[(0.91893853320467 + N[(N[(N[Log[x], $MachinePrecision] * N[(x - 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision] + N[(N[(0.083333333333333 + t$95$1), $MachinePrecision] / x), $MachinePrecision]), $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) - 0.0027777777777778\\
t_1 := z \cdot t\_0\\
\mathbf{if}\;t\_1 \leq 10^{+307}:\\
\;\;\;\;\left(0.91893853320467 + \left(\log x \cdot \left(x - 0.5\right) - x\right)\right) + \frac{0.083333333333333 + t\_1}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{t\_0}{x}\\
\end{array}
\end{array}
if (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) < 9.99999999999999986e306Initial program 98.1%
if 9.99999999999999986e306 < (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) Initial program 84.2%
remove-double-neg84.2%
distribute-frac-neg284.2%
sub-neg84.2%
associate-+l+84.2%
fma-define84.2%
sub-neg84.2%
metadata-eval84.2%
+-commutative84.2%
unsub-neg84.2%
distribute-frac-neg284.2%
remove-double-neg84.2%
Simplified84.2%
Taylor expanded in z around inf 84.2%
Taylor expanded in z around 0 94.8%
sub-neg94.8%
distribute-rgt-in91.4%
associate-*r/91.4%
metadata-eval91.4%
associate-*l/91.4%
associate-*r/91.3%
associate-*l/91.4%
associate-/l*79.3%
distribute-rgt-out94.8%
associate-*r/94.8%
metadata-eval94.8%
distribute-neg-frac94.8%
metadata-eval94.8%
Simplified94.8%
Taylor expanded in x around 0 94.9%
Final simplification97.4%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (* z (+ y 0.0007936500793651)) 0.0027777777777778))
(t_1 (* z t_0)))
(if (<= t_1 1e+307)
(+ (/ (+ 0.083333333333333 t_1) x) (* x (+ (log x) -1.0)))
(* z (/ t_0 x)))))
double code(double x, double y, double z) {
double t_0 = (z * (y + 0.0007936500793651)) - 0.0027777777777778;
double t_1 = z * t_0;
double tmp;
if (t_1 <= 1e+307) {
tmp = ((0.083333333333333 + t_1) / x) + (x * (log(x) + -1.0));
} 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) :: t_1
real(8) :: tmp
t_0 = (z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0
t_1 = z * t_0
if (t_1 <= 1d+307) then
tmp = ((0.083333333333333d0 + t_1) / x) + (x * (log(x) + (-1.0d0)))
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)) - 0.0027777777777778;
double t_1 = z * t_0;
double tmp;
if (t_1 <= 1e+307) {
tmp = ((0.083333333333333 + t_1) / x) + (x * (Math.log(x) + -1.0));
} else {
tmp = z * (t_0 / x);
}
return tmp;
}
def code(x, y, z): t_0 = (z * (y + 0.0007936500793651)) - 0.0027777777777778 t_1 = z * t_0 tmp = 0 if t_1 <= 1e+307: tmp = ((0.083333333333333 + t_1) / x) + (x * (math.log(x) + -1.0)) else: tmp = z * (t_0 / x) return tmp
function code(x, y, z) t_0 = Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) t_1 = Float64(z * t_0) tmp = 0.0 if (t_1 <= 1e+307) tmp = Float64(Float64(Float64(0.083333333333333 + t_1) / x) + Float64(x * Float64(log(x) + -1.0))); 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)) - 0.0027777777777778; t_1 = z * t_0; tmp = 0.0; if (t_1 <= 1e+307) tmp = ((0.083333333333333 + t_1) / x) + (x * (log(x) + -1.0)); else tmp = z * (t_0 / x); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]}, Block[{t$95$1 = N[(z * t$95$0), $MachinePrecision]}, If[LessEqual[t$95$1, 1e+307], N[(N[(N[(0.083333333333333 + t$95$1), $MachinePrecision] / x), $MachinePrecision] + N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $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) - 0.0027777777777778\\
t_1 := z \cdot t\_0\\
\mathbf{if}\;t\_1 \leq 10^{+307}:\\
\;\;\;\;\frac{0.083333333333333 + t\_1}{x} + x \cdot \left(\log x + -1\right)\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{t\_0}{x}\\
\end{array}
\end{array}
if (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) < 9.99999999999999986e306Initial program 98.1%
flip--79.4%
metadata-eval79.4%
metadata-eval79.4%
clear-num79.3%
fmm-def79.3%
metadata-eval79.3%
metadata-eval79.3%
Applied egg-rr79.3%
associate-+l-79.3%
associate-*l/79.3%
*-un-lft-identity79.3%
clear-num79.4%
metadata-eval79.4%
metadata-eval79.4%
fmm-def79.4%
*-un-lft-identity79.4%
fma-define79.4%
metadata-eval79.4%
fmm-def79.4%
*-un-lft-identity79.4%
flip-+98.1%
Applied egg-rr98.1%
Taylor expanded in x around inf 97.5%
sub-neg97.5%
mul-1-neg97.5%
log-rec97.5%
remove-double-neg97.5%
metadata-eval97.5%
Simplified97.5%
if 9.99999999999999986e306 < (*.f64 (-.f64 (*.f64 (+.f64 y #s(literal 7936500793651/10000000000000000 binary64)) z) #s(literal 13888888888889/5000000000000000 binary64)) z) Initial program 84.2%
remove-double-neg84.2%
distribute-frac-neg284.2%
sub-neg84.2%
associate-+l+84.2%
fma-define84.2%
sub-neg84.2%
metadata-eval84.2%
+-commutative84.2%
unsub-neg84.2%
distribute-frac-neg284.2%
remove-double-neg84.2%
Simplified84.2%
Taylor expanded in z around inf 84.2%
Taylor expanded in z around 0 94.8%
sub-neg94.8%
distribute-rgt-in91.4%
associate-*r/91.4%
metadata-eval91.4%
associate-*l/91.4%
associate-*r/91.3%
associate-*l/91.4%
associate-/l*79.3%
distribute-rgt-out94.8%
associate-*r/94.8%
metadata-eval94.8%
distribute-neg-frac94.8%
metadata-eval94.8%
Simplified94.8%
Taylor expanded in x around 0 94.9%
Final simplification96.9%
(FPCore (x y z)
:precision binary64
(if (<= x 5.6e+18)
(/
(+
0.083333333333333
(* z (- (* z (* y (+ 1.0 (/ 0.0007936500793651 y)))) 0.0027777777777778)))
x)
(* x (+ (log x) -1.0))))
double code(double x, double y, double z) {
double tmp;
if (x <= 5.6e+18) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = x * (log(x) + -1.0);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= 5.6d+18) then
tmp = (0.083333333333333d0 + (z * ((z * (y * (1.0d0 + (0.0007936500793651d0 / y)))) - 0.0027777777777778d0))) / x
else
tmp = x * (log(x) + (-1.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 5.6e+18) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = x * (Math.log(x) + -1.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 5.6e+18: tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x else: tmp = x * (math.log(x) + -1.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 5.6e+18) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y * Float64(1.0 + Float64(0.0007936500793651 / y)))) - 0.0027777777777778))) / x); else tmp = Float64(x * Float64(log(x) + -1.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 5.6e+18) tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x; else tmp = x * (log(x) + -1.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 5.6e+18], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y * N[(1.0 + N[(0.0007936500793651 / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(x * N[(N[Log[x], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.6 \cdot 10^{+18}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y \cdot \left(1 + \frac{0.0007936500793651}{y}\right)\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log x + -1\right)\\
\end{array}
\end{array}
if x < 5.6e18Initial program 99.8%
remove-double-neg99.8%
distribute-frac-neg299.8%
sub-neg99.8%
associate-+l+99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 98.7%
Taylor expanded in y around inf 98.7%
associate-*r/98.8%
metadata-eval98.8%
Simplified98.8%
if 5.6e18 < x Initial program 89.4%
remove-double-neg89.4%
distribute-frac-neg289.4%
sub-neg89.4%
associate-+l+89.4%
fma-define89.6%
sub-neg89.6%
metadata-eval89.6%
+-commutative89.6%
unsub-neg89.6%
distribute-frac-neg289.6%
remove-double-neg89.6%
Simplified89.6%
Taylor expanded in x around inf 71.2%
sub-neg71.2%
mul-1-neg71.2%
log-rec71.2%
remove-double-neg71.2%
metadata-eval71.2%
Simplified71.2%
(FPCore (x y z)
:precision binary64
(if (<= x 9.2e-58)
(/
(+
0.083333333333333
(* z (- (* z (* y (+ 1.0 (/ 0.0007936500793651 y)))) 0.0027777777777778)))
x)
(+
(* 0.083333333333333 (/ 1.0 x))
(*
z
(+
(* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x)))
(* 0.0027777777777778 (/ -1.0 x)))))))
double code(double x, double y, double z) {
double tmp;
if (x <= 9.2e-58) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = (0.083333333333333 * (1.0 / x)) + (z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.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) :: tmp
if (x <= 9.2d-58) then
tmp = (0.083333333333333d0 + (z * ((z * (y * (1.0d0 + (0.0007936500793651d0 / y)))) - 0.0027777777777778d0))) / x
else
tmp = (0.083333333333333d0 * (1.0d0 / x)) + (z * ((z * ((0.0007936500793651d0 * (1.0d0 / x)) + (y / x))) + (0.0027777777777778d0 * ((-1.0d0) / x))))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 9.2e-58) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = (0.083333333333333 * (1.0 / x)) + (z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x))));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 9.2e-58: tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x else: tmp = (0.083333333333333 * (1.0 / x)) + (z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 9.2e-58) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y * Float64(1.0 + Float64(0.0007936500793651 / y)))) - 0.0027777777777778))) / x); else tmp = Float64(Float64(0.083333333333333 * Float64(1.0 / x)) + Float64(z * Float64(Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x))) + Float64(0.0027777777777778 * Float64(-1.0 / x))))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 9.2e-58) tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x; else tmp = (0.083333333333333 * (1.0 / x)) + (z * ((z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) + (0.0027777777777778 * (-1.0 / x)))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 9.2e-58], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y * N[(1.0 + N[(0.0007936500793651 / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(0.083333333333333 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(z * N[(N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0027777777777778 * N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.2 \cdot 10^{-58}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y \cdot \left(1 + \frac{0.0007936500793651}{y}\right)\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;0.083333333333333 \cdot \frac{1}{x} + z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right) + 0.0027777777777778 \cdot \frac{-1}{x}\right)\\
\end{array}
\end{array}
if x < 9.1999999999999995e-58Initial program 99.8%
remove-double-neg99.8%
distribute-frac-neg299.8%
sub-neg99.8%
associate-+l+99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 99.8%
Taylor expanded in y around inf 99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
if 9.1999999999999995e-58 < x Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.9%
Taylor expanded in x around 0 39.9%
Taylor expanded in z around 0 45.2%
Final simplification67.4%
(FPCore (x y z) :precision binary64 (if (or (<= z -13500000.0) (not (<= z 11.5))) (* z (+ (* (+ y 0.0007936500793651) (/ z x)) (/ -0.0027777777777778 x))) (/ (+ 0.083333333333333 (* z (- (* z y) 0.0027777777777778))) x)))
double code(double x, double y, double z) {
double tmp;
if ((z <= -13500000.0) || !(z <= 11.5)) {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x));
} else {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / 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 <= (-13500000.0d0)) .or. (.not. (z <= 11.5d0))) then
tmp = z * (((y + 0.0007936500793651d0) * (z / x)) + ((-0.0027777777777778d0) / x))
else
tmp = (0.083333333333333d0 + (z * ((z * y) - 0.0027777777777778d0))) / x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -13500000.0) || !(z <= 11.5)) {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x));
} else {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -13500000.0) or not (z <= 11.5): tmp = z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x)) else: tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -13500000.0) || !(z <= 11.5)) tmp = Float64(z * Float64(Float64(Float64(y + 0.0007936500793651) * Float64(z / x)) + Float64(-0.0027777777777778 / x))); else tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * y) - 0.0027777777777778))) / x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -13500000.0) || ~((z <= 11.5))) tmp = z * (((y + 0.0007936500793651) * (z / x)) + (-0.0027777777777778 / x)); else tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -13500000.0], N[Not[LessEqual[z, 11.5]], $MachinePrecision]], N[(z * N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision] + N[(-0.0027777777777778 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.083333333333333 + N[(z * N[(N[(z * y), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -13500000 \lor \neg \left(z \leq 11.5\right):\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x} + \frac{-0.0027777777777778}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot y - 0.0027777777777778\right)}{x}\\
\end{array}
\end{array}
if z < -1.35e7 or 11.5 < z Initial program 91.2%
remove-double-neg91.2%
distribute-frac-neg291.2%
sub-neg91.2%
associate-+l+91.2%
fma-define91.3%
sub-neg91.3%
metadata-eval91.3%
+-commutative91.3%
unsub-neg91.3%
distribute-frac-neg291.3%
remove-double-neg91.3%
Simplified91.3%
Taylor expanded in z around inf 78.5%
Taylor expanded in z around 0 83.5%
sub-neg83.5%
distribute-rgt-in75.7%
associate-*r/75.7%
metadata-eval75.7%
associate-*l/75.7%
associate-*r/75.7%
associate-*l/75.1%
associate-/l*70.0%
distribute-rgt-out83.5%
associate-*r/83.5%
metadata-eval83.5%
distribute-neg-frac83.5%
metadata-eval83.5%
Simplified83.5%
if -1.35e7 < z < 11.5Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.9%
Taylor expanded in y around inf 46.9%
*-commutative46.9%
Simplified46.9%
Final simplification67.1%
(FPCore (x y z)
:precision binary64
(if (<= z -13500000.0)
(* z (/ (- (* z (+ y 0.0007936500793651)) 0.0027777777777778) x))
(if (<= z 0.108)
(/ (+ 0.083333333333333 (* z (- (* z y) 0.0027777777777778))) x)
(* z (* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x)))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -13500000.0) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 0.108) {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x)));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (z <= (-13500000.0d0)) then
tmp = z * (((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0) / x)
else if (z <= 0.108d0) then
tmp = (0.083333333333333d0 + (z * ((z * y) - 0.0027777777777778d0))) / x
else
tmp = z * (z * ((0.0007936500793651d0 * (1.0d0 / x)) + (y / x)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -13500000.0) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 0.108) {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -13500000.0: tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x) elif z <= 0.108: tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x else: tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -13500000.0) tmp = Float64(z * Float64(Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) / x)); elseif (z <= 0.108) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * y) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -13500000.0) tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x); elseif (z <= 0.108) tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x; else tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -13500000.0], N[(z * N[(N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 0.108], N[(N[(0.083333333333333 + N[(z * N[(N[(z * y), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -13500000:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778}{x}\\
\mathbf{elif}\;z \leq 0.108:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot y - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if z < -1.35e7Initial program 90.8%
remove-double-neg90.8%
distribute-frac-neg290.8%
sub-neg90.8%
associate-+l+90.8%
fma-define90.9%
sub-neg90.9%
metadata-eval90.9%
+-commutative90.9%
unsub-neg90.9%
distribute-frac-neg290.9%
remove-double-neg90.9%
Simplified90.9%
Taylor expanded in z around inf 78.6%
Taylor expanded in z around 0 82.9%
sub-neg82.9%
distribute-rgt-in78.1%
associate-*r/78.0%
metadata-eval78.0%
associate-*l/78.0%
associate-*r/78.1%
associate-*l/78.1%
associate-/l*71.6%
distribute-rgt-out82.9%
associate-*r/82.9%
metadata-eval82.9%
distribute-neg-frac82.9%
metadata-eval82.9%
Simplified82.9%
Taylor expanded in x around 0 82.8%
if -1.35e7 < z < 0.107999999999999999Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.8%
Taylor expanded in y around inf 46.8%
*-commutative46.8%
Simplified46.8%
if 0.107999999999999999 < z Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.8%
Taylor expanded in z around inf 77.7%
Taylor expanded in z around 0 83.2%
sub-neg83.2%
distribute-rgt-in73.3%
associate-*r/73.4%
metadata-eval73.4%
associate-*l/73.3%
associate-*r/73.3%
associate-*l/72.1%
associate-/l*68.4%
distribute-rgt-out83.2%
associate-*r/83.2%
metadata-eval83.2%
distribute-neg-frac83.2%
metadata-eval83.2%
Simplified83.2%
Taylor expanded in z around inf 82.9%
Final simplification66.9%
(FPCore (x y z)
:precision binary64
(if (<= x 9.2e+76)
(/
(+
0.083333333333333
(* z (- (* z (* y (+ 1.0 (/ 0.0007936500793651 y)))) 0.0027777777777778)))
x)
(*
z
(+
(* (+ y 0.0007936500793651) (/ z x))
(* (/ 1.0 x) 0.0027777777777778)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 9.2e+76) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778));
}
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 <= 9.2d+76) then
tmp = (0.083333333333333d0 + (z * ((z * (y * (1.0d0 + (0.0007936500793651d0 / y)))) - 0.0027777777777778d0))) / x
else
tmp = z * (((y + 0.0007936500793651d0) * (z / x)) + ((1.0d0 / x) * 0.0027777777777778d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 9.2e+76) {
tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x;
} else {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 9.2e+76: tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x else: tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 9.2e+76) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y * Float64(1.0 + Float64(0.0007936500793651 / y)))) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(Float64(Float64(y + 0.0007936500793651) * Float64(z / x)) + Float64(Float64(1.0 / x) * 0.0027777777777778))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 9.2e+76) tmp = (0.083333333333333 + (z * ((z * (y * (1.0 + (0.0007936500793651 / y)))) - 0.0027777777777778))) / x; else tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 9.2e+76], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y * N[(1.0 + N[(0.0007936500793651 / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / x), $MachinePrecision] * 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.2 \cdot 10^{+76}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y \cdot \left(1 + \frac{0.0007936500793651}{y}\right)\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x} + \frac{1}{x} \cdot 0.0027777777777778\right)\\
\end{array}
\end{array}
if x < 9.20000000000000005e76Initial program 99.7%
remove-double-neg99.7%
distribute-frac-neg299.7%
sub-neg99.7%
associate-+l+99.7%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 91.2%
Taylor expanded in y around inf 91.2%
associate-*r/91.2%
metadata-eval91.2%
Simplified91.2%
if 9.20000000000000005e76 < x Initial program 87.2%
remove-double-neg87.2%
distribute-frac-neg287.2%
sub-neg87.2%
associate-+l+87.2%
fma-define87.4%
sub-neg87.4%
metadata-eval87.4%
+-commutative87.4%
unsub-neg87.4%
distribute-frac-neg287.4%
remove-double-neg87.4%
Simplified87.4%
Taylor expanded in z around inf 20.7%
Taylor expanded in z around 0 28.0%
sub-neg28.0%
distribute-rgt-in28.0%
associate-*r/28.0%
metadata-eval28.0%
associate-*l/28.0%
associate-*r/27.9%
associate-*l/27.0%
associate-/l*28.0%
distribute-rgt-out28.0%
associate-*r/28.0%
metadata-eval28.0%
distribute-neg-frac28.0%
metadata-eval28.0%
Simplified28.0%
add-sqr-sqrt0.0%
sqrt-unprod28.0%
frac-times28.0%
metadata-eval28.0%
metadata-eval28.0%
frac-times28.0%
un-div-inv28.0%
un-div-inv28.0%
sqrt-unprod28.0%
add-sqr-sqrt28.0%
*-commutative28.0%
Applied egg-rr28.0%
Final simplification67.3%
(FPCore (x y z)
:precision binary64
(if (<= z -13500000.0)
(* z (/ (- (* z (+ y 0.0007936500793651)) 0.0027777777777778) x))
(if (<= z 0.108)
(/ (+ 0.083333333333333 (* z (- (* z y) 0.0027777777777778))) x)
(* z (* (+ y 0.0007936500793651) (/ z x))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -13500000.0) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 0.108) {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / 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 <= (-13500000.0d0)) then
tmp = z * (((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0) / x)
else if (z <= 0.108d0) then
tmp = (0.083333333333333d0 + (z * ((z * y) - 0.0027777777777778d0))) / x
else
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -13500000.0) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 0.108) {
tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / x));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -13500000.0: tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x) elif z <= 0.108: tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x else: tmp = z * ((y + 0.0007936500793651) * (z / x)) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -13500000.0) tmp = Float64(z * Float64(Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) / x)); elseif (z <= 0.108) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * y) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -13500000.0) tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x); elseif (z <= 0.108) tmp = (0.083333333333333 + (z * ((z * y) - 0.0027777777777778))) / x; else tmp = z * ((y + 0.0007936500793651) * (z / x)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -13500000.0], N[(z * N[(N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 0.108], N[(N[(0.083333333333333 + N[(z * N[(N[(z * y), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -13500000:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778}{x}\\
\mathbf{elif}\;z \leq 0.108:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot y - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\end{array}
\end{array}
if z < -1.35e7Initial program 90.8%
remove-double-neg90.8%
distribute-frac-neg290.8%
sub-neg90.8%
associate-+l+90.8%
fma-define90.9%
sub-neg90.9%
metadata-eval90.9%
+-commutative90.9%
unsub-neg90.9%
distribute-frac-neg290.9%
remove-double-neg90.9%
Simplified90.9%
Taylor expanded in z around inf 78.6%
Taylor expanded in z around 0 82.9%
sub-neg82.9%
distribute-rgt-in78.1%
associate-*r/78.0%
metadata-eval78.0%
associate-*l/78.0%
associate-*r/78.1%
associate-*l/78.1%
associate-/l*71.6%
distribute-rgt-out82.9%
associate-*r/82.9%
metadata-eval82.9%
distribute-neg-frac82.9%
metadata-eval82.9%
Simplified82.9%
Taylor expanded in x around 0 82.8%
if -1.35e7 < z < 0.107999999999999999Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.8%
Taylor expanded in y around inf 46.8%
*-commutative46.8%
Simplified46.8%
if 0.107999999999999999 < z Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.8%
Taylor expanded in z around inf 77.7%
Taylor expanded in z around inf 77.3%
*-commutative77.3%
associate-*r/77.4%
metadata-eval77.4%
unpow277.4%
associate-*r*82.9%
*-commutative82.9%
distribute-rgt-in73.0%
associate-*l/73.0%
associate-*r/72.9%
associate-*l/71.8%
associate-/l*68.0%
distribute-rgt-out82.8%
Simplified82.8%
Final simplification66.9%
(FPCore (x y z)
:precision binary64
(if (<= z -9e-80)
(* z (/ (- (* z (+ y 0.0007936500793651)) 0.0027777777777778) x))
(if (<= z 1e-5)
(/
(+
0.083333333333333
(* z (- (* z 0.0007936500793651) 0.0027777777777778)))
x)
(* z (* (+ y 0.0007936500793651) (/ z x))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -9e-80) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 1e-5) {
tmp = (0.083333333333333 + (z * ((z * 0.0007936500793651) - 0.0027777777777778))) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / 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 <= (-9d-80)) then
tmp = z * (((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0) / x)
else if (z <= 1d-5) then
tmp = (0.083333333333333d0 + (z * ((z * 0.0007936500793651d0) - 0.0027777777777778d0))) / x
else
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -9e-80) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 1e-5) {
tmp = (0.083333333333333 + (z * ((z * 0.0007936500793651) - 0.0027777777777778))) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / x));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -9e-80: tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x) elif z <= 1e-5: tmp = (0.083333333333333 + (z * ((z * 0.0007936500793651) - 0.0027777777777778))) / x else: tmp = z * ((y + 0.0007936500793651) * (z / x)) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -9e-80) tmp = Float64(z * Float64(Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) / x)); elseif (z <= 1e-5) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * 0.0007936500793651) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -9e-80) tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x); elseif (z <= 1e-5) tmp = (0.083333333333333 + (z * ((z * 0.0007936500793651) - 0.0027777777777778))) / x; else tmp = z * ((y + 0.0007936500793651) * (z / x)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -9e-80], N[(z * N[(N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1e-5], N[(N[(0.083333333333333 + N[(z * N[(N[(z * 0.0007936500793651), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9 \cdot 10^{-80}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778}{x}\\
\mathbf{elif}\;z \leq 10^{-5}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot 0.0007936500793651 - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\end{array}
\end{array}
if z < -9.0000000000000006e-80Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.7%
sub-neg91.7%
metadata-eval91.7%
+-commutative91.7%
unsub-neg91.7%
distribute-frac-neg291.7%
remove-double-neg91.7%
Simplified91.8%
Taylor expanded in z around inf 72.6%
Taylor expanded in z around 0 76.4%
sub-neg76.4%
distribute-rgt-in72.1%
associate-*r/72.0%
metadata-eval72.0%
associate-*l/72.0%
associate-*r/72.1%
associate-*l/73.4%
associate-/l*67.6%
distribute-rgt-out77.8%
associate-*r/77.8%
metadata-eval77.8%
distribute-neg-frac77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in x around 0 77.7%
if -9.0000000000000006e-80 < z < 1.00000000000000008e-5Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.0%
Taylor expanded in y around 0 40.3%
*-commutative40.3%
Simplified40.3%
if 1.00000000000000008e-5 < z Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.8%
Taylor expanded in z around inf 77.7%
Taylor expanded in z around inf 77.3%
*-commutative77.3%
associate-*r/77.4%
metadata-eval77.4%
unpow277.4%
associate-*r*82.9%
*-commutative82.9%
distribute-rgt-in73.0%
associate-*l/73.0%
associate-*r/72.9%
associate-*l/71.8%
associate-/l*68.0%
distribute-rgt-out82.8%
Simplified82.8%
Final simplification63.8%
(FPCore (x y z)
:precision binary64
(if (<= x 5e+76)
(/
(+
0.083333333333333
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))
x)
(*
z
(+
(* (+ y 0.0007936500793651) (/ z x))
(* (/ 1.0 x) 0.0027777777777778)))))
double code(double x, double y, double z) {
double tmp;
if (x <= 5e+76) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x;
} else {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778));
}
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 <= 5d+76) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0))) / x
else
tmp = z * (((y + 0.0007936500793651d0) * (z / x)) + ((1.0d0 / x) * 0.0027777777777778d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 5e+76) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x;
} else {
tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 5e+76: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x else: tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 5e+76) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(Float64(Float64(y + 0.0007936500793651) * Float64(z / x)) + Float64(Float64(1.0 / x) * 0.0027777777777778))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 5e+76) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x; else tmp = z * (((y + 0.0007936500793651) * (z / x)) + ((1.0 / x) * 0.0027777777777778)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 5e+76], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / x), $MachinePrecision] * 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5 \cdot 10^{+76}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x} + \frac{1}{x} \cdot 0.0027777777777778\right)\\
\end{array}
\end{array}
if x < 4.99999999999999991e76Initial program 99.7%
remove-double-neg99.7%
distribute-frac-neg299.7%
sub-neg99.7%
associate-+l+99.7%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 91.2%
if 4.99999999999999991e76 < x Initial program 87.2%
remove-double-neg87.2%
distribute-frac-neg287.2%
sub-neg87.2%
associate-+l+87.2%
fma-define87.4%
sub-neg87.4%
metadata-eval87.4%
+-commutative87.4%
unsub-neg87.4%
distribute-frac-neg287.4%
remove-double-neg87.4%
Simplified87.4%
Taylor expanded in z around inf 20.7%
Taylor expanded in z around 0 28.0%
sub-neg28.0%
distribute-rgt-in28.0%
associate-*r/28.0%
metadata-eval28.0%
associate-*l/28.0%
associate-*r/27.9%
associate-*l/27.0%
associate-/l*28.0%
distribute-rgt-out28.0%
associate-*r/28.0%
metadata-eval28.0%
distribute-neg-frac28.0%
metadata-eval28.0%
Simplified28.0%
add-sqr-sqrt0.0%
sqrt-unprod28.0%
frac-times28.0%
metadata-eval28.0%
metadata-eval28.0%
frac-times28.0%
un-div-inv28.0%
un-div-inv28.0%
sqrt-unprod28.0%
add-sqr-sqrt28.0%
*-commutative28.0%
Applied egg-rr28.0%
Final simplification67.3%
(FPCore (x y z) :precision binary64 (if (or (<= z -9e-80) (not (<= z 3.1e-6))) (* z (* (+ y 0.0007936500793651) (/ z x))) (/ (+ 0.083333333333333 (* z -0.0027777777777778)) x)))
double code(double x, double y, double z) {
double tmp;
if ((z <= -9e-80) || !(z <= 3.1e-6)) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / 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 <= (-9d-80)) .or. (.not. (z <= 3.1d-6))) then
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
else
tmp = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -9e-80) || !(z <= 3.1e-6)) {
tmp = z * ((y + 0.0007936500793651) * (z / x));
} else {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -9e-80) or not (z <= 3.1e-6): tmp = z * ((y + 0.0007936500793651) * (z / x)) else: tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -9e-80) || !(z <= 3.1e-6)) tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); else tmp = Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -9e-80) || ~((z <= 3.1e-6))) tmp = z * ((y + 0.0007936500793651) * (z / x)); else tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -9e-80], N[Not[LessEqual[z, 3.1e-6]], $MachinePrecision]], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9 \cdot 10^{-80} \lor \neg \left(z \leq 3.1 \cdot 10^{-6}\right):\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}\\
\end{array}
\end{array}
if z < -9.0000000000000006e-80 or 3.1e-6 < z Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.8%
Taylor expanded in z around inf 75.3%
Taylor expanded in z around inf 74.9%
*-commutative74.9%
associate-*r/74.9%
metadata-eval74.9%
unpow274.9%
associate-*r*79.6%
*-commutative79.6%
distribute-rgt-in72.3%
associate-*l/72.3%
associate-*r/72.2%
associate-*l/72.2%
associate-/l*67.5%
distribute-rgt-out80.2%
Simplified80.2%
if -9.0000000000000006e-80 < z < 3.1e-6Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.0%
Taylor expanded in z around 0 40.3%
*-commutative40.3%
Simplified40.3%
Final simplification63.7%
(FPCore (x y z)
:precision binary64
(if (<= z -3.8e-80)
(* z (/ (- (* z (+ y 0.0007936500793651)) 0.0027777777777778) x))
(if (<= z 9.5e-6)
(/ (+ 0.083333333333333 (* z -0.0027777777777778)) x)
(* z (* (+ y 0.0007936500793651) (/ z x))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -3.8e-80) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 9.5e-6) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / 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 <= (-3.8d-80)) then
tmp = z * (((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0) / x)
else if (z <= 9.5d-6) then
tmp = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / x
else
tmp = z * ((y + 0.0007936500793651d0) * (z / x))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -3.8e-80) {
tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x);
} else if (z <= 9.5e-6) {
tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x;
} else {
tmp = z * ((y + 0.0007936500793651) * (z / x));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -3.8e-80: tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x) elif z <= 9.5e-6: tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x else: tmp = z * ((y + 0.0007936500793651) * (z / x)) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -3.8e-80) tmp = Float64(z * Float64(Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778) / x)); elseif (z <= 9.5e-6) tmp = Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / x); else tmp = Float64(z * Float64(Float64(y + 0.0007936500793651) * Float64(z / x))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -3.8e-80) tmp = z * (((z * (y + 0.0007936500793651)) - 0.0027777777777778) / x); elseif (z <= 9.5e-6) tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; else tmp = z * ((y + 0.0007936500793651) * (z / x)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -3.8e-80], N[(z * N[(N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 9.5e-6], N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(N[(y + 0.0007936500793651), $MachinePrecision] * N[(z / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.8 \cdot 10^{-80}:\\
\;\;\;\;z \cdot \frac{z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778}{x}\\
\mathbf{elif}\;z \leq 9.5 \cdot 10^{-6}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(\left(y + 0.0007936500793651\right) \cdot \frac{z}{x}\right)\\
\end{array}
\end{array}
if z < -3.79999999999999967e-80Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.7%
sub-neg91.7%
metadata-eval91.7%
+-commutative91.7%
unsub-neg91.7%
distribute-frac-neg291.7%
remove-double-neg91.7%
Simplified91.8%
Taylor expanded in z around inf 72.6%
Taylor expanded in z around 0 76.4%
sub-neg76.4%
distribute-rgt-in72.1%
associate-*r/72.0%
metadata-eval72.0%
associate-*l/72.0%
associate-*r/72.1%
associate-*l/73.4%
associate-/l*67.6%
distribute-rgt-out77.8%
associate-*r/77.8%
metadata-eval77.8%
distribute-neg-frac77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in x around 0 77.7%
if -3.79999999999999967e-80 < z < 9.5000000000000005e-6Initial program 99.6%
remove-double-neg99.6%
distribute-frac-neg299.6%
sub-neg99.6%
associate-+l+99.6%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
+-commutative99.7%
unsub-neg99.7%
distribute-frac-neg299.7%
remove-double-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 46.0%
Taylor expanded in z around 0 40.3%
*-commutative40.3%
Simplified40.3%
if 9.5000000000000005e-6 < z Initial program 91.7%
remove-double-neg91.7%
distribute-frac-neg291.7%
sub-neg91.7%
associate-+l+91.7%
fma-define91.8%
sub-neg91.8%
metadata-eval91.8%
+-commutative91.8%
unsub-neg91.8%
distribute-frac-neg291.8%
remove-double-neg91.8%
Simplified91.8%
Taylor expanded in z around inf 77.7%
Taylor expanded in z around inf 77.3%
*-commutative77.3%
associate-*r/77.4%
metadata-eval77.4%
unpow277.4%
associate-*r*82.9%
*-commutative82.9%
distribute-rgt-in73.0%
associate-*l/73.0%
associate-*r/72.9%
associate-*l/71.8%
associate-/l*68.0%
distribute-rgt-out82.8%
Simplified82.8%
Final simplification63.8%
(FPCore (x y z)
:precision binary64
(if (<= x 6e+76)
(/
(+
0.083333333333333
(* z (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))
x)
(* z (* z (+ (* 0.0007936500793651 (/ 1.0 x)) (/ y x))))))
double code(double x, double y, double z) {
double tmp;
if (x <= 6e+76) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x)));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= 6d+76) then
tmp = (0.083333333333333d0 + (z * ((z * (y + 0.0007936500793651d0)) - 0.0027777777777778d0))) / x
else
tmp = z * (z * ((0.0007936500793651d0 * (1.0d0 / x)) + (y / x)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= 6e+76) {
tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x;
} else {
tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= 6e+76: tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x else: tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x))) return tmp
function code(x, y, z) tmp = 0.0 if (x <= 6e+76) tmp = Float64(Float64(0.083333333333333 + Float64(z * Float64(Float64(z * Float64(y + 0.0007936500793651)) - 0.0027777777777778))) / x); else tmp = Float64(z * Float64(z * Float64(Float64(0.0007936500793651 * Float64(1.0 / x)) + Float64(y / x)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= 6e+76) tmp = (0.083333333333333 + (z * ((z * (y + 0.0007936500793651)) - 0.0027777777777778))) / x; else tmp = z * (z * ((0.0007936500793651 * (1.0 / x)) + (y / x))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, 6e+76], N[(N[(0.083333333333333 + N[(z * N[(N[(z * N[(y + 0.0007936500793651), $MachinePrecision]), $MachinePrecision] - 0.0027777777777778), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(z * N[(z * N[(N[(0.0007936500793651 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision] + N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6 \cdot 10^{+76}:\\
\;\;\;\;\frac{0.083333333333333 + z \cdot \left(z \cdot \left(y + 0.0007936500793651\right) - 0.0027777777777778\right)}{x}\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(z \cdot \left(0.0007936500793651 \cdot \frac{1}{x} + \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if x < 5.9999999999999996e76Initial program 99.7%
remove-double-neg99.7%
distribute-frac-neg299.7%
sub-neg99.7%
associate-+l+99.7%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
+-commutative99.8%
unsub-neg99.8%
distribute-frac-neg299.8%
remove-double-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 91.2%
if 5.9999999999999996e76 < x Initial program 87.2%
remove-double-neg87.2%
distribute-frac-neg287.2%
sub-neg87.2%
associate-+l+87.2%
fma-define87.4%
sub-neg87.4%
metadata-eval87.4%
+-commutative87.4%
unsub-neg87.4%
distribute-frac-neg287.4%
remove-double-neg87.4%
Simplified87.4%
Taylor expanded in z around inf 20.7%
Taylor expanded in z around 0 28.0%
sub-neg28.0%
distribute-rgt-in28.0%
associate-*r/28.0%
metadata-eval28.0%
associate-*l/28.0%
associate-*r/27.9%
associate-*l/27.0%
associate-/l*28.0%
distribute-rgt-out28.0%
associate-*r/28.0%
metadata-eval28.0%
distribute-neg-frac28.0%
metadata-eval28.0%
Simplified28.0%
Taylor expanded in z around inf 28.0%
Final simplification67.3%
(FPCore (x y z) :precision binary64 (if (or (<= z -30.0) (not (<= z 1.9e+71))) (* -0.0027777777777778 (/ z x)) (/ 0.083333333333333 x)))
double code(double x, double y, double z) {
double tmp;
if ((z <= -30.0) || !(z <= 1.9e+71)) {
tmp = -0.0027777777777778 * (z / x);
} else {
tmp = 0.083333333333333 / x;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((z <= (-30.0d0)) .or. (.not. (z <= 1.9d+71))) then
tmp = (-0.0027777777777778d0) * (z / x)
else
tmp = 0.083333333333333d0 / x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -30.0) || !(z <= 1.9e+71)) {
tmp = -0.0027777777777778 * (z / x);
} else {
tmp = 0.083333333333333 / x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -30.0) or not (z <= 1.9e+71): tmp = -0.0027777777777778 * (z / x) else: tmp = 0.083333333333333 / x return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -30.0) || !(z <= 1.9e+71)) tmp = Float64(-0.0027777777777778 * Float64(z / x)); else tmp = Float64(0.083333333333333 / x); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -30.0) || ~((z <= 1.9e+71))) tmp = -0.0027777777777778 * (z / x); else tmp = 0.083333333333333 / x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -30.0], N[Not[LessEqual[z, 1.9e+71]], $MachinePrecision]], N[(-0.0027777777777778 * N[(z / x), $MachinePrecision]), $MachinePrecision], N[(0.083333333333333 / x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -30 \lor \neg \left(z \leq 1.9 \cdot 10^{+71}\right):\\
\;\;\;\;-0.0027777777777778 \cdot \frac{z}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.083333333333333}{x}\\
\end{array}
\end{array}
if z < -30 or 1.9e71 < z Initial program 90.6%
remove-double-neg90.6%
distribute-frac-neg290.6%
sub-neg90.6%
associate-+l+90.6%
fma-define90.6%
sub-neg90.6%
metadata-eval90.6%
+-commutative90.6%
unsub-neg90.6%
distribute-frac-neg290.6%
remove-double-neg90.6%
Simplified90.6%
Taylor expanded in z around inf 81.2%
Taylor expanded in z around 0 18.8%
if -30 < z < 1.9e71Initial program 98.8%
remove-double-neg98.8%
distribute-frac-neg298.8%
sub-neg98.8%
associate-+l+98.8%
fma-define99.0%
sub-neg99.0%
metadata-eval99.0%
+-commutative99.0%
unsub-neg99.0%
distribute-frac-neg299.0%
remove-double-neg99.0%
Simplified99.0%
Taylor expanded in x around 0 49.4%
Taylor expanded in y around inf 49.4%
Taylor expanded in z around 0 32.5%
Final simplification26.1%
(FPCore (x y z) :precision binary64 (/ (+ 0.083333333333333 (* z -0.0027777777777778)) x))
double code(double x, double y, double z) {
return (0.083333333333333 + (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 = (0.083333333333333d0 + (z * (-0.0027777777777778d0))) / x
end function
public static double code(double x, double y, double z) {
return (0.083333333333333 + (z * -0.0027777777777778)) / x;
}
def code(x, y, z): return (0.083333333333333 + (z * -0.0027777777777778)) / x
function code(x, y, z) return Float64(Float64(0.083333333333333 + Float64(z * -0.0027777777777778)) / x) end
function tmp = code(x, y, z) tmp = (0.083333333333333 + (z * -0.0027777777777778)) / x; end
code[x_, y_, z_] := N[(N[(0.083333333333333 + N[(z * -0.0027777777777778), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333 + z \cdot -0.0027777777777778}{x}
\end{array}
Initial program 95.0%
remove-double-neg95.0%
distribute-frac-neg295.0%
sub-neg95.0%
associate-+l+95.0%
fma-define95.1%
sub-neg95.1%
metadata-eval95.1%
+-commutative95.1%
unsub-neg95.1%
distribute-frac-neg295.1%
remove-double-neg95.1%
Simplified95.1%
Taylor expanded in x around 0 64.3%
Taylor expanded in z around 0 26.1%
*-commutative26.1%
Simplified26.1%
(FPCore (x y z) :precision binary64 (/ 0.083333333333333 x))
double code(double x, double y, double z) {
return 0.083333333333333 / x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.083333333333333d0 / x
end function
public static double code(double x, double y, double z) {
return 0.083333333333333 / x;
}
def code(x, y, z): return 0.083333333333333 / x
function code(x, y, z) return Float64(0.083333333333333 / x) end
function tmp = code(x, y, z) tmp = 0.083333333333333 / x; end
code[x_, y_, z_] := N[(0.083333333333333 / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.083333333333333}{x}
\end{array}
Initial program 95.0%
remove-double-neg95.0%
distribute-frac-neg295.0%
sub-neg95.0%
associate-+l+95.0%
fma-define95.1%
sub-neg95.1%
metadata-eval95.1%
+-commutative95.1%
unsub-neg95.1%
distribute-frac-neg295.1%
remove-double-neg95.1%
Simplified95.1%
Taylor expanded in x around 0 64.3%
Taylor expanded in y around inf 62.8%
Taylor expanded in z around 0 18.9%
(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 2024163
(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)))