Alternative 1 | |
---|---|
Error | 0.3 |
Cost | 13376 |
\[\frac{1}{\frac{1}{\sqrt{1 + x} + \sqrt{x}}}
\]
(FPCore (x) :precision binary64 (/ 1.0 (- (sqrt (+ x 1.0)) (sqrt x))))
(FPCore (x) :precision binary64 (/ (+ (sqrt (+ 1.0 x)) (sqrt x)) (+ x (- 1.0 x))))
double code(double x) { return 1.0 / (sqrt((x + 1.0)) - sqrt(x)); }
double code(double x) { return (sqrt((1.0 + x)) + sqrt(x)) / (x + (1.0 - x)); }
real(8) function code(x) real(8), intent (in) :: x code = 1.0d0 / (sqrt((x + 1.0d0)) - sqrt(x)) end function
real(8) function code(x) real(8), intent (in) :: x code = (sqrt((1.0d0 + x)) + sqrt(x)) / (x + (1.0d0 - x)) end function
public static double code(double x) { return 1.0 / (Math.sqrt((x + 1.0)) - Math.sqrt(x)); }
public static double code(double x) { return (Math.sqrt((1.0 + x)) + Math.sqrt(x)) / (x + (1.0 - x)); }
def code(x): return 1.0 / (math.sqrt((x + 1.0)) - math.sqrt(x))
def code(x): return (math.sqrt((1.0 + x)) + math.sqrt(x)) / (x + (1.0 - x))
function code(x) return Float64(1.0 / Float64(sqrt(Float64(x + 1.0)) - sqrt(x))) end
function code(x) return Float64(Float64(sqrt(Float64(1.0 + x)) + sqrt(x)) / Float64(x + Float64(1.0 - x))) end
function tmp = code(x) tmp = 1.0 / (sqrt((x + 1.0)) - sqrt(x)); end
function tmp = code(x) tmp = (sqrt((1.0 + x)) + sqrt(x)) / (x + (1.0 - x)); end
code[x_] := N[(1.0 / N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_] := N[(N[(N[Sqrt[N[(1.0 + x), $MachinePrecision]], $MachinePrecision] + N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / N[(x + N[(1.0 - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{1}{\sqrt{x + 1} - \sqrt{x}}
\frac{\sqrt{1 + x} + \sqrt{x}}{x + \left(1 - x\right)}
Results
Initial program 4.7
Applied egg-rr0.3
Simplified0.3
[Start]0.3 | \[ \frac{1}{x + \left(1 - x\right)} \cdot \sqrt{x + 1} + \frac{1}{x + \left(1 - x\right)} \cdot \sqrt{x}
\] |
---|---|
distribute-lft-in [<=]0.3 | \[ \color{blue}{\frac{1}{x + \left(1 - x\right)} \cdot \left(\sqrt{x + 1} + \sqrt{x}\right)}
\] |
*-commutative [<=]0.3 | \[ \color{blue}{\left(\sqrt{x + 1} + \sqrt{x}\right) \cdot \frac{1}{x + \left(1 - x\right)}}
\] |
/-rgt-identity [<=]0.3 | \[ \color{blue}{\frac{\sqrt{x + 1} + \sqrt{x}}{1}} \cdot \frac{1}{x + \left(1 - x\right)}
\] |
associate-/r/ [<=]0.3 | \[ \color{blue}{\frac{\sqrt{x + 1} + \sqrt{x}}{\frac{1}{\frac{1}{x + \left(1 - x\right)}}}}
\] |
+-commutative [=>]0.3 | \[ \frac{\sqrt{\color{blue}{1 + x}} + \sqrt{x}}{\frac{1}{\frac{1}{x + \left(1 - x\right)}}}
\] |
remove-double-div [=>]0.3 | \[ \frac{\sqrt{1 + x} + \sqrt{x}}{\color{blue}{x + \left(1 - x\right)}}
\] |
Final simplification0.3
Alternative 1 | |
---|---|
Error | 0.3 |
Cost | 13376 |
Alternative 2 | |
---|---|
Error | 4.7 |
Cost | 13248 |
Alternative 3 | |
---|---|
Error | 41.1 |
Cost | 7232 |
Alternative 4 | |
---|---|
Error | 42.5 |
Cost | 6976 |
Alternative 5 | |
---|---|
Error | 50.8 |
Cost | 64 |
herbie shell --seed 1
(FPCore (x)
:name "1 / (sqrt(x+1) - sqrt(x))"
:precision binary64
:pre (and (<= 1e-6 x) (<= x 1000000.0))
(/ 1.0 (- (sqrt (+ x 1.0)) (sqrt x))))