Average Error: 0.1 → 0.1
Time: 43.6s
Precision: 64
Internal Precision: 128
\[\frac{1 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}{2 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}\]
\[\frac{1}{\frac{(\left(\frac{2 \cdot t}{t + 1}\right) \cdot \left(\frac{2 \cdot t}{t + 1}\right) + 2)_*}{(\left(\frac{2 \cdot t}{t + 1}\right) \cdot \left(\frac{2 \cdot t}{t + 1}\right) + 1)_*}}\]

Error

Bits error versus t

Derivation

  1. Initial program 0.1

    \[\frac{1 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}{2 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}\]
  2. Simplified0.1

    \[\leadsto \color{blue}{\frac{(\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 1)_*}{(\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 2)_*}}\]
  3. Using strategy rm
  4. Applied *-un-lft-identity0.1

    \[\leadsto \frac{\color{blue}{1 \cdot (\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 1)_*}}{(\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 2)_*}\]
  5. Applied associate-/l*0.1

    \[\leadsto \color{blue}{\frac{1}{\frac{(\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 2)_*}{(\left(\frac{t \cdot 2}{1 + t}\right) \cdot \left(\frac{t \cdot 2}{1 + t}\right) + 1)_*}}}\]
  6. Final simplification0.1

    \[\leadsto \frac{1}{\frac{(\left(\frac{2 \cdot t}{t + 1}\right) \cdot \left(\frac{2 \cdot t}{t + 1}\right) + 2)_*}{(\left(\frac{2 \cdot t}{t + 1}\right) \cdot \left(\frac{2 \cdot t}{t + 1}\right) + 1)_*}}\]

Reproduce

herbie shell --seed 2019005 +o rules:numerics
(FPCore (t)
  :name "Kahan p13 Example 1"
  (/ (+ 1 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t)))) (+ 2 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t))))))

Details

Time bar (total: 43.3s)Debug log

sample53.0ms

Algorithm
intervals

simplify518.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

518.0ms
(/ (+ 1 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t)))) (+ 2 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t)))))

prune17.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 0.0b

localize40.0ms

Local error

Found 4 expressions with local error:

0.0b
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1) (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2))
0.0b
(/ (* t 2) (+ 1 t))
0.0b
(/ (* t 2) (+ 1 t))
0.0b
(/ (* t 2) (+ 1 t))

rewrite28.0ms

Algorithm
rewrite-expression-head
Counts
4 → 98
Calls

4 calls. Slowest were:

8.0ms
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1) (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2))
6.0ms
(/ (* t 2) (+ 1 t))
6.0ms
(/ (* t 2) (+ 1 t))

series115.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

70.0ms
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1) (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2))
18.0ms
(/ (* t 2) (+ 1 t))
14.0ms
(/ (* t 2) (+ 1 t))
13.0ms
(/ (* t 2) (+ 1 t))

simplify11.8s

Counts
76 → 110
Calls

76 calls. Slowest were:

712.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))
672.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))
648.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))

prune1.1s

Pruning

3 alts after pruning (2 fresh and 1 done)

Merged error: 0b

localize10.0ms

Local error

Found 4 expressions with local error:

0.0b
(/ 1 (/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2) (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))
0.0b
(/ (* t 2) (+ 1 t))
0.0b
(/ (* t 2) (+ 1 t))
0.0b
(/ (* t 2) (+ 1 t))

rewrite12.0ms

Algorithm
rewrite-expression-head
Counts
4 → 99
Calls

4 calls. Slowest were:

3.0ms
(/ (* t 2) (+ 1 t))
3.0ms
(/ (* t 2) (+ 1 t))
3.0ms
(/ (* t 2) (+ 1 t))

series112.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

72.0ms
(/ 1 (/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2) (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))
18.0ms
(/ (* t 2) (+ 1 t))
11.0ms
(/ (* t 2) (+ 1 t))
11.0ms
(/ (* t 2) (+ 1 t))

simplify8.5s

Counts
66 → 111
Calls

66 calls. Slowest were:

749.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))
728.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))
632.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))

prune1.2s

Pruning

3 alts after pruning (1 fresh and 2 done)

Merged error: 0b

localize12.0ms

Local error

Found 4 expressions with local error:

0.0b
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2) (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))
0.0b
(sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1))
0.0b
(hypot 1 (/ (* t 2) (+ 1 t)))
0.0b
(/ (* t 2) (+ 1 t))

rewrite9.0ms

Algorithm
rewrite-expression-head
Counts
4 → 88
Calls

4 calls. Slowest were:

4.0ms
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2) (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))
3.0ms
(/ (* t 2) (+ 1 t))
1.0ms
(sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1))

series214.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

95.0ms
(/ (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2) (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))
55.0ms
(sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1))
46.0ms
(hypot 1 (/ (* t 2) (+ 1 t)))
17.0ms
(/ (* t 2) (+ 1 t))

simplify17.4s

Counts
75 → 100
Calls

75 calls. Slowest were:

695.0ms
(/ (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2)) (sqrt (* (cbrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)) (cbrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1)))))
654.0ms
(- (+ (* 2 t) (* 2 (pow t 3))) (* 2 (pow t 2)))
420.0ms
(/ (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 2)) (cbrt (sqrt (fma (/ (* t 2) (+ 1 t)) (/ (* t 2) (+ 1 t)) 1))))

prune1.3s

Pruning

3 alts after pruning (0 fresh and 3 done)

Merged error: 0b

regimes36.0ms

Accuracy

0% (0.0b remaining)

Error of 0.1b against oracle of 0.1b and baseline of 0.1b

bsearch1.0ms

end0.0ms

sample993.0ms

Algorithm
intervals