Average Error: 10.9 → 0.1
Time: 47.6s
Precision: 64
Internal Precision: 1152
\[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
↓
\[\begin{array}{l}
\mathbf{if}\;x \le -1.7387217572195604 \cdot 10^{-16}:\\
\;\;\;\;\left(\frac{2}{{x}^{7}} + \frac{2}{{x}^{5}}\right) + \frac{\frac{\frac{2}{x}}{x}}{x}\\
\mathbf{if}\;x \le 2500116.352981572:\\
\;\;\;\;-\left(2 \cdot \frac{1}{x} + \left(2 \cdot {x}^{3} + 2 \cdot x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{2}{{x}^{7}} + \frac{2}{{x}^{5}}\right) + \frac{\frac{2}{x}}{x \cdot x}\\
\end{array}\]
Target
| Original | 10.9 |
|---|
| Target | 1.6 |
|---|
| Herbie | 0.1 |
|---|
\[\frac{2}{x \cdot \left(x \cdot x - 1\right)}\]
Derivation
- Split input into 3 regimes
if x < -1.7387217572195604e-16
Initial program 21.9
\[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
Taylor expanded around inf 0.5
\[\leadsto \color{blue}{2 \cdot \frac{1}{{x}^{3}} + \left(2 \cdot \frac{1}{{x}^{5}} + 2 \cdot \frac{1}{{x}^{7}}\right)}\]
Applied simplify0.1
\[\leadsto \color{blue}{\left(\frac{2}{{x}^{7}} + \frac{2}{{x}^{5}}\right) + \frac{\frac{2}{x}}{x \cdot x}}\]
- Using strategy
rm Applied associate-/r*0.1
\[\leadsto \left(\frac{2}{{x}^{7}} + \frac{2}{{x}^{5}}\right) + \color{blue}{\frac{\frac{\frac{2}{x}}{x}}{x}}\]
if -1.7387217572195604e-16 < x < 2500116.352981572
Initial program 0.7
\[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
Taylor expanded around 0 0.0
\[\leadsto \color{blue}{-\left(2 \cdot \frac{1}{x} + \left(2 \cdot {x}^{3} + 2 \cdot x\right)\right)}\]
if 2500116.352981572 < x
Initial program 19.2
\[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
Taylor expanded around inf 0.5
\[\leadsto \color{blue}{2 \cdot \frac{1}{{x}^{3}} + \left(2 \cdot \frac{1}{{x}^{5}} + 2 \cdot \frac{1}{{x}^{7}}\right)}\]
Applied simplify0.1
\[\leadsto \color{blue}{\left(\frac{2}{{x}^{7}} + \frac{2}{{x}^{5}}\right) + \frac{\frac{2}{x}}{x \cdot x}}\]
- Recombined 3 regimes into one program.
- Removed slow
pow expressions.
Runtime
herbie shell --seed '#(1062803647 245428163 493620569 3595423923 1908391097 2390014376)'
(FPCore (x)
:name "3frac (problem 3.3.3)"
:herbie-target
(/ 2 (* x (- (* x x) 1)))
(+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1))))