Average Error: 52.4 → 12.6
Time: 2.4m
Precision: 64
Internal Precision: 128
\[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
\[\begin{array}{l} \mathbf{if}\;i \le 8.925236338593809 \cdot 10^{+74}:\\ \;\;\;\;\frac{\frac{\beta + \left(i + \alpha\right)}{(2 \cdot i + \alpha)_* + \beta} \cdot \frac{i}{(2 \cdot i + \alpha)_* + \beta}}{\frac{\left(\left(\beta + \alpha\right) + 2 \cdot i\right) \cdot \left(\left(\beta + \alpha\right) + 2 \cdot i\right) - 1.0}{i \cdot \left(\left(\beta + \alpha\right) + i\right) + \beta \cdot \alpha}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\log \left(e^{\frac{i}{(2 \cdot i + \alpha)_* + \beta}}\right) \cdot \frac{\beta + \left(i + \alpha\right)}{(2 \cdot i + \alpha)_* + \beta}}{4}\\ \end{array}\]

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Derivation

  1. Split input into 2 regimes
  2. if i < 8.925236338593809e+74

    1. Initial program 23.6

      \[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    2. Using strategy rm
    3. Applied associate-/l*11.5

      \[\leadsto \frac{\color{blue}{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    4. Using strategy rm
    5. Applied associate-/r/11.5

      \[\leadsto \frac{\color{blue}{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)} \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    6. Applied associate-/l*11.5

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

      \[\leadsto \frac{\color{blue}{\frac{i}{(2 \cdot i + \alpha)_* + \beta} \cdot \frac{\left(\alpha + i\right) + \beta}{(2 \cdot i + \alpha)_* + \beta}}}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}\]

    if 8.925236338593809e+74 < i

    1. Initial program 61.8

      \[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    2. Using strategy rm
    3. Applied associate-/l*46.9

      \[\leadsto \frac{\color{blue}{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    4. Using strategy rm
    5. Applied associate-/r/46.9

      \[\leadsto \frac{\color{blue}{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)} \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
    6. Applied associate-/l*46.9

      \[\leadsto \color{blue}{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}}\]
    7. Simplified46.9

      \[\leadsto \frac{\color{blue}{\frac{i}{(2 \cdot i + \alpha)_* + \beta} \cdot \frac{\left(\alpha + i\right) + \beta}{(2 \cdot i + \alpha)_* + \beta}}}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}\]
    8. Taylor expanded around 0 14.2

      \[\leadsto \frac{\frac{i}{(2 \cdot i + \alpha)_* + \beta} \cdot \frac{\left(\alpha + i\right) + \beta}{(2 \cdot i + \alpha)_* + \beta}}{\color{blue}{4}}\]
    9. Using strategy rm
    10. Applied add-log-exp13.0

      \[\leadsto \frac{\color{blue}{\log \left(e^{\frac{i}{(2 \cdot i + \alpha)_* + \beta}}\right)} \cdot \frac{\left(\alpha + i\right) + \beta}{(2 \cdot i + \alpha)_* + \beta}}{4}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification12.6

    \[\leadsto \begin{array}{l} \mathbf{if}\;i \le 8.925236338593809 \cdot 10^{+74}:\\ \;\;\;\;\frac{\frac{\beta + \left(i + \alpha\right)}{(2 \cdot i + \alpha)_* + \beta} \cdot \frac{i}{(2 \cdot i + \alpha)_* + \beta}}{\frac{\left(\left(\beta + \alpha\right) + 2 \cdot i\right) \cdot \left(\left(\beta + \alpha\right) + 2 \cdot i\right) - 1.0}{i \cdot \left(\left(\beta + \alpha\right) + i\right) + \beta \cdot \alpha}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\log \left(e^{\frac{i}{(2 \cdot i + \alpha)_* + \beta}}\right) \cdot \frac{\beta + \left(i + \alpha\right)}{(2 \cdot i + \alpha)_* + \beta}}{4}\\ \end{array}\]

Reproduce

herbie shell --seed 2019002 +o rules:numerics
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/4"
  :pre (and (> alpha -1) (> beta -1) (> i 1))
  (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))

Details

Time bar (total: 2.4m)Debug log

sample45.0ms

Algorithm
intervals

simplify308.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

308.0ms
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))

prune14.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 52.3b

localize33.0ms

Local error

Found 4 expressions with local error:

46.1b
(/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))))
35.3b
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
0.1b
(* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
0.0b
(* i (+ (+ alpha beta) i))

rewrite119.0ms

Algorithm
rewrite-expression-head
Counts
4 → 163
Calls

4 calls. Slowest were:

53.0ms
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
27.0ms
(/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))))
24.0ms
(* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))

series1.0s

Counts
4 → 12
Calls

4 calls. Slowest were:

772.0ms
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
152.0ms
(/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))))
88.0ms
(* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
25.0ms
(* i (+ (+ alpha beta) i))

simplify43.2s

Counts
173 → 175
Calls

173 calls. Slowest were:

1.0s
(* (* i (+ (+ alpha beta) i)) (* i (+ (+ alpha beta) i)))
965.0ms
(+ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)
943.0ms
(* (* i (+ (+ alpha beta) i)) (* i (+ (+ alpha beta) i)))

prune2.6s

Pruning

5 alts after pruning (5 fresh and 0 done)

Merged error: 35.2b

localize8.0ms

Local error

Found 4 expressions with local error:

40.5b
(/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
35.3b
(/ (/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
4.3b
(/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
0.0b
(* i (+ (+ alpha beta) i))

rewrite74.0ms

Algorithm
rewrite-expression-head
Counts
4 → 172
Calls

4 calls. Slowest were:

27.0ms
(/ (/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
23.0ms
(/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
16.0ms
(/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))

series1.1s

Counts
4 → 12
Calls

4 calls. Slowest were:

840.0ms
(/ (/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
142.0ms
(/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
135.0ms
(/ (* i (+ (+ alpha beta) i)) (/ (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
24.0ms
(* i (+ (+ alpha beta) i))

simplify41.2s

Counts
190 → 184
Calls

190 calls. Slowest were:

1.0s
(/ (/ i (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) 1)
930.0ms
(* (+ (+ alpha beta) (* 2 i)) (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i))))
829.0ms
(* (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i))) (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i))))

prune2.7s

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 35.2b

localize46.0ms

Local error

Found 4 expressions with local error:

40.5b
(/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
0.1b
(/ (* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
0.0b
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))
0.0b
(* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta)))

rewrite110.0ms

Algorithm
rewrite-expression-head
Counts
4 → 165
Calls

4 calls. Slowest were:

43.0ms
(/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
38.0ms
(/ (* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
16.0ms
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))

series1.0s

Counts
4 → 12
Calls

4 calls. Slowest were:

745.0ms
(/ (* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
191.0ms
(/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
68.0ms
(* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta)))
42.0ms
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))

simplify35.9s

Counts
161 → 177
Calls

161 calls. Slowest were:

976.0ms
(/ (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
940.0ms
(/ (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (cbrt (+ (* beta alpha) (* i (+ (+ alpha beta) i)))))
895.0ms
(/ (/ i (+ (fma 2 i alpha) beta)) (/ (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) 1))

prune2.7s

Pruning

6 alts after pruning (4 fresh and 2 done)

Merged error: 9.3b

localize9.0ms

Local error

Found 3 expressions with local error:

0.0b
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))
0.0b
(* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta)))
0.0b
(/ i (+ (fma 2 i alpha) beta))

rewrite24.0ms

Algorithm
rewrite-expression-head
Counts
3 → 89
Calls

3 calls. Slowest were:

11.0ms
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))
11.0ms
(* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta)))
2.0ms
(/ i (+ (fma 2 i alpha) beta))

series183.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

84.0ms
(* (/ i (+ (fma 2 i alpha) beta)) (/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta)))
61.0ms
(/ i (+ (fma 2 i alpha) beta))
39.0ms
(/ (+ (+ alpha i) beta) (+ (fma 2 i alpha) beta))

simplify4.8s

Counts
71 → 98
Calls

71 calls. Slowest were:

720.0ms
(/ i (- (* (fma 2 i alpha) (fma 2 i alpha)) (* beta beta)))
604.0ms
(/ i (+ (pow (fma 2 i alpha) 3) (pow beta 3)))
547.0ms
(+ (pow (+ alpha i) 3) (pow beta 3))

prune1.1s

Pruning

13 alts after pruning (11 fresh and 2 done)

Merged error: 9.3b

regimes318.0ms

Accuracy

31.5% (4.0b remaining)

Error of 12.6b against oracle of 8.6b and baseline of 14.5b

bsearch94.0ms

end0.0ms

sample2.3s

Algorithm
intervals