| Alternative 1 | |
|---|---|
| Error | 2.6 |
| Cost | 52800 |
\[\left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right) + \left(\left(\sqrt{1 + y} - \sqrt{y}\right) + \frac{1}{\sqrt{x} + \sqrt{1 + x}}\right)
\]
(FPCore (x y z t) :precision binary64 (+ (+ (+ (- (sqrt (+ x 1.0)) (sqrt x)) (- (sqrt (+ y 1.0)) (sqrt y))) (- (sqrt (+ z 1.0)) (sqrt z))) (- (sqrt (+ t 1.0)) (sqrt t))))
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (sqrt y) (sqrt (+ 1.0 y)))))
(+
(/ (+ 1.0 (* t_1 (/ 1.0 (+ (sqrt x) (sqrt (+ 1.0 x)))))) t_1)
(+ (- (sqrt (+ 1.0 z)) (sqrt z)) (- (sqrt (+ 1.0 t)) (sqrt t))))))double code(double x, double y, double z, double t) {
return (((sqrt((x + 1.0)) - sqrt(x)) + (sqrt((y + 1.0)) - sqrt(y))) + (sqrt((z + 1.0)) - sqrt(z))) + (sqrt((t + 1.0)) - sqrt(t));
}
double code(double x, double y, double z, double t) {
double t_1 = sqrt(y) + sqrt((1.0 + y));
return ((1.0 + (t_1 * (1.0 / (sqrt(x) + sqrt((1.0 + x)))))) / t_1) + ((sqrt((1.0 + z)) - sqrt(z)) + (sqrt((1.0 + t)) - sqrt(t)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((sqrt((x + 1.0d0)) - sqrt(x)) + (sqrt((y + 1.0d0)) - sqrt(y))) + (sqrt((z + 1.0d0)) - sqrt(z))) + (sqrt((t + 1.0d0)) - sqrt(t))
end function
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
t_1 = sqrt(y) + sqrt((1.0d0 + y))
code = ((1.0d0 + (t_1 * (1.0d0 / (sqrt(x) + sqrt((1.0d0 + x)))))) / t_1) + ((sqrt((1.0d0 + z)) - sqrt(z)) + (sqrt((1.0d0 + t)) - sqrt(t)))
end function
public static double code(double x, double y, double z, double t) {
return (((Math.sqrt((x + 1.0)) - Math.sqrt(x)) + (Math.sqrt((y + 1.0)) - Math.sqrt(y))) + (Math.sqrt((z + 1.0)) - Math.sqrt(z))) + (Math.sqrt((t + 1.0)) - Math.sqrt(t));
}
public static double code(double x, double y, double z, double t) {
double t_1 = Math.sqrt(y) + Math.sqrt((1.0 + y));
return ((1.0 + (t_1 * (1.0 / (Math.sqrt(x) + Math.sqrt((1.0 + x)))))) / t_1) + ((Math.sqrt((1.0 + z)) - Math.sqrt(z)) + (Math.sqrt((1.0 + t)) - Math.sqrt(t)));
}
def code(x, y, z, t): return (((math.sqrt((x + 1.0)) - math.sqrt(x)) + (math.sqrt((y + 1.0)) - math.sqrt(y))) + (math.sqrt((z + 1.0)) - math.sqrt(z))) + (math.sqrt((t + 1.0)) - math.sqrt(t))
def code(x, y, z, t): t_1 = math.sqrt(y) + math.sqrt((1.0 + y)) return ((1.0 + (t_1 * (1.0 / (math.sqrt(x) + math.sqrt((1.0 + x)))))) / t_1) + ((math.sqrt((1.0 + z)) - math.sqrt(z)) + (math.sqrt((1.0 + t)) - math.sqrt(t)))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(sqrt(Float64(x + 1.0)) - sqrt(x)) + Float64(sqrt(Float64(y + 1.0)) - sqrt(y))) + Float64(sqrt(Float64(z + 1.0)) - sqrt(z))) + Float64(sqrt(Float64(t + 1.0)) - sqrt(t))) end
function code(x, y, z, t) t_1 = Float64(sqrt(y) + sqrt(Float64(1.0 + y))) return Float64(Float64(Float64(1.0 + Float64(t_1 * Float64(1.0 / Float64(sqrt(x) + sqrt(Float64(1.0 + x)))))) / t_1) + Float64(Float64(sqrt(Float64(1.0 + z)) - sqrt(z)) + Float64(sqrt(Float64(1.0 + t)) - sqrt(t)))) end
function tmp = code(x, y, z, t) tmp = (((sqrt((x + 1.0)) - sqrt(x)) + (sqrt((y + 1.0)) - sqrt(y))) + (sqrt((z + 1.0)) - sqrt(z))) + (sqrt((t + 1.0)) - sqrt(t)); end
function tmp = code(x, y, z, t) t_1 = sqrt(y) + sqrt((1.0 + y)); tmp = ((1.0 + (t_1 * (1.0 / (sqrt(x) + sqrt((1.0 + x)))))) / t_1) + ((sqrt((1.0 + z)) - sqrt(z)) + (sqrt((1.0 + t)) - sqrt(t))); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[x], $MachinePrecision]), $MachinePrecision] + N[(N[Sqrt[N[(y + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Sqrt[N[(z + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Sqrt[N[(t + 1.0), $MachinePrecision]], $MachinePrecision] - N[Sqrt[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[Sqrt[y], $MachinePrecision] + N[Sqrt[N[(1.0 + y), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(1.0 + N[(t$95$1 * N[(1.0 / N[(N[Sqrt[x], $MachinePrecision] + N[Sqrt[N[(1.0 + x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(N[(N[Sqrt[N[(1.0 + z), $MachinePrecision]], $MachinePrecision] - N[Sqrt[z], $MachinePrecision]), $MachinePrecision] + N[(N[Sqrt[N[(1.0 + t), $MachinePrecision]], $MachinePrecision] - N[Sqrt[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\left(\left(\left(\sqrt{x + 1} - \sqrt{x}\right) + \left(\sqrt{y + 1} - \sqrt{y}\right)\right) + \left(\sqrt{z + 1} - \sqrt{z}\right)\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)
\begin{array}{l}
t_1 := \sqrt{y} + \sqrt{1 + y}\\
\frac{1 + t_1 \cdot \frac{1}{\sqrt{x} + \sqrt{1 + x}}}{t_1} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\end{array}
Results
| Original | 5.2 |
|---|---|
| Target | 0.4 |
| Herbie | 1.3 |
Initial program 5.2
Simplified5.2
[Start]5.2 | \[ \left(\left(\left(\sqrt{x + 1} - \sqrt{x}\right) + \left(\sqrt{y + 1} - \sqrt{y}\right)\right) + \left(\sqrt{z + 1} - \sqrt{z}\right)\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)
\] |
|---|---|
associate-+l+ [=>]5.2 | \[ \color{blue}{\left(\left(\sqrt{x + 1} - \sqrt{x}\right) + \left(\sqrt{y + 1} - \sqrt{y}\right)\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)}
\] |
associate-+l- [=>]5.2 | \[ \color{blue}{\left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{y + 1} - \sqrt{y}\right)\right)\right)} + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
associate--r- [=>]5.9 | \[ \left(\sqrt{x + 1} - \color{blue}{\left(\left(\sqrt{x} - \sqrt{y + 1}\right) + \sqrt{y}\right)}\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
remove-double-neg [<=]5.9 | \[ \left(\sqrt{x + 1} - \left(\left(\sqrt{x} - \sqrt{y + 1}\right) + \color{blue}{\left(-\left(-\sqrt{y}\right)\right)}\right)\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
associate-+l- [=>]5.2 | \[ \left(\sqrt{x + 1} - \color{blue}{\left(\sqrt{x} - \left(\sqrt{y + 1} - \left(-\left(-\sqrt{y}\right)\right)\right)\right)}\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
+-commutative [=>]5.2 | \[ \left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{\color{blue}{1 + y}} - \left(-\left(-\sqrt{y}\right)\right)\right)\right)\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
remove-double-neg [=>]5.2 | \[ \left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{1 + y} - \color{blue}{\sqrt{y}}\right)\right)\right) + \left(\left(\sqrt{z + 1} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
sub-neg [=>]5.2 | \[ \left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{1 + y} - \sqrt{y}\right)\right)\right) + \left(\color{blue}{\left(\sqrt{z + 1} + \left(-\sqrt{z}\right)\right)} + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
sub-neg [<=]5.2 | \[ \left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{1 + y} - \sqrt{y}\right)\right)\right) + \left(\color{blue}{\left(\sqrt{z + 1} - \sqrt{z}\right)} + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
+-commutative [=>]5.2 | \[ \left(\sqrt{x + 1} - \left(\sqrt{x} - \left(\sqrt{1 + y} - \sqrt{y}\right)\right)\right) + \left(\left(\sqrt{\color{blue}{1 + z}} - \sqrt{z}\right) + \left(\sqrt{t + 1} - \sqrt{t}\right)\right)
\] |
Applied egg-rr5.2
Applied egg-rr1.3
Applied egg-rr1.3
Simplified1.3
[Start]1.3 | \[ \frac{\left(\left(-\left(\sqrt{1 + x} + \sqrt{x}\right)\right) - \left(\sqrt{1 + y} + \sqrt{y}\right)\right) \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
|---|---|
*-commutative [=>]1.3 | \[ \frac{\color{blue}{\frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)} \cdot \left(\left(-\left(\sqrt{1 + x} + \sqrt{x}\right)\right) - \left(\sqrt{1 + y} + \sqrt{y}\right)\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
sub-neg [=>]1.3 | \[ \frac{\frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)} \cdot \color{blue}{\left(\left(-\left(\sqrt{1 + x} + \sqrt{x}\right)\right) + \left(-\left(\sqrt{1 + y} + \sqrt{y}\right)\right)\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
distribute-rgt-in [=>]1.3 | \[ \frac{\color{blue}{\left(-\left(\sqrt{1 + x} + \sqrt{x}\right)\right) \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)} + \left(-\left(\sqrt{1 + y} + \sqrt{y}\right)\right) \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
rgt-mult-inverse [=>]1.3 | \[ \frac{\color{blue}{1} + \left(-\left(\sqrt{1 + y} + \sqrt{y}\right)\right) \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
+-commutative [=>]1.3 | \[ \frac{1 + \left(-\color{blue}{\left(\sqrt{y} + \sqrt{1 + y}\right)}\right) \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
distribute-neg-in [=>]1.3 | \[ \frac{1 + \color{blue}{\left(\left(-\sqrt{y}\right) + \left(-\sqrt{1 + y}\right)\right)} \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
sub-neg [<=]1.3 | \[ \frac{1 + \color{blue}{\left(\left(-\sqrt{y}\right) - \sqrt{1 + y}\right)} \cdot \frac{1}{-\left(\sqrt{1 + x} + \sqrt{x}\right)}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
neg-mul-1 [=>]1.3 | \[ \frac{1 + \left(\left(-\sqrt{y}\right) - \sqrt{1 + y}\right) \cdot \frac{1}{\color{blue}{-1 \cdot \left(\sqrt{1 + x} + \sqrt{x}\right)}}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
associate-/r* [=>]1.3 | \[ \frac{1 + \left(\left(-\sqrt{y}\right) - \sqrt{1 + y}\right) \cdot \color{blue}{\frac{\frac{1}{-1}}{\sqrt{1 + x} + \sqrt{x}}}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
metadata-eval [=>]1.3 | \[ \frac{1 + \left(\left(-\sqrt{y}\right) - \sqrt{1 + y}\right) \cdot \frac{\color{blue}{-1}}{\sqrt{1 + x} + \sqrt{x}}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
+-commutative [=>]1.3 | \[ \frac{1 + \left(\left(-\sqrt{y}\right) - \sqrt{1 + y}\right) \cdot \frac{-1}{\color{blue}{\sqrt{x} + \sqrt{1 + x}}}}{\sqrt{1 + y} + \sqrt{y}} + \left(\left(\sqrt{1 + z} - \sqrt{z}\right) + \left(\sqrt{1 + t} - \sqrt{t}\right)\right)
\] |
Final simplification1.3
| Alternative 1 | |
|---|---|
| Error | 2.6 |
| Cost | 52800 |
| Alternative 2 | |
|---|---|
| Error | 4.4 |
| Cost | 40008 |
| Alternative 3 | |
|---|---|
| Error | 4.4 |
| Cost | 39876 |
| Alternative 4 | |
|---|---|
| Error | 4.5 |
| Cost | 39872 |
| Alternative 5 | |
|---|---|
| Error | 5.9 |
| Cost | 39752 |
| Alternative 6 | |
|---|---|
| Error | 5.4 |
| Cost | 39748 |
| Alternative 7 | |
|---|---|
| Error | 6.2 |
| Cost | 39620 |
| Alternative 8 | |
|---|---|
| Error | 6.0 |
| Cost | 26696 |
| Alternative 9 | |
|---|---|
| Error | 9.0 |
| Cost | 26564 |
| Alternative 10 | |
|---|---|
| Error | 11.5 |
| Cost | 13508 |
| Alternative 11 | |
|---|---|
| Error | 25.3 |
| Cost | 13380 |
| Alternative 12 | |
|---|---|
| Error | 22.5 |
| Cost | 13380 |
| Alternative 13 | |
|---|---|
| Error | 21.8 |
| Cost | 13376 |
| Alternative 14 | |
|---|---|
| Error | 41.3 |
| Cost | 13248 |
| Alternative 15 | |
|---|---|
| Error | 41.8 |
| Cost | 64 |
herbie shell --seed 2023034
(FPCore (x y z t)
:name "Main:z from "
:precision binary64
:herbie-target
(+ (+ (+ (/ 1.0 (+ (sqrt (+ x 1.0)) (sqrt x))) (/ 1.0 (+ (sqrt (+ y 1.0)) (sqrt y)))) (/ 1.0 (+ (sqrt (+ z 1.0)) (sqrt z)))) (- (sqrt (+ t 1.0)) (sqrt t)))
(+ (+ (+ (- (sqrt (+ x 1.0)) (sqrt x)) (- (sqrt (+ y 1.0)) (sqrt y))) (- (sqrt (+ z 1.0)) (sqrt z))) (- (sqrt (+ t 1.0)) (sqrt t))))