Average Error: 3.5 → 2.6
Time: 2.8m
Precision: 64
Internal Precision: 128
\[\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}\]
\[\begin{array}{l} \mathbf{if}\;\alpha \le 2.0115520703248476 \cdot 10^{+107}:\\ \;\;\;\;\frac{\frac{1}{\left(\beta + \alpha\right) + 2} \cdot \frac{1.0 + \left(\beta \cdot \alpha + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2}}{\left(\alpha + 2\right) + \left(\beta + 1.0\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\frac{\frac{2.0}{\alpha}}{\alpha} + \left(1 - \frac{1.0}{\alpha}\right)}{\left(\beta + \alpha\right) + 2}}{1.0 + \left(\left(\beta + \alpha\right) + 2\right)}\\ \end{array}\]

Error

Bits error versus alpha

Bits error versus beta

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if alpha < 2.0115520703248476e+107

    1. Initial program 0.7

      \[\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}\]
    2. Using strategy rm
    3. Applied *-un-lft-identity0.7

      \[\leadsto \frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\color{blue}{1 \cdot \left(\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0\right)}}\]
    4. Applied associate-/r*0.7

      \[\leadsto \color{blue}{\frac{\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}}\]
    5. Simplified0.7

      \[\leadsto \frac{\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{1}}{\color{blue}{\left(2 + \alpha\right) + \left(1.0 + \beta\right)}}\]
    6. Using strategy rm
    7. Applied div-inv0.7

      \[\leadsto \frac{\frac{\color{blue}{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1} \cdot \frac{1}{\left(\alpha + \beta\right) + 2 \cdot 1}}}{1}}{\left(2 + \alpha\right) + \left(1.0 + \beta\right)}\]

    if 2.0115520703248476e+107 < alpha

    1. Initial program 13.4

      \[\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1.0}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}\]
    2. Taylor expanded around inf 9.4

      \[\leadsto \frac{\frac{\color{blue}{\left(2.0 \cdot \frac{1}{{\alpha}^{2}} + 1\right) - 1.0 \cdot \frac{1}{\alpha}}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}\]
    3. Simplified9.4

      \[\leadsto \frac{\frac{\color{blue}{\frac{\frac{2.0}{\alpha}}{\alpha} + \left(1 - \frac{1.0}{\alpha}\right)}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1.0}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification2.6

    \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \le 2.0115520703248476 \cdot 10^{+107}:\\ \;\;\;\;\frac{\frac{1}{\left(\beta + \alpha\right) + 2} \cdot \frac{1.0 + \left(\beta \cdot \alpha + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2}}{\left(\alpha + 2\right) + \left(\beta + 1.0\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\frac{\frac{2.0}{\alpha}}{\alpha} + \left(1 - \frac{1.0}{\alpha}\right)}{\left(\beta + \alpha\right) + 2}}{1.0 + \left(\left(\beta + \alpha\right) + 2\right)}\\ \end{array}\]

Reproduce

herbie shell --seed 2019004 
(FPCore (alpha beta)
  :name "Octave 3.8, jcobi/3"
  :pre (and (> alpha -1) (> beta -1))
  (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0)))

Details

Time bar (total: 2.8m)Debug log

sample48.0ms

Algorithm
intervals

simplify780.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

780.0ms
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0))

prune22.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 3.5b

localize35.0ms

Local error

Found 4 expressions with local error:

3.4b
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
0.1b
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0))
0.0b
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
0.0b
(+ (+ (+ alpha beta) (* 2 1)) 1.0)

rewrite225.0ms

Algorithm
rewrite-expression-head
Counts
4 → 348
Calls

4 calls. Slowest were:

76.0ms
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
71.0ms
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0))
41.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))

series236.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

92.0ms
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
80.0ms
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0))
56.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
8.0ms
(+ (+ (+ alpha beta) (* 2 1)) 1.0)

simplify1.4m

Counts
524 → 360
Calls

524 calls. Slowest were:

794.0ms
(/ (/ (* (cbrt (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))) (cbrt (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))))) 1) (* (cbrt (+ (+ (+ alpha beta) (* 2 1)) 1.0)) (cbrt (+ (+ (+ alpha beta) (* 2 1)) 1.0))))
733.0ms
(/ (/ (/ (cbrt (+ (+ (+ alpha beta) (* beta alpha)) 1.0)) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1.0))
700.0ms
(+ (pow (+ (+ alpha beta) (* beta alpha)) 3) (pow 1.0 3))

prune5.9s

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 0.4b

localize16.0ms

Local error

Found 4 expressions with local error:

3.4b
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
0.1b
(/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
0.0b
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
0.0b
(+ (+ 2 alpha) (+ 1.0 beta))

rewrite150.0ms

Algorithm
rewrite-expression-head
Counts
4 → 172
Calls

4 calls. Slowest were:

72.0ms
(/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
41.0ms
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
31.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))

series265.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

115.0ms
(/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
96.0ms
(/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1)))
39.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
14.0ms
(+ (+ 2 alpha) (+ 1.0 beta))

simplify27.4s

Counts
213 → 184
Calls

213 calls. Slowest were:

849.0ms
(* (+ (* 2 2) (- (* alpha alpha) (* 2 alpha))) (+ (* 1.0 1.0) (- (* beta beta) (* 1.0 beta))))
651.0ms
(+ (pow (+ (+ alpha beta) (* beta alpha)) 3) (pow 1.0 3))
531.0ms
(+ (* (+ (pow 2 3) (pow alpha 3)) (- 1.0 beta)) (* (+ (* 2 2) (- (* alpha alpha) (* 2 alpha))) (- (* 1.0 1.0) (* beta beta))))

prune2.6s

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.4b

localize20.0ms

Local error

Found 4 expressions with local error:

3.4b
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
0.1b
(/ (/ (* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1)))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
0.1b
(* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1))))
0.0b
(/ 1 (+ (+ alpha beta) (* 2 1)))

rewrite128.0ms

Algorithm
rewrite-expression-head
Counts
4 → 132
Calls

4 calls. Slowest were:

59.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
38.0ms
(/ (/ (* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1)))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
26.0ms
(* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1))))

series244.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

95.0ms
(/ (/ (* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1)))) 1) (+ (+ 2 alpha) (+ 1.0 beta)))
83.0ms
(* (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (/ 1 (+ (+ alpha beta) (* 2 1))))
47.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
18.0ms
(/ 1 (+ (+ alpha beta) (* 2 1)))

simplify16.6s

Counts
127 → 144
Calls

127 calls. Slowest were:

773.0ms
(+ (pow (+ (+ alpha beta) (* beta alpha)) 3) (pow 1.0 3))
579.0ms
(- (* (+ (+ alpha beta) (* beta alpha)) (+ (+ alpha beta) (* beta alpha))) (* 1.0 1.0))
541.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (pow (+ alpha beta) 3) (pow (* 2 1) 3)))

prune2.4s

Pruning

8 alts after pruning (7 fresh and 1 done)

Merged error: 0.4b

localize8.0ms

Local error

Found 4 expressions with local error:

13.3b
(sqrt (/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta))))
13.3b
(sqrt (/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta))))
3.4b
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
3.4b
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))

rewrite147.0ms

Algorithm
rewrite-expression-head
Counts
4 → 136
Calls

4 calls. Slowest were:

41.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
38.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
38.0ms
(sqrt (/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta))))

series411.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

154.0ms
(sqrt (/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta))))
147.0ms
(sqrt (/ (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) 1) (+ (+ 2 alpha) (+ 1.0 beta))))
60.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))
50.0ms
(/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2 1)))

simplify18.9s

Counts
128 → 148
Calls

128 calls. Slowest were:

669.0ms
(+ (pow (+ (+ alpha beta) (* beta alpha)) 3) (pow 1.0 3))
660.0ms
(- (* (+ (+ alpha beta) (* beta alpha)) (+ (+ alpha beta) (* beta alpha))) (* 1.0 1.0))
583.0ms
(+ (pow (+ (+ alpha beta) (* beta alpha)) 3) (pow 1.0 3))

prune3.1s

Pruning

8 alts after pruning (6 fresh and 2 done)

Merged error: 0.4b

regimes182.0ms

Accuracy

36.3% (1.5b remaining)

Error of 2.6b against oracle of 1.1b and baseline of 3.5b

bsearch115.0ms

end0.0ms

sample1.7s

Algorithm
intervals