Average Error: 28.1 → 28.1
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}\]
\[\frac{t + \left(y \cdot \left(y \cdot \left(z + x \cdot y\right) + 27464.7644705\right) + 230661.510616\right) \cdot y}{(\left((\left(y \cdot y\right) \cdot \left(y + a\right) + \left((y \cdot b + c)_*\right))_*\right) \cdot y + i)_*}\]

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.1

    \[\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.1

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

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

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

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

Reproduce

herbie shell --seed 2019002 +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.0m)Debug log

sample98.0ms

Algorithm
intervals

simplify172.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune48.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 28.6b

localize76.0ms

Local error

Found 4 expressions with local error:

28.5b
(/ (+ (* (+ (* (+ (* (+ (* 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)
0.1b
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
0.1b
(* (+ (* (+ y a) y) b) y)

rewrite119.0ms

Algorithm
rewrite-expression-head
Counts
4 → 78
Calls

4 calls. Slowest were:

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

series437.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

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

simplify10.9s

Counts
53 → 90
Calls

53 calls. Slowest were:

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

prune2.1s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 27.9b

localize41.0ms

Local error

Found 4 expressions with local error:

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

rewrite41.0ms

Algorithm
rewrite-expression-head
Counts
4 → 74
Calls

4 calls. Slowest were:

26.0ms
(/ (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1) (fma (fma (* y y) (+ a y) (fma y b c)) y i))
10.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
4.0ms
(* (+ (* x y) z) y)

series434.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

287.0ms
(/ (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1) (fma (fma (* y y) (+ a y) (fma y b c)) y i))
57.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
52.0ms
(fma (* y y) (+ a y) (fma y b c))
38.0ms
(* (+ (* x y) z) y)

simplify8.9s

Counts
53 → 86
Calls

53 calls. Slowest were:

493.0ms
(+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)
480.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1)
440.0ms
(/ (sqrt (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1)) 1)

prune2.1s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 27.9b

localize14.0ms

Local error

Found 4 expressions with local error:

28.7b
(/ (+ (* (+ (* (+ (* (+ 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
(* (+ (* x y) z) y)
0.1b
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)

rewrite98.0ms

Algorithm
rewrite-expression-head
Counts
4 → 95
Calls

4 calls. Slowest were:

49.0ms
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
27.0ms
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
12.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)

series959.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

490.0ms
(/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))
385.0ms
(/ 1 (/ (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i) (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
45.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
39.0ms
(* (+ (* x y) z) y)

simplify14.0s

Counts
66 → 107
Calls

66 calls. Slowest were:

592.0ms
(- (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)
581.0ms
(- (log (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (log (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
564.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)))

prune2.4s

Pruning

16 alts after pruning (14 fresh and 2 done)

Merged error: 27.9b

localize28.0ms

Local error

Found 4 expressions with local error:

28.5b
(/ (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1) (fma (* (* (cbrt (fma (* y y) (+ a y) (fma y b c))) (cbrt (fma (* y y) (+ a y) (fma y b c)))) (cbrt (fma (* y y) (+ a y) (fma y b c)))) y i))
23.3b
(cbrt (fma (* y y) (+ a y) (fma y b c)))
23.3b
(cbrt (fma (* y y) (+ a y) (fma y b c)))
23.3b
(cbrt (fma (* y y) (+ a y) (fma y b c)))

rewrite26.0ms

Algorithm
rewrite-expression-head
Counts
4 → 73
Calls

4 calls. Slowest were:

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

series694.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

309.0ms
(/ (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) 1) (fma (* (* (cbrt (fma (* y y) (+ a y) (fma y b c))) (cbrt (fma (* y y) (+ a y) (fma y b c)))) (cbrt (fma (* y y) (+ a y) (fma y b c)))) y i))
146.0ms
(cbrt (fma (* y y) (+ a y) (fma y b c)))
120.0ms
(cbrt (fma (* y y) (+ a y) (fma y b c)))
118.0ms
(cbrt (fma (* y y) (+ a y) (fma y b c)))

simplify11.6s

Counts
54 → 85
Calls

54 calls. Slowest were:

1.2s
(exp (* 1/3 (- (log -1) (* 3 (log (/ -1 y))))))
1.2s
(exp (* 1/3 (- (log -1) (* 3 (log (/ -1 y))))))
1.1s
(exp (* 1/3 (- (log -1) (* 3 (log (/ -1 y))))))

prune2.2s

Pruning

16 alts after pruning (14 fresh and 2 done)

Merged error: 27.9b

regimes1.8s

Accuracy

0% (0.7b remaining)

Error of 28.1b against oracle of 27.4b and baseline of 28.1b

bsearch3.0ms

end0.0ms

sample3.0s

Algorithm
intervals