Average Error: 15.9 → 6.4
Time: 22.7s
Precision: 64
Internal Precision: 128
\[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
\[\begin{array}{l} \mathbf{if}\;\alpha \le 1.7503839507741062 \cdot 10^{+25}:\\ \;\;\;\;{e}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - (\left(\frac{1}{\alpha \cdot \alpha}\right) \cdot \left(4.0 - \frac{8.0}{\alpha}\right) + \left(-\frac{2.0}{\alpha}\right))_*}{2.0}\\ \end{array}\]

Error

Bits error versus alpha

Bits error versus beta

Derivation

  1. Split input into 2 regimes
  2. if alpha < 1.7503839507741062e+25

    1. Initial program 1.0

      \[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
    2. Using strategy rm
    3. Applied div-sub1.0

      \[\leadsto \frac{\color{blue}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right)} + 1.0}{2.0}\]
    4. Applied associate-+l-1.0

      \[\leadsto \frac{\color{blue}{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}}{2.0}\]
    5. Using strategy rm
    6. Applied add-exp-log1.0

      \[\leadsto \color{blue}{e^{\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)}}\]
    7. Using strategy rm
    8. Applied pow11.0

      \[\leadsto e^{\log \color{blue}{\left({\left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)}^{1}\right)}}\]
    9. Applied log-pow1.0

      \[\leadsto e^{\color{blue}{1 \cdot \log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)}}\]
    10. Applied exp-prod1.0

      \[\leadsto \color{blue}{{\left(e^{1}\right)}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}}\]
    11. Simplified1.0

      \[\leadsto {\color{blue}{e}}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}\]

    if 1.7503839507741062e+25 < alpha

    1. Initial program 50.4

      \[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
    2. Using strategy rm
    3. Applied div-sub50.4

      \[\leadsto \frac{\color{blue}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right)} + 1.0}{2.0}\]
    4. Applied associate-+l-48.7

      \[\leadsto \frac{\color{blue}{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}}{2.0}\]
    5. Taylor expanded around -inf 19.0

      \[\leadsto \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{blue}{\left(4.0 \cdot \frac{1}{{\alpha}^{2}} - \left(2.0 \cdot \frac{1}{\alpha} + 8.0 \cdot \frac{1}{{\alpha}^{3}}\right)\right)}}{2.0}\]
    6. Simplified19.0

      \[\leadsto \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{blue}{(\left(\frac{1}{\alpha \cdot \alpha}\right) \cdot \left(4.0 - \frac{8.0}{\alpha}\right) + \left(-\frac{2.0}{\alpha}\right))_*}}{2.0}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification6.4

    \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \le 1.7503839507741062 \cdot 10^{+25}:\\ \;\;\;\;{e}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - (\left(\frac{1}{\alpha \cdot \alpha}\right) \cdot \left(4.0 - \frac{8.0}{\alpha}\right) + \left(-\frac{2.0}{\alpha}\right))_*}{2.0}\\ \end{array}\]

Reproduce

herbie shell --seed 2019008 +o rules:numerics
(FPCore (alpha beta)
  :name "Octave 3.8, jcobi/1"
  :pre (and (> alpha -1) (> beta -1))
  (/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))

Details

Time bar (total: 21.5s)Debug log

sample80.0ms

Algorithm
intervals

simplify103.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

103.0ms
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0)

prune8.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 17.0b

localize32.0ms

Local error

Found 2 expressions with local error:

1.4b
(+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0)
0.0b
(/ (- beta alpha) (+ (+ alpha beta) 2.0))

rewrite46.0ms

Algorithm
rewrite-expression-head
Counts
2 → 62
Calls

2 calls. Slowest were:

38.0ms
(+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0)
7.0ms
(/ (- beta alpha) (+ (+ alpha beta) 2.0))

series74.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

49.0ms
(+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0)
25.0ms
(/ (- beta alpha) (+ (+ alpha beta) 2.0))

simplify2.6s

Counts
51 → 68
Calls

51 calls. Slowest were:

518.0ms
(/ (- beta alpha) (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
235.0ms
(- (+ (* 4.0 (/ 1 (pow beta 2))) 2.0) (* 2.0 (/ 1 beta)))
202.0ms
(- (+ (* 4.0 (/ 1 (pow beta 2))) 2.0) (* 2.0 (/ 1 beta)))

prune699.0ms

Pruning

7 alts after pruning (7 fresh and 0 done)

Merged error: 16.2b

localize16.0ms

Local error

Found 4 expressions with local error:

1.4b
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
0.0b
(/ alpha (+ (+ alpha beta) 2.0))
0.0b
(/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)
0.0b
(/ beta (+ (+ alpha beta) 2.0))

rewrite90.0ms

Algorithm
rewrite-expression-head
Counts
4 → 75
Calls

4 calls. Slowest were:

76.0ms
(/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)
7.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
4.0ms
(/ beta (+ (+ alpha beta) 2.0))

series140.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

53.0ms
(/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)
39.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
24.0ms
(/ beta (+ (+ alpha beta) 2.0))
23.0ms
(/ alpha (+ (+ alpha beta) 2.0))

simplify4.2s

Counts
34 → 87
Calls

34 calls. Slowest were:

653.0ms
(- (* beta (+ (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) (* (+ (+ alpha beta) 2.0) (- (* (/ alpha (+ (+ alpha beta) 2.0)) (/ alpha (+ (+ alpha beta) 2.0))) (* 1.0 1.0))))
487.0ms
(/ alpha (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
455.0ms
(- (* 0.5 beta) (+ (* 0.25 (pow beta 2)) (* 0.25 (* beta alpha))))

prune834.0ms

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 3.6b

localize13.0ms

Local error

Found 4 expressions with local error:

8.2b
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
1.4b
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
0.7b
(exp (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
0.0b
(/ alpha (+ (+ alpha beta) 2.0))

rewrite43.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

16.0ms
(exp (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
15.0ms
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
7.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)

series270.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

122.0ms
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
63.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
45.0ms
(exp (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
41.0ms
(/ alpha (+ (+ alpha beta) 2.0))

simplify3.4s

Counts
37 → 83
Calls

37 calls. Slowest were:

472.0ms
(/ alpha (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
227.0ms
(- (* 4.0 (/ 1 (pow alpha 2))) (+ (* 2.0 (/ 1 alpha)) (* 8.0 (/ 1 (pow alpha 3)))))
224.0ms
(- (* 0.5 alpha) (+ (* 0.25 (pow alpha 2)) (* 0.25 (* beta alpha))))

prune998.0ms

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 3.6b

localize14.0ms

Local error

Found 4 expressions with local error:

8.2b
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
1.4b
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
0.8b
(pow E (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
0.0b
(/ alpha (+ (+ alpha beta) 2.0))

rewrite47.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

25.0ms
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
11.0ms
(pow E (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
8.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)

series281.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

110.0ms
(log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0))
104.0ms
(pow E (log (/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)))
44.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
23.0ms
(/ alpha (+ (+ alpha beta) 2.0))

simplify4.3s

Counts
36 → 83
Calls

36 calls. Slowest were:

476.0ms
(/ alpha (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
470.0ms
(- (* 0.5 alpha) (+ 1.0 (* 0.25 (pow alpha 2))))
306.0ms
(- (* 0.5 alpha) (+ (* 0.25 (pow alpha 2)) (* 0.25 (* beta alpha))))

prune902.0ms

Pruning

10 alts after pruning (9 fresh and 1 done)

Merged error: 3.6b

regimes162.0ms

Accuracy

73.3% (3.3b remaining)

Error of 6.4b against oracle of 3.1b and baseline of 15.5b

bsearch56.0ms

end0.0ms

sample2.1s

Algorithm
intervals