
(FPCore (x) :precision binary64 (- (sqrt (+ x 1.0)) (sqrt x)))
double code(double x) {
return sqrt((x + 1.0)) - sqrt(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt((x + 1.0d0)) - sqrt(x)
end function
public static double code(double x) {
return Math.sqrt((x + 1.0)) - Math.sqrt(x);
}
def code(x): return math.sqrt((x + 1.0)) - math.sqrt(x)
function code(x) return Float64(sqrt(Float64(x + 1.0)) - sqrt(x)) end
function tmp = code(x) tmp = sqrt((x + 1.0)) - sqrt(x); end
code[x_] := N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sqrt{x + 1} - \sqrt{x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- (sqrt (+ x 1.0)) (sqrt x)))
double code(double x) {
return sqrt((x + 1.0)) - sqrt(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt((x + 1.0d0)) - sqrt(x)
end function
public static double code(double x) {
return Math.sqrt((x + 1.0)) - Math.sqrt(x);
}
def code(x): return math.sqrt((x + 1.0)) - math.sqrt(x)
function code(x) return Float64(sqrt(Float64(x + 1.0)) - sqrt(x)) end
function tmp = code(x) tmp = sqrt((x + 1.0)) - sqrt(x); end
code[x_] := N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sqrt{x + 1} - \sqrt{x}
\end{array}
(FPCore (x) :precision binary64 (/ 1.0 (+ (sqrt x) (sqrt (+ 1.0 x)))))
double code(double x) {
return 1.0 / (sqrt(x) + sqrt((1.0 + x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (sqrt(x) + sqrt((1.0d0 + x)))
end function
public static double code(double x) {
return 1.0 / (Math.sqrt(x) + Math.sqrt((1.0 + x)));
}
def code(x): return 1.0 / (math.sqrt(x) + math.sqrt((1.0 + x)))
function code(x) return Float64(1.0 / Float64(sqrt(x) + sqrt(Float64(1.0 + x)))) end
function tmp = code(x) tmp = 1.0 / (sqrt(x) + sqrt((1.0 + x))); end
code[x_] := N[(1.0 / N[(N[Sqrt[x], $MachinePrecision] + N[Sqrt[N[(1.0 + x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\sqrt{x} + \sqrt{1 + x}}
\end{array}
Initial program 50.0%
flip--50.5%
div-inv50.5%
add-sqr-sqrt50.3%
add-sqr-sqrt51.1%
Applied egg-rr51.1%
*-commutative51.1%
associate-/r/51.1%
+-commutative51.1%
associate--l+99.7%
+-inverses99.7%
metadata-eval99.7%
/-rgt-identity99.7%
+-commutative99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x) :precision binary64 (let* ((t_0 (- (sqrt (+ 1.0 x)) (sqrt x)))) (if (<= t_0 2e-5) (* (pow x -0.5) 0.5) t_0)))
double code(double x) {
double t_0 = sqrt((1.0 + x)) - sqrt(x);
double tmp;
if (t_0 <= 2e-5) {
tmp = pow(x, -0.5) * 0.5;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = sqrt((1.0d0 + x)) - sqrt(x)
if (t_0 <= 2d-5) then
tmp = (x ** (-0.5d0)) * 0.5d0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x) {
double t_0 = Math.sqrt((1.0 + x)) - Math.sqrt(x);
double tmp;
if (t_0 <= 2e-5) {
tmp = Math.pow(x, -0.5) * 0.5;
} else {
tmp = t_0;
}
return tmp;
}
def code(x): t_0 = math.sqrt((1.0 + x)) - math.sqrt(x) tmp = 0 if t_0 <= 2e-5: tmp = math.pow(x, -0.5) * 0.5 else: tmp = t_0 return tmp
function code(x) t_0 = Float64(sqrt(Float64(1.0 + x)) - sqrt(x)) tmp = 0.0 if (t_0 <= 2e-5) tmp = Float64((x ^ -0.5) * 0.5); else tmp = t_0; end return tmp end
function tmp_2 = code(x) t_0 = sqrt((1.0 + x)) - sqrt(x); tmp = 0.0; if (t_0 <= 2e-5) tmp = (x ^ -0.5) * 0.5; else tmp = t_0; end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(N[Sqrt[N[(1.0 + x), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 2e-5], N[(N[Power[x, -0.5], $MachinePrecision] * 0.5), $MachinePrecision], t$95$0]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{1 + x} - \sqrt{x}\\
\mathbf{if}\;t_0 \leq 2 \cdot 10^{-5}:\\
\;\;\;\;{x}^{-0.5} \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if (-.f64 (sqrt.f64 (+.f64 x 1)) (sqrt.f64 x)) < 2.00000000000000016e-5Initial program 5.0%
flip--5.7%
div-inv5.7%
add-sqr-sqrt5.7%
add-sqr-sqrt6.0%
Applied egg-rr6.0%
*-commutative6.0%
associate-/r/6.0%
+-commutative6.0%
associate--l+99.4%
+-inverses99.4%
metadata-eval99.4%
/-rgt-identity99.4%
+-commutative99.4%
Simplified99.4%
flip3-+67.2%
div-inv67.1%
+-commutative67.1%
+-commutative67.1%
sqrt-pow267.0%
metadata-eval67.0%
sqrt-pow266.9%
metadata-eval66.9%
add-sqr-sqrt67.1%
+-commutative67.1%
associate-+l+67.1%
Applied egg-rr49.3%
+-commutative49.3%
+-commutative49.3%
associate-+r-49.3%
*-commutative49.3%
+-commutative49.3%
Simplified49.3%
Taylor expanded in x around inf 67.1%
Taylor expanded in x around inf 99.2%
unpow1/299.2%
unpow-199.2%
exp-to-pow91.8%
*-commutative91.8%
neg-mul-191.8%
exp-prod91.8%
distribute-lft-neg-out91.8%
distribute-rgt-neg-in91.8%
metadata-eval91.8%
exp-to-pow99.4%
*-commutative99.4%
Simplified99.4%
if 2.00000000000000016e-5 < (-.f64 (sqrt.f64 (+.f64 x 1)) (sqrt.f64 x)) Initial program 98.6%
Final simplification99.0%
(FPCore (x) :precision binary64 (if (<= x 1.0) (/ 1.0 (+ 1.0 (sqrt x))) (* (pow x -0.5) 0.5)))
double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 1.0 / (1.0 + sqrt(x));
} else {
tmp = pow(x, -0.5) * 0.5;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.0d0) then
tmp = 1.0d0 / (1.0d0 + sqrt(x))
else
tmp = (x ** (-0.5d0)) * 0.5d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 1.0 / (1.0 + Math.sqrt(x));
} else {
tmp = Math.pow(x, -0.5) * 0.5;
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.0: tmp = 1.0 / (1.0 + math.sqrt(x)) else: tmp = math.pow(x, -0.5) * 0.5 return tmp
function code(x) tmp = 0.0 if (x <= 1.0) tmp = Float64(1.0 / Float64(1.0 + sqrt(x))); else tmp = Float64((x ^ -0.5) * 0.5); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.0) tmp = 1.0 / (1.0 + sqrt(x)); else tmp = (x ^ -0.5) * 0.5; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.0], N[(1.0 / N[(1.0 + N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[x, -0.5], $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\frac{1}{1 + \sqrt{x}}\\
\mathbf{else}:\\
\;\;\;\;{x}^{-0.5} \cdot 0.5\\
\end{array}
\end{array}
if x < 1Initial program 100.0%
flip--99.9%
div-inv99.9%
add-sqr-sqrt99.9%
add-sqr-sqrt99.9%
Applied egg-rr99.9%
*-commutative99.9%
associate-/r/99.9%
+-commutative99.9%
associate--l+99.9%
+-inverses99.9%
metadata-eval99.9%
/-rgt-identity99.9%
+-commutative99.9%
Simplified99.9%
inv-pow99.9%
add-sqr-sqrt99.9%
unpow-prod-down99.9%
pow-prod-up99.9%
+-commutative99.9%
add-sqr-sqrt99.9%
+-commutative99.9%
add-sqr-sqrt99.9%
hypot-def99.9%
pow1/299.9%
sqrt-pow199.9%
metadata-eval99.9%
pow1/299.9%
sqrt-pow199.9%
metadata-eval99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0 99.0%
if 1 < x Initial program 7.9%
flip--8.9%
div-inv8.9%
add-sqr-sqrt8.6%
add-sqr-sqrt10.1%
Applied egg-rr10.1%
*-commutative10.1%
associate-/r/10.1%
+-commutative10.1%
associate--l+99.4%
+-inverses99.4%
metadata-eval99.4%
/-rgt-identity99.4%
+-commutative99.4%
Simplified99.4%
flip3-+68.5%
div-inv68.5%
+-commutative68.5%
+-commutative68.5%
sqrt-pow268.4%
metadata-eval68.4%
sqrt-pow268.3%
metadata-eval68.3%
add-sqr-sqrt68.5%
+-commutative68.5%
associate-+l+68.5%
Applied egg-rr51.5%
+-commutative51.5%
+-commutative51.5%
associate-+r-51.5%
*-commutative51.5%
+-commutative51.5%
Simplified51.5%
Taylor expanded in x around inf 67.6%
Taylor expanded in x around inf 97.0%
unpow1/297.0%
unpow-197.0%
exp-to-pow90.0%
*-commutative90.0%
neg-mul-190.0%
exp-prod90.0%
distribute-lft-neg-out90.0%
distribute-rgt-neg-in90.0%
metadata-eval90.0%
exp-to-pow97.2%
*-commutative97.2%
Simplified97.2%
Final simplification98.0%
(FPCore (x) :precision binary64 (if (<= x 0.25) 1.0 (* (pow x -0.5) 0.5)))
double code(double x) {
double tmp;
if (x <= 0.25) {
tmp = 1.0;
} else {
tmp = pow(x, -0.5) * 0.5;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 0.25d0) then
tmp = 1.0d0
else
tmp = (x ** (-0.5d0)) * 0.5d0
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 0.25) {
tmp = 1.0;
} else {
tmp = Math.pow(x, -0.5) * 0.5;
}
return tmp;
}
def code(x): tmp = 0 if x <= 0.25: tmp = 1.0 else: tmp = math.pow(x, -0.5) * 0.5 return tmp
function code(x) tmp = 0.0 if (x <= 0.25) tmp = 1.0; else tmp = Float64((x ^ -0.5) * 0.5); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 0.25) tmp = 1.0; else tmp = (x ^ -0.5) * 0.5; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 0.25], 1.0, N[(N[Power[x, -0.5], $MachinePrecision] * 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.25:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;{x}^{-0.5} \cdot 0.5\\
\end{array}
\end{array}
if x < 0.25Initial program 100.0%
Taylor expanded in x around 0 97.1%
if 0.25 < x Initial program 7.9%
flip--8.9%
div-inv8.9%
add-sqr-sqrt8.6%
add-sqr-sqrt10.1%
Applied egg-rr10.1%
*-commutative10.1%
associate-/r/10.1%
+-commutative10.1%
associate--l+99.4%
+-inverses99.4%
metadata-eval99.4%
/-rgt-identity99.4%
+-commutative99.4%
Simplified99.4%
flip3-+68.5%
div-inv68.5%
+-commutative68.5%
+-commutative68.5%
sqrt-pow268.4%
metadata-eval68.4%
sqrt-pow268.3%
metadata-eval68.3%
add-sqr-sqrt68.5%
+-commutative68.5%
associate-+l+68.5%
Applied egg-rr51.5%
+-commutative51.5%
+-commutative51.5%
associate-+r-51.5%
*-commutative51.5%
+-commutative51.5%
Simplified51.5%
Taylor expanded in x around inf 67.6%
Taylor expanded in x around inf 97.0%
unpow1/297.0%
unpow-197.0%
exp-to-pow90.0%
*-commutative90.0%
neg-mul-190.0%
exp-prod90.0%
distribute-lft-neg-out90.0%
distribute-rgt-neg-in90.0%
metadata-eval90.0%
exp-to-pow97.2%
*-commutative97.2%
Simplified97.2%
Final simplification97.2%
(FPCore (x) :precision binary64 (if (<= x 0.25) 1.0 (sqrt (/ 0.25 x))))
double code(double x) {
double tmp;
if (x <= 0.25) {
tmp = 1.0;
} else {
tmp = sqrt((0.25 / x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 0.25d0) then
tmp = 1.0d0
else
tmp = sqrt((0.25d0 / x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 0.25) {
tmp = 1.0;
} else {
tmp = Math.sqrt((0.25 / x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 0.25: tmp = 1.0 else: tmp = math.sqrt((0.25 / x)) return tmp
function code(x) tmp = 0.0 if (x <= 0.25) tmp = 1.0; else tmp = sqrt(Float64(0.25 / x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 0.25) tmp = 1.0; else tmp = sqrt((0.25 / x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 0.25], 1.0, N[Sqrt[N[(0.25 / x), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.25:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{0.25}{x}}\\
\end{array}
\end{array}
if x < 0.25Initial program 100.0%
Taylor expanded in x around 0 97.1%
if 0.25 < x Initial program 7.9%
flip--8.9%
div-inv8.9%
add-sqr-sqrt8.6%
add-sqr-sqrt10.1%
Applied egg-rr10.1%
*-commutative10.1%
associate-/r/10.1%
+-commutative10.1%
associate--l+99.4%
+-inverses99.4%
metadata-eval99.4%
/-rgt-identity99.4%
+-commutative99.4%
Simplified99.4%
flip3-+68.5%
div-inv68.5%
+-commutative68.5%
+-commutative68.5%
sqrt-pow268.4%
metadata-eval68.4%
sqrt-pow268.3%
metadata-eval68.3%
add-sqr-sqrt68.5%
+-commutative68.5%
associate-+l+68.5%
Applied egg-rr51.5%
+-commutative51.5%
+-commutative51.5%
associate-+r-51.5%
*-commutative51.5%
+-commutative51.5%
Simplified51.5%
Taylor expanded in x around inf 96.7%
add-sqr-sqrt96.2%
sqrt-unprod96.7%
associate-/r*96.7%
metadata-eval96.7%
associate-/r*96.7%
metadata-eval96.7%
frac-times96.7%
metadata-eval96.7%
add-sqr-sqrt97.0%
Applied egg-rr97.0%
Final simplification97.1%
(FPCore (x) :precision binary64 1.0)
double code(double x) {
return 1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double x) {
return 1.0;
}
def code(x): return 1.0
function code(x) return 1.0 end
function tmp = code(x) tmp = 1.0; end
code[x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 50.0%
Taylor expanded in x around 0 48.1%
Final simplification48.1%
(FPCore (x) :precision binary64 (/ 1.0 (+ (sqrt (+ x 1.0)) (sqrt x))))
double code(double x) {
return 1.0 / (sqrt((x + 1.0)) + sqrt(x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 / (sqrt((x + 1.0d0)) + sqrt(x))
end function
public static double code(double x) {
return 1.0 / (Math.sqrt((x + 1.0)) + Math.sqrt(x));
}
def code(x): return 1.0 / (math.sqrt((x + 1.0)) + math.sqrt(x))
function code(x) return Float64(1.0 / Float64(sqrt(Float64(x + 1.0)) + sqrt(x))) end
function tmp = code(x) tmp = 1.0 / (sqrt((x + 1.0)) + sqrt(x)); end
code[x_] := N[(1.0 / N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] + N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\sqrt{x + 1} + \sqrt{x}}
\end{array}
herbie shell --seed 2023290
(FPCore (x)
:name "Main:bigenough3 from C"
:precision binary64
:herbie-target
(/ 1.0 (+ (sqrt (+ x 1.0)) (sqrt x)))
(- (sqrt (+ x 1.0)) (sqrt x)))