Average Error: 3.7 → 1.7
Time: 34.3s
Precision: 64
Internal Precision: 128
\[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
\[\frac{x}{x + y \cdot e^{(\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{3.0}}{t}\right) \cdot \left(c - b\right) + \left(\sqrt{t + a} \cdot \frac{z}{t}\right))_* \cdot 2.0}}\]

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

Derivation

  1. Initial program 3.7

    \[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
  2. Using strategy rm
  3. Applied div-inv3.7

    \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\color{blue}{\left(z \cdot \sqrt{t + a}\right) \cdot \frac{1}{t}} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
  4. Applied prod-diff21.9

    \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \color{blue}{\left((\left(z \cdot \sqrt{t + a}\right) \cdot \left(\frac{1}{t}\right) + \left(-\left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right) \cdot \left(b - c\right)\right))_* + (\left(-\left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right) \cdot \left(b - c\right) + \left(\left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right) \cdot \left(b - c\right)\right))_*\right)}}}\]
  5. Simplified20.1

    \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\color{blue}{(\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{3.0}}{t}\right) \cdot \left(c - b\right) + \left(\sqrt{t + a} \cdot \frac{z}{t}\right))_*} + (\left(-\left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right) \cdot \left(b - c\right) + \left(\left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right) \cdot \left(b - c\right)\right))_*\right)}}\]
  6. Simplified1.7

    \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left((\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{3.0}}{t}\right) \cdot \left(c - b\right) + \left(\sqrt{t + a} \cdot \frac{z}{t}\right))_* + \color{blue}{0}\right)}}\]
  7. Final simplification1.7

    \[\leadsto \frac{x}{x + y \cdot e^{(\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{3.0}}{t}\right) \cdot \left(c - b\right) + \left(\sqrt{t + a} \cdot \frac{z}{t}\right))_* \cdot 2.0}}\]

Reproduce

herbie shell --seed 2019005 +o rules:numerics
(FPCore (x y z t a b c)
  :name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
  (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))

Details

Time bar (total: 33.8s)Debug log

sample228.0ms

Algorithm
intervals

simplify244.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

244.0ms
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))

prune13.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 1.7b

localize110.0ms

Local error

Found 4 expressions with local error:

3.2b
(/ (* z (sqrt (+ t a))) t)
0.9b
(- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))
0.2b
(/ 2.0 (* t 3.0))
0.2b
(* z (sqrt (+ t a)))

rewrite25.0ms

Algorithm
rewrite-expression-head
Counts
4 → 112
Calls

4 calls. Slowest were:

13.0ms
(- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))
8.0ms
(/ (* z (sqrt (+ t a))) t)
2.0ms
(* z (sqrt (+ t a)))

series312.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

189.0ms
(- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))
56.0ms
(* z (sqrt (+ t a)))
53.0ms
(/ (* z (sqrt (+ t a))) t)
15.0ms
(/ 2.0 (* t 3.0))

simplify16.7s

Counts
92 → 124
Calls

92 calls. Slowest were:

631.0ms
(* t (* (+ b c) (* (- a (/ 5.0 6.0)) (* t 3.0))))
611.0ms
(* (fma (- (sqrt (/ 2.0 (* t 3.0)))) (sqrt (/ 2.0 (* t 3.0))) (* (sqrt (/ 2.0 (* t 3.0))) (sqrt (/ 2.0 (* t 3.0))))) (- b c))
576.0ms
(- (* (* z (sqrt (+ t a))) (* (+ (* b b) (+ (* c c) (* b c))) (+ (* (+ a (/ 5.0 6.0)) (+ a (/ 5.0 6.0))) (+ (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (* (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))) (* t (* (- (pow b 3) (pow c 3)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3)))))

prune3.4s

Pruning

2 alts after pruning (1 fresh and 1 done)

Merged error: 0.0b

localize62.0ms

Local error

Found 4 expressions with local error:

5.0b
(* (sqrt (+ t a)) (/ z t))
0.3b
(/ (/ 2.0 3.0) t)
0.1b
(fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t)))
0.0b
(* y (exp (* 2.0 (+ (fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t))) 0))))

rewrite12.0ms

Algorithm
rewrite-expression-head
Counts
4 → 70
Calls

4 calls. Slowest were:

5.0ms
(* (sqrt (+ t a)) (/ z t))
4.0ms
(* y (exp (* 2.0 (+ (fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t))) 0))))
3.0ms
(/ (/ 2.0 3.0) t)

series528.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

335.0ms
(* y (exp (* 2.0 (+ (fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t))) 0))))
129.0ms
(fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t)))
57.0ms
(* (sqrt (+ t a)) (/ z t))
8.0ms
(/ (/ 2.0 3.0) t)

simplify2.1s

Counts
37 → 82
Calls

37 calls. Slowest were:

289.0ms
(* (* (* (sqrt (+ t a)) (sqrt (+ t a))) (sqrt (+ t a))) (* (* (/ z t) (/ z t)) (/ z t)))
217.0ms
(sqrt (exp (* 2.0 (+ (fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t))) 0))))
203.0ms
(cbrt (exp (* 2.0 (+ (fma (- (+ (/ 5.0 6.0) a) (/ (/ 2.0 3.0) t)) (- c b) (* (sqrt (+ t a)) (/ z t))) 0))))

prune1.8s

Pruning

2 alts after pruning (0 fresh and 2 done)

Merged error: 0.0b

regimes1.1s

Accuracy

0% (0.4b remaining)

Error of 1.7b against oracle of 1.3b and baseline of 1.7b

bsearch9.0ms

end0.0ms

sample7.2s

Algorithm
intervals