Average Error: 2.0 → 2.2
Time: 46.0s
Precision: 64
Internal Precision: 128
\[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]
\[\begin{array}{l} \mathbf{if}\;y \le -4.3882314950956804 \cdot 10^{-250} \lor \neg \left(y \le -3.9719064524468603 \cdot 10^{-283}\right):\\ \;\;\;\;\frac{x \cdot {e}^{\left(\left(\log a \cdot \left(t - 1.0\right) + y \cdot \log z\right) - b\right)}}{y}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{\frac{\frac{y}{\left(e^{-b} \cdot {z}^{y}\right) \cdot \left({a}^{\left(-1.0\right)} \cdot {a}^{t}\right)}}{x}}\\ \end{array}\]

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

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if y < -4.3882314950956804e-250 or -3.9719064524468603e-283 < y

    1. Initial program 1.9

      \[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]
    2. Using strategy rm
    3. Applied *-un-lft-identity1.9

      \[\leadsto \frac{x \cdot e^{\color{blue}{1 \cdot \left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}}{y}\]
    4. Applied exp-prod2.0

      \[\leadsto \frac{x \cdot \color{blue}{{\left(e^{1}\right)}^{\left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}}{y}\]
    5. Simplified2.0

      \[\leadsto \frac{x \cdot {\color{blue}{e}}^{\left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}{y}\]

    if -4.3882314950956804e-250 < y < -3.9719064524468603e-283

    1. Initial program 5.0

      \[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]
    2. Using strategy rm
    3. Applied *-un-lft-identity5.0

      \[\leadsto \frac{x \cdot e^{\color{blue}{1 \cdot \left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}}{y}\]
    4. Applied exp-prod5.1

      \[\leadsto \frac{x \cdot \color{blue}{{\left(e^{1}\right)}^{\left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}}{y}\]
    5. Simplified5.1

      \[\leadsto \frac{x \cdot {\color{blue}{e}}^{\left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}}{y}\]
    6. Taylor expanded around inf 5.0

      \[\leadsto \frac{x \cdot \color{blue}{e^{1.0 \cdot \log \left(\frac{1}{a}\right) - \left(\log \left(\frac{1}{z}\right) \cdot y + \left(b + t \cdot \log \left(\frac{1}{a}\right)\right)\right)}}}{y}\]
    7. Simplified10.2

      \[\leadsto \frac{x \cdot \color{blue}{\left(\left({z}^{y} \cdot e^{-b}\right) \cdot \left({a}^{t} \cdot {a}^{\left(-1.0\right)}\right)\right)}}{y}\]
    8. Using strategy rm
    9. Applied associate-/l*10.1

      \[\leadsto \color{blue}{\frac{x}{\frac{y}{\left({z}^{y} \cdot e^{-b}\right) \cdot \left({a}^{t} \cdot {a}^{\left(-1.0\right)}\right)}}}\]
    10. Using strategy rm
    11. Applied clear-num10.1

      \[\leadsto \color{blue}{\frac{1}{\frac{\frac{y}{\left({z}^{y} \cdot e^{-b}\right) \cdot \left({a}^{t} \cdot {a}^{\left(-1.0\right)}\right)}}{x}}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification2.2

    \[\leadsto \begin{array}{l} \mathbf{if}\;y \le -4.3882314950956804 \cdot 10^{-250} \lor \neg \left(y \le -3.9719064524468603 \cdot 10^{-283}\right):\\ \;\;\;\;\frac{x \cdot {e}^{\left(\left(\log a \cdot \left(t - 1.0\right) + y \cdot \log z\right) - b\right)}}{y}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{\frac{\frac{y}{\left(e^{-b} \cdot {z}^{y}\right) \cdot \left({a}^{\left(-1.0\right)} \cdot {a}^{t}\right)}}{x}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019005 
(FPCore (x y z t a b)
  :name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
  (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))

Details

Time bar (total: 44.7s)Debug log

sample679.0ms

Algorithm
intervals

simplify340.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

339.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)

prune35.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.3b

localize33.0ms

Local error

Found 4 expressions with local error:

2.4b
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
0.9b
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
0.2b
(* y (log z))
0.1b
(* (- t 1.0) (log a))

rewrite51.0ms

Algorithm
rewrite-expression-head
Counts
4 → 57
Calls

4 calls. Slowest were:

32.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
11.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
5.0ms
(* (- t 1.0) (log a))

series439.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

175.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
145.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
67.0ms
(* (- t 1.0) (log a))
52.0ms
(* y (log z))

simplify2.7s

Counts
34 → 69
Calls

34 calls. Slowest were:

349.0ms
(+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y))))
220.0ms
(- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1))))
197.0ms
(* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a)))

prune1.2s

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.0b

localize21.0ms

Local error

Found 4 expressions with local error:

2.4b
(/ (* x (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
0.9b
(pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
0.2b
(* y (log z))
0.1b
(* (- t 1.0) (log a))

rewrite25.0ms

Algorithm
rewrite-expression-head
Counts
4 → 57
Calls

4 calls. Slowest were:

13.0ms
(/ (* x (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
6.0ms
(pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
4.0ms
(* (- t 1.0) (log a))

series469.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

218.0ms
(/ (* x (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
137.0ms
(pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
58.0ms
(* (- t 1.0) (log a))
55.0ms
(* y (log z))

simplify2.5s

Counts
33 → 69
Calls

33 calls. Slowest were:

348.0ms
(+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y))))
198.0ms
(* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a)))
192.0ms
(- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1))))

prune1.2s

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 0.0b

localize34.0ms

Local error

Found 4 expressions with local error:

2.4b
(/ (* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))) y)
0.4b
(* (pow z y) (exp (- b)))
0.0b
(* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
0.0b
(* (pow a t) (pow a (- 1.0)))

rewrite74.0ms

Algorithm
rewrite-expression-head
Counts
4 → 93
Calls

4 calls. Slowest were:

37.0ms
(* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
21.0ms
(/ (* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))) y)
7.0ms
(* (pow a t) (pow a (- 1.0)))

series1.2s

Counts
4 → 12
Calls

4 calls. Slowest were:

593.0ms
(/ (* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))) y)
347.0ms
(* x (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
195.0ms
(* (pow a t) (pow a (- 1.0)))
83.0ms
(* (pow z y) (exp (- b)))

simplify4.0s

Counts
68 → 105
Calls

68 calls. Slowest were:

331.0ms
(+ (* 1/2 (* (pow (/ 1 (pow a 1.0)) 1.0) (* (pow t 2) (pow (log a) 2)))) (+ (* (pow (/ 1 (pow a 1.0)) 1.0) (* t (log a))) (pow (/ 1 (pow a 1.0)) 1.0)))
321.0ms
(* x (* (* (pow z y) (exp (- b))) (pow a t)))
310.0ms
(* (pow (/ 1 (pow a 1.0)) 1.0) (* x (* (exp (* -1 (* (log (/ 1 z)) y))) (* (exp (* -1 (* t (log (/ 1 a))))) (exp (- b))))))

prune1.6s

Pruning

6 alts after pruning (4 fresh and 2 done)

Merged error: 0b

localize12.0ms

Local error

Found 4 expressions with local error:

1.0b
(/ x (/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))))
0.4b
(* (pow z y) (exp (- b)))
0.0b
(/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
0.0b
(* (pow a t) (pow a (- 1.0)))

rewrite32.0ms

Algorithm
rewrite-expression-head
Counts
4 → 99
Calls

4 calls. Slowest were:

14.0ms
(/ x (/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))))
6.0ms
(/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
6.0ms
(* (pow z y) (exp (- b)))

series1.2s

Counts
4 → 12
Calls

4 calls. Slowest were:

513.0ms
(/ x (/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0))))))
409.0ms
(/ y (* (* (pow z y) (exp (- b))) (* (pow a t) (pow a (- 1.0)))))
198.0ms
(* (pow a t) (pow a (- 1.0)))
78.0ms
(* (pow z y) (exp (- b)))

simplify3.8s

Counts
73 → 111
Calls

73 calls. Slowest were:

373.0ms
(- (* a y) (+ (* a (* (log z) (pow y 2))) (* a (* t (* y (log a))))))
329.0ms
(* (* (* (pow z y) (pow z y)) (pow z y)) (* (* (exp (- b)) (exp (- b))) (exp (- b))))
294.0ms
(- (+ (* (pow (/ 1 (pow a 1.0)) 1.0) (* x (log z))) (* (pow (/ 1 (pow a 1.0)) 1.0) (/ x y))) (* (pow (/ 1 (pow a 1.0)) 1.0) (/ (* x b) y)))

prune1.9s

Pruning

5 alts after pruning (3 fresh and 2 done)

Merged error: 0b

regimes953.0ms

Accuracy

-6.6% (2.2b remaining)

Error of 2.2b against oracle of 0.0b and baseline of 2.1b

bsearch2.3s

end0.0ms

sample18.0s

Algorithm
intervals