
(FPCore (x y) :precision binary64 (- (- 1.0 (/ 1.0 (* x 9.0))) (/ y (* 3.0 (sqrt x)))))
double code(double x, double y) {
return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * sqrt(x)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (1.0d0 / (x * 9.0d0))) - (y / (3.0d0 * sqrt(x)))
end function
public static double code(double x, double y) {
return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * Math.sqrt(x)));
}
def code(x, y): return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * math.sqrt(x)))
function code(x, y) return Float64(Float64(1.0 - Float64(1.0 / Float64(x * 9.0))) - Float64(y / Float64(3.0 * sqrt(x)))) end
function tmp = code(x, y) tmp = (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * sqrt(x))); end
code[x_, y_] := N[(N[(1.0 - N[(1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y / N[(3.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{1}{x \cdot 9}\right) - \frac{y}{3 \cdot \sqrt{x}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (- (- 1.0 (/ 1.0 (* x 9.0))) (/ y (* 3.0 (sqrt x)))))
double code(double x, double y) {
return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * sqrt(x)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (1.0d0 / (x * 9.0d0))) - (y / (3.0d0 * sqrt(x)))
end function
public static double code(double x, double y) {
return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * Math.sqrt(x)));
}
def code(x, y): return (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * math.sqrt(x)))
function code(x, y) return Float64(Float64(1.0 - Float64(1.0 / Float64(x * 9.0))) - Float64(y / Float64(3.0 * sqrt(x)))) end
function tmp = code(x, y) tmp = (1.0 - (1.0 / (x * 9.0))) - (y / (3.0 * sqrt(x))); end
code[x_, y_] := N[(N[(1.0 - N[(1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y / N[(3.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{1}{x \cdot 9}\right) - \frac{y}{3 \cdot \sqrt{x}}
\end{array}
(FPCore (x y) :precision binary64 (- (+ 1.0 (/ -1.0 (* x 9.0))) (/ y (sqrt (* x 9.0)))))
double code(double x, double y) {
return (1.0 + (-1.0 / (x * 9.0))) - (y / sqrt((x * 9.0)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 + ((-1.0d0) / (x * 9.0d0))) - (y / sqrt((x * 9.0d0)))
end function
public static double code(double x, double y) {
return (1.0 + (-1.0 / (x * 9.0))) - (y / Math.sqrt((x * 9.0)));
}
def code(x, y): return (1.0 + (-1.0 / (x * 9.0))) - (y / math.sqrt((x * 9.0)))
function code(x, y) return Float64(Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))) - Float64(y / sqrt(Float64(x * 9.0)))) end
function tmp = code(x, y) tmp = (1.0 + (-1.0 / (x * 9.0))) - (y / sqrt((x * 9.0))); end
code[x_, y_] := N[(N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y / N[Sqrt[N[(x * 9.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 + \frac{-1}{x \cdot 9}\right) - \frac{y}{\sqrt{x \cdot 9}}
\end{array}
Initial program 99.7%
*-commutative99.7%
metadata-eval99.7%
sqrt-prod99.7%
pow1/299.7%
Applied egg-rr99.7%
unpow1/299.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y) :precision binary64 (if (or (<= y -5.1e+42) (not (<= y 2.6e+16))) (- 1.0 (/ y (sqrt (* x 9.0)))) (- 1.0 (/ 0.1111111111111111 x))))
double code(double x, double y) {
double tmp;
if ((y <= -5.1e+42) || !(y <= 2.6e+16)) {
tmp = 1.0 - (y / sqrt((x * 9.0)));
} else {
tmp = 1.0 - (0.1111111111111111 / x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-5.1d+42)) .or. (.not. (y <= 2.6d+16))) then
tmp = 1.0d0 - (y / sqrt((x * 9.0d0)))
else
tmp = 1.0d0 - (0.1111111111111111d0 / x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -5.1e+42) || !(y <= 2.6e+16)) {
tmp = 1.0 - (y / Math.sqrt((x * 9.0)));
} else {
tmp = 1.0 - (0.1111111111111111 / x);
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -5.1e+42) or not (y <= 2.6e+16): tmp = 1.0 - (y / math.sqrt((x * 9.0))) else: tmp = 1.0 - (0.1111111111111111 / x) return tmp
function code(x, y) tmp = 0.0 if ((y <= -5.1e+42) || !(y <= 2.6e+16)) tmp = Float64(1.0 - Float64(y / sqrt(Float64(x * 9.0)))); else tmp = Float64(1.0 - Float64(0.1111111111111111 / x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -5.1e+42) || ~((y <= 2.6e+16))) tmp = 1.0 - (y / sqrt((x * 9.0))); else tmp = 1.0 - (0.1111111111111111 / x); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -5.1e+42], N[Not[LessEqual[y, 2.6e+16]], $MachinePrecision]], N[(1.0 - N[(y / N[Sqrt[N[(x * 9.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.1 \cdot 10^{+42} \lor \neg \left(y \leq 2.6 \cdot 10^{+16}\right):\\
\;\;\;\;1 - \frac{y}{\sqrt{x \cdot 9}}\\
\mathbf{else}:\\
\;\;\;\;1 - \frac{0.1111111111111111}{x}\\
\end{array}
\end{array}
if y < -5.0999999999999999e42 or 2.6e16 < y Initial program 99.6%
*-commutative99.6%
metadata-eval99.6%
sqrt-prod99.8%
pow1/299.8%
Applied egg-rr99.8%
unpow1/299.8%
Simplified99.8%
Taylor expanded in x around inf 92.0%
if -5.0999999999999999e42 < y < 2.6e16Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification96.5%
(FPCore (x y) :precision binary64 (if (or (<= y -6.5e+42) (not (<= y 2.1e+16))) (+ 1.0 (* -0.3333333333333333 (/ y (sqrt x)))) (- 1.0 (/ 0.1111111111111111 x))))
double code(double x, double y) {
double tmp;
if ((y <= -6.5e+42) || !(y <= 2.1e+16)) {
tmp = 1.0 + (-0.3333333333333333 * (y / sqrt(x)));
} else {
tmp = 1.0 - (0.1111111111111111 / x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-6.5d+42)) .or. (.not. (y <= 2.1d+16))) then
tmp = 1.0d0 + ((-0.3333333333333333d0) * (y / sqrt(x)))
else
tmp = 1.0d0 - (0.1111111111111111d0 / x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -6.5e+42) || !(y <= 2.1e+16)) {
tmp = 1.0 + (-0.3333333333333333 * (y / Math.sqrt(x)));
} else {
tmp = 1.0 - (0.1111111111111111 / x);
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -6.5e+42) or not (y <= 2.1e+16): tmp = 1.0 + (-0.3333333333333333 * (y / math.sqrt(x))) else: tmp = 1.0 - (0.1111111111111111 / x) return tmp
function code(x, y) tmp = 0.0 if ((y <= -6.5e+42) || !(y <= 2.1e+16)) tmp = Float64(1.0 + Float64(-0.3333333333333333 * Float64(y / sqrt(x)))); else tmp = Float64(1.0 - Float64(0.1111111111111111 / x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -6.5e+42) || ~((y <= 2.1e+16))) tmp = 1.0 + (-0.3333333333333333 * (y / sqrt(x))); else tmp = 1.0 - (0.1111111111111111 / x); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -6.5e+42], N[Not[LessEqual[y, 2.1e+16]], $MachinePrecision]], N[(1.0 + N[(-0.3333333333333333 * N[(y / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.5 \cdot 10^{+42} \lor \neg \left(y \leq 2.1 \cdot 10^{+16}\right):\\
\;\;\;\;1 + -0.3333333333333333 \cdot \frac{y}{\sqrt{x}}\\
\mathbf{else}:\\
\;\;\;\;1 - \frac{0.1111111111111111}{x}\\
\end{array}
\end{array}
if y < -6.50000000000000052e42 or 2.1e16 < y Initial program 99.6%
sub-neg99.6%
*-commutative99.6%
associate-/r*99.6%
metadata-eval99.6%
distribute-frac-neg99.6%
neg-mul-199.6%
times-frac99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 91.8%
if -6.50000000000000052e42 < y < 2.1e16Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification96.4%
(FPCore (x y) :precision binary64 (if (or (<= y -2.8e+71) (not (<= y 4.5e+81))) (/ -0.3333333333333333 (/ (sqrt x) y)) (+ 1.0 (/ (/ -1.0 x) 9.0))))
double code(double x, double y) {
double tmp;
if ((y <= -2.8e+71) || !(y <= 4.5e+81)) {
tmp = -0.3333333333333333 / (sqrt(x) / y);
} else {
tmp = 1.0 + ((-1.0 / x) / 9.0);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-2.8d+71)) .or. (.not. (y <= 4.5d+81))) then
tmp = (-0.3333333333333333d0) / (sqrt(x) / y)
else
tmp = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -2.8e+71) || !(y <= 4.5e+81)) {
tmp = -0.3333333333333333 / (Math.sqrt(x) / y);
} else {
tmp = 1.0 + ((-1.0 / x) / 9.0);
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -2.8e+71) or not (y <= 4.5e+81): tmp = -0.3333333333333333 / (math.sqrt(x) / y) else: tmp = 1.0 + ((-1.0 / x) / 9.0) return tmp
function code(x, y) tmp = 0.0 if ((y <= -2.8e+71) || !(y <= 4.5e+81)) tmp = Float64(-0.3333333333333333 / Float64(sqrt(x) / y)); else tmp = Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -2.8e+71) || ~((y <= 4.5e+81))) tmp = -0.3333333333333333 / (sqrt(x) / y); else tmp = 1.0 + ((-1.0 / x) / 9.0); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -2.8e+71], N[Not[LessEqual[y, 4.5e+81]], $MachinePrecision]], N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.8 \cdot 10^{+71} \lor \neg \left(y \leq 4.5 \cdot 10^{+81}\right):\\
\;\;\;\;\frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{\frac{-1}{x}}{9}\\
\end{array}
\end{array}
if y < -2.80000000000000002e71 or 4.50000000000000017e81 < y Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.5%
fma-neg99.5%
associate-/r*99.5%
metadata-eval99.5%
*-commutative99.5%
associate-/r*99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf 95.8%
associate-*r*95.7%
Simplified95.7%
associate-*l*95.8%
sqrt-div95.8%
metadata-eval95.8%
associate-*l/95.9%
*-un-lft-identity95.9%
clear-num95.9%
un-div-inv96.1%
Applied egg-rr96.1%
if -2.80000000000000002e71 < y < 4.50000000000000017e81Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 96.1%
associate-*r/96.2%
metadata-eval96.2%
Simplified96.2%
metadata-eval96.2%
associate-/r*96.2%
*-commutative96.2%
inv-pow96.2%
Applied egg-rr96.2%
unpow-196.2%
associate-/r*96.2%
Applied egg-rr96.2%
Final simplification96.1%
(FPCore (x y) :precision binary64 (if (or (<= y -4.7e+76) (not (<= y 1.35e+80))) (* y (/ -0.3333333333333333 (sqrt x))) (+ 1.0 (/ (/ -1.0 x) 9.0))))
double code(double x, double y) {
double tmp;
if ((y <= -4.7e+76) || !(y <= 1.35e+80)) {
tmp = y * (-0.3333333333333333 / sqrt(x));
} else {
tmp = 1.0 + ((-1.0 / x) / 9.0);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-4.7d+76)) .or. (.not. (y <= 1.35d+80))) then
tmp = y * ((-0.3333333333333333d0) / sqrt(x))
else
tmp = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -4.7e+76) || !(y <= 1.35e+80)) {
tmp = y * (-0.3333333333333333 / Math.sqrt(x));
} else {
tmp = 1.0 + ((-1.0 / x) / 9.0);
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -4.7e+76) or not (y <= 1.35e+80): tmp = y * (-0.3333333333333333 / math.sqrt(x)) else: tmp = 1.0 + ((-1.0 / x) / 9.0) return tmp
function code(x, y) tmp = 0.0 if ((y <= -4.7e+76) || !(y <= 1.35e+80)) tmp = Float64(y * Float64(-0.3333333333333333 / sqrt(x))); else tmp = Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -4.7e+76) || ~((y <= 1.35e+80))) tmp = y * (-0.3333333333333333 / sqrt(x)); else tmp = 1.0 + ((-1.0 / x) / 9.0); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -4.7e+76], N[Not[LessEqual[y, 1.35e+80]], $MachinePrecision]], N[(y * N[(-0.3333333333333333 / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.7 \cdot 10^{+76} \lor \neg \left(y \leq 1.35 \cdot 10^{+80}\right):\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{\sqrt{x}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{\frac{-1}{x}}{9}\\
\end{array}
\end{array}
if y < -4.7000000000000003e76 or 1.34999999999999991e80 < y Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.5%
fma-neg99.5%
associate-/r*99.5%
metadata-eval99.5%
*-commutative99.5%
associate-/r*99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf 95.8%
associate-*r*95.7%
Simplified95.7%
sqrt-div95.8%
metadata-eval95.8%
un-div-inv95.9%
Applied egg-rr95.9%
if -4.7000000000000003e76 < y < 1.34999999999999991e80Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 96.1%
associate-*r/96.2%
metadata-eval96.2%
Simplified96.2%
metadata-eval96.2%
associate-/r*96.2%
*-commutative96.2%
inv-pow96.2%
Applied egg-rr96.2%
unpow-196.2%
associate-/r*96.2%
Applied egg-rr96.2%
Final simplification96.1%
(FPCore (x y) :precision binary64 (if (<= x 0.11) (+ (/ -0.3333333333333333 (/ (sqrt x) y)) (/ -0.1111111111111111 x)) (- 1.0 (/ y (sqrt (* x 9.0))))))
double code(double x, double y) {
double tmp;
if (x <= 0.11) {
tmp = (-0.3333333333333333 / (sqrt(x) / y)) + (-0.1111111111111111 / x);
} else {
tmp = 1.0 - (y / sqrt((x * 9.0)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 0.11d0) then
tmp = ((-0.3333333333333333d0) / (sqrt(x) / y)) + ((-0.1111111111111111d0) / x)
else
tmp = 1.0d0 - (y / sqrt((x * 9.0d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 0.11) {
tmp = (-0.3333333333333333 / (Math.sqrt(x) / y)) + (-0.1111111111111111 / x);
} else {
tmp = 1.0 - (y / Math.sqrt((x * 9.0)));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 0.11: tmp = (-0.3333333333333333 / (math.sqrt(x) / y)) + (-0.1111111111111111 / x) else: tmp = 1.0 - (y / math.sqrt((x * 9.0))) return tmp
function code(x, y) tmp = 0.0 if (x <= 0.11) tmp = Float64(Float64(-0.3333333333333333 / Float64(sqrt(x) / y)) + Float64(-0.1111111111111111 / x)); else tmp = Float64(1.0 - Float64(y / sqrt(Float64(x * 9.0)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 0.11) tmp = (-0.3333333333333333 / (sqrt(x) / y)) + (-0.1111111111111111 / x); else tmp = 1.0 - (y / sqrt((x * 9.0))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 0.11], N[(N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision] + N[(-0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision], N[(1.0 - N[(y / N[Sqrt[N[(x * 9.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.11:\\
\;\;\;\;\frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}} + \frac{-0.1111111111111111}{x}\\
\mathbf{else}:\\
\;\;\;\;1 - \frac{y}{\sqrt{x \cdot 9}}\\
\end{array}
\end{array}
if x < 0.110000000000000001Initial program 99.6%
*-commutative99.6%
metadata-eval99.6%
sqrt-prod99.6%
pow1/299.6%
Applied egg-rr99.6%
unpow1/299.6%
Simplified99.6%
Taylor expanded in x around 0 98.4%
sub-neg98.4%
distribute-neg-frac298.4%
sqrt-prod98.4%
distribute-rgt-neg-in98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
*-rgt-identity98.4%
times-frac98.3%
metadata-eval98.3%
metadata-eval98.3%
times-frac98.3%
*-commutative98.3%
times-frac98.3%
/-rgt-identity98.3%
associate-/r/98.4%
Simplified98.4%
if 0.110000000000000001 < x Initial program 99.8%
*-commutative99.8%
metadata-eval99.8%
sqrt-prod99.9%
pow1/299.9%
Applied egg-rr99.9%
unpow1/299.9%
Simplified99.9%
Taylor expanded in x around inf 99.5%
Final simplification98.9%
(FPCore (x y) :precision binary64 (+ (- 1.0 (/ 0.1111111111111111 x)) (/ -0.3333333333333333 (/ (sqrt x) y))))
double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 / (sqrt(x) / y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (0.1111111111111111d0 / x)) + ((-0.3333333333333333d0) / (sqrt(x) / y))
end function
public static double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 / (Math.sqrt(x) / y));
}
def code(x, y): return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 / (math.sqrt(x) / y))
function code(x, y) return Float64(Float64(1.0 - Float64(0.1111111111111111 / x)) + Float64(-0.3333333333333333 / Float64(sqrt(x) / y))) end
function tmp = code(x, y) tmp = (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 / (sqrt(x) / y)); end
code[x_, y_] := N[(N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision] + N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{0.1111111111111111}{x}\right) + \frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}}
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
associate-/r*99.7%
metadata-eval99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
times-frac99.6%
metadata-eval99.6%
Simplified99.6%
clear-num99.6%
un-div-inv99.7%
Applied egg-rr99.7%
(FPCore (x y) :precision binary64 (+ (- 1.0 (/ 0.1111111111111111 x)) (* -0.3333333333333333 (/ y (sqrt x)))))
double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 * (y / sqrt(x)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (0.1111111111111111d0 / x)) + ((-0.3333333333333333d0) * (y / sqrt(x)))
end function
public static double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 * (y / Math.sqrt(x)));
}
def code(x, y): return (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 * (y / math.sqrt(x)))
function code(x, y) return Float64(Float64(1.0 - Float64(0.1111111111111111 / x)) + Float64(-0.3333333333333333 * Float64(y / sqrt(x)))) end
function tmp = code(x, y) tmp = (1.0 - (0.1111111111111111 / x)) + (-0.3333333333333333 * (y / sqrt(x))); end
code[x_, y_] := N[(N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision] + N[(-0.3333333333333333 * N[(y / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{0.1111111111111111}{x}\right) + -0.3333333333333333 \cdot \frac{y}{\sqrt{x}}
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
associate-/r*99.7%
metadata-eval99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
times-frac99.6%
metadata-eval99.6%
Simplified99.6%
(FPCore (x y)
:precision binary64
(if (<= y -8.8e+80)
(+
1.0
(/
(+
0.1111111111111111
(/ (+ 0.024691358024691357 (* (/ 1.0 x) 0.0027434842249657062)) x))
x))
(if (<= y 1.4e+154)
(+ 1.0 (/ (/ -1.0 x) 9.0))
(/
(+ 1.0 (/ 0.012345679012345678 (* x x)))
(+ 1.0 (/ -0.1111111111111111 x))))))
double code(double x, double y) {
double tmp;
if (y <= -8.8e+80) {
tmp = 1.0 + ((0.1111111111111111 + ((0.024691358024691357 + ((1.0 / x) * 0.0027434842249657062)) / x)) / x);
} else if (y <= 1.4e+154) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-8.8d+80)) then
tmp = 1.0d0 + ((0.1111111111111111d0 + ((0.024691358024691357d0 + ((1.0d0 / x) * 0.0027434842249657062d0)) / x)) / x)
else if (y <= 1.4d+154) then
tmp = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
else
tmp = (1.0d0 + (0.012345679012345678d0 / (x * x))) / (1.0d0 + ((-0.1111111111111111d0) / x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -8.8e+80) {
tmp = 1.0 + ((0.1111111111111111 + ((0.024691358024691357 + ((1.0 / x) * 0.0027434842249657062)) / x)) / x);
} else if (y <= 1.4e+154) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -8.8e+80: tmp = 1.0 + ((0.1111111111111111 + ((0.024691358024691357 + ((1.0 / x) * 0.0027434842249657062)) / x)) / x) elif y <= 1.4e+154: tmp = 1.0 + ((-1.0 / x) / 9.0) else: tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x)) return tmp
function code(x, y) tmp = 0.0 if (y <= -8.8e+80) tmp = Float64(1.0 + Float64(Float64(0.1111111111111111 + Float64(Float64(0.024691358024691357 + Float64(Float64(1.0 / x) * 0.0027434842249657062)) / x)) / x)); elseif (y <= 1.4e+154) tmp = Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)); else tmp = Float64(Float64(1.0 + Float64(0.012345679012345678 / Float64(x * x))) / Float64(1.0 + Float64(-0.1111111111111111 / x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -8.8e+80) tmp = 1.0 + ((0.1111111111111111 + ((0.024691358024691357 + ((1.0 / x) * 0.0027434842249657062)) / x)) / x); elseif (y <= 1.4e+154) tmp = 1.0 + ((-1.0 / x) / 9.0); else tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -8.8e+80], N[(1.0 + N[(N[(0.1111111111111111 + N[(N[(0.024691358024691357 + N[(N[(1.0 / x), $MachinePrecision] * 0.0027434842249657062), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.4e+154], N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(0.012345679012345678 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(-0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8.8 \cdot 10^{+80}:\\
\;\;\;\;1 + \frac{0.1111111111111111 + \frac{0.024691358024691357 + \frac{1}{x} \cdot 0.0027434842249657062}{x}}{x}\\
\mathbf{elif}\;y \leq 1.4 \cdot 10^{+154}:\\
\;\;\;\;1 + \frac{\frac{-1}{x}}{9}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \frac{0.012345679012345678}{x \cdot x}}{1 + \frac{-0.1111111111111111}{x}}\\
\end{array}
\end{array}
if y < -8.80000000000000011e80Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.5%
fma-neg99.5%
associate-/r*99.4%
metadata-eval99.4%
*-commutative99.4%
associate-/r*99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around 0 1.9%
associate-*r/1.9%
metadata-eval1.9%
Simplified1.9%
Applied egg-rr1.8%
associate-*r/1.8%
*-rgt-identity1.8%
Simplified1.8%
Taylor expanded in x around -inf 33.4%
if -8.80000000000000011e80 < y < 1.4e154Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 90.9%
associate-*r/91.0%
metadata-eval91.0%
Simplified91.0%
metadata-eval91.0%
associate-/r*90.9%
*-commutative90.9%
inv-pow90.9%
Applied egg-rr90.9%
unpow-190.9%
associate-/r*91.0%
Applied egg-rr91.0%
if 1.4e154 < y Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 3.2%
associate-*r/3.2%
metadata-eval3.2%
Simplified3.2%
Applied egg-rr21.4%
associate-*r/21.4%
*-rgt-identity21.4%
Simplified21.4%
unpow221.4%
Applied egg-rr21.4%
Final simplification72.6%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ 1.0 (/ -0.1111111111111111 x))))
(if (<= y -6.5e+88)
(/ (+ 1.0 (* (/ 1.0 x) (/ -0.012345679012345678 x))) t_0)
(if (<= y 1.06e+155)
(+ 1.0 (/ (/ -1.0 x) 9.0))
(/ (+ 1.0 (/ 0.012345679012345678 (* x x))) t_0)))))
double code(double x, double y) {
double t_0 = 1.0 + (-0.1111111111111111 / x);
double tmp;
if (y <= -6.5e+88) {
tmp = (1.0 + ((1.0 / x) * (-0.012345679012345678 / x))) / t_0;
} else if (y <= 1.06e+155) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + ((-0.1111111111111111d0) / x)
if (y <= (-6.5d+88)) then
tmp = (1.0d0 + ((1.0d0 / x) * ((-0.012345679012345678d0) / x))) / t_0
else if (y <= 1.06d+155) then
tmp = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
else
tmp = (1.0d0 + (0.012345679012345678d0 / (x * x))) / t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 1.0 + (-0.1111111111111111 / x);
double tmp;
if (y <= -6.5e+88) {
tmp = (1.0 + ((1.0 / x) * (-0.012345679012345678 / x))) / t_0;
} else if (y <= 1.06e+155) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / t_0;
}
return tmp;
}
def code(x, y): t_0 = 1.0 + (-0.1111111111111111 / x) tmp = 0 if y <= -6.5e+88: tmp = (1.0 + ((1.0 / x) * (-0.012345679012345678 / x))) / t_0 elif y <= 1.06e+155: tmp = 1.0 + ((-1.0 / x) / 9.0) else: tmp = (1.0 + (0.012345679012345678 / (x * x))) / t_0 return tmp
function code(x, y) t_0 = Float64(1.0 + Float64(-0.1111111111111111 / x)) tmp = 0.0 if (y <= -6.5e+88) tmp = Float64(Float64(1.0 + Float64(Float64(1.0 / x) * Float64(-0.012345679012345678 / x))) / t_0); elseif (y <= 1.06e+155) tmp = Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)); else tmp = Float64(Float64(1.0 + Float64(0.012345679012345678 / Float64(x * x))) / t_0); end return tmp end
function tmp_2 = code(x, y) t_0 = 1.0 + (-0.1111111111111111 / x); tmp = 0.0; if (y <= -6.5e+88) tmp = (1.0 + ((1.0 / x) * (-0.012345679012345678 / x))) / t_0; elseif (y <= 1.06e+155) tmp = 1.0 + ((-1.0 / x) / 9.0); else tmp = (1.0 + (0.012345679012345678 / (x * x))) / t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(1.0 + N[(-0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.5e+88], N[(N[(1.0 + N[(N[(1.0 / x), $MachinePrecision] * N[(-0.012345679012345678 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision], If[LessEqual[y, 1.06e+155], N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(0.012345679012345678 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{-0.1111111111111111}{x}\\
\mathbf{if}\;y \leq -6.5 \cdot 10^{+88}:\\
\;\;\;\;\frac{1 + \frac{1}{x} \cdot \frac{-0.012345679012345678}{x}}{t\_0}\\
\mathbf{elif}\;y \leq 1.06 \cdot 10^{+155}:\\
\;\;\;\;1 + \frac{\frac{-1}{x}}{9}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \frac{0.012345679012345678}{x \cdot x}}{t\_0}\\
\end{array}
\end{array}
if y < -6.5000000000000002e88Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.5%
fma-neg99.5%
associate-/r*99.3%
metadata-eval99.3%
*-commutative99.3%
associate-/r*99.3%
distribute-neg-frac99.3%
metadata-eval99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in y around 0 1.7%
associate-*r/1.7%
metadata-eval1.7%
Simplified1.7%
Applied egg-rr1.6%
associate-*r/1.6%
*-rgt-identity1.6%
Simplified1.6%
clear-num1.6%
inv-pow1.6%
div-inv1.6%
unpow21.6%
metadata-eval1.6%
metadata-eval1.6%
swap-sqr1.6%
pow-prod-down1.6%
unpow-prod-down1.6%
inv-pow1.6%
associate-*l*1.6%
metadata-eval1.6%
unpow-11.6%
add-sqr-sqrt1.6%
metadata-eval1.6%
sqrt-unprod1.6%
swap-sqr1.6%
unpow21.6%
metadata-eval1.6%
metadata-eval1.6%
div-inv1.6%
sqrt-div1.6%
clear-num1.6%
metadata-eval1.6%
unpow21.6%
frac-times1.6%
Applied egg-rr28.9%
associate-*r/28.9%
metadata-eval28.9%
Simplified28.9%
if -6.5000000000000002e88 < y < 1.06000000000000005e155Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
fma-neg99.7%
associate-/r*99.7%
metadata-eval99.7%
*-commutative99.7%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 89.9%
associate-*r/90.0%
metadata-eval90.0%
Simplified90.0%
metadata-eval90.0%
associate-/r*90.0%
*-commutative90.0%
inv-pow90.0%
Applied egg-rr90.0%
unpow-190.0%
associate-/r*90.0%
Applied egg-rr90.0%
if 1.06000000000000005e155 < y Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 3.3%
associate-*r/3.3%
metadata-eval3.3%
Simplified3.3%
Applied egg-rr22.0%
associate-*r/22.0%
*-rgt-identity22.0%
Simplified22.0%
unpow222.0%
Applied egg-rr22.0%
Final simplification71.9%
(FPCore (x y)
:precision binary64
(if (<= y 1.06e+155)
(+ 1.0 (/ (/ -1.0 x) 9.0))
(/
(+ 1.0 (/ 0.012345679012345678 (* x x)))
(+ 1.0 (/ -0.1111111111111111 x)))))
double code(double x, double y) {
double tmp;
if (y <= 1.06e+155) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 1.06d+155) then
tmp = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
else
tmp = (1.0d0 + (0.012345679012345678d0 / (x * x))) / (1.0d0 + ((-0.1111111111111111d0) / x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 1.06e+155) {
tmp = 1.0 + ((-1.0 / x) / 9.0);
} else {
tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 1.06e+155: tmp = 1.0 + ((-1.0 / x) / 9.0) else: tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x)) return tmp
function code(x, y) tmp = 0.0 if (y <= 1.06e+155) tmp = Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)); else tmp = Float64(Float64(1.0 + Float64(0.012345679012345678 / Float64(x * x))) / Float64(1.0 + Float64(-0.1111111111111111 / x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 1.06e+155) tmp = 1.0 + ((-1.0 / x) / 9.0); else tmp = (1.0 + (0.012345679012345678 / (x * x))) / (1.0 + (-0.1111111111111111 / x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 1.06e+155], N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(0.012345679012345678 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(-0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.06 \cdot 10^{+155}:\\
\;\;\;\;1 + \frac{\frac{-1}{x}}{9}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \frac{0.012345679012345678}{x \cdot x}}{1 + \frac{-0.1111111111111111}{x}}\\
\end{array}
\end{array}
if y < 1.06000000000000005e155Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 75.1%
associate-*r/75.2%
metadata-eval75.2%
Simplified75.2%
metadata-eval75.2%
associate-/r*75.2%
*-commutative75.2%
inv-pow75.2%
Applied egg-rr75.2%
unpow-175.2%
associate-/r*75.2%
Applied egg-rr75.2%
if 1.06000000000000005e155 < y Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 3.3%
associate-*r/3.3%
metadata-eval3.3%
Simplified3.3%
Applied egg-rr22.0%
associate-*r/22.0%
*-rgt-identity22.0%
Simplified22.0%
unpow222.0%
Applied egg-rr22.0%
Final simplification68.0%
(FPCore (x y) :precision binary64 (if (<= x 0.11) (/ -0.1111111111111111 x) 1.0))
double code(double x, double y) {
double tmp;
if (x <= 0.11) {
tmp = -0.1111111111111111 / x;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 0.11d0) then
tmp = (-0.1111111111111111d0) / x
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 0.11) {
tmp = -0.1111111111111111 / x;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 0.11: tmp = -0.1111111111111111 / x else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (x <= 0.11) tmp = Float64(-0.1111111111111111 / x); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 0.11) tmp = -0.1111111111111111 / x; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 0.11], N[(-0.1111111111111111 / x), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.11:\\
\;\;\;\;\frac{-0.1111111111111111}{x}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < 0.110000000000000001Initial program 99.6%
associate--l-99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
distribute-frac-neg99.6%
sub-neg99.6%
neg-mul-199.6%
*-commutative99.6%
associate-/l*99.5%
fma-neg99.5%
associate-/r*99.5%
metadata-eval99.5%
*-commutative99.5%
associate-/r*99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around 0 98.3%
Taylor expanded in y around 0 60.1%
if 0.110000000000000001 < x Initial program 99.8%
associate--l-99.8%
sub-neg99.8%
+-commutative99.8%
distribute-neg-in99.8%
distribute-frac-neg99.8%
sub-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around 0 70.3%
associate-*r/70.3%
metadata-eval70.3%
Simplified70.3%
Taylor expanded in x around inf 69.9%
(FPCore (x y) :precision binary64 (+ 1.0 (/ (/ -1.0 x) 9.0)))
double code(double x, double y) {
return 1.0 + ((-1.0 / x) / 9.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 + (((-1.0d0) / x) / 9.0d0)
end function
public static double code(double x, double y) {
return 1.0 + ((-1.0 / x) / 9.0);
}
def code(x, y): return 1.0 + ((-1.0 / x) / 9.0)
function code(x, y) return Float64(1.0 + Float64(Float64(-1.0 / x) / 9.0)) end
function tmp = code(x, y) tmp = 1.0 + ((-1.0 / x) / 9.0); end
code[x_, y_] := N[(1.0 + N[(N[(-1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \frac{\frac{-1}{x}}{9}
\end{array}
Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 65.3%
associate-*r/65.4%
metadata-eval65.4%
Simplified65.4%
metadata-eval65.4%
associate-/r*65.4%
*-commutative65.4%
inv-pow65.4%
Applied egg-rr65.4%
unpow-165.4%
associate-/r*65.4%
Applied egg-rr65.4%
Final simplification65.4%
(FPCore (x y) :precision binary64 (- 1.0 (/ 0.1111111111111111 x)))
double code(double x, double y) {
return 1.0 - (0.1111111111111111 / x);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 - (0.1111111111111111d0 / x)
end function
public static double code(double x, double y) {
return 1.0 - (0.1111111111111111 / x);
}
def code(x, y): return 1.0 - (0.1111111111111111 / x)
function code(x, y) return Float64(1.0 - Float64(0.1111111111111111 / x)) end
function tmp = code(x, y) tmp = 1.0 - (0.1111111111111111 / x); end
code[x_, y_] := N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \frac{0.1111111111111111}{x}
\end{array}
Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 65.3%
associate-*r/65.4%
metadata-eval65.4%
Simplified65.4%
(FPCore (x y) :precision binary64 1.0)
double code(double x, double y) {
return 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0
end function
public static double code(double x, double y) {
return 1.0;
}
def code(x, y): return 1.0
function code(x, y) return 1.0 end
function tmp = code(x, y) tmp = 1.0; end
code[x_, y_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 99.7%
associate--l-99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
distribute-frac-neg99.7%
sub-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.6%
fma-neg99.6%
associate-/r*99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around 0 65.3%
associate-*r/65.4%
metadata-eval65.4%
Simplified65.4%
Taylor expanded in x around inf 33.0%
(FPCore (x y) :precision binary64 (- (- 1.0 (/ (/ 1.0 x) 9.0)) (/ y (* 3.0 (sqrt x)))))
double code(double x, double y) {
return (1.0 - ((1.0 / x) / 9.0)) - (y / (3.0 * sqrt(x)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - ((1.0d0 / x) / 9.0d0)) - (y / (3.0d0 * sqrt(x)))
end function
public static double code(double x, double y) {
return (1.0 - ((1.0 / x) / 9.0)) - (y / (3.0 * Math.sqrt(x)));
}
def code(x, y): return (1.0 - ((1.0 / x) / 9.0)) - (y / (3.0 * math.sqrt(x)))
function code(x, y) return Float64(Float64(1.0 - Float64(Float64(1.0 / x) / 9.0)) - Float64(y / Float64(3.0 * sqrt(x)))) end
function tmp = code(x, y) tmp = (1.0 - ((1.0 / x) / 9.0)) - (y / (3.0 * sqrt(x))); end
code[x_, y_] := N[(N[(1.0 - N[(N[(1.0 / x), $MachinePrecision] / 9.0), $MachinePrecision]), $MachinePrecision] - N[(y / N[(3.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{\frac{1}{x}}{9}\right) - \frac{y}{3 \cdot \sqrt{x}}
\end{array}
herbie shell --seed 2024137
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, D"
:precision binary64
:alt
(! :herbie-platform default (- (- 1 (/ (/ 1 x) 9)) (/ y (* 3 (sqrt x)))))
(- (- 1.0 (/ 1.0 (* x 9.0))) (/ y (* 3.0 (sqrt x)))))