Average Error: 9.2 → 0.1
Time: 1.5m
Precision: 64
Internal Precision: 128
\[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
\[\begin{array}{l} \mathbf{if}\;x \le -114.64678370673498 \lor \neg \left(x \le 104.33422414858215\right):\\ \;\;\;\;\frac{2}{{x}^{5}} + \left(\frac{\frac{\frac{2}{x}}{x}}{x} + \frac{2}{{x}^{7}}\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{-2}{x} - \left(\frac{-1}{-1 + x} - \frac{1}{1 + x}\right)\\ \end{array}\]

Error

Bits error versus x

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original9.2
Target0.3
Herbie0.1
\[\frac{2}{x \cdot \left(x \cdot x - 1\right)}\]

Derivation

  1. Split input into 2 regimes
  2. if x < -114.64678370673498 or 104.33422414858215 < x

    1. Initial program 18.7

      \[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    2. Using strategy rm
    3. Applied add-sqr-sqrt48.3

      \[\leadsto \left(\frac{1}{\color{blue}{\sqrt{x + 1} \cdot \sqrt{x + 1}}} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    4. Applied associate-/r*50.2

      \[\leadsto \left(\color{blue}{\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x + 1}}} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    5. Taylor expanded around inf 0.5

      \[\leadsto \color{blue}{2 \cdot \frac{1}{{x}^{7}} + \left(2 \cdot \frac{1}{{x}^{3}} + 2 \cdot \frac{1}{{x}^{5}}\right)}\]
    6. Simplified0.1

      \[\leadsto \color{blue}{\left(\frac{\frac{2}{x}}{x \cdot x} + \frac{2}{{x}^{7}}\right) + \frac{2}{{x}^{5}}}\]
    7. Using strategy rm
    8. Applied associate-/r*0.1

      \[\leadsto \left(\color{blue}{\frac{\frac{\frac{2}{x}}{x}}{x}} + \frac{2}{{x}^{7}}\right) + \frac{2}{{x}^{5}}\]

    if -114.64678370673498 < x < 104.33422414858215

    1. Initial program 0.0

      \[\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    2. Using strategy rm
    3. Applied add-sqr-sqrt0.2

      \[\leadsto \left(\frac{1}{\color{blue}{\sqrt{x + 1} \cdot \sqrt{x + 1}}} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    4. Applied associate-/r*0.2

      \[\leadsto \left(\color{blue}{\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x + 1}}} - \frac{2}{x}\right) + \frac{1}{x - 1}\]
    5. Using strategy rm
    6. Applied *-un-lft-identity0.2

      \[\leadsto \left(\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x + 1}} - \frac{2}{x}\right) + \color{blue}{1 \cdot \frac{1}{x - 1}}\]
    7. Applied *-un-lft-identity0.2

      \[\leadsto \color{blue}{1 \cdot \left(\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x + 1}} - \frac{2}{x}\right)} + 1 \cdot \frac{1}{x - 1}\]
    8. Applied distribute-lft-out0.2

      \[\leadsto \color{blue}{1 \cdot \left(\left(\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x + 1}} - \frac{2}{x}\right) + \frac{1}{x - 1}\right)}\]
    9. Simplified0.0

      \[\leadsto 1 \cdot \color{blue}{\left(\frac{-2}{x} - \left(\frac{-1}{x + -1} - \frac{1}{x + 1}\right)\right)}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.1

    \[\leadsto \begin{array}{l} \mathbf{if}\;x \le -114.64678370673498 \lor \neg \left(x \le 104.33422414858215\right):\\ \;\;\;\;\frac{2}{{x}^{5}} + \left(\frac{\frac{\frac{2}{x}}{x}}{x} + \frac{2}{{x}^{7}}\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{-2}{x} - \left(\frac{-1}{-1 + x} - \frac{1}{1 + x}\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2018360 +o rules:numerics
(FPCore (x)
  :name "3frac (problem 3.3.3)"

  :herbie-target
  (/ 2 (* x (- (* x x) 1)))

  (+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1))))

Details

Time bar (total: 1.4m)Debug log

start181.0ms

Algorithm
intervals

setup93.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 8.2b

localize34.0ms

Local error

Found 4 expressions with local error:

2.3b
(+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1)))
0.0b
(- (/ 1 (+ x 1)) (/ 2 x))
0.0b
(/ 1 (+ x 1))
0.0b
(/ 1 (- x 1))

rewrite34.0ms

Algorithm
rewrite-expression-head
Counts
4 → 144
Calls

4 calls. Slowest were:

24.0ms
(+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1)))
6.0ms
(- (/ 1 (+ x 1)) (/ 2 x))
1.0ms
(/ 1 (- x 1))

series72.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

44.0ms
(+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1)))
11.0ms
(- (/ 1 (+ x 1)) (/ 2 x))
9.0ms
(/ 1 (- x 1))
7.0ms
(/ 1 (+ x 1))

simplify14.6s

Counts
127 → 156
Calls

127 calls. Slowest were:

1.1s
(+ (* (- (* 1 x) (* (+ x 1) 2)) (- x 1)) (* (* (+ x 1) x) 1))
959.0ms
(fma (/ 1 (- (* x x) (* 1 1))) (- x 1) (- (* (/ 1 x) 2)))
947.0ms
(fma (/ 1 (- (* x x) (* 1 1))) (- x 1) (- (* (/ 2 x) 1)))

prune1.2s

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.1b

localize13.0ms

Local error

Found 4 expressions with local error:

6.5b
(+ (- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x)) (/ 1 (- x 1)))
0.1b
(/ 1 (sqrt (+ x 1)))
0.1b
(/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1)))
0.0b
(- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x))

rewrite104.0ms

Algorithm
rewrite-expression-head
Counts
4 → 523
Calls

4 calls. Slowest were:

66.0ms
(+ (- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x)) (/ 1 (- x 1)))
19.0ms
(- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x))
5.0ms
(/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1)))

series74.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

37.0ms
(+ (- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x)) (/ 1 (- x 1)))
18.0ms
(/ 1 (sqrt (+ x 1)))
11.0ms
(- (/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1))) (/ 2 x))
7.0ms
(/ (/ 1 (sqrt (+ x 1))) (sqrt (+ x 1)))

simplify43.2s

Counts
693 → 535
Calls

693 calls. Slowest were:

747.0ms
(+ (* (- (* (/ 1 (sqrt (+ x 1))) x) (* (sqrt (+ x 1)) 2)) (- x 1)) (* (* (sqrt (+ x 1)) x) 1))
618.0ms
(fma (/ (/ 1 (sqrt (+ (pow x 3) (pow 1 3)))) (* (cbrt (sqrt (+ x 1))) (cbrt (sqrt (+ x 1))))) (/ (sqrt (+ (* x x) (- (* 1 1) (* x 1)))) (cbrt (sqrt (+ x 1)))) (- (* (/ 2 x) 1)))
599.0ms
(fma (/ (/ 1 (sqrt (- (* x x) (* 1 1)))) (* (cbrt (sqrt (+ x 1))) (cbrt (sqrt (+ x 1))))) (/ (sqrt (- x 1)) (cbrt (sqrt (+ x 1)))) (- (* (cbrt (/ 2 x)) (* (cbrt (/ 2 x)) (cbrt (/ 2 x))))))

prune4.8s

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.1b

localize15.0ms

Local error

Found 4 expressions with local error:

3.0b
(- (/ -1 (+ x -1)) (/ 1 (+ x 1)))
1.8b
(- (/ -2 x) (- (/ -1 (+ x -1)) (/ 1 (+ x 1))))
0.0b
(/ 1 (+ x 1))
0.0b
(/ -1 (+ x -1))

rewrite25.0ms

Algorithm
rewrite-expression-head
Counts
4 → 186
Calls

4 calls. Slowest were:

15.0ms
(- (/ -2 x) (- (/ -1 (+ x -1)) (/ 1 (+ x 1))))
6.0ms
(- (/ -1 (+ x -1)) (/ 1 (+ x 1)))
1.0ms
(/ 1 (+ x 1))

series76.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

36.0ms
(- (/ -2 x) (- (/ -1 (+ x -1)) (/ 1 (+ x 1))))
27.0ms
(- (/ -1 (+ x -1)) (/ 1 (+ x 1)))
7.0ms
(/ -1 (+ x -1))
6.0ms
(/ 1 (+ x 1))

simplify18.2s

Counts
208 → 198
Calls

208 calls. Slowest were:

537.0ms
(fma (/ -1 (- (* x x) (* -1 -1))) (- x -1) (- (* (- x 1) (/ 1 (- (* x x) (* 1 1))))))
385.0ms
(fma -1 (/ 1 (+ x -1)) (- (* (- x 1) (/ 1 (- (* x x) (* 1 1))))))
370.0ms
(fma 1 (/ -1 (+ x -1)) (- (* (- x 1) (/ 1 (- (* x x) (* 1 1))))))

prune1.6s

Pruning

2 alts after pruning (1 fresh and 1 done)

Merged error: 0.1b

localize20.0ms

Local error

Found 4 expressions with local error:

0.3b
(/ 2 (pow x 7))
0.1b
(/ 2 (pow x 5))
0.1b
(/ (/ 2 x) (* x x))
0.0b
(+ (+ (/ (/ 2 x) (* x x)) (/ 2 (pow x 7))) (/ 2 (pow x 5)))

rewrite39.0ms

Algorithm
rewrite-expression-head
Counts
4 → 78
Calls

4 calls. Slowest were:

35.0ms
(+ (+ (/ (/ 2 x) (* x x)) (/ 2 (pow x 7))) (/ 2 (pow x 5)))
3.0ms
(/ (/ 2 x) (* x x))
1.0ms
(/ 2 (pow x 7))

series86.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

26.0ms
(+ (+ (/ (/ 2 x) (* x x)) (/ 2 (pow x 7))) (/ 2 (pow x 5)))
24.0ms
(/ 2 (pow x 7))
20.0ms
(/ 2 (pow x 5))
16.0ms
(/ (/ 2 x) (* x x))

simplify1.2s

Counts
43 → 90
Calls

43 calls. Slowest were:

275.0ms
(/ (* (* (/ 2 x) (/ 2 x)) (/ 2 x)) (* (* (* x x) (* x x)) (* x x)))
160.0ms
(+ (* (- (* (/ (/ 2 x) (* x x)) (/ (/ 2 x) (* x x))) (* (/ 2 (pow x 7)) (/ 2 (pow x 7)))) (pow x 5)) (* (- (/ (/ 2 x) (* x x)) (/ 2 (pow x 7))) 2))
147.0ms
(* (+ (* (/ (/ 2 x) (* x x)) (/ (/ 2 x) (* x x))) (- (* (/ 2 (pow x 7)) (/ 2 (pow x 7))) (* (/ (/ 2 x) (* x x)) (/ 2 (pow x 7))))) (pow x 5))

prune787.0ms

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 0b

regimes49.0ms

Accuracy

99.5% (0.0b remaining)

Error of 0.1b against oracle of 0.0b and baseline of 9.2b

bsearch146.0ms