
(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 14 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 -4.6e+81) (not (<= y 9.5e+91))) (* -0.3333333333333333 (* y (sqrt (/ 1.0 x)))) (+ 1.0 (/ -1.0 (* x 9.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -4.6e+81) || !(y <= 9.5e+91)) {
tmp = -0.3333333333333333 * (y * sqrt((1.0 / 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.6d+81)) .or. (.not. (y <= 9.5d+91))) then
tmp = (-0.3333333333333333d0) * (y * sqrt((1.0d0 / 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.6e+81) || !(y <= 9.5e+91)) {
tmp = -0.3333333333333333 * (y * Math.sqrt((1.0 / x)));
} else {
tmp = 1.0 + (-1.0 / (x * 9.0));
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -4.6e+81) or not (y <= 9.5e+91): tmp = -0.3333333333333333 * (y * math.sqrt((1.0 / x))) else: tmp = 1.0 + (-1.0 / (x * 9.0)) return tmp
function code(x, y) tmp = 0.0 if ((y <= -4.6e+81) || !(y <= 9.5e+91)) tmp = Float64(-0.3333333333333333 * Float64(y * sqrt(Float64(1.0 / x)))); else tmp = Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -4.6e+81) || ~((y <= 9.5e+91))) tmp = -0.3333333333333333 * (y * sqrt((1.0 / x))); else tmp = 1.0 + (-1.0 / (x * 9.0)); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -4.6e+81], N[Not[LessEqual[y, 9.5e+91]], $MachinePrecision]], N[(-0.3333333333333333 * N[(y * N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.6 \cdot 10^{+81} \lor \neg \left(y \leq 9.5 \cdot 10^{+91}\right):\\
\;\;\;\;-0.3333333333333333 \cdot \left(y \cdot \sqrt{\frac{1}{x}}\right)\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{x \cdot 9}\\
\end{array}
\end{array}
if y < -4.5999999999999998e81 or 9.5000000000000001e91 < y Initial program 99.6%
*-commutative99.6%
metadata-eval99.6%
sqrt-prod99.7%
pow1/299.7%
Applied egg-rr99.7%
unpow1/299.7%
Simplified99.7%
Taylor expanded in y around inf 94.3%
if -4.5999999999999998e81 < y < 9.5000000000000001e91Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.7%
fma-neg99.7%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 93.7%
cancel-sign-sub-inv93.7%
metadata-eval93.7%
associate-*r/93.7%
metadata-eval93.7%
+-commutative93.7%
Simplified93.7%
add-sqr-sqrt0.0%
sqrt-unprod43.0%
frac-times43.0%
metadata-eval43.0%
metadata-eval43.0%
frac-times43.0%
sqrt-unprod43.0%
add-sqr-sqrt43.0%
clear-num43.0%
inv-pow43.0%
div-inv43.0%
metadata-eval43.0%
Applied egg-rr43.0%
*-commutative43.0%
metadata-eval43.0%
add-sqr-sqrt43.0%
swap-sqr43.0%
pow-prod-down43.0%
inv-pow43.0%
frac-2neg43.0%
metadata-eval43.0%
inv-pow43.0%
clear-num43.0%
frac-times43.0%
metadata-eval43.0%
/-rgt-identity43.0%
distribute-lft-neg-in43.0%
metadata-eval43.0%
metadata-eval43.0%
sqrt-prod43.0%
*-commutative43.0%
add-sqr-sqrt0.0%
sqrt-unprod93.6%
swap-sqr93.6%
Applied egg-rr93.8%
Final simplification94.0%
(FPCore (x y) :precision binary64 (if (or (<= y -2.4e+81) (not (<= y 9.2e+91))) (* (sqrt (/ 1.0 x)) (* y -0.3333333333333333)) (+ 1.0 (/ -1.0 (* x 9.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -2.4e+81) || !(y <= 9.2e+91)) {
tmp = sqrt((1.0 / x)) * (y * -0.3333333333333333);
} 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.4d+81)) .or. (.not. (y <= 9.2d+91))) then
tmp = sqrt((1.0d0 / x)) * (y * (-0.3333333333333333d0))
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.4e+81) || !(y <= 9.2e+91)) {
tmp = Math.sqrt((1.0 / x)) * (y * -0.3333333333333333);
} else {
tmp = 1.0 + (-1.0 / (x * 9.0));
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -2.4e+81) or not (y <= 9.2e+91): tmp = math.sqrt((1.0 / x)) * (y * -0.3333333333333333) else: tmp = 1.0 + (-1.0 / (x * 9.0)) return tmp
function code(x, y) tmp = 0.0 if ((y <= -2.4e+81) || !(y <= 9.2e+91)) tmp = Float64(sqrt(Float64(1.0 / x)) * Float64(y * -0.3333333333333333)); else tmp = Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -2.4e+81) || ~((y <= 9.2e+91))) tmp = sqrt((1.0 / x)) * (y * -0.3333333333333333); else tmp = 1.0 + (-1.0 / (x * 9.0)); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -2.4e+81], N[Not[LessEqual[y, 9.2e+91]], $MachinePrecision]], N[(N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision] * N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.4 \cdot 10^{+81} \lor \neg \left(y \leq 9.2 \cdot 10^{+91}\right):\\
\;\;\;\;\sqrt{\frac{1}{x}} \cdot \left(y \cdot -0.3333333333333333\right)\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{x \cdot 9}\\
\end{array}
\end{array}
if y < -2.3999999999999999e81 or 9.19999999999999965e91 < y Initial program 99.6%
*-commutative99.6%
metadata-eval99.6%
sqrt-prod99.7%
pow1/299.7%
Applied egg-rr99.7%
unpow1/299.7%
Simplified99.7%
Taylor expanded in y around inf 94.3%
*-commutative94.3%
associate-*l*94.5%
Simplified94.5%
if -2.3999999999999999e81 < y < 9.19999999999999965e91Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.7%
fma-neg99.7%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 93.7%
cancel-sign-sub-inv93.7%
metadata-eval93.7%
associate-*r/93.7%
metadata-eval93.7%
+-commutative93.7%
Simplified93.7%
add-sqr-sqrt0.0%
sqrt-unprod43.0%
frac-times43.0%
metadata-eval43.0%
metadata-eval43.0%
frac-times43.0%
sqrt-unprod43.0%
add-sqr-sqrt43.0%
clear-num43.0%
inv-pow43.0%
div-inv43.0%
metadata-eval43.0%
Applied egg-rr43.0%
*-commutative43.0%
metadata-eval43.0%
add-sqr-sqrt43.0%
swap-sqr43.0%
pow-prod-down43.0%
inv-pow43.0%
frac-2neg43.0%
metadata-eval43.0%
inv-pow43.0%
clear-num43.0%
frac-times43.0%
metadata-eval43.0%
/-rgt-identity43.0%
distribute-lft-neg-in43.0%
metadata-eval43.0%
metadata-eval43.0%
sqrt-prod43.0%
*-commutative43.0%
add-sqr-sqrt0.0%
sqrt-unprod93.6%
swap-sqr93.6%
Applied egg-rr93.8%
Final simplification94.1%
(FPCore (x y) :precision binary64 (if (or (<= y -4.8e+57) (not (<= y 2.7e+40))) (+ 1.0 (* y (/ -0.3333333333333333 (sqrt x)))) (+ 1.0 (/ -1.0 (* x 9.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -4.8e+57) || !(y <= 2.7e+40)) {
tmp = 1.0 + (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.8d+57)) .or. (.not. (y <= 2.7d+40))) then
tmp = 1.0d0 + (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.8e+57) || !(y <= 2.7e+40)) {
tmp = 1.0 + (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.8e+57) or not (y <= 2.7e+40): tmp = 1.0 + (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.8e+57) || !(y <= 2.7e+40)) tmp = Float64(1.0 + Float64(y * Float64(-0.3333333333333333 / sqrt(x)))); else tmp = Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -4.8e+57) || ~((y <= 2.7e+40))) tmp = 1.0 + (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.8e+57], N[Not[LessEqual[y, 2.7e+40]], $MachinePrecision]], N[(1.0 + N[(y * N[(-0.3333333333333333 / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.8 \cdot 10^{+57} \lor \neg \left(y \leq 2.7 \cdot 10^{+40}\right):\\
\;\;\;\;1 + y \cdot \frac{-0.3333333333333333}{\sqrt{x}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{x \cdot 9}\\
\end{array}
\end{array}
if y < -4.80000000000000009e57 or 2.70000000000000009e40 < y Initial program 99.6%
associate--l-99.6%
+-commutative99.6%
associate--r+99.6%
sub-neg99.6%
distribute-frac-neg99.6%
associate-+r-99.6%
neg-mul-199.6%
associate-*l/99.5%
fma-neg99.5%
associate-/r*99.4%
metadata-eval99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
*-commutative99.4%
associate-/r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 96.1%
*-commutative96.1%
associate-*l*96.2%
*-commutative96.2%
Simplified96.2%
pow196.2%
*-commutative96.2%
*-commutative96.2%
inv-pow96.2%
sqrt-pow196.2%
metadata-eval96.2%
Applied egg-rr96.2%
unpow196.2%
associate-*l*96.1%
Simplified96.1%
add-sqr-sqrt0.0%
sqrt-unprod5.3%
swap-sqr5.3%
metadata-eval5.3%
metadata-eval5.3%
pow-prod-up5.3%
metadata-eval5.3%
unpow-prod-down5.3%
metadata-eval5.3%
add-sqr-sqrt5.3%
swap-sqr5.3%
pow-prod-down5.3%
inv-pow5.3%
inv-pow5.3%
sqrt-unprod0.0%
add-sqr-sqrt96.2%
expm1-log1p-u45.5%
expm1-udef7.2%
associate-/r*7.2%
metadata-eval7.2%
Applied egg-rr7.2%
expm1-def45.4%
expm1-log1p96.1%
Simplified96.1%
if -4.80000000000000009e57 < y < 2.70000000000000009e40Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
associate-*r/97.8%
metadata-eval97.8%
+-commutative97.8%
Simplified97.8%
add-sqr-sqrt0.0%
sqrt-unprod43.7%
frac-times43.7%
metadata-eval43.7%
metadata-eval43.7%
frac-times43.7%
sqrt-unprod43.8%
add-sqr-sqrt43.8%
clear-num43.8%
inv-pow43.8%
div-inv43.8%
metadata-eval43.8%
Applied egg-rr43.8%
*-commutative43.8%
metadata-eval43.8%
add-sqr-sqrt43.8%
swap-sqr43.8%
pow-prod-down43.8%
inv-pow43.8%
frac-2neg43.8%
metadata-eval43.8%
inv-pow43.8%
clear-num43.8%
frac-times43.8%
metadata-eval43.8%
/-rgt-identity43.8%
distribute-lft-neg-in43.8%
metadata-eval43.8%
metadata-eval43.8%
sqrt-prod43.8%
*-commutative43.8%
add-sqr-sqrt0.0%
sqrt-unprod97.7%
swap-sqr97.6%
Applied egg-rr97.9%
Final simplification97.1%
(FPCore (x y) :precision binary64 (if (or (<= y -4.8e+57) (not (<= y 3.1e+45))) (+ 1.0 (/ y (* (sqrt x) -3.0))) (+ 1.0 (/ -1.0 (* x 9.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -4.8e+57) || !(y <= 3.1e+45)) {
tmp = 1.0 + (y / (sqrt(x) * -3.0));
} 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.8d+57)) .or. (.not. (y <= 3.1d+45))) then
tmp = 1.0d0 + (y / (sqrt(x) * (-3.0d0)))
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.8e+57) || !(y <= 3.1e+45)) {
tmp = 1.0 + (y / (Math.sqrt(x) * -3.0));
} else {
tmp = 1.0 + (-1.0 / (x * 9.0));
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -4.8e+57) or not (y <= 3.1e+45): tmp = 1.0 + (y / (math.sqrt(x) * -3.0)) else: tmp = 1.0 + (-1.0 / (x * 9.0)) return tmp
function code(x, y) tmp = 0.0 if ((y <= -4.8e+57) || !(y <= 3.1e+45)) tmp = Float64(1.0 + Float64(y / Float64(sqrt(x) * -3.0))); else tmp = Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -4.8e+57) || ~((y <= 3.1e+45))) tmp = 1.0 + (y / (sqrt(x) * -3.0)); else tmp = 1.0 + (-1.0 / (x * 9.0)); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -4.8e+57], N[Not[LessEqual[y, 3.1e+45]], $MachinePrecision]], N[(1.0 + N[(y / N[(N[Sqrt[x], $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.8 \cdot 10^{+57} \lor \neg \left(y \leq 3.1 \cdot 10^{+45}\right):\\
\;\;\;\;1 + \frac{y}{\sqrt{x} \cdot -3}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{x \cdot 9}\\
\end{array}
\end{array}
if y < -4.80000000000000009e57 or 3.09999999999999988e45 < y Initial program 99.6%
associate--l-99.6%
+-commutative99.6%
associate--r+99.6%
sub-neg99.6%
distribute-frac-neg99.6%
associate-+r-99.6%
neg-mul-199.6%
associate-*l/99.5%
fma-neg99.5%
associate-/r*99.4%
metadata-eval99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
*-commutative99.4%
associate-/r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 96.1%
*-commutative96.1%
associate-*l*96.2%
*-commutative96.2%
Simplified96.2%
pow196.2%
*-commutative96.2%
*-commutative96.2%
inv-pow96.2%
sqrt-pow196.2%
metadata-eval96.2%
Applied egg-rr96.2%
unpow196.2%
associate-*l*96.1%
Simplified96.1%
add-sqr-sqrt44.7%
sqrt-unprod30.5%
swap-sqr25.9%
swap-sqr25.9%
metadata-eval25.9%
metadata-eval25.9%
pow-prod-up26.0%
metadata-eval26.0%
unpow-prod-down25.9%
*-commutative25.9%
add-sqr-sqrt25.9%
unpow-prod-down25.9%
inv-pow25.9%
inv-pow25.9%
swap-sqr30.5%
div-inv30.6%
div-inv30.6%
Applied egg-rr96.3%
*-commutative96.3%
Simplified96.3%
if -4.80000000000000009e57 < y < 3.09999999999999988e45Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
associate-*r/97.8%
metadata-eval97.8%
+-commutative97.8%
Simplified97.8%
add-sqr-sqrt0.0%
sqrt-unprod43.7%
frac-times43.7%
metadata-eval43.7%
metadata-eval43.7%
frac-times43.7%
sqrt-unprod43.8%
add-sqr-sqrt43.8%
clear-num43.8%
inv-pow43.8%
div-inv43.8%
metadata-eval43.8%
Applied egg-rr43.8%
*-commutative43.8%
metadata-eval43.8%
add-sqr-sqrt43.8%
swap-sqr43.8%
pow-prod-down43.8%
inv-pow43.8%
frac-2neg43.8%
metadata-eval43.8%
inv-pow43.8%
clear-num43.8%
frac-times43.8%
metadata-eval43.8%
/-rgt-identity43.8%
distribute-lft-neg-in43.8%
metadata-eval43.8%
metadata-eval43.8%
sqrt-prod43.8%
*-commutative43.8%
add-sqr-sqrt0.0%
sqrt-unprod97.7%
swap-sqr97.6%
Applied egg-rr97.9%
Final simplification97.2%
(FPCore (x y)
:precision binary64
(if (<= y -4.6e+57)
(+ 1.0 (/ (/ y -3.0) (sqrt x)))
(if (<= y 1.1e+45)
(+ 1.0 (/ -1.0 (* x 9.0)))
(+ 1.0 (/ y (* (sqrt x) -3.0))))))
double code(double x, double y) {
double tmp;
if (y <= -4.6e+57) {
tmp = 1.0 + ((y / -3.0) / sqrt(x));
} else if (y <= 1.1e+45) {
tmp = 1.0 + (-1.0 / (x * 9.0));
} else {
tmp = 1.0 + (y / (sqrt(x) * -3.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.6d+57)) then
tmp = 1.0d0 + ((y / (-3.0d0)) / sqrt(x))
else if (y <= 1.1d+45) then
tmp = 1.0d0 + ((-1.0d0) / (x * 9.0d0))
else
tmp = 1.0d0 + (y / (sqrt(x) * (-3.0d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -4.6e+57) {
tmp = 1.0 + ((y / -3.0) / Math.sqrt(x));
} else if (y <= 1.1e+45) {
tmp = 1.0 + (-1.0 / (x * 9.0));
} else {
tmp = 1.0 + (y / (Math.sqrt(x) * -3.0));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -4.6e+57: tmp = 1.0 + ((y / -3.0) / math.sqrt(x)) elif y <= 1.1e+45: tmp = 1.0 + (-1.0 / (x * 9.0)) else: tmp = 1.0 + (y / (math.sqrt(x) * -3.0)) return tmp
function code(x, y) tmp = 0.0 if (y <= -4.6e+57) tmp = Float64(1.0 + Float64(Float64(y / -3.0) / sqrt(x))); elseif (y <= 1.1e+45) tmp = Float64(1.0 + Float64(-1.0 / Float64(x * 9.0))); else tmp = Float64(1.0 + Float64(y / Float64(sqrt(x) * -3.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -4.6e+57) tmp = 1.0 + ((y / -3.0) / sqrt(x)); elseif (y <= 1.1e+45) tmp = 1.0 + (-1.0 / (x * 9.0)); else tmp = 1.0 + (y / (sqrt(x) * -3.0)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -4.6e+57], N[(1.0 + N[(N[(y / -3.0), $MachinePrecision] / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.1e+45], N[(1.0 + N[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(y / N[(N[Sqrt[x], $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.6 \cdot 10^{+57}:\\
\;\;\;\;1 + \frac{\frac{y}{-3}}{\sqrt{x}}\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{+45}:\\
\;\;\;\;1 + \frac{-1}{x \cdot 9}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{y}{\sqrt{x} \cdot -3}\\
\end{array}
\end{array}
if y < -4.5999999999999998e57Initial program 99.5%
associate--l-99.5%
+-commutative99.5%
associate--r+99.5%
sub-neg99.5%
distribute-frac-neg99.5%
associate-+r-99.5%
neg-mul-199.5%
associate-*l/99.5%
fma-neg99.5%
associate-/r*99.4%
metadata-eval99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
*-commutative99.4%
associate-/r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 96.0%
*-commutative96.0%
associate-*l*96.1%
*-commutative96.1%
Simplified96.1%
pow196.1%
*-commutative96.1%
*-commutative96.1%
inv-pow96.1%
sqrt-pow196.0%
metadata-eval96.0%
Applied egg-rr96.0%
unpow196.0%
associate-*l*96.1%
Simplified96.1%
add-sqr-sqrt95.8%
sqrt-unprod61.7%
swap-sqr52.0%
swap-sqr52.0%
metadata-eval52.0%
metadata-eval52.0%
pow-prod-up52.1%
metadata-eval52.1%
unpow-prod-down52.0%
*-commutative52.0%
add-sqr-sqrt52.0%
unpow-prod-down52.0%
inv-pow52.0%
inv-pow52.0%
swap-sqr61.7%
div-inv61.8%
div-inv61.8%
Applied egg-rr96.3%
if -4.5999999999999998e57 < y < 1.1e45Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
associate-*r/97.8%
metadata-eval97.8%
+-commutative97.8%
Simplified97.8%
add-sqr-sqrt0.0%
sqrt-unprod43.7%
frac-times43.7%
metadata-eval43.7%
metadata-eval43.7%
frac-times43.7%
sqrt-unprod43.8%
add-sqr-sqrt43.8%
clear-num43.8%
inv-pow43.8%
div-inv43.8%
metadata-eval43.8%
Applied egg-rr43.8%
*-commutative43.8%
metadata-eval43.8%
add-sqr-sqrt43.8%
swap-sqr43.8%
pow-prod-down43.8%
inv-pow43.8%
frac-2neg43.8%
metadata-eval43.8%
inv-pow43.8%
clear-num43.8%
frac-times43.8%
metadata-eval43.8%
/-rgt-identity43.8%
distribute-lft-neg-in43.8%
metadata-eval43.8%
metadata-eval43.8%
sqrt-prod43.8%
*-commutative43.8%
add-sqr-sqrt0.0%
sqrt-unprod97.7%
swap-sqr97.6%
Applied egg-rr97.9%
if 1.1e45 < y Initial program 99.6%
associate--l-99.6%
+-commutative99.6%
associate--r+99.6%
sub-neg99.6%
distribute-frac-neg99.6%
associate-+r-99.6%
neg-mul-199.6%
associate-*l/99.4%
fma-neg99.5%
associate-/r*99.3%
metadata-eval99.3%
distribute-neg-frac99.3%
metadata-eval99.3%
*-commutative99.3%
associate-/r*99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in y around inf 96.2%
*-commutative96.2%
associate-*l*96.3%
*-commutative96.3%
Simplified96.3%
pow196.3%
*-commutative96.3%
*-commutative96.3%
inv-pow96.3%
sqrt-pow196.3%
metadata-eval96.3%
Applied egg-rr96.3%
unpow196.3%
associate-*l*96.2%
Simplified96.2%
add-sqr-sqrt0.0%
sqrt-unprod3.2%
swap-sqr3.1%
swap-sqr3.1%
metadata-eval3.1%
metadata-eval3.1%
pow-prod-up3.1%
metadata-eval3.1%
unpow-prod-down3.1%
*-commutative3.1%
add-sqr-sqrt3.1%
unpow-prod-down3.1%
inv-pow3.1%
inv-pow3.1%
swap-sqr3.2%
div-inv3.2%
div-inv3.2%
Applied egg-rr96.5%
*-commutative96.5%
Simplified96.5%
Final simplification97.3%
(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%
distribute-frac-neg99.7%
*-commutative99.7%
associate-/r*99.6%
metadata-eval99.6%
neg-mul-199.6%
times-frac99.6%
metadata-eval99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (x y) :precision binary64 (+ (- 1.0 (/ 0.1111111111111111 x)) (/ (* y -0.3333333333333333) (sqrt x))))
double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + ((y * -0.3333333333333333) / sqrt(x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (0.1111111111111111d0 / x)) + ((y * (-0.3333333333333333d0)) / sqrt(x))
end function
public static double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) + ((y * -0.3333333333333333) / Math.sqrt(x));
}
def code(x, y): return (1.0 - (0.1111111111111111 / x)) + ((y * -0.3333333333333333) / math.sqrt(x))
function code(x, y) return Float64(Float64(1.0 - Float64(0.1111111111111111 / x)) + Float64(Float64(y * -0.3333333333333333) / sqrt(x))) end
function tmp = code(x, y) tmp = (1.0 - (0.1111111111111111 / x)) + ((y * -0.3333333333333333) / sqrt(x)); end
code[x_, y_] := N[(N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y * -0.3333333333333333), $MachinePrecision] / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{0.1111111111111111}{x}\right) + \frac{y \cdot -0.3333333333333333}{\sqrt{x}}
\end{array}
Initial program 99.7%
sub-neg99.7%
distribute-frac-neg99.7%
*-commutative99.7%
associate-/r*99.6%
metadata-eval99.6%
neg-mul-199.6%
times-frac99.6%
metadata-eval99.6%
Simplified99.6%
associate-*r/99.6%
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (x y) :precision binary64 (- (- 1.0 (/ 0.1111111111111111 x)) (/ y (sqrt (* x 9.0)))))
double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) - (y / sqrt((x * 9.0)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - (0.1111111111111111d0 / x)) - (y / sqrt((x * 9.0d0)))
end function
public static double code(double x, double y) {
return (1.0 - (0.1111111111111111 / x)) - (y / Math.sqrt((x * 9.0)));
}
def code(x, y): return (1.0 - (0.1111111111111111 / x)) - (y / math.sqrt((x * 9.0)))
function code(x, y) return Float64(Float64(1.0 - Float64(0.1111111111111111 / x)) - Float64(y / sqrt(Float64(x * 9.0)))) end
function tmp = code(x, y) tmp = (1.0 - (0.1111111111111111 / x)) - (y / sqrt((x * 9.0))); end
code[x_, y_] := N[(N[(1.0 - N[(0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision] - N[(y / N[Sqrt[N[(x * 9.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - \frac{0.1111111111111111}{x}\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%
Taylor expanded in x around 0 99.7%
Final simplification99.7%
(FPCore (x y) :precision binary64 (if (<= x 0.11) (* (/ 1.0 x) -0.1111111111111111) 1.0))
double code(double x, double y) {
double tmp;
if (x <= 0.11) {
tmp = (1.0 / x) * -0.1111111111111111;
} 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 = (1.0d0 / x) * (-0.1111111111111111d0)
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 = (1.0 / x) * -0.1111111111111111;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 0.11: tmp = (1.0 / x) * -0.1111111111111111 else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (x <= 0.11) tmp = Float64(Float64(1.0 / x) * -0.1111111111111111); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 0.11) tmp = (1.0 / x) * -0.1111111111111111; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 0.11], N[(N[(1.0 / x), $MachinePrecision] * -0.1111111111111111), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.11:\\
\;\;\;\;\frac{1}{x} \cdot -0.1111111111111111\\
\mathbf{else}:\\
\;\;\;\;1\\
\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 57.6%
clear-num57.6%
associate-/r/57.6%
Applied egg-rr57.6%
if 0.110000000000000001 < x Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 62.4%
Final simplification59.7%
(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(-1.0 / Float64(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[(-1.0 / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \frac{-1}{x \cdot 9}
\end{array}
Initial program 99.7%
associate--l-99.7%
+-commutative99.7%
associate--r+99.7%
sub-neg99.7%
distribute-frac-neg99.7%
associate-+r-99.7%
neg-mul-199.7%
associate-*l/99.6%
fma-neg99.7%
associate-/r*99.6%
metadata-eval99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around 0 61.7%
cancel-sign-sub-inv61.7%
metadata-eval61.7%
associate-*r/61.7%
metadata-eval61.7%
+-commutative61.7%
Simplified61.7%
add-sqr-sqrt0.0%
sqrt-unprod31.6%
frac-times31.6%
metadata-eval31.6%
metadata-eval31.6%
frac-times31.6%
sqrt-unprod29.0%
add-sqr-sqrt29.0%
clear-num29.0%
inv-pow29.0%
div-inv29.0%
metadata-eval29.0%
Applied egg-rr29.0%
*-commutative29.0%
metadata-eval29.0%
add-sqr-sqrt29.0%
swap-sqr29.0%
pow-prod-down29.0%
inv-pow29.0%
frac-2neg29.0%
metadata-eval29.0%
inv-pow29.0%
clear-num29.0%
frac-times29.0%
metadata-eval29.0%
/-rgt-identity29.0%
distribute-lft-neg-in29.0%
metadata-eval29.0%
metadata-eval29.0%
sqrt-prod29.0%
*-commutative29.0%
add-sqr-sqrt0.0%
sqrt-unprod61.6%
swap-sqr61.6%
Applied egg-rr61.8%
Final simplification61.8%
(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%
*-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 57.6%
if 0.110000000000000001 < x Initial program 99.8%
associate--l-99.8%
+-commutative99.8%
associate--r+99.8%
sub-neg99.8%
distribute-frac-neg99.8%
associate-+r-99.8%
neg-mul-199.8%
associate-*l/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 62.4%
Final simplification59.7%
(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%
+-commutative99.7%
associate--r+99.7%
sub-neg99.7%
distribute-frac-neg99.7%
associate-+r-99.7%
neg-mul-199.7%
associate-*l/99.6%
fma-neg99.7%
associate-/r*99.6%
metadata-eval99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around 0 61.7%
cancel-sign-sub-inv61.7%
metadata-eval61.7%
associate-*r/61.7%
metadata-eval61.7%
+-commutative61.7%
Simplified61.7%
Final simplification61.7%
(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%
+-commutative99.7%
associate--r+99.7%
sub-neg99.7%
distribute-frac-neg99.7%
associate-+r-99.7%
neg-mul-199.7%
associate-*l/99.6%
fma-neg99.7%
associate-/r*99.6%
metadata-eval99.6%
distribute-neg-frac99.6%
metadata-eval99.6%
*-commutative99.6%
associate-/r*99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 28.9%
Final simplification28.9%
(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 2023320
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, D"
:precision binary64
:herbie-target
(- (- 1.0 (/ (/ 1.0 x) 9.0)) (/ y (* 3.0 (sqrt x))))
(- (- 1.0 (/ 1.0 (* x 9.0))) (/ y (* 3.0 (sqrt x)))))