Average Error: 28.7 → 28.8
Time: 1.2m
Precision: 64
Internal Precision: 128
\[\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\]
\[(y \cdot \left((\left(y \cdot y\right) \cdot \left((y \cdot x + z)_*\right) + \left((27464.7644705 \cdot y + 230661.510616)_*\right))_*\right) + t)_* \cdot \frac{1}{i + y \cdot \left(\left(b + \left(y + a\right) \cdot y\right) \cdot y + c\right)}\]

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Bits error versus b

Bits error versus c

Bits error versus i

Derivation

  1. Initial program 28.7

    \[\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\]
  2. Using strategy rm
  3. Applied clear-num28.9

    \[\leadsto \color{blue}{\frac{1}{\frac{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}}}\]
  4. Using strategy rm
  5. Applied div-inv29.0

    \[\leadsto \frac{1}{\color{blue}{\left(\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i\right) \cdot \frac{1}{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}}}\]
  6. Applied associate-/r*28.8

    \[\leadsto \color{blue}{\frac{\frac{1}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}}{\frac{1}{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}}}\]
  7. Simplified28.8

    \[\leadsto \frac{\frac{1}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}}{\color{blue}{\frac{1}{(y \cdot \left((\left(y \cdot y\right) \cdot \left((y \cdot x + z)_*\right) + \left((27464.7644705 \cdot y + 230661.510616)_*\right))_*\right) + t)_*}}}\]
  8. Using strategy rm
  9. Applied associate-/r/28.8

    \[\leadsto \color{blue}{\frac{\frac{1}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}}{1} \cdot (y \cdot \left((\left(y \cdot y\right) \cdot \left((y \cdot x + z)_*\right) + \left((27464.7644705 \cdot y + 230661.510616)_*\right))_*\right) + t)_*}\]
  10. Final simplification28.8

    \[\leadsto (y \cdot \left((\left(y \cdot y\right) \cdot \left((y \cdot x + z)_*\right) + \left((27464.7644705 \cdot y + 230661.510616)_*\right))_*\right) + t)_* \cdot \frac{1}{i + y \cdot \left(\left(b + \left(y + a\right) \cdot y\right) \cdot y + c\right)}\]

Reproduce

herbie shell --seed 2019008 +o rules:numerics
(FPCore (x y z t a b c i)
  :name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
  (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))

Details

Time bar (total: 1.1m)Debug log

sample90.0ms

Algorithm
intervals

simplify187.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

187.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

prune42.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 28.3b

localize97.0ms

Local error

Found 4 expressions with local error:

28.3b
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
0.1b
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
0.1b
(* (+ (* (+ y a) y) b) y)
0.1b
(* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)

rewrite125.0ms

Algorithm
rewrite-expression-head
Counts
4 → 78
Calls

4 calls. Slowest were:

76.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
22.0ms
(* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)
18.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)

series501.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

360.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
63.0ms
(* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)
57.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
21.0ms
(* (+ (* (+ y a) y) b) y)

simplify12.2s

Counts
53 → 90
Calls

53 calls. Slowest were:

725.0ms
(* (- (* (* (+ (* (+ (* x y) z) y) 27464.7644705) y) (* (+ (* (+ (* x y) z) y) 27464.7644705) y)) (* 230661.510616 230661.510616)) y)
529.0ms
(* (+ (pow (* (+ (* (+ y a) y) b) y) 3) (pow c 3)) y)
520.0ms
(* (+ (pow (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 3) (pow 230661.510616 3)) y)

prune2.2s

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 27.7b

localize73.0ms

Local error

Found 4 expressions with local error:

28.3b
(/ (+ (* (* (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616))) (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))) (* (cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616)) y)) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
21.8b
(cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))
21.8b
(cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))
21.8b
(cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616))

rewrite100.0ms

Algorithm
rewrite-expression-head
Counts
4 → 74
Calls

4 calls. Slowest were:

92.0ms
(/ (+ (* (* (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616))) (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))) (* (cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616)) y)) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
6.0ms
(cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616))
1.0ms
(cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))

series1.0s

Counts
4 → 12
Calls

4 calls. Slowest were:

571.0ms
(/ (+ (* (* (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616))) (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))) (* (cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616)) y)) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
175.0ms
(cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))
157.0ms
(cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))
142.0ms
(cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616))

simplify11.7s

Counts
52 → 86
Calls

52 calls. Slowest were:

514.0ms
(+ (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) (* (+ (* (+ (* (+ y a) y) b) y) c) y)) (- (* i i) (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
444.0ms
(- (* (* (* (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616))) (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))) (* (cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616)) y)) (* (* (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616))) (cbrt (fma (* y y) (fma y x z) (fma y 27464.7644705 230661.510616)))) (* (cbrt (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616)) y))) (* t t))
428.0ms
(- (+ (pow 230661.510616 1/3) (* 0.0396898531780084 (* y (pow 230661.510616 1/3)))) (* 0.0015752844452918634 (* (pow y 2) (pow 230661.510616 1/3))))

prune2.6s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 27.7b

localize13.0ms

Local error

Found 4 expressions with local error:

28.4b
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
0.4b
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
0.1b
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
0.1b
(* (+ (* (+ y a) y) b) y)

rewrite97.0ms

Algorithm
rewrite-expression-head
Counts
4 → 95
Calls

4 calls. Slowest were:

54.0ms
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
16.0ms
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
15.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)

series1.0s

Counts
4 → 12
Calls

4 calls. Slowest were:

489.0ms
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
455.0ms
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
56.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
37.0ms
(* (+ (* (+ y a) y) b) y)

simplify14.4s

Counts
66 → 107
Calls

66 calls. Slowest were:

742.0ms
(+ (pow (* (+ (* (+ (* (+ y a) y) b) y) c) y) 3) (pow i 3))
566.0ms
(+ (* (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)) (- (* t t) (* (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
484.0ms
(/ 1 (* (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)) (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))))

prune2.5s

Pruning

13 alts after pruning (11 fresh and 2 done)

Merged error: 27.7b

localize21.0ms

Local error

Found 4 expressions with local error:

28.2b
(/ (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t)))
0.3b
(fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616))
0.3b
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
0.1b
(/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t))

rewrite63.0ms

Algorithm
rewrite-expression-head
Counts
4 → 101
Calls

4 calls. Slowest were:

41.0ms
(/ (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t)))
20.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
1.0ms
(/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t))

series600.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

394.0ms
(/ (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t)))
86.0ms
(/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t))
71.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
48.0ms
(fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616))

simplify13.4s

Counts
84 → 113
Calls

84 calls. Slowest were:

497.0ms
(/ (/ 1 (+ (pow (* (+ (* (+ (* (+ y a) y) b) y) c) y) 3) (pow i 3))) (sqrt (/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t))))
467.0ms
(/ (- (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (/ 1 (fma y (fma (* y y) (fma y x z) (fma 27464.7644705 y 230661.510616)) t)))
416.0ms
(/ (/ 1 (+ (pow (* (+ (* (+ (* (+ y a) y) b) y) c) y) 3) (pow i 3))) 1)

prune1.8s

Pruning

12 alts after pruning (10 fresh and 2 done)

Merged error: 27.7b

regimes938.0ms

Accuracy

0% (0.9b remaining)

Error of 28.8b against oracle of 27.9b and baseline of 28.8b

bsearch3.0ms

end0.0ms

sample2.1s

Algorithm
intervals