Average Error: 28.2 → 28.3
Time: 1.1m
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}\]
\[(\left(y \cdot y\right) \cdot \left((y \cdot \left((y \cdot x + z)_*\right) + 27464.7644705)_*\right) + \left((y \cdot 230661.510616 + t)_*\right))_* \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.2

    \[\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 *-un-lft-identity28.2

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

    \[\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}}}\]
  5. Using strategy rm
  6. Applied associate-/r/28.3

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

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

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

Reproduce

herbie shell --seed 2018362 +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: 58.6s)Debug log

start81.0ms

Algorithm
intervals

setup238.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 27.7b

localize55.0ms

Local error

Found 4 expressions with local error:

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

rewrite98.0ms

Algorithm
rewrite-expression-head
Counts
4 → 78
Calls

4 calls. Slowest were:

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

series478.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

361.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
59.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
36.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
22.0ms
(* (+ (* x y) z) y)

simplify11.7s

Counts
53 → 90
Calls

53 calls. Slowest were:

624.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
591.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)))
488.0ms
(+ (pow (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) 3) (pow t 3))

prune2.3s

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 27.2b

localize10.0ms

Local error

Found 4 expressions with local error:

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

rewrite150.0ms

Algorithm
rewrite-expression-head
Counts
4 → 79
Calls

4 calls. Slowest were:

68.0ms
(* (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
31.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
29.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)

series574.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

395.0ms
(* (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
75.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
58.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
46.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

simplify11.1s

Counts
46 → 91
Calls

46 calls. Slowest were:

737.0ms
(* (+ (pow (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) 3) (pow t 3)) 1)
600.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)))
598.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)))

prune2.7s

Pruning

12 alts after pruning (11 fresh and 1 done)

Merged error: 27.2b

localize8.0ms

Local error

Found 4 expressions with local error:

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

rewrite79.0ms

Algorithm
rewrite-expression-head
Counts
4 → 95
Calls

4 calls. Slowest were:

38.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:

461.0ms
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
430.0ms
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
66.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
46.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)

simplify14.8s

Counts
66 → 107
Calls

66 calls. Slowest were:

546.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)))
544.0ms
(* (+ (pow (* (+ (* (+ y a) y) b) y) 3) (pow c 3)) y)
521.0ms
(* (- (* (* (+ (* (+ y a) y) b) y) (* (+ (* (+ y a) y) b) y)) (* c c)) y)

prune2.9s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 27.2b

localize39.0ms

Local error

Found 4 expressions with local error:

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

rewrite77.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

28.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
27.0ms
(* (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (fma (* y y) (fma y (fma y x z) 27464.7644705) (fma y 230661.510616 t)))
21.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)

series530.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

375.0ms
(* (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (fma (* y y) (fma y (fma y x z) 27464.7644705) (fma y 230661.510616 t)))
61.0ms
(fma (* y y) (fma y (fma y x z) 27464.7644705) (fma y 230661.510616 t))
49.0ms
(* (+ (* (+ (* (+ y a) y) b) y) c) y)
44.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

simplify6.7s

Counts
36 → 83
Calls

36 calls. Slowest were:

665.0ms
(sqrt (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
484.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)))
480.0ms
(* (- (* (* (+ (* (+ y a) y) b) y) (* (+ (* (+ y a) y) b) y)) (* c c)) y)

prune1.7s

Pruning

10 alts after pruning (9 fresh and 1 done)

Merged error: 27.2b

regimes1.2s

Accuracy

0% (0.8b remaining)

Error of 28.3b against oracle of 27.5b and baseline of 28.3b

bsearch3.0ms