| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 1224 |
(FPCore (x y) :precision binary64 (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))
(FPCore (x y)
:precision binary64
(if (<= y -2.05e+20)
(- x (/ -1.0 y))
(if (<= y 12500.0)
(+ 1.0 (* y (/ (+ x -1.0) (+ y 1.0))))
(+
x
(+
(* (/ (+ x -1.0) y) (+ -1.0 (/ -1.0 (* y y))))
(/ (+ x -1.0) (* y y)))))))double code(double x, double y) {
return 1.0 - (((1.0 - x) * y) / (y + 1.0));
}
double code(double x, double y) {
double tmp;
if (y <= -2.05e+20) {
tmp = x - (-1.0 / y);
} else if (y <= 12500.0) {
tmp = 1.0 + (y * ((x + -1.0) / (y + 1.0)));
} else {
tmp = x + ((((x + -1.0) / y) * (-1.0 + (-1.0 / (y * y)))) + ((x + -1.0) / (y * y)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 - (((1.0d0 - x) * y) / (y + 1.0d0))
end function
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-2.05d+20)) then
tmp = x - ((-1.0d0) / y)
else if (y <= 12500.0d0) then
tmp = 1.0d0 + (y * ((x + (-1.0d0)) / (y + 1.0d0)))
else
tmp = x + ((((x + (-1.0d0)) / y) * ((-1.0d0) + ((-1.0d0) / (y * y)))) + ((x + (-1.0d0)) / (y * y)))
end if
code = tmp
end function
public static double code(double x, double y) {
return 1.0 - (((1.0 - x) * y) / (y + 1.0));
}
public static double code(double x, double y) {
double tmp;
if (y <= -2.05e+20) {
tmp = x - (-1.0 / y);
} else if (y <= 12500.0) {
tmp = 1.0 + (y * ((x + -1.0) / (y + 1.0)));
} else {
tmp = x + ((((x + -1.0) / y) * (-1.0 + (-1.0 / (y * y)))) + ((x + -1.0) / (y * y)));
}
return tmp;
}
def code(x, y): return 1.0 - (((1.0 - x) * y) / (y + 1.0))
def code(x, y): tmp = 0 if y <= -2.05e+20: tmp = x - (-1.0 / y) elif y <= 12500.0: tmp = 1.0 + (y * ((x + -1.0) / (y + 1.0))) else: tmp = x + ((((x + -1.0) / y) * (-1.0 + (-1.0 / (y * y)))) + ((x + -1.0) / (y * y))) return tmp
function code(x, y) return Float64(1.0 - Float64(Float64(Float64(1.0 - x) * y) / Float64(y + 1.0))) end
function code(x, y) tmp = 0.0 if (y <= -2.05e+20) tmp = Float64(x - Float64(-1.0 / y)); elseif (y <= 12500.0) tmp = Float64(1.0 + Float64(y * Float64(Float64(x + -1.0) / Float64(y + 1.0)))); else tmp = Float64(x + Float64(Float64(Float64(Float64(x + -1.0) / y) * Float64(-1.0 + Float64(-1.0 / Float64(y * y)))) + Float64(Float64(x + -1.0) / Float64(y * y)))); end return tmp end
function tmp = code(x, y) tmp = 1.0 - (((1.0 - x) * y) / (y + 1.0)); end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -2.05e+20) tmp = x - (-1.0 / y); elseif (y <= 12500.0) tmp = 1.0 + (y * ((x + -1.0) / (y + 1.0))); else tmp = x + ((((x + -1.0) / y) * (-1.0 + (-1.0 / (y * y)))) + ((x + -1.0) / (y * y))); end tmp_2 = tmp; end
code[x_, y_] := N[(1.0 - N[(N[(N[(1.0 - x), $MachinePrecision] * y), $MachinePrecision] / N[(y + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_] := If[LessEqual[y, -2.05e+20], N[(x - N[(-1.0 / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 12500.0], N[(1.0 + N[(y * N[(N[(x + -1.0), $MachinePrecision] / N[(y + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(N[(N[(x + -1.0), $MachinePrecision] / y), $MachinePrecision] * N[(-1.0 + N[(-1.0 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x + -1.0), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
1 - \frac{\left(1 - x\right) \cdot y}{y + 1}
\begin{array}{l}
\mathbf{if}\;y \leq -2.05 \cdot 10^{+20}:\\
\;\;\;\;x - \frac{-1}{y}\\
\mathbf{elif}\;y \leq 12500:\\
\;\;\;\;1 + y \cdot \frac{x + -1}{y + 1}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{x + -1}{y} \cdot \left(-1 + \frac{-1}{y \cdot y}\right) + \frac{x + -1}{y \cdot y}\right)\\
\end{array}
Results
| Original | 63.6% |
|---|---|
| Target | 99.6% |
| Herbie | 99.5% |
if y < -2.05e20Initial program 23.9%
Simplified53.5%
[Start]23.9 | \[ 1 - \frac{\left(1 - x\right) \cdot y}{y + 1}
\] |
|---|---|
sub-neg [=>]23.9 | \[ \color{blue}{1 + \left(-\frac{\left(1 - x\right) \cdot y}{y + 1}\right)}
\] |
distribute-neg-frac [=>]23.9 | \[ 1 + \color{blue}{\frac{-\left(1 - x\right) \cdot y}{y + 1}}
\] |
neg-mul-1 [=>]23.9 | \[ 1 + \frac{\color{blue}{-1 \cdot \left(\left(1 - x\right) \cdot y\right)}}{y + 1}
\] |
associate-*l/ [<=]23.8 | \[ 1 + \color{blue}{\frac{-1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
metadata-eval [<=]23.8 | \[ 1 + \frac{\color{blue}{1 \cdot -1}}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
associate-*l/ [<=]23.8 | \[ 1 + \color{blue}{\left(\frac{1}{y + 1} \cdot -1\right)} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
associate-/r/ [<=]23.8 | \[ 1 + \color{blue}{\frac{1}{\frac{y + 1}{-1}}} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
metadata-eval [<=]23.8 | \[ 1 + \frac{\color{blue}{--1}}{\frac{y + 1}{-1}} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
distribute-neg-frac [<=]23.8 | \[ 1 + \color{blue}{\left(-\frac{-1}{\frac{y + 1}{-1}}\right)} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
cancel-sign-sub-inv [<=]23.8 | \[ \color{blue}{1 - \frac{-1}{\frac{y + 1}{-1}} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
associate-/r/ [<=]23.8 | \[ 1 - \color{blue}{\frac{-1}{\frac{\frac{y + 1}{-1}}{\left(1 - x\right) \cdot y}}}
\] |
associate-/r* [<=]23.8 | \[ 1 - \frac{-1}{\color{blue}{\frac{y + 1}{-1 \cdot \left(\left(1 - x\right) \cdot y\right)}}}
\] |
neg-mul-1 [<=]23.8 | \[ 1 - \frac{-1}{\frac{y + 1}{\color{blue}{-\left(1 - x\right) \cdot y}}}
\] |
associate-/r/ [=>]23.8 | \[ 1 - \color{blue}{\frac{-1}{y + 1} \cdot \left(-\left(1 - x\right) \cdot y\right)}
\] |
distribute-rgt-neg-in [<=]23.8 | \[ 1 - \color{blue}{\left(-\frac{-1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)\right)}
\] |
associate-/r/ [<=]23.8 | \[ 1 - \left(-\color{blue}{\frac{-1}{\frac{y + 1}{\left(1 - x\right) \cdot y}}}\right)
\] |
distribute-neg-frac [=>]23.8 | \[ 1 - \color{blue}{\frac{--1}{\frac{y + 1}{\left(1 - x\right) \cdot y}}}
\] |
metadata-eval [=>]23.8 | \[ 1 - \frac{\color{blue}{1}}{\frac{y + 1}{\left(1 - x\right) \cdot y}}
\] |
associate-/r/ [=>]23.8 | \[ 1 - \color{blue}{\frac{1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
Taylor expanded in y around inf 100.0%
Simplified100.0%
[Start]100.0 | \[ \left(\frac{1}{y} + x\right) - \frac{x}{y}
\] |
|---|---|
+-commutative [=>]100.0 | \[ \color{blue}{\left(x + \frac{1}{y}\right)} - \frac{x}{y}
\] |
associate--l+ [=>]100.0 | \[ \color{blue}{x + \left(\frac{1}{y} - \frac{x}{y}\right)}
\] |
div-sub [<=]100.0 | \[ x + \color{blue}{\frac{1 - x}{y}}
\] |
sub-neg [=>]100.0 | \[ x + \frac{\color{blue}{1 + \left(-x\right)}}{y}
\] |
+-commutative [=>]100.0 | \[ x + \frac{\color{blue}{\left(-x\right) + 1}}{y}
\] |
neg-sub0 [=>]100.0 | \[ x + \frac{\color{blue}{\left(0 - x\right)} + 1}{y}
\] |
associate-+l- [=>]100.0 | \[ x + \frac{\color{blue}{0 - \left(x - 1\right)}}{y}
\] |
neg-sub0 [<=]100.0 | \[ x + \frac{\color{blue}{-\left(x - 1\right)}}{y}
\] |
mul-1-neg [<=]100.0 | \[ x + \frac{\color{blue}{-1 \cdot \left(x - 1\right)}}{y}
\] |
associate-*r/ [<=]100.0 | \[ x + \color{blue}{-1 \cdot \frac{x - 1}{y}}
\] |
mul-1-neg [=>]100.0 | \[ x + \color{blue}{\left(-\frac{x - 1}{y}\right)}
\] |
unsub-neg [=>]100.0 | \[ \color{blue}{x - \frac{x - 1}{y}}
\] |
sub-neg [=>]100.0 | \[ x - \frac{\color{blue}{x + \left(-1\right)}}{y}
\] |
metadata-eval [=>]100.0 | \[ x - \frac{x + \color{blue}{-1}}{y}
\] |
Taylor expanded in x around 0 100.0%
if -2.05e20 < y < 12500Initial program 100.0%
Simplified100.0%
[Start]100.0 | \[ 1 - \frac{\left(1 - x\right) \cdot y}{y + 1}
\] |
|---|---|
sub-neg [=>]100.0 | \[ \color{blue}{1 + \left(-\frac{\left(1 - x\right) \cdot y}{y + 1}\right)}
\] |
distribute-neg-frac [=>]100.0 | \[ 1 + \color{blue}{\frac{-\left(1 - x\right) \cdot y}{y + 1}}
\] |
neg-mul-1 [=>]100.0 | \[ 1 + \frac{\color{blue}{-1 \cdot \left(\left(1 - x\right) \cdot y\right)}}{y + 1}
\] |
associate-*l/ [<=]100.0 | \[ 1 + \color{blue}{\frac{-1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
metadata-eval [<=]100.0 | \[ 1 + \frac{\color{blue}{1 \cdot -1}}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
associate-*l/ [<=]100.0 | \[ 1 + \color{blue}{\left(\frac{1}{y + 1} \cdot -1\right)} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
associate-/r/ [<=]100.0 | \[ 1 + \color{blue}{\frac{1}{\frac{y + 1}{-1}}} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
metadata-eval [<=]100.0 | \[ 1 + \frac{\color{blue}{--1}}{\frac{y + 1}{-1}} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
distribute-neg-frac [<=]100.0 | \[ 1 + \color{blue}{\left(-\frac{-1}{\frac{y + 1}{-1}}\right)} \cdot \left(\left(1 - x\right) \cdot y\right)
\] |
cancel-sign-sub-inv [<=]100.0 | \[ \color{blue}{1 - \frac{-1}{\frac{y + 1}{-1}} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
associate-/r/ [<=]99.9 | \[ 1 - \color{blue}{\frac{-1}{\frac{\frac{y + 1}{-1}}{\left(1 - x\right) \cdot y}}}
\] |
associate-/r* [<=]99.9 | \[ 1 - \frac{-1}{\color{blue}{\frac{y + 1}{-1 \cdot \left(\left(1 - x\right) \cdot y\right)}}}
\] |
neg-mul-1 [<=]99.9 | \[ 1 - \frac{-1}{\frac{y + 1}{\color{blue}{-\left(1 - x\right) \cdot y}}}
\] |
associate-/r/ [=>]100.0 | \[ 1 - \color{blue}{\frac{-1}{y + 1} \cdot \left(-\left(1 - x\right) \cdot y\right)}
\] |
distribute-rgt-neg-in [<=]100.0 | \[ 1 - \color{blue}{\left(-\frac{-1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)\right)}
\] |
associate-/r/ [<=]99.9 | \[ 1 - \left(-\color{blue}{\frac{-1}{\frac{y + 1}{\left(1 - x\right) \cdot y}}}\right)
\] |
distribute-neg-frac [=>]99.9 | \[ 1 - \color{blue}{\frac{--1}{\frac{y + 1}{\left(1 - x\right) \cdot y}}}
\] |
metadata-eval [=>]99.9 | \[ 1 - \frac{\color{blue}{1}}{\frac{y + 1}{\left(1 - x\right) \cdot y}}
\] |
associate-/r/ [=>]100.0 | \[ 1 - \color{blue}{\frac{1}{y + 1} \cdot \left(\left(1 - x\right) \cdot y\right)}
\] |
if 12500 < y Initial program 30.8%
Simplified55.5%
[Start]30.8 | \[ 1 - \frac{\left(1 - x\right) \cdot y}{y + 1}
\] |
|---|---|
sub-neg [=>]30.8 | \[ \color{blue}{1 + \left(-\frac{\left(1 - x\right) \cdot y}{y + 1}\right)}
\] |
+-commutative [=>]30.8 | \[ \color{blue}{\left(-\frac{\left(1 - x\right) \cdot y}{y + 1}\right) + 1}
\] |
associate-/l* [=>]55.3 | \[ \left(-\color{blue}{\frac{1 - x}{\frac{y + 1}{y}}}\right) + 1
\] |
distribute-neg-frac [=>]55.3 | \[ \color{blue}{\frac{-\left(1 - x\right)}{\frac{y + 1}{y}}} + 1
\] |
*-lft-identity [<=]55.3 | \[ \frac{\color{blue}{1 \cdot \left(-\left(1 - x\right)\right)}}{\frac{y + 1}{y}} + 1
\] |
associate-*l/ [<=]55.3 | \[ \color{blue}{\frac{1}{\frac{y + 1}{y}} \cdot \left(-\left(1 - x\right)\right)} + 1
\] |
fma-def [=>]55.3 | \[ \color{blue}{\mathsf{fma}\left(\frac{1}{\frac{y + 1}{y}}, -\left(1 - x\right), 1\right)}
\] |
associate-/l* [<=]55.5 | \[ \mathsf{fma}\left(\color{blue}{\frac{1 \cdot y}{y + 1}}, -\left(1 - x\right), 1\right)
\] |
*-lft-identity [=>]55.5 | \[ \mathsf{fma}\left(\frac{\color{blue}{y}}{y + 1}, -\left(1 - x\right), 1\right)
\] |
+-commutative [=>]55.5 | \[ \mathsf{fma}\left(\frac{y}{\color{blue}{1 + y}}, -\left(1 - x\right), 1\right)
\] |
neg-sub0 [=>]55.5 | \[ \mathsf{fma}\left(\frac{y}{1 + y}, \color{blue}{0 - \left(1 - x\right)}, 1\right)
\] |
associate--r- [=>]55.5 | \[ \mathsf{fma}\left(\frac{y}{1 + y}, \color{blue}{\left(0 - 1\right) + x}, 1\right)
\] |
metadata-eval [=>]55.5 | \[ \mathsf{fma}\left(\frac{y}{1 + y}, \color{blue}{-1} + x, 1\right)
\] |
+-commutative [=>]55.5 | \[ \mathsf{fma}\left(\frac{y}{1 + y}, \color{blue}{x + -1}, 1\right)
\] |
Taylor expanded in y around inf 99.7%
Simplified99.7%
[Start]99.7 | \[ \left(\frac{x}{{y}^{2}} + \left(-1 \cdot \frac{x - 1}{y} + \left(-1 \cdot \frac{x - 1}{{y}^{3}} + x\right)\right)\right) - \frac{1}{{y}^{2}}
\] |
|---|---|
+-commutative [=>]99.7 | \[ \color{blue}{\left(\left(-1 \cdot \frac{x - 1}{y} + \left(-1 \cdot \frac{x - 1}{{y}^{3}} + x\right)\right) + \frac{x}{{y}^{2}}\right)} - \frac{1}{{y}^{2}}
\] |
associate--l+ [=>]99.7 | \[ \color{blue}{\left(-1 \cdot \frac{x - 1}{y} + \left(-1 \cdot \frac{x - 1}{{y}^{3}} + x\right)\right) + \left(\frac{x}{{y}^{2}} - \frac{1}{{y}^{2}}\right)}
\] |
div-sub [<=]99.7 | \[ \left(-1 \cdot \frac{x - 1}{y} + \left(-1 \cdot \frac{x - 1}{{y}^{3}} + x\right)\right) + \color{blue}{\frac{x - 1}{{y}^{2}}}
\] |
associate-+r+ [=>]99.7 | \[ \color{blue}{\left(\left(-1 \cdot \frac{x - 1}{y} + -1 \cdot \frac{x - 1}{{y}^{3}}\right) + x\right)} + \frac{x - 1}{{y}^{2}}
\] |
+-commutative [=>]99.7 | \[ \color{blue}{\left(x + \left(-1 \cdot \frac{x - 1}{y} + -1 \cdot \frac{x - 1}{{y}^{3}}\right)\right)} + \frac{x - 1}{{y}^{2}}
\] |
associate-+l+ [=>]99.7 | \[ \color{blue}{x + \left(\left(-1 \cdot \frac{x - 1}{y} + -1 \cdot \frac{x - 1}{{y}^{3}}\right) + \frac{x - 1}{{y}^{2}}\right)}
\] |
Final simplification99.9%
| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 1224 |
| Alternative 2 | |
|---|---|
| Accuracy | 99.4% |
| Cost | 968 |
| Alternative 3 | |
|---|---|
| Accuracy | 98.3% |
| Cost | 841 |
| Alternative 4 | |
|---|---|
| Accuracy | 98.2% |
| Cost | 840 |
| Alternative 5 | |
|---|---|
| Accuracy | 98.0% |
| Cost | 713 |
| Alternative 6 | |
|---|---|
| Accuracy | 98.3% |
| Cost | 713 |
| Alternative 7 | |
|---|---|
| Accuracy | 74.1% |
| Cost | 588 |
| Alternative 8 | |
|---|---|
| Accuracy | 74.2% |
| Cost | 588 |
| Alternative 9 | |
|---|---|
| Accuracy | 97.8% |
| Cost | 585 |
| Alternative 10 | |
|---|---|
| Accuracy | 85.5% |
| Cost | 584 |
| Alternative 11 | |
|---|---|
| Accuracy | 74.1% |
| Cost | 328 |
| Alternative 12 | |
|---|---|
| Accuracy | 38.5% |
| Cost | 64 |
herbie shell --seed 2023157
(FPCore (x y)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, D"
:precision binary64
:herbie-target
(if (< y -3693.8482788297247) (- (/ 1.0 y) (- (/ x y) x)) (if (< y 6799310503.41891) (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))) (- (/ 1.0 y) (- (/ x y) x))))
(- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))