Average Error: 2.2 → 0.3
Time: 27.1s
Precision: 64
Internal Precision: 128
\[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}\]
\[\begin{array}{l} \mathbf{if}\;k \le 7.573246734597647 \cdot 10^{+136}:\\ \;\;\;\;\frac{\frac{{k}^{m} \cdot a}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{{k}^{4}}{a}}\right) \cdot 99 + \left((\left(\frac{-10}{k}\right) \cdot \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right) + \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right))_*\right))_*\\ \end{array}\]

Error

Bits error versus a

Bits error versus k

Bits error versus m

Derivation

  1. Split input into 2 regimes
  2. if k < 7.573246734597647e+136

    1. Initial program 0.1

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}\]
    2. Simplified0.1

      \[\leadsto \color{blue}{\frac{{k}^{m} \cdot a}{(k \cdot \left(k + 10\right) + 1)_*}}\]
    3. Using strategy rm
    4. Applied add-sqr-sqrt0.2

      \[\leadsto \frac{{k}^{m} \cdot a}{\color{blue}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*} \cdot \sqrt{(k \cdot \left(k + 10\right) + 1)_*}}}\]
    5. Applied associate-/r*0.2

      \[\leadsto \color{blue}{\frac{\frac{{k}^{m} \cdot a}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}}\]

    if 7.573246734597647e+136 < k

    1. Initial program 9.4

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}\]
    2. Simplified9.4

      \[\leadsto \color{blue}{\frac{{k}^{m} \cdot a}{(k \cdot \left(k + 10\right) + 1)_*}}\]
    3. Taylor expanded around -inf 63.0

      \[\leadsto \color{blue}{\left(99 \cdot \frac{a \cdot e^{m \cdot \left(\log -1 - \log \left(\frac{-1}{k}\right)\right)}}{{k}^{4}} + \frac{a \cdot e^{m \cdot \left(\log -1 - \log \left(\frac{-1}{k}\right)\right)}}{{k}^{2}}\right) - 10 \cdot \frac{a \cdot e^{m \cdot \left(\log -1 - \log \left(\frac{-1}{k}\right)\right)}}{{k}^{3}}}\]
    4. Simplified0.7

      \[\leadsto \color{blue}{(\left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{{k}^{4}}{a}}\right) \cdot 99 + \left((\left(\frac{-10}{k}\right) \cdot \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right) + \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right))_*\right))_*}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;k \le 7.573246734597647 \cdot 10^{+136}:\\ \;\;\;\;\frac{\frac{{k}^{m} \cdot a}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}}{\sqrt{(k \cdot \left(k + 10\right) + 1)_*}}\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{{k}^{4}}{a}}\right) \cdot 99 + \left((\left(\frac{-10}{k}\right) \cdot \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right) + \left(\frac{{\left(e^{m}\right)}^{\left(\log k\right)}}{\frac{k}{\frac{a}{k}}}\right))_*\right))_*\\ \end{array}\]

Reproduce

herbie shell --seed 2019010 +o rules:numerics
(FPCore (a k m)
  :name "Falkner and Boettcher, Appendix A"
  (/ (* a (pow k m)) (+ (+ 1 (* 10 k)) (* k k))))

Details

Time bar (total: 26.2s)Debug log

sample138.0ms

Algorithm
intervals

simplify81.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
81.0ms
(/ (* a (pow k m)) (+ (+ 1 (* 10 k)) (* k k)))

prune16.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 2.8b

localize16.0ms

Local error

Found 2 expressions with local error:

2.8b
(/ (* (pow k m) a) (fma k (+ k 10) 1))
0.0b
(* (pow k m) a)

rewrite7.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
add-exp-log
add-cbrt-cube
associate-*l*
times-frac
associate-/r*
add-log-exp
log1p-expm1-u
pow1
expm1-log1p-u
associate-/l*
div-inv
*-commutative
div-exp
frac-2neg
clear-num
cbrt-undiv
Counts
2 → 34
Calls
2 calls:
Slowest
5.0ms
(/ (* (pow k m) a) (fma k (+ k 10) 1))
2.0ms
(* (pow k m) a)

series234.0ms

Counts
2 → 6
Calls
2 calls:
Slowest
131.0ms
(/ (* (pow k m) a) (fma k (+ k 10) 1))
103.0ms
(* (pow k m) a)

simplify2.1s

Counts
20 → 40
Calls
20 calls:
Slowest
831.0ms
(/ (* (* (* (pow k m) a) (* (pow k m) a)) (* (pow k m) a)) (* (* (fma k (+ k 10) 1) (fma k (+ k 10) 1)) (fma k (+ k 10) 1)))
442.0ms
(- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3))))
212.0ms
(- (+ a (* (log k) (* m a))) (* 10 (* a k)))
190.0ms
(- (+ (* 99 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 4))) (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 2))) (* 10 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 3))))
134.0ms
(* a (exp (* m (- (log -1) (log (/ -1 k))))))

prune479.0ms

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 0.1b

localize1.0ms

Local error

Found 2 expressions with local error:

2.8b
(/ (* (pow k m) a) (fma k (+ k 10) 1))
0.0b
(* (pow k m) a)

rewrite13.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
add-exp-log
add-cbrt-cube
associate-*l*
times-frac
associate-/r*
add-log-exp
log1p-expm1-u
pow1
expm1-log1p-u
associate-/l*
div-inv
*-commutative
div-exp
frac-2neg
clear-num
cbrt-undiv
Counts
2 → 34
Calls
2 calls:
Slowest
9.0ms
(/ (* (pow k m) a) (fma k (+ k 10) 1))
3.0ms
(* (pow k m) a)

series179.0ms

Counts
2 → 6
Calls
2 calls:
Slowest
93.0ms
(* (pow k m) a)
86.0ms
(/ (* (pow k m) a) (fma k (+ k 10) 1))

simplify2.4s

Counts
20 → 40
Calls
20 calls:
Slowest
693.0ms
(/ (* (* (* (pow k m) a) (* (pow k m) a)) (* (pow k m) a)) (* (* (fma k (+ k 10) 1) (fma k (+ k 10) 1)) (fma k (+ k 10) 1)))
605.0ms
(- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3))))
322.0ms
(- (+ (* 99 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 4))) (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 2))) (* 10 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 3))))
259.0ms
(- (+ a (* (log k) (* m a))) (* 10 (* a k)))
147.0ms
(* a (exp (* m (- (log -1) (log (/ -1 k))))))

prune435.0ms

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 0.1b

localize17.0ms

Local error

Found 4 expressions with local error:

12.3b
(sqrt (fma k (+ k 10) 1))
12.3b
(sqrt (fma k (+ k 10) 1))
0.1b
(/ (/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1))) (sqrt (fma k (+ k 10) 1)))
0.0b
(/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1)))

rewrite39.0ms

Algorithm
rewrite-expression-head
Rules
108×times-frac
66×sqrt-prod
55×add-sqr-sqrt
53×add-cube-cbrt
53×*-un-lft-identity
12×associate-/r*
11×associate-/l*
11×add-exp-log
11×add-cbrt-cube
div-inv
add-log-exp
log1p-expm1-u
div-exp
cbrt-undiv
pow1
expm1-log1p-u
pow1/2
frac-2neg
clear-num
rem-sqrt-square
associate-/l/
Counts
4 → 148
Calls
4 calls:
Slowest
17.0ms
(/ (/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1))) (sqrt (fma k (+ k 10) 1)))
15.0ms
(/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1)))
2.0ms
(sqrt (fma k (+ k 10) 1))
2.0ms
(sqrt (fma k (+ k 10) 1))

series251.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
116.0ms
(/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1)))
103.0ms
(/ (/ (* (pow k m) a) (sqrt (fma k (+ k 10) 1))) (sqrt (fma k (+ k 10) 1)))
17.0ms
(sqrt (fma k (+ k 10) 1))
15.0ms
(sqrt (fma k (+ k 10) 1))

simplify10.3s

Counts
181 → 160
Calls
181 calls:
Slowest
685.0ms
(- (+ (* 37 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) k)) (* 5 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))))
661.0ms
(/ (* (* (* (pow k m) a) (* (pow k m) a)) (* (pow k m) a)) (* (* (sqrt (fma k (+ k 10) 1)) (sqrt (fma k (+ k 10) 1))) (sqrt (fma k (+ k 10) 1))))
628.0ms
(- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3))))
338.0ms
(/ 1 (sqrt (sqrt (fma k (+ k 10) 1))))
338.0ms
(- (* 5 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 2))) (+ (* 37 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 3))) (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) k)))

prune2.3s

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 0.1b

localize13.0ms

Local error

Found 4 expressions with local error:

12.3b
(sqrt (fma k (+ k 10) 1))
12.3b
(sqrt (fma k (+ k 10) 1))
12.3b
(sqrt (fma k (+ k 10) 1))
0.1b
(/ a (sqrt (sqrt (fma k (+ k 10) 1))))

rewrite4.0ms

Algorithm
rewrite-expression-head
Rules
18×sqrt-prod
13×add-sqr-sqrt
10×add-cube-cbrt
10×*-un-lft-identity
associate-/r*
add-log-exp
log1p-expm1-u
add-exp-log
add-cbrt-cube
pow1
expm1-log1p-u
pow1/2
rem-sqrt-square
div-inv
frac-2neg
clear-num
Counts
4 → 63
Calls
4 calls:
Slowest
1.0ms
(/ a (sqrt (sqrt (fma k (+ k 10) 1))))
1.0ms
(sqrt (fma k (+ k 10) 1))
1.0ms
(sqrt (fma k (+ k 10) 1))
1.0ms
(sqrt (fma k (+ k 10) 1))

series183.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
119.0ms
(/ a (sqrt (sqrt (fma k (+ k 10) 1))))
25.0ms
(sqrt (fma k (+ k 10) 1))
24.0ms
(sqrt (fma k (+ k 10) 1))
15.0ms
(sqrt (fma k (+ k 10) 1))

simplify2.0s

Counts
33 → 75
Calls
33 calls:
Slowest
337.0ms
(- (+ (* 123/8 (* a (pow k 2))) a) (* 5/2 (* a k)))
266.0ms
(- (+ (* 123/8 (* a (sqrt (/ 1 (pow k 5))))) (* a (sqrt (/ 1 k)))) (* 5/2 (* a (sqrt (/ 1 (pow k 3))))))
235.0ms
(- (+ (* 123/8 (/ (* a (sqrt (/ -1 k))) (pow k 2))) (* a (sqrt (/ -1 k)))) (* 5/2 (/ (* a (sqrt (/ -1 k))) k)))
105.0ms
(- (+ 5 k) (* 12 (/ 1 k)))
96.0ms
(- (+ 1 (* 5 k)) (* 12 (pow k 2)))

prune1.2s

Pruning

5 alts after pruning (3 fresh and 2 done)

Merged error: 0.1b

regimes239.0ms

Accuracy

94.1% (0.1b remaining)

Error of 0.3b against oracle of 0.2b and baseline of 2.2b

bsearch109.0ms

end0.0ms

sample3.5s

Algorithm
intervals