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

Error

Bits error versus alpha

Bits error versus beta

Derivation

  1. Split input into 2 regimes
  2. if (/ (- beta alpha) (+ (+ alpha beta) 2.0)) < -0.9999999999943141

    1. Initial program 60.2

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

      \[\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-58.4

      \[\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 clear-num58.4

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

      \[\leadsto \frac{\frac{1}{\frac{\left(\alpha + \beta\right) + 2.0}{\beta}} - \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}\]
    8. Simplified10.8

      \[\leadsto \frac{\frac{1}{\frac{\left(\alpha + \beta\right) + 2.0}{\beta}} - \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}\]

    if -0.9999999999943141 < (/ (- beta alpha) (+ (+ alpha beta) 2.0))

    1. Initial program 0.3

      \[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
    2. Using strategy rm
    3. Applied add-log-exp0.3

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

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

Reproduce

herbie shell --seed 2019005 +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: 25.3s)Debug log

sample53.0ms

Algorithm
intervals

simplify75.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune4.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 14.3b

localize19.0ms

Local error

Found 2 expressions with local error:

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

rewrite47.0ms

Algorithm
rewrite-expression-head
Counts
2 → 62
Calls

2 calls. Slowest were:

34.0ms
(+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0)
11.0ms
(/ (- beta alpha) (+ (+ alpha beta) 2.0))

series98.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

64.0ms
(+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0)
33.0ms
(/ (- beta alpha) (+ (+ alpha beta) 2.0))

simplify2.7s

Counts
51 → 68
Calls

51 calls. Slowest were:

609.0ms
(/ (- beta alpha) (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
199.0ms
(- (+ (* 4.0 (/ 1 (pow beta 2))) 2.0) (* 2.0 (/ 1 beta)))
189.0ms
(/ (* (* (- beta alpha) (- beta alpha)) (- beta alpha)) (* (* (+ (+ alpha beta) 2.0) (+ (+ alpha beta) 2.0)) (+ (+ alpha beta) 2.0)))

prune685.0ms

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 13.6b

localize26.0ms

Local error

Found 4 expressions with local error:

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

rewrite57.0ms

Algorithm
rewrite-expression-head
Counts
4 → 98
Calls

4 calls. Slowest were:

34.0ms
(- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0))
13.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
4.0ms
(/ beta (+ (+ alpha beta) 2.0))

series165.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

65.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
39.0ms
(/ beta (+ (+ alpha beta) 2.0))
35.0ms
(- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0))
26.0ms
(/ alpha (+ (+ alpha beta) 2.0))

simplify9.3s

Counts
76 → 110
Calls

76 calls. Slowest were:

702.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))))
484.0ms
(/ beta (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
482.0ms
(fma (/ beta (- (* (+ alpha beta) (+ alpha beta)) (* 2.0 2.0))) (- (+ alpha beta) 2.0) (- (* (cbrt (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) (* (cbrt (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) (cbrt (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0))))))

prune1.1s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 3.8b

localize12.0ms

Local error

Found 4 expressions with local error:

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

rewrite34.0ms

Algorithm
rewrite-expression-head
Counts
4 → 73
Calls

4 calls. Slowest were:

19.0ms
(log (exp (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)))
7.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
3.0ms
(/ beta (+ (+ alpha beta) 2.0))

series177.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

64.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
50.0ms
(log (exp (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)))
39.0ms
(/ beta (+ (+ alpha beta) 2.0))
24.0ms
(/ alpha (+ (+ alpha beta) 2.0))

simplify3.0s

Counts
36 → 85
Calls

36 calls. Slowest were:

533.0ms
(/ alpha (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
414.0ms
(- (* 0.5 beta) (+ (* 0.25 (pow beta 2)) (* 0.25 (* beta alpha))))
262.0ms
(- (* 4.0 (/ 1 (pow alpha 2))) (+ (* 2.0 (/ 1 alpha)) (* 8.0 (/ 1 (pow alpha 3)))))

prune1.1s

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 3.8b

localize11.0ms

Local error

Found 4 expressions with local error:

2.0b
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
0.4b
(/ 1 (/ (+ (+ alpha beta) 2.0) beta))
0.0b
(/ alpha (+ (+ alpha beta) 2.0))
0.0b
(/ (+ (+ alpha beta) 2.0) beta)

rewrite19.0ms

Algorithm
rewrite-expression-head
Counts
4 → 73
Calls

4 calls. Slowest were:

8.0ms
(/ (+ (+ alpha beta) 2.0) beta)
6.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
2.0ms
(/ 1 (/ (+ (+ alpha beta) 2.0) beta))

series180.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

67.0ms
(- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)
51.0ms
(/ 1 (/ (+ (+ alpha beta) 2.0) beta))
39.0ms
(/ alpha (+ (+ alpha beta) 2.0))
23.0ms
(/ (+ (+ alpha beta) 2.0) beta)

simplify3.1s

Counts
30 → 85
Calls

30 calls. Slowest were:

884.0ms
(/ alpha (+ (pow (+ alpha beta) 3) (pow 2.0 3)))
490.0ms
(+ (pow (+ alpha beta) 3) (pow 2.0 3))
456.0ms
(- (* 0.5 beta) (+ (* 0.25 (pow beta 2)) (* 0.25 (* beta alpha))))

prune958.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 3.8b

regimes414.0ms

Accuracy

98.6% (0.2b remaining)

Error of 3.1b against oracle of 2.9b and baseline of 15.6b

bsearch6.0ms

end0.0ms

sample1.9s

Algorithm
intervals