| Alternative 1 | |
|---|---|
| Accuracy | 99.7% |
| Cost | 19584 |
\[\frac{1}{\sqrt{x} + \mathsf{hypot}\left(1, \sqrt{x}\right)}
\]

(FPCore (x) :precision binary64 (- (sqrt (+ x 1.0)) (sqrt x)))
(FPCore (x) :precision binary64 (/ 1.0 (+ (sqrt x) (hypot 1.0 (sqrt x)))))
double code(double x) {
return sqrt((x + 1.0)) - sqrt(x);
}
double code(double x) {
return 1.0 / (sqrt(x) + hypot(1.0, sqrt(x)));
}
public static double code(double x) {
return Math.sqrt((x + 1.0)) - Math.sqrt(x);
}
public static double code(double x) {
return 1.0 / (Math.sqrt(x) + Math.hypot(1.0, Math.sqrt(x)));
}
def code(x): return math.sqrt((x + 1.0)) - math.sqrt(x)
def code(x): return 1.0 / (math.sqrt(x) + math.hypot(1.0, math.sqrt(x)))
function code(x) return Float64(sqrt(Float64(x + 1.0)) - sqrt(x)) end
function code(x) return Float64(1.0 / Float64(sqrt(x) + hypot(1.0, sqrt(x)))) end
function tmp = code(x) tmp = sqrt((x + 1.0)) - sqrt(x); end
function tmp = code(x) tmp = 1.0 / (sqrt(x) + hypot(1.0, sqrt(x))); end
code[x_] := N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision]
code[x_] := N[(1.0 / N[(N[Sqrt[x], $MachinePrecision] + N[Sqrt[1.0 ^ 2 + N[Sqrt[x], $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\sqrt{x + 1} - \sqrt{x}
\begin{array}{l}
\\
\frac{1}{\sqrt{x} + \mathsf{hypot}\left(1, \sqrt{x}\right)}
\end{array}
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
Results
| Original | 52.1% |
|---|---|
| Target | 99.7% |
| Herbie | 99.7% |
Initial program 56.5%
Applied egg-rr57.1%
[Start]56.5% | \[ \sqrt{x + 1} - \sqrt{x}
\] |
|---|---|
flip-- [=>]56.7% | \[ \color{blue}{\frac{\sqrt{x + 1} \cdot \sqrt{x + 1} - \sqrt{x} \cdot \sqrt{x}}{\sqrt{x + 1} + \sqrt{x}}}
\] |
div-inv [=>]56.7% | \[ \color{blue}{\left(\sqrt{x + 1} \cdot \sqrt{x + 1} - \sqrt{x} \cdot \sqrt{x}\right) \cdot \frac{1}{\sqrt{x + 1} + \sqrt{x}}}
\] |
add-sqr-sqrt [<=]56.8% | \[ \left(\color{blue}{\left(x + 1\right)} - \sqrt{x} \cdot \sqrt{x}\right) \cdot \frac{1}{\sqrt{x + 1} + \sqrt{x}}
\] |
add-sqr-sqrt [<=]57.1% | \[ \left(\left(x + 1\right) - \color{blue}{x}\right) \cdot \frac{1}{\sqrt{x + 1} + \sqrt{x}}
\] |
Simplified99.7%
[Start]57.1% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\sqrt{x + 1} + \sqrt{x}}
\] |
|---|---|
associate-*r/ [=>]57.1% | \[ \color{blue}{\frac{\left(\left(x + 1\right) - x\right) \cdot 1}{\sqrt{x + 1} + \sqrt{x}}}
\] |
*-rgt-identity [=>]57.1% | \[ \frac{\color{blue}{\left(x + 1\right) - x}}{\sqrt{x + 1} + \sqrt{x}}
\] |
remove-double-neg [<=]57.1% | \[ \frac{\left(x + 1\right) - x}{\sqrt{x + 1} + \color{blue}{\left(-\left(-\sqrt{x}\right)\right)}}
\] |
sub-neg [<=]57.1% | \[ \frac{\left(x + 1\right) - x}{\color{blue}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}}
\] |
div-sub [=>]56.5% | \[ \color{blue}{\frac{x + 1}{\sqrt{x + 1} - \left(-\sqrt{x}\right)} - \frac{x}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}}
\] |
rem-square-sqrt [<=]56.4% | \[ \frac{x + 1}{\sqrt{x + 1} - \left(-\sqrt{x}\right)} - \frac{\color{blue}{\sqrt{x} \cdot \sqrt{x}}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
sqr-neg [<=]56.4% | \[ \frac{x + 1}{\sqrt{x + 1} - \left(-\sqrt{x}\right)} - \frac{\color{blue}{\left(-\sqrt{x}\right) \cdot \left(-\sqrt{x}\right)}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
div-sub [<=]56.8% | \[ \color{blue}{\frac{\left(x + 1\right) - \left(-\sqrt{x}\right) \cdot \left(-\sqrt{x}\right)}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}}
\] |
sqr-neg [=>]56.8% | \[ \frac{\left(x + 1\right) - \color{blue}{\sqrt{x} \cdot \sqrt{x}}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
+-commutative [=>]56.8% | \[ \frac{\color{blue}{\left(1 + x\right)} - \sqrt{x} \cdot \sqrt{x}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
rem-square-sqrt [=>]57.1% | \[ \frac{\left(1 + x\right) - \color{blue}{x}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
associate--l+ [=>]99.7% | \[ \frac{\color{blue}{1 + \left(x - x\right)}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
+-inverses [=>]99.7% | \[ \frac{1 + \color{blue}{0}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
metadata-eval [=>]99.7% | \[ \frac{\color{blue}{1}}{\sqrt{x + 1} - \left(-\sqrt{x}\right)}
\] |
sub-neg [=>]99.7% | \[ \frac{1}{\color{blue}{\sqrt{x + 1} + \left(-\left(-\sqrt{x}\right)\right)}}
\] |
Applied egg-rr99.7%
[Start]99.7% | \[ \frac{1}{\sqrt{1 + x} + \sqrt{x}}
\] |
|---|---|
add-sqr-sqrt [=>]99.7% | \[ \frac{1}{\sqrt{1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}} + \sqrt{x}}
\] |
hypot-1-def [=>]99.7% | \[ \frac{1}{\color{blue}{\mathsf{hypot}\left(1, \sqrt{x}\right)} + \sqrt{x}}
\] |
Final simplification99.7%
| Alternative 1 | |
|---|---|
| Accuracy | 99.7% |
| Cost | 19584 |
| Alternative 2 | |
|---|---|
| Accuracy | 99.4% |
| Cost | 26308 |
| Alternative 3 | |
|---|---|
| Accuracy | 99.7% |
| Cost | 13248 |
| Alternative 4 | |
|---|---|
| Accuracy | 98.6% |
| Cost | 6980 |
| Alternative 5 | |
|---|---|
| Accuracy | 96.9% |
| Cost | 6852 |
| Alternative 6 | |
|---|---|
| Accuracy | 98.0% |
| Cost | 6852 |
| Alternative 7 | |
|---|---|
| Accuracy | 50.1% |
| Cost | 64 |
herbie shell --seed 2023167
(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)))