
(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 20 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 (<= y -9.5e+48)
(+ 1.0 (* y (* -0.3333333333333333 (pow x -0.5))))
(if (<= y 1.5e+46)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(+ 1.0 (* -0.3333333333333333 (* y (sqrt (/ 1.0 x))))))))
double code(double x, double y) {
double tmp;
if (y <= -9.5e+48) {
tmp = 1.0 + (y * (-0.3333333333333333 * pow(x, -0.5)));
} else if (y <= 1.5e+46) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = 1.0 + (-0.3333333333333333 * (y * sqrt((1.0 / x))));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-9.5d+48)) then
tmp = 1.0d0 + (y * ((-0.3333333333333333d0) * (x ** (-0.5d0))))
else if (y <= 1.5d+46) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = 1.0d0 + ((-0.3333333333333333d0) * (y * sqrt((1.0d0 / x))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -9.5e+48) {
tmp = 1.0 + (y * (-0.3333333333333333 * Math.pow(x, -0.5)));
} else if (y <= 1.5e+46) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = 1.0 + (-0.3333333333333333 * (y * Math.sqrt((1.0 / x))));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -9.5e+48: tmp = 1.0 + (y * (-0.3333333333333333 * math.pow(x, -0.5))) elif y <= 1.5e+46: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = 1.0 + (-0.3333333333333333 * (y * math.sqrt((1.0 / x)))) return tmp
function code(x, y) tmp = 0.0 if (y <= -9.5e+48) tmp = Float64(1.0 + Float64(y * Float64(-0.3333333333333333 * (x ^ -0.5)))); elseif (y <= 1.5e+46) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(1.0 + Float64(-0.3333333333333333 * Float64(y * sqrt(Float64(1.0 / x))))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -9.5e+48) tmp = 1.0 + (y * (-0.3333333333333333 * (x ^ -0.5))); elseif (y <= 1.5e+46) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = 1.0 + (-0.3333333333333333 * (y * sqrt((1.0 / x)))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -9.5e+48], N[(1.0 + N[(y * N[(-0.3333333333333333 * N[Power[x, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.5e+46], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.3333333333333333 * N[(y * N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.5 \cdot 10^{+48}:\\
\;\;\;\;1 + y \cdot \left(-0.3333333333333333 \cdot {x}^{-0.5}\right)\\
\mathbf{elif}\;y \leq 1.5 \cdot 10^{+46}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;1 + -0.3333333333333333 \cdot \left(y \cdot \sqrt{\frac{1}{x}}\right)\\
\end{array}
\end{array}
if y < -9.4999999999999997e48Initial program 99.6%
sub-neg99.6%
distribute-frac-neg99.6%
+-commutative99.6%
associate-+r-99.6%
+-commutative99.6%
associate-+r-99.6%
neg-mul-199.6%
*-commutative99.6%
associate-*r/99.5%
fma-neg99.5%
associate-/r*99.5%
metadata-eval99.5%
*-commutative99.5%
associate-/r*99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 96.1%
associate-*r*96.1%
Simplified96.1%
expm1-log1p-u50.5%
expm1-udef14.0%
inv-pow14.0%
sqrt-pow114.0%
metadata-eval14.0%
Applied egg-rr14.0%
expm1-def50.5%
expm1-log1p96.1%
Simplified96.1%
if -9.4999999999999997e48 < y < 1.50000000000000012e46Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.5%
cancel-sign-sub-inv97.5%
metadata-eval97.5%
associate-*r/97.4%
metadata-eval97.4%
+-commutative97.4%
Simplified97.4%
add-sqr-sqrt0.0%
sqrt-unprod56.4%
frac-times56.4%
metadata-eval56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-unprod56.4%
add-sqr-sqrt56.4%
clear-num56.4%
div-inv56.4%
metadata-eval56.4%
inv-pow56.4%
Applied egg-rr56.4%
inv-pow56.4%
add-sqr-sqrt56.4%
associate-/r*56.4%
metadata-eval56.4%
sqrt-div56.4%
inv-pow56.4%
inv-pow56.4%
metadata-eval56.4%
div-inv56.4%
clear-num56.4%
add-sqr-sqrt56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-prod0.0%
add-sqr-sqrt97.4%
*-commutative97.4%
sqrt-prod97.3%
metadata-eval97.3%
Applied egg-rr97.3%
associate-/l/97.4%
*-commutative97.4%
metadata-eval97.4%
associate-*r*97.4%
neg-mul-197.4%
associate-*l*97.4%
distribute-rgt-neg-out97.4%
*-commutative97.4%
associate-*r*97.4%
rem-square-sqrt97.5%
distribute-rgt-neg-in97.5%
metadata-eval97.5%
Simplified97.5%
if 1.50000000000000012e46 < y Initial program 99.5%
sub-neg99.5%
distribute-frac-neg99.5%
+-commutative99.5%
associate-+r-99.5%
+-commutative99.5%
associate-+r-99.5%
neg-mul-199.5%
*-commutative99.5%
associate-*r/99.3%
fma-neg99.3%
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 inf 91.1%
Final simplification95.9%
(FPCore (x y)
:precision binary64
(let* ((t_0 (sqrt (/ 1.0 x))))
(if (<= y -8e+48)
(+ 1.0 (* t_0 (* y -0.3333333333333333)))
(if (<= y 6.2e+44)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(+ 1.0 (* -0.3333333333333333 (* y t_0)))))))
double code(double x, double y) {
double t_0 = sqrt((1.0 / x));
double tmp;
if (y <= -8e+48) {
tmp = 1.0 + (t_0 * (y * -0.3333333333333333));
} else if (y <= 6.2e+44) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = 1.0 + (-0.3333333333333333 * (y * 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 = sqrt((1.0d0 / x))
if (y <= (-8d+48)) then
tmp = 1.0d0 + (t_0 * (y * (-0.3333333333333333d0)))
else if (y <= 6.2d+44) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = 1.0d0 + ((-0.3333333333333333d0) * (y * t_0))
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = Math.sqrt((1.0 / x));
double tmp;
if (y <= -8e+48) {
tmp = 1.0 + (t_0 * (y * -0.3333333333333333));
} else if (y <= 6.2e+44) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = 1.0 + (-0.3333333333333333 * (y * t_0));
}
return tmp;
}
def code(x, y): t_0 = math.sqrt((1.0 / x)) tmp = 0 if y <= -8e+48: tmp = 1.0 + (t_0 * (y * -0.3333333333333333)) elif y <= 6.2e+44: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = 1.0 + (-0.3333333333333333 * (y * t_0)) return tmp
function code(x, y) t_0 = sqrt(Float64(1.0 / x)) tmp = 0.0 if (y <= -8e+48) tmp = Float64(1.0 + Float64(t_0 * Float64(y * -0.3333333333333333))); elseif (y <= 6.2e+44) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(1.0 + Float64(-0.3333333333333333 * Float64(y * t_0))); end return tmp end
function tmp_2 = code(x, y) t_0 = sqrt((1.0 / x)); tmp = 0.0; if (y <= -8e+48) tmp = 1.0 + (t_0 * (y * -0.3333333333333333)); elseif (y <= 6.2e+44) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = 1.0 + (-0.3333333333333333 * (y * t_0)); end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[y, -8e+48], N[(1.0 + N[(t$95$0 * N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.2e+44], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.3333333333333333 * N[(y * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\frac{1}{x}}\\
\mathbf{if}\;y \leq -8 \cdot 10^{+48}:\\
\;\;\;\;1 + t_0 \cdot \left(y \cdot -0.3333333333333333\right)\\
\mathbf{elif}\;y \leq 6.2 \cdot 10^{+44}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;1 + -0.3333333333333333 \cdot \left(y \cdot t_0\right)\\
\end{array}
\end{array}
if y < -8.00000000000000035e48Initial program 99.6%
sub-neg99.6%
distribute-frac-neg99.6%
+-commutative99.6%
associate-+r-99.6%
+-commutative99.6%
associate-+r-99.6%
neg-mul-199.6%
*-commutative99.6%
associate-*r/99.5%
fma-neg99.5%
associate-/r*99.5%
metadata-eval99.5%
*-commutative99.5%
associate-/r*99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 96.1%
*-commutative96.1%
associate-*l*96.1%
*-commutative96.1%
Simplified96.1%
if -8.00000000000000035e48 < y < 6.19999999999999991e44Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.5%
cancel-sign-sub-inv97.5%
metadata-eval97.5%
associate-*r/97.4%
metadata-eval97.4%
+-commutative97.4%
Simplified97.4%
add-sqr-sqrt0.0%
sqrt-unprod56.4%
frac-times56.4%
metadata-eval56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-unprod56.4%
add-sqr-sqrt56.4%
clear-num56.4%
div-inv56.4%
metadata-eval56.4%
inv-pow56.4%
Applied egg-rr56.4%
inv-pow56.4%
add-sqr-sqrt56.4%
associate-/r*56.4%
metadata-eval56.4%
sqrt-div56.4%
inv-pow56.4%
inv-pow56.4%
metadata-eval56.4%
div-inv56.4%
clear-num56.4%
add-sqr-sqrt56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-prod0.0%
add-sqr-sqrt97.4%
*-commutative97.4%
sqrt-prod97.3%
metadata-eval97.3%
Applied egg-rr97.3%
associate-/l/97.4%
*-commutative97.4%
metadata-eval97.4%
associate-*r*97.4%
neg-mul-197.4%
associate-*l*97.4%
distribute-rgt-neg-out97.4%
*-commutative97.4%
associate-*r*97.4%
rem-square-sqrt97.5%
distribute-rgt-neg-in97.5%
metadata-eval97.5%
Simplified97.5%
if 6.19999999999999991e44 < y Initial program 99.5%
sub-neg99.5%
distribute-frac-neg99.5%
+-commutative99.5%
associate-+r-99.5%
+-commutative99.5%
associate-+r-99.5%
neg-mul-199.5%
*-commutative99.5%
associate-*r/99.3%
fma-neg99.3%
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 inf 91.1%
Final simplification95.9%
(FPCore (x y) :precision binary64 (if (or (<= y -4.2e+48) (not (<= y 3.4e+45))) (+ 1.0 (* y (* -0.3333333333333333 (pow x -0.5)))) (+ 1.0 (/ -0.3333333333333333 (* x 3.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -4.2e+48) || !(y <= 3.4e+45)) {
tmp = 1.0 + (y * (-0.3333333333333333 * pow(x, -0.5)));
} else {
tmp = 1.0 + (-0.3333333333333333 / (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.2d+48)) .or. (.not. (y <= 3.4d+45))) then
tmp = 1.0d0 + (y * ((-0.3333333333333333d0) * (x ** (-0.5d0))))
else
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -4.2e+48) || !(y <= 3.4e+45)) {
tmp = 1.0 + (y * (-0.3333333333333333 * Math.pow(x, -0.5)));
} else {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -4.2e+48) or not (y <= 3.4e+45): tmp = 1.0 + (y * (-0.3333333333333333 * math.pow(x, -0.5))) else: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) return tmp
function code(x, y) tmp = 0.0 if ((y <= -4.2e+48) || !(y <= 3.4e+45)) tmp = Float64(1.0 + Float64(y * Float64(-0.3333333333333333 * (x ^ -0.5)))); else tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -4.2e+48) || ~((y <= 3.4e+45))) tmp = 1.0 + (y * (-0.3333333333333333 * (x ^ -0.5))); else tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -4.2e+48], N[Not[LessEqual[y, 3.4e+45]], $MachinePrecision]], N[(1.0 + N[(y * N[(-0.3333333333333333 * N[Power[x, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.2 \cdot 10^{+48} \lor \neg \left(y \leq 3.4 \cdot 10^{+45}\right):\\
\;\;\;\;1 + y \cdot \left(-0.3333333333333333 \cdot {x}^{-0.5}\right)\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\end{array}
\end{array}
if y < -4.1999999999999997e48 or 3.4e45 < y Initial program 99.5%
sub-neg99.5%
distribute-frac-neg99.5%
+-commutative99.5%
associate-+r-99.5%
+-commutative99.5%
associate-+r-99.5%
neg-mul-199.5%
*-commutative99.5%
associate-*r/99.4%
fma-neg99.4%
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 inf 93.8%
associate-*r*93.7%
Simplified93.7%
expm1-log1p-u51.9%
expm1-udef11.6%
inv-pow11.6%
sqrt-pow111.6%
metadata-eval11.6%
Applied egg-rr11.6%
expm1-def51.9%
expm1-log1p93.7%
Simplified93.7%
if -4.1999999999999997e48 < y < 3.4e45Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 97.5%
cancel-sign-sub-inv97.5%
metadata-eval97.5%
associate-*r/97.4%
metadata-eval97.4%
+-commutative97.4%
Simplified97.4%
add-sqr-sqrt0.0%
sqrt-unprod56.4%
frac-times56.4%
metadata-eval56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-unprod56.4%
add-sqr-sqrt56.4%
clear-num56.4%
div-inv56.4%
metadata-eval56.4%
inv-pow56.4%
Applied egg-rr56.4%
inv-pow56.4%
add-sqr-sqrt56.4%
associate-/r*56.4%
metadata-eval56.4%
sqrt-div56.4%
inv-pow56.4%
inv-pow56.4%
metadata-eval56.4%
div-inv56.4%
clear-num56.4%
add-sqr-sqrt56.4%
metadata-eval56.4%
frac-times56.4%
sqrt-prod0.0%
add-sqr-sqrt97.4%
*-commutative97.4%
sqrt-prod97.3%
metadata-eval97.3%
Applied egg-rr97.3%
associate-/l/97.4%
*-commutative97.4%
metadata-eval97.4%
associate-*r*97.4%
neg-mul-197.4%
associate-*l*97.4%
distribute-rgt-neg-out97.4%
*-commutative97.4%
associate-*r*97.4%
rem-square-sqrt97.5%
distribute-rgt-neg-in97.5%
metadata-eval97.5%
Simplified97.5%
Final simplification95.8%
(FPCore (x y)
:precision binary64
(if (<= y -7.6e+48)
(/ -0.3333333333333333 (/ (sqrt x) y))
(if (<= y 1.52e+66)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(* -0.3333333333333333 (* y (pow x -0.5))))))
double code(double x, double y) {
double tmp;
if (y <= -7.6e+48) {
tmp = -0.3333333333333333 / (sqrt(x) / y);
} else if (y <= 1.52e+66) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * pow(x, -0.5));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-7.6d+48)) then
tmp = (-0.3333333333333333d0) / (sqrt(x) / y)
else if (y <= 1.52d+66) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = (-0.3333333333333333d0) * (y * (x ** (-0.5d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -7.6e+48) {
tmp = -0.3333333333333333 / (Math.sqrt(x) / y);
} else if (y <= 1.52e+66) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * Math.pow(x, -0.5));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -7.6e+48: tmp = -0.3333333333333333 / (math.sqrt(x) / y) elif y <= 1.52e+66: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = -0.3333333333333333 * (y * math.pow(x, -0.5)) return tmp
function code(x, y) tmp = 0.0 if (y <= -7.6e+48) tmp = Float64(-0.3333333333333333 / Float64(sqrt(x) / y)); elseif (y <= 1.52e+66) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(-0.3333333333333333 * Float64(y * (x ^ -0.5))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -7.6e+48) tmp = -0.3333333333333333 / (sqrt(x) / y); elseif (y <= 1.52e+66) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = -0.3333333333333333 * (y * (x ^ -0.5)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -7.6e+48], N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.52e+66], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(y * N[Power[x, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.6 \cdot 10^{+48}:\\
\;\;\;\;\frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}}\\
\mathbf{elif}\;y \leq 1.52 \cdot 10^{+66}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \left(y \cdot {x}^{-0.5}\right)\\
\end{array}
\end{array}
if y < -7.60000000000000001e48Initial 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 y around inf 87.3%
*-commutative87.3%
Simplified87.3%
*-commutative87.3%
sqrt-div87.3%
metadata-eval87.3%
div-inv87.4%
*-commutative87.4%
clear-num87.3%
un-div-inv87.4%
Applied egg-rr87.4%
if -7.60000000000000001e48 < y < 1.52000000000000004e66Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
if 1.52000000000000004e66 < y Initial program 99.4%
*-commutative99.4%
metadata-eval99.4%
sqrt-prod99.5%
pow1/299.5%
Applied egg-rr99.5%
unpow1/299.5%
Simplified99.5%
Taylor expanded in y around inf 88.8%
*-commutative88.8%
Simplified88.8%
inv-pow88.8%
sqrt-pow188.9%
metadata-eval88.9%
expm1-log1p-u87.9%
expm1-udef41.7%
Applied egg-rr41.7%
expm1-def87.9%
expm1-log1p88.9%
Simplified88.9%
Final simplification92.2%
(FPCore (x y)
:precision binary64
(if (<= y -1.45e+49)
(* y (* -0.3333333333333333 (sqrt (/ 1.0 x))))
(if (<= y 2.1e+66)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(* -0.3333333333333333 (* y (pow x -0.5))))))
double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = y * (-0.3333333333333333 * sqrt((1.0 / x)));
} else if (y <= 2.1e+66) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * pow(x, -0.5));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-1.45d+49)) then
tmp = y * ((-0.3333333333333333d0) * sqrt((1.0d0 / x)))
else if (y <= 2.1d+66) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = (-0.3333333333333333d0) * (y * (x ** (-0.5d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = y * (-0.3333333333333333 * Math.sqrt((1.0 / x)));
} else if (y <= 2.1e+66) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * Math.pow(x, -0.5));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.45e+49: tmp = y * (-0.3333333333333333 * math.sqrt((1.0 / x))) elif y <= 2.1e+66: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = -0.3333333333333333 * (y * math.pow(x, -0.5)) return tmp
function code(x, y) tmp = 0.0 if (y <= -1.45e+49) tmp = Float64(y * Float64(-0.3333333333333333 * sqrt(Float64(1.0 / x)))); elseif (y <= 2.1e+66) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(-0.3333333333333333 * Float64(y * (x ^ -0.5))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.45e+49) tmp = y * (-0.3333333333333333 * sqrt((1.0 / x))); elseif (y <= 2.1e+66) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = -0.3333333333333333 * (y * (x ^ -0.5)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.45e+49], N[(y * N[(-0.3333333333333333 * N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.1e+66], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(y * N[Power[x, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+49}:\\
\;\;\;\;y \cdot \left(-0.3333333333333333 \cdot \sqrt{\frac{1}{x}}\right)\\
\mathbf{elif}\;y \leq 2.1 \cdot 10^{+66}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \left(y \cdot {x}^{-0.5}\right)\\
\end{array}
\end{array}
if y < -1.45e49Initial 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 y around inf 87.3%
*-commutative87.3%
Simplified87.3%
Taylor expanded in y around 0 87.3%
*-commutative87.3%
*-commutative87.3%
associate-*l*87.4%
Simplified87.4%
if -1.45e49 < y < 2.10000000000000005e66Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
if 2.10000000000000005e66 < y Initial program 99.4%
*-commutative99.4%
metadata-eval99.4%
sqrt-prod99.5%
pow1/299.5%
Applied egg-rr99.5%
unpow1/299.5%
Simplified99.5%
Taylor expanded in y around inf 88.8%
*-commutative88.8%
Simplified88.8%
inv-pow88.8%
sqrt-pow188.9%
metadata-eval88.9%
expm1-log1p-u87.9%
expm1-udef41.7%
Applied egg-rr41.7%
expm1-def87.9%
expm1-log1p88.9%
Simplified88.9%
Final simplification92.2%
(FPCore (x y)
:precision binary64
(if (<= y -1.45e+49)
(* (sqrt (/ 1.0 x)) (* y -0.3333333333333333))
(if (<= y 4.6e+69)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(* -0.3333333333333333 (* y (pow x -0.5))))))
double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = sqrt((1.0 / x)) * (y * -0.3333333333333333);
} else if (y <= 4.6e+69) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * pow(x, -0.5));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-1.45d+49)) then
tmp = sqrt((1.0d0 / x)) * (y * (-0.3333333333333333d0))
else if (y <= 4.6d+69) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = (-0.3333333333333333d0) * (y * (x ** (-0.5d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = Math.sqrt((1.0 / x)) * (y * -0.3333333333333333);
} else if (y <= 4.6e+69) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y * Math.pow(x, -0.5));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.45e+49: tmp = math.sqrt((1.0 / x)) * (y * -0.3333333333333333) elif y <= 4.6e+69: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = -0.3333333333333333 * (y * math.pow(x, -0.5)) return tmp
function code(x, y) tmp = 0.0 if (y <= -1.45e+49) tmp = Float64(sqrt(Float64(1.0 / x)) * Float64(y * -0.3333333333333333)); elseif (y <= 4.6e+69) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(-0.3333333333333333 * Float64(y * (x ^ -0.5))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.45e+49) tmp = sqrt((1.0 / x)) * (y * -0.3333333333333333); elseif (y <= 4.6e+69) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = -0.3333333333333333 * (y * (x ^ -0.5)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.45e+49], N[(N[Sqrt[N[(1.0 / x), $MachinePrecision]], $MachinePrecision] * N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.6e+69], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(y * N[Power[x, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+49}:\\
\;\;\;\;\sqrt{\frac{1}{x}} \cdot \left(y \cdot -0.3333333333333333\right)\\
\mathbf{elif}\;y \leq 4.6 \cdot 10^{+69}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \left(y \cdot {x}^{-0.5}\right)\\
\end{array}
\end{array}
if y < -1.45e49Initial 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 y around inf 87.3%
*-commutative87.3%
associate-*l*87.4%
Simplified87.4%
if -1.45e49 < y < 4.60000000000000033e69Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
if 4.60000000000000033e69 < y Initial program 99.4%
*-commutative99.4%
metadata-eval99.4%
sqrt-prod99.5%
pow1/299.5%
Applied egg-rr99.5%
unpow1/299.5%
Simplified99.5%
Taylor expanded in y around inf 88.8%
*-commutative88.8%
Simplified88.8%
inv-pow88.8%
sqrt-pow188.9%
metadata-eval88.9%
expm1-log1p-u87.9%
expm1-udef41.7%
Applied egg-rr41.7%
expm1-def87.9%
expm1-log1p88.9%
Simplified88.9%
Final simplification92.2%
(FPCore (x y) :precision binary64 (if (or (<= y -1.45e+49) (not (<= y 1.3e+66))) (* -0.3333333333333333 (/ y (sqrt x))) (+ 1.0 (/ -0.3333333333333333 (* x 3.0)))))
double code(double x, double y) {
double tmp;
if ((y <= -1.45e+49) || !(y <= 1.3e+66)) {
tmp = -0.3333333333333333 * (y / sqrt(x));
} else {
tmp = 1.0 + (-0.3333333333333333 / (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 <= (-1.45d+49)) .or. (.not. (y <= 1.3d+66))) then
tmp = (-0.3333333333333333d0) * (y / sqrt(x))
else
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -1.45e+49) || !(y <= 1.3e+66)) {
tmp = -0.3333333333333333 * (y / Math.sqrt(x));
} else {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -1.45e+49) or not (y <= 1.3e+66): tmp = -0.3333333333333333 * (y / math.sqrt(x)) else: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) return tmp
function code(x, y) tmp = 0.0 if ((y <= -1.45e+49) || !(y <= 1.3e+66)) tmp = Float64(-0.3333333333333333 * Float64(y / sqrt(x))); else tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -1.45e+49) || ~((y <= 1.3e+66))) tmp = -0.3333333333333333 * (y / sqrt(x)); else tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -1.45e+49], N[Not[LessEqual[y, 1.3e+66]], $MachinePrecision]], N[(-0.3333333333333333 * N[(y / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+49} \lor \neg \left(y \leq 1.3 \cdot 10^{+66}\right):\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{\sqrt{x}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\end{array}
\end{array}
if y < -1.45e49 or 1.30000000000000006e66 < y Initial program 99.5%
*-commutative99.5%
metadata-eval99.5%
sqrt-prod99.6%
pow1/299.6%
Applied egg-rr99.6%
unpow1/299.6%
Simplified99.6%
Taylor expanded in y around inf 87.9%
*-commutative87.9%
Simplified87.9%
*-commutative87.9%
sqrt-div87.9%
metadata-eval87.9%
div-inv87.9%
expm1-log1p-u34.5%
expm1-udef34.4%
Applied egg-rr34.4%
expm1-def34.5%
expm1-log1p87.9%
Simplified87.9%
if -1.45e49 < y < 1.30000000000000006e66Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
Final simplification92.2%
(FPCore (x y)
:precision binary64
(if (<= y -1.45e+49)
(/ -0.3333333333333333 (/ (sqrt x) y))
(if (<= y 9.2e+69)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(* -0.3333333333333333 (/ y (sqrt x))))))
double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = -0.3333333333333333 / (sqrt(x) / y);
} else if (y <= 9.2e+69) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y / sqrt(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.45d+49)) then
tmp = (-0.3333333333333333d0) / (sqrt(x) / y)
else if (y <= 9.2d+69) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = (-0.3333333333333333d0) * (y / sqrt(x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = -0.3333333333333333 / (Math.sqrt(x) / y);
} else if (y <= 9.2e+69) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = -0.3333333333333333 * (y / Math.sqrt(x));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.45e+49: tmp = -0.3333333333333333 / (math.sqrt(x) / y) elif y <= 9.2e+69: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = -0.3333333333333333 * (y / math.sqrt(x)) return tmp
function code(x, y) tmp = 0.0 if (y <= -1.45e+49) tmp = Float64(-0.3333333333333333 / Float64(sqrt(x) / y)); elseif (y <= 9.2e+69) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(-0.3333333333333333 * Float64(y / sqrt(x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.45e+49) tmp = -0.3333333333333333 / (sqrt(x) / y); elseif (y <= 9.2e+69) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = -0.3333333333333333 * (y / sqrt(x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.45e+49], N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9.2e+69], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(y / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+49}:\\
\;\;\;\;\frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}}\\
\mathbf{elif}\;y \leq 9.2 \cdot 10^{+69}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{\sqrt{x}}\\
\end{array}
\end{array}
if y < -1.45e49Initial 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 y around inf 87.3%
*-commutative87.3%
Simplified87.3%
*-commutative87.3%
sqrt-div87.3%
metadata-eval87.3%
div-inv87.4%
*-commutative87.4%
clear-num87.3%
un-div-inv87.4%
Applied egg-rr87.4%
if -1.45e49 < y < 9.20000000000000067e69Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
if 9.20000000000000067e69 < y Initial program 99.4%
*-commutative99.4%
metadata-eval99.4%
sqrt-prod99.5%
pow1/299.5%
Applied egg-rr99.5%
unpow1/299.5%
Simplified99.5%
Taylor expanded in y around inf 88.8%
*-commutative88.8%
Simplified88.8%
*-commutative88.8%
sqrt-div88.8%
metadata-eval88.8%
div-inv88.7%
expm1-log1p-u82.5%
expm1-udef82.5%
Applied egg-rr82.5%
expm1-def82.5%
expm1-log1p88.7%
Simplified88.7%
Final simplification92.2%
(FPCore (x y)
:precision binary64
(if (<= y -1.45e+49)
(/ -0.3333333333333333 (/ (sqrt x) y))
(if (<= y 8.9e+67)
(+ 1.0 (/ -0.3333333333333333 (* x 3.0)))
(/ (* y -0.3333333333333333) (sqrt x)))))
double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = -0.3333333333333333 / (sqrt(x) / y);
} else if (y <= 8.9e+67) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = (y * -0.3333333333333333) / sqrt(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.45d+49)) then
tmp = (-0.3333333333333333d0) / (sqrt(x) / y)
else if (y <= 8.9d+67) then
tmp = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
else
tmp = (y * (-0.3333333333333333d0)) / sqrt(x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.45e+49) {
tmp = -0.3333333333333333 / (Math.sqrt(x) / y);
} else if (y <= 8.9e+67) {
tmp = 1.0 + (-0.3333333333333333 / (x * 3.0));
} else {
tmp = (y * -0.3333333333333333) / Math.sqrt(x);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.45e+49: tmp = -0.3333333333333333 / (math.sqrt(x) / y) elif y <= 8.9e+67: tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)) else: tmp = (y * -0.3333333333333333) / math.sqrt(x) return tmp
function code(x, y) tmp = 0.0 if (y <= -1.45e+49) tmp = Float64(-0.3333333333333333 / Float64(sqrt(x) / y)); elseif (y <= 8.9e+67) tmp = Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))); else tmp = Float64(Float64(y * -0.3333333333333333) / sqrt(x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.45e+49) tmp = -0.3333333333333333 / (sqrt(x) / y); elseif (y <= 8.9e+67) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); else tmp = (y * -0.3333333333333333) / sqrt(x); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.45e+49], N[(-0.3333333333333333 / N[(N[Sqrt[x], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.9e+67], N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * -0.3333333333333333), $MachinePrecision] / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+49}:\\
\;\;\;\;\frac{-0.3333333333333333}{\frac{\sqrt{x}}{y}}\\
\mathbf{elif}\;y \leq 8.9 \cdot 10^{+67}:\\
\;\;\;\;1 + \frac{-0.3333333333333333}{x \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot -0.3333333333333333}{\sqrt{x}}\\
\end{array}
\end{array}
if y < -1.45e49Initial 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 y around inf 87.3%
*-commutative87.3%
Simplified87.3%
*-commutative87.3%
sqrt-div87.3%
metadata-eval87.3%
div-inv87.4%
*-commutative87.4%
clear-num87.3%
un-div-inv87.4%
Applied egg-rr87.4%
if -1.45e49 < y < 8.89999999999999983e67Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/99.8%
fma-neg99.8%
associate-/r*99.8%
metadata-eval99.8%
*-commutative99.8%
associate-/r*99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around 0 95.1%
cancel-sign-sub-inv95.1%
metadata-eval95.1%
associate-*r/95.1%
metadata-eval95.1%
+-commutative95.1%
Simplified95.1%
add-sqr-sqrt0.0%
sqrt-unprod54.4%
frac-times54.4%
metadata-eval54.4%
metadata-eval54.4%
frac-times54.4%
sqrt-unprod54.5%
add-sqr-sqrt54.5%
clear-num54.5%
div-inv54.5%
metadata-eval54.5%
inv-pow54.5%
Applied egg-rr54.5%
inv-pow54.5%
add-sqr-sqrt54.5%
associate-/r*54.5%
metadata-eval54.5%
sqrt-div54.5%
inv-pow54.5%
inv-pow54.5%
metadata-eval54.5%
div-inv54.5%
clear-num54.5%
add-sqr-sqrt54.5%
metadata-eval54.5%
frac-times54.5%
sqrt-prod0.0%
add-sqr-sqrt95.1%
*-commutative95.1%
sqrt-prod95.0%
metadata-eval95.0%
Applied egg-rr95.0%
associate-/l/95.0%
*-commutative95.0%
metadata-eval95.0%
associate-*r*95.0%
neg-mul-195.0%
associate-*l*95.0%
distribute-rgt-neg-out95.0%
*-commutative95.0%
associate-*r*95.0%
rem-square-sqrt95.1%
distribute-rgt-neg-in95.1%
metadata-eval95.1%
Simplified95.1%
if 8.89999999999999983e67 < y Initial program 99.4%
*-commutative99.4%
metadata-eval99.4%
sqrt-prod99.5%
pow1/299.5%
Applied egg-rr99.5%
unpow1/299.5%
Simplified99.5%
Taylor expanded in y around inf 88.8%
*-commutative88.8%
Simplified88.8%
*-commutative88.8%
sqrt-div88.8%
metadata-eval88.8%
div-inv88.7%
associate-*l/88.7%
Applied egg-rr88.7%
Final simplification92.2%
(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.6%
metadata-eval99.6%
distribute-frac-neg99.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%
*-commutative99.7%
associate-/r*99.6%
metadata-eval99.6%
distribute-frac-neg99.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 (/ 0.1111111111111111 x)))))
double code(double x, double y) {
return (1.0 + (-0.1111111111111111 / x)) - (y * sqrt((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)) - (y * sqrt((0.1111111111111111d0 / x)))
end function
public static double code(double x, double y) {
return (1.0 + (-0.1111111111111111 / x)) - (y * Math.sqrt((0.1111111111111111 / x)));
}
def code(x, y): return (1.0 + (-0.1111111111111111 / x)) - (y * math.sqrt((0.1111111111111111 / x)))
function code(x, y) return Float64(Float64(1.0 + Float64(-0.1111111111111111 / x)) - Float64(y * sqrt(Float64(0.1111111111111111 / x)))) end
function tmp = code(x, y) tmp = (1.0 + (-0.1111111111111111 / x)) - (y * sqrt((0.1111111111111111 / x))); end
code[x_, y_] := N[(N[(1.0 + N[(-0.1111111111111111 / x), $MachinePrecision]), $MachinePrecision] - N[(y * N[Sqrt[N[(0.1111111111111111 / x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 + \frac{-0.1111111111111111}{x}\right) - y \cdot \sqrt{\frac{0.1111111111111111}{x}}
\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%
*-un-lft-identity99.7%
sqrt-prod99.7%
metadata-eval99.7%
times-frac99.7%
pow1/299.7%
pow-flip99.7%
metadata-eval99.7%
Applied egg-rr99.7%
Taylor expanded in x around 0 99.7%
associate--l-99.7%
sub-neg99.7%
*-commutative99.7%
div-inv99.6%
metadata-eval99.6%
metadata-eval99.6%
associate-*l*99.6%
metadata-eval99.6%
metadata-eval99.6%
metadata-eval99.6%
sqrt-pow199.6%
inv-pow99.6%
sqrt-prod99.6%
div-inv99.7%
Applied egg-rr99.7%
distribute-neg-in99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
associate-+r+99.7%
sub-neg99.7%
Simplified99.7%
Final simplification99.7%
(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.024) (* -0.1111111111111111 (/ 1.0 x)) 1.0))
double code(double x, double y) {
double tmp;
if (x <= 0.024) {
tmp = -0.1111111111111111 * (1.0 / 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.024d0) then
tmp = (-0.1111111111111111d0) * (1.0d0 / x)
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 0.024) {
tmp = -0.1111111111111111 * (1.0 / x);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 0.024: tmp = -0.1111111111111111 * (1.0 / x) else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (x <= 0.024) tmp = Float64(-0.1111111111111111 * Float64(1.0 / x)); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 0.024) tmp = -0.1111111111111111 * (1.0 / x); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 0.024], N[(-0.1111111111111111 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.024:\\
\;\;\;\;-0.1111111111111111 \cdot \frac{1}{x}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < 0.024Initial 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 56.6%
clear-num56.5%
associate-/r/56.6%
Applied egg-rr56.6%
if 0.024 < x Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/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 x around inf 63.0%
Final simplification60.2%
(FPCore (x y) :precision binary64 (if (<= x 0.024) (/ -0.1111111111111111 x) 1.0))
double code(double x, double y) {
double tmp;
if (x <= 0.024) {
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.024d0) 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.024) {
tmp = -0.1111111111111111 / x;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 0.024: tmp = -0.1111111111111111 / x else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (x <= 0.024) 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.024) tmp = -0.1111111111111111 / x; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 0.024], N[(-0.1111111111111111 / x), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.024:\\
\;\;\;\;\frac{-0.1111111111111111}{x}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if x < 0.024Initial 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 56.6%
if 0.024 < x Initial program 99.8%
sub-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
associate-+r-99.8%
+-commutative99.8%
associate-+r-99.8%
neg-mul-199.8%
*-commutative99.8%
associate-*r/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 x around inf 63.0%
Final simplification60.2%
(FPCore (x y) :precision binary64 (+ 1.0 (/ -0.3333333333333333 (* x 3.0))))
double code(double x, double y) {
return 1.0 + (-0.3333333333333333 / (x * 3.0));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 + ((-0.3333333333333333d0) / (x * 3.0d0))
end function
public static double code(double x, double y) {
return 1.0 + (-0.3333333333333333 / (x * 3.0));
}
def code(x, y): return 1.0 + (-0.3333333333333333 / (x * 3.0))
function code(x, y) return Float64(1.0 + Float64(-0.3333333333333333 / Float64(x * 3.0))) end
function tmp = code(x, y) tmp = 1.0 + (-0.3333333333333333 / (x * 3.0)); end
code[x_, y_] := N[(1.0 + N[(-0.3333333333333333 / N[(x * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \frac{-0.3333333333333333}{x \cdot 3}
\end{array}
Initial program 99.7%
sub-neg99.7%
distribute-frac-neg99.7%
+-commutative99.7%
associate-+r-99.7%
+-commutative99.7%
associate-+r-99.7%
neg-mul-199.7%
*-commutative99.7%
associate-*r/99.7%
fma-neg99.7%
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 60.9%
cancel-sign-sub-inv60.9%
metadata-eval60.9%
associate-*r/60.8%
metadata-eval60.8%
+-commutative60.8%
Simplified60.8%
add-sqr-sqrt0.0%
sqrt-unprod38.3%
frac-times38.3%
metadata-eval38.3%
metadata-eval38.3%
frac-times38.3%
sqrt-unprod36.1%
add-sqr-sqrt36.1%
clear-num36.1%
div-inv36.1%
metadata-eval36.1%
inv-pow36.1%
Applied egg-rr36.1%
inv-pow36.1%
add-sqr-sqrt36.1%
associate-/r*36.1%
metadata-eval36.1%
sqrt-div36.1%
inv-pow36.1%
inv-pow36.1%
metadata-eval36.1%
div-inv36.1%
clear-num36.1%
add-sqr-sqrt36.1%
metadata-eval36.1%
frac-times36.1%
sqrt-prod0.0%
add-sqr-sqrt60.8%
*-commutative60.8%
sqrt-prod60.8%
metadata-eval60.8%
Applied egg-rr60.8%
associate-/l/60.8%
*-commutative60.8%
metadata-eval60.8%
associate-*r*60.8%
neg-mul-160.8%
associate-*l*60.8%
distribute-rgt-neg-out60.8%
*-commutative60.8%
associate-*r*60.8%
rem-square-sqrt60.9%
distribute-rgt-neg-in60.9%
metadata-eval60.9%
Simplified60.9%
Final simplification60.9%
(FPCore (x y) :precision binary64 (+ 1.0 (* 0.1111111111111111 (/ -1.0 x))))
double code(double x, double y) {
return 1.0 + (0.1111111111111111 * (-1.0 / x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 + (0.1111111111111111d0 * ((-1.0d0) / x))
end function
public static double code(double x, double y) {
return 1.0 + (0.1111111111111111 * (-1.0 / x));
}
def code(x, y): return 1.0 + (0.1111111111111111 * (-1.0 / x))
function code(x, y) return Float64(1.0 + Float64(0.1111111111111111 * Float64(-1.0 / x))) end
function tmp = code(x, y) tmp = 1.0 + (0.1111111111111111 * (-1.0 / x)); end
code[x_, y_] := N[(1.0 + N[(0.1111111111111111 * N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + 0.1111111111111111 \cdot \frac{-1}{x}
\end{array}
Initial program 99.7%
sub-neg99.7%
distribute-frac-neg99.7%
+-commutative99.7%
associate-+r-99.7%
+-commutative99.7%
associate-+r-99.7%
neg-mul-199.7%
*-commutative99.7%
associate-*r/99.7%
fma-neg99.7%
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 60.9%
Final simplification60.9%
(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%
sub-neg99.7%
distribute-frac-neg99.7%
+-commutative99.7%
associate-+r-99.7%
+-commutative99.7%
associate-+r-99.7%
neg-mul-199.7%
*-commutative99.7%
associate-*r/99.7%
fma-neg99.7%
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 60.9%
cancel-sign-sub-inv60.9%
metadata-eval60.9%
associate-*r/60.8%
metadata-eval60.8%
+-commutative60.8%
Simplified60.8%
Final simplification60.8%
(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%
sub-neg99.7%
distribute-frac-neg99.7%
+-commutative99.7%
associate-+r-99.7%
+-commutative99.7%
associate-+r-99.7%
neg-mul-199.7%
*-commutative99.7%
associate-*r/99.7%
fma-neg99.7%
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 x around inf 36.0%
Final simplification36.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 2024021
(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)))))