Average Error: 45.7 → 0.0
Time: 7.5s
Precision: 64
Internal Precision: 128
\[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
\[\begin{array}{l} \mathbf{if}\;i \le 233.33044365925295:\\ \;\;\;\;\frac{\frac{i}{4} \cdot i}{(\left(4 \cdot i\right) \cdot i + \left(-1.0\right))_*}\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{0.015625}{i}\right) \cdot \left(\frac{1}{i}\right) + \frac{1}{16})_* + \frac{0.00390625}{{i}^{4}}\\ \end{array}\]

Error

Bits error versus i

Derivation

  1. Split input into 2 regimes
  2. if i < 233.33044365925295

    1. Initial program 44.3

      \[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
    2. Simplified0.0

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

    if 233.33044365925295 < i

    1. Initial program 47.2

      \[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
    2. Simplified30.8

      \[\leadsto \color{blue}{\frac{i \cdot \frac{i}{4}}{(\left(i \cdot 4\right) \cdot i + \left(-1.0\right))_*}}\]
    3. Taylor expanded around inf 0.0

      \[\leadsto \color{blue}{0.015625 \cdot \frac{1}{{i}^{2}} + \left(\frac{1}{16} + 0.00390625 \cdot \frac{1}{{i}^{4}}\right)}\]
    4. Simplified0.0

      \[\leadsto \color{blue}{(\left(\frac{0.015625}{i}\right) \cdot \left(\frac{1}{i}\right) + \frac{1}{16})_* + \frac{0.00390625}{{i}^{4}}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.0

    \[\leadsto \begin{array}{l} \mathbf{if}\;i \le 233.33044365925295:\\ \;\;\;\;\frac{\frac{i}{4} \cdot i}{(\left(4 \cdot i\right) \cdot i + \left(-1.0\right))_*}\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{0.015625}{i}\right) \cdot \left(\frac{1}{i}\right) + \frac{1}{16})_* + \frac{0.00390625}{{i}^{4}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019005 +o rules:numerics
(FPCore (i)
  :name "Octave 3.8, jcobi/4, as called"
  :pre (and (> i 0))
  (/ (/ (* (* i i) (* i i)) (* (* 2 i) (* 2 i))) (- (* (* 2 i) (* 2 i)) 1.0)))

Details

Time bar (total: 7.1s)Debug log

sample55.0ms

Algorithm
intervals

simplify177.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

177.0ms
(/ (/ (* (* i i) (* i i)) (* (* 2 i) (* 2 i))) (- (* (* 2 i) (* 2 i)) 1.0))

prune9.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 16.0b

localize18.0ms

Local error

Found 1 expressions with local error:

16.0b
(/ (* i (/ i 4)) (fma (* i 4) i (- 1.0)))

rewrite8.0ms

Algorithm
rewrite-expression-head
Counts
1 → 22
Calls

1 calls. Slowest were:

7.0ms
(/ (* i (/ i 4)) (fma (* i 4) i (- 1.0)))

series56.0ms

Counts
1 → 3
Calls

1 calls. Slowest were:

56.0ms
(/ (* i (/ i 4)) (fma (* i 4) i (- 1.0)))

simplify633.0ms

Counts
15 → 25
Calls

15 calls. Slowest were:

153.0ms
(- (log (* i (/ i 4))) (log (fma (* i 4) i (- 1.0))))
146.0ms
(- (+ (* 0.25 (pow i 2)) (+ (* 1.0 (pow i 4)) (* 4.0 (pow i 6)))))
76.0ms
(/ (* (* (* i (/ i 4)) (* i (/ i 4))) (* i (/ i 4))) (* (* (fma (* i 4) i (- 1.0)) (fma (* i 4) i (- 1.0))) (fma (* i 4) i (- 1.0))))

prune238.0ms

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 0b

localize14.0ms

Local error

Found 2 expressions with local error:

15.6b
(/ i (fma 16 (* i i) (* -4 1.0)))
0.1b
(* i (/ i (fma 16 (* i i) (* -4 1.0))))

rewrite2.0ms

Algorithm
rewrite-expression-head
Counts
2 → 30
Calls

2 calls. Slowest were:

1.0ms
(* i (/ i (fma 16 (* i i) (* -4 1.0))))
1.0ms
(/ i (fma 16 (* i i) (* -4 1.0)))

series76.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

45.0ms
(/ i (fma 16 (* i i) (* -4 1.0)))
31.0ms
(* i (/ i (fma 16 (* i i) (* -4 1.0))))

simplify958.0ms

Counts
14 → 36
Calls

14 calls. Slowest were:

318.0ms
(- (+ (* 0.25 i) (+ (* 1.0 (pow i 3)) (* 4.0 (pow i 5)))))
150.0ms
(- (+ (* 0.25 (pow i 2)) (+ (* 1.0 (pow i 4)) (* 4.0 (pow i 6)))))
133.0ms
(+ (* 0.00390625 (/ 1 (pow i 5))) (+ (* 0.015625 (/ 1 (pow i 3))) (* 1/16 (/ 1 i))))

prune264.0ms

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 0b

localize8.0ms

Local error

Found 2 expressions with local error:

16.1b
(* (* i i) (/ 1 (fma 16 (* i i) (* -4 1.0))))
0.4b
(/ 1 (fma 16 (* i i) (* -4 1.0)))

rewrite5.0ms

Algorithm
rewrite-expression-head
Counts
2 → 39
Calls

2 calls. Slowest were:

4.0ms
(* (* i i) (/ 1 (fma 16 (* i i) (* -4 1.0))))
1.0ms
(/ 1 (fma 16 (* i i) (* -4 1.0)))

series53.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

33.0ms
(* (* i i) (/ 1 (fma 16 (* i i) (* -4 1.0))))
20.0ms
(/ 1 (fma 16 (* i i) (* -4 1.0)))

simplify1.6s

Counts
20 → 45
Calls

20 calls. Slowest were:

730.0ms
(* (* (* (* i i) (* i i)) (* i i)) (* (* (/ 1 (fma 16 (* i i) (* -4 1.0))) (/ 1 (fma 16 (* i i) (* -4 1.0)))) (/ 1 (fma 16 (* i i) (* -4 1.0)))))
167.0ms
(- (+ (* 0.25 (pow i 2)) (+ (* 1.0 (pow i 4)) (* 4.0 (pow i 6)))))
132.0ms
(- (+ (* 1.0 (pow i 2)) (+ (* 4.0 (pow i 4)) 0.25)))

prune374.0ms

Pruning

4 alts after pruning (1 fresh and 3 done)

Merged error: 0b

localize18.0ms

Local error

Found 3 expressions with local error:

0.2b
(/ 0.00390625 (pow i 4))
0.1b
(fma (/ 0.015625 i) (/ 1 i) 1/16)
0.0b
(+ (fma (/ 0.015625 i) (/ 1 i) 1/16) (/ 0.00390625 (pow i 4)))

rewrite9.0ms

Algorithm
rewrite-expression-head
Counts
3 → 43
Calls

3 calls. Slowest were:

7.0ms
(+ (fma (/ 0.015625 i) (/ 1 i) 1/16) (/ 0.00390625 (pow i 4)))
1.0ms
(/ 0.00390625 (pow i 4))
0.0ms
(fma (/ 0.015625 i) (/ 1 i) 1/16)

series73.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

40.0ms
(fma (/ 0.015625 i) (/ 1 i) 1/16)
18.0ms
(+ (fma (/ 0.015625 i) (/ 1 i) 1/16) (/ 0.00390625 (pow i 4)))
15.0ms
(/ 0.00390625 (pow i 4))

simplify509.0ms

Counts
17 → 52
Calls

17 calls. Slowest were:

106.0ms
(* (exp (fma (/ 0.015625 i) (/ 1 i) 1/16)) (exp (/ 0.00390625 (pow i 4))))
77.0ms
(+ (* 0.015625 (/ 1 (pow i 2))) (+ 1/16 (* 0.00390625 (/ 1 (pow i 4)))))
75.0ms
(+ (* 0.015625 (/ 1 (pow i 2))) (+ 1/16 (* 0.00390625 (/ 1 (pow i 4)))))

prune501.0ms

Pruning

4 alts after pruning (0 fresh and 4 done)

Merged error: 0b

regimes56.0ms

Accuracy

100% (0.0b remaining)

Error of 0.0b against oracle of 0.0b and baseline of 15.2b

bsearch82.0ms

end0.0ms

sample1.2s

Algorithm
intervals