Average Error: 28.7 → 28.8
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(t + \left(y \cdot \left(y \cdot \left(z + x \cdot y\right) + 27464.7644705\right) + 230661.510616\right) \cdot y\right) \cdot \frac{1}{i + y \cdot \left(\left(\left(a + y\right) \cdot y + b\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

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

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 div-inv28.8

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

  4. Final simplification28.8

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

Reproduce

herbie shell --seed 2019016 
(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

sample136.0ms

Algorithm
intervals
Results
264×(pre true 80)
264×(body real 80)

simplify165.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
165.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

prune26.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 28.0b

localize52.0ms

Local error

Found 4 expressions with local error:

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

rewrite108.0ms

Algorithm
rewrite-expression-head
Rules
15×add-cube-cbrt
15×*-un-lft-identity
15×add-sqr-sqrt
9×associate-*l*
9×times-frac
6×add-exp-log
6×associate-*l/
6×add-cbrt-cube
5×flip-+
5×flip3-+
4×add-log-exp
4×pow1
3×associate-/l*
3×*-commutative
3×associate-/r*
2×associate-/l/
2×associate-/r/
1×div-inv
1×div-exp
1×frac-2neg
1×clear-num
1×cbrt-undiv
Counts
4 → 70
Calls
4 calls:
Slowest
64.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
19.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
15.0ms
(* (+ (* (+ y a) y) b) y)
8.0ms
(* (+ (* x y) z) y)

series443.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
347.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
42.0ms
(* (+ (* (+ (* x y) z) y) 27464.7644705) y)
32.0ms
(* (+ (* (+ y a) y) b) y)
22.0ms
(* (+ (* x y) z) y)

simplify8.5s

Counts
53 → 82
Calls
53 calls:
Slowest
514.0ms
(/ (sqrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)) (sqrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
430.0ms
(+ (pow (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) 3) (pow t 3))
389.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)))
367.0ms
(* (+ (pow (* (+ (* x y) z) y) 3) (pow 27464.7644705 3)) y)
349.0ms
(/ (sqrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)) 1)

prune1.6s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 27.5b

localize17.0ms

Local error

Found 4 expressions with local error:

27.6b
(* (+ (* (+ (* (+ (* (+ (* 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)
0.1b
(* (+ (* (+ y a) y) b) y)

rewrite141.0ms

Algorithm
rewrite-expression-head
Rules
9×add-cube-cbrt
9×add-exp-log
9×associate-*l*
9×*-un-lft-identity
9×add-sqr-sqrt
7×pow1
6×flip-+
6×associate-*l/
6×associate-*r*
6×add-cbrt-cube
6×flip3-+
4×add-log-exp
4×associate-/r/
3×*-commutative
3×associate-/r*
2×div-inv
2×prod-exp
2×rec-exp
2×frac-times
1×inv-pow
1×pow-flip
1×cbrt-unprod
1×un-div-inv
1×associate-*r/
1×pow-prod-down
1×frac-2neg
1×clear-num
Counts
4 → 71
Calls
4 calls:
Slowest
88.0ms
(* (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
27.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
15.0ms
(* (+ (* (+ y a) y) b) y)
8.0ms
(* (+ (* x y) z) y)

series422.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
321.0ms
(* (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
43.0ms
(/ 1 (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
37.0ms
(* (+ (* (+ y a) y) b) y)
21.0ms
(* (+ (* x y) z) y)

simplify7.4s

Counts
46 → 83
Calls
46 calls:
Slowest
613.0ms
(* (+ (pow (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) 3) (pow t 3)) 1)
444.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)))
420.0ms
(* (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)) (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)))
398.0ms
(* (- (* (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)) (* t t)) 1)
371.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)))

prune1.8s

Pruning

10 alts after pruning (9 fresh and 1 done)

Merged error: 27.5b

localize27.0ms

Local error

Found 4 expressions with local error:

28.0b
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
28.0b
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
28.0b
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
2.7b
(cbrt (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))

rewrite198.0ms

Algorithm
rewrite-expression-head
Rules
36×times-frac
35×add-cube-cbrt
35×*-un-lft-identity
35×add-sqr-sqrt
15×cbrt-prod
10×add-exp-log
10×add-cbrt-cube
9×associate-/l*
9×associate-/r*
8×associate-/r/
7×flip-+
7×flip3-+
6×associate-/l/
4×add-log-exp
4×div-inv
4×pow1
3×div-exp
3×frac-2neg
3×clear-num
3×cbrt-undiv
1×pow1/3
1×cbrt-div
Counts
4 → 117
Calls
4 calls:
Slowest
72.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
62.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
36.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
25.0ms
(cbrt (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))

series1.6s

Counts
4 → 12
Calls
4 calls:
Slowest
684.0ms
(cbrt (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
313.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
311.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
284.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

simplify21.7s

Counts
115 → 129
Calls
115 calls:
Slowest
668.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)))
439.0ms
(+ (pow (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) 3) (pow t 3))
385.0ms
(- (* (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)) (* t t))
380.0ms
(cbrt (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (pow (* (+ (* (+ (* (+ y a) y) b) y) c) y) 3) (pow i 3))))
369.0ms
(/ (* (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t)) (cbrt (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t))) 1)

prune2.8s

Pruning

10 alts after pruning (8 fresh and 2 done)

Merged error: 27.5b

localize29.0ms

Local error

Found 4 expressions with local error:

24.6b
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (* (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))))
18.5b
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
18.5b
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
18.5b
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

rewrite160.0ms

Algorithm
rewrite-expression-head
Rules
18×cbrt-div
10×flip-+
10×flip3-+
9×add-cube-cbrt
9×add-exp-log
9×*-un-lft-identity
9×add-sqr-sqrt
9×cbrt-prod
8×associate-/r/
7×add-cbrt-cube
4×add-log-exp
4×pow1
4×frac-times
3×associate-/l*
3×pow1/3
3×times-frac
2×associate-/l/
2×associate-*r/
2×associate-*l/
2×div-exp
2×cbrt-undiv
1×div-inv
1×cbrt-unprod
1×prod-exp
1×frac-2neg
1×associate-/r*
1×clear-num
Counts
4 → 70
Calls
4 calls:
Slowest
108.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (* (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))))
17.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
17.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
16.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

series847.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
502.0ms
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (* (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)) (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))))
123.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
112.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))
110.0ms
(cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))

simplify9.2s

Counts
50 → 82
Calls
50 calls:
Slowest
609.0ms
(/ 1 (cbrt (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
430.0ms
(* (cbrt (+ (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) (* (+ (* (+ (* (+ y a) y) b) y) c) y)) (- (* i i) (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))) (cbrt (+ (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) (* (+ (* (+ (* (+ y a) y) b) y) c) y)) (- (* i i) (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))))
426.0ms
(cbrt (+ (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) (* (+ (* (+ (* (+ y a) y) b) y) c) y)) (- (* i i) (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) i))))
406.0ms
(cbrt (+ (pow (* (+ (* (+ (* (+ y a) y) b) y) c) y) 3) (pow i 3)))
403.0ms
(* (cbrt (+ (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) (* (+ (* (+ (* (+ y a) y) b) y) c) y)) (- (* i i) (* (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))) (cbrt (- (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))

prune1.8s

Pruning

11 alts after pruning (8 fresh and 3 done)

Merged error: 27.0b

regimes1.2s

Accuracy

0% (1.2b remaining)

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

bsearch4.0ms

end0.0ms

sample2.7s

Algorithm
intervals
Results
8252×(pre true 80)
8251×(body real 80)
1×(body real 160)