Average Error: 1.9 → 2.8
Time: 45.2s
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}\;t \le 8.105866356025745 \cdot 10^{-274} \lor \neg \left(t \le 3.6483414341945934 \cdot 10^{-83}\right):\\ \;\;\;\;\frac{x}{\frac{y}{e^{\left(\log z \cdot y + \log a \cdot \left(t - 1.0\right)\right) - b}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{{a}^{\left(t - 1.0\right)} \cdot \left({z}^{y} \cdot x\right)}{y \cdot e^{b}}\\ \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 t < 8.105866356025745e-274 or 3.6483414341945934e-83 < t

    1. Initial program 1.6

      \[\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 associate-/l*1.4

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

    if 8.105866356025745e-274 < t < 3.6483414341945934e-83

    1. Initial program 3.7

      \[\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 exp-diff10.4

      \[\leadsto \frac{x \cdot \color{blue}{\frac{e^{y \cdot \log z + \left(t - 1.0\right) \cdot \log a}}{e^{b}}}}{y}\]
    4. Applied associate-*r/11.5

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

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

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

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

Reproduce

herbie shell --seed 2019010 +o rules:numerics
(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.1s)Debug log

sample628.0ms

Algorithm
intervals

simplify365.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
365.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)

prune33.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.9b

localize45.0ms

Local error

Found 4 expressions with local error:

1.5b
(/ (* 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.3b
(* y (log z))
0.1b
(* (- t 1.0) (log a))

rewrite63.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
add-exp-log
associate-*r*
add-cbrt-cube
pow1
add-log-exp
log1p-expm1-u
expm1-log1p-u
exp-sum
associate-*l*
exp-prod
*-commutative
associate-*l/
exp-diff
associate-/l/
rem-exp-log
associate-/l*
flip--
div-inv
cbrt-unprod
associate-*r/
associate--l+
prod-exp
pow-prod-down
flip3--
frac-2neg
sub-neg
clear-num
Counts
4 → 65
Calls
4 calls:
Slowest
32.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
25.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
4.0ms
(* (- t 1.0) (log a))
1.0ms
(* y (log z))

series443.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
216.0ms
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y)
143.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
51.0ms
(* (- t 1.0) (log a))
32.0ms
(* y (log z))

simplify3.5s

Counts
34 → 77
Calls
34 calls:
Slowest
417.0ms
(- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1))))
335.0ms
(+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y))))
234.0ms
(* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a)))
219.0ms
(exp (* (cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) (cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))))
210.0ms
(- (* 1.0 (log (/ 1 a))) (* t (log (/ 1 a))))

prune1.3s

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.1b

localize13.0ms

Local error

Found 4 expressions with local error:

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

rewrite46.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
add-exp-log
associate-*r*
add-cbrt-cube
pow1
add-log-exp
associate-/r*
log1p-expm1-u
expm1-log1p-u
exp-sum
associate-*l*
exp-prod
div-inv
*-commutative
associate-/r/
associate-*l/
exp-diff
rem-exp-log
flip--
cbrt-unprod
associate--l+
prod-exp
pow-prod-down
flip3--
frac-2neg
sub-neg
clear-num
Counts
4 → 69
Calls
4 calls:
Slowest
35.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
5.0ms
(/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))))
4.0ms
(* (- t 1.0) (log a))
1.0ms
(* y (log z))

series413.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
197.0ms
(/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))))
111.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
55.0ms
(* (- t 1.0) (log a))
49.0ms
(* y (log z))

simplify4.0s

Counts
39 → 81
Calls
39 calls:
Slowest
324.0ms
(exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
301.0ms
(- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1))))
278.0ms
(+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y))))
251.0ms
(* (- (* t t) (* 1.0 1.0)) (log a))
204.0ms
(- (+ (* y (log z)) (* (- t 1.0) (log a))) b)

prune1.4s

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 0.1b

localize33.0ms

Local error

Found 4 expressions with local error:

1.4b
(/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))))
0.1b
(fma t (log a) (- b))
0.1b
(/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))
0.0b
(* (pow a (- 1.0)) (pow z y))

rewrite34.0ms

Algorithm
rewrite-expression-head
Rules
10×associate-*l/
add-exp-log
add-cube-cbrt
*-un-lft-identity
associate-/r*
add-sqr-sqrt
add-cbrt-cube
associate-/r/
add-log-exp
log1p-expm1-u
prod-exp
pow1
pow-to-exp
expm1-log1p-u
pow-sub
div-inv
associate-*l*
associate-*r*
pow-neg
neg-sub0
frac-2neg
clear-num
fma-udef
cbrt-unprod
*-commutative
Counts
4 → 67
Calls
4 calls:
Slowest
10.0ms
(/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))))
9.0ms
(/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))
8.0ms
(* (pow a (- 1.0)) (pow z y))
0.0ms
(fma t (log a) (- b))

series882.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
419.0ms
(/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))))
230.0ms
(/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))
170.0ms
(* (pow a (- 1.0)) (pow z y))
62.0ms
(fma t (log a) (- b))

simplify3.0s

Counts
34 → 79
Calls
34 calls:
Slowest
311.0ms
(+ (* (pow (/ 1 (pow a 1.0)) 1.0) (* (log z) y)) (+ (* 1/2 (* (pow (/ 1 (pow a 1.0)) 1.0) (* (pow (log z) 2) (pow y 2)))) (pow (/ 1 (pow a 1.0)) 1.0)))
296.0ms
(/ (* (exp (* -1 (* (log (/ 1 z)) y))) (* x (exp (- (+ b (* t (log (/ 1 a)))))))) (* a y))
221.0ms
(- (* t (log -1)) (+ (* t (log (/ -1 a))) b))
213.0ms
(* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))
213.0ms
(cbrt (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))))

prune1.5s

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 0b

localize18.0ms

Local error

Found 4 expressions with local error:

1.4b
(/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)))
0.1b
(fma t (log a) (- b))
0.0b
(* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))
0.0b
(/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))))

rewrite20.0ms

Algorithm
rewrite-expression-head
Rules
add-exp-log
add-cube-cbrt
add-cbrt-cube
*-un-lft-identity
add-sqr-sqrt
add-log-exp
log1p-expm1-u
associate-*l*
pow1
expm1-log1p-u
div-inv
associate-*r*
prod-exp
associate-*l/
frac-2neg
associate-/r*
clear-num
fma-udef
cbrt-unprod
*-commutative
associate-/r/
pow-to-exp
Counts
4 → 58
Calls
4 calls:
Slowest
14.0ms
(* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))
4.0ms
(/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))))
1.0ms
(/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)))
0.0ms
(fma t (log a) (- b))

series1.5s

Counts
4 → 12
Calls
4 calls:
Slowest
573.0ms
(/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))))
462.0ms
(/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)))
412.0ms
(* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))
81.0ms
(fma t (log a) (- b))

simplify3.6s

Counts
24 → 70
Calls
24 calls:
Slowest
749.0ms
(- (+ (* a (* y b)) (* a y)) (* a (* (log z) (pow y 2))))
346.0ms
(/ (* (exp (* -1 (* (log (/ 1 z)) y))) (* x (exp (- (+ b (* t (log (/ 1 a)))))))) (* a y))
300.0ms
(/ y (* (exp (* -1 (* (log (/ 1 z)) y))) (exp (- (+ b (* t (log (/ 1 a))))))))
282.0ms
(- y (+ (* t (* y (log a))) (* (log z) (pow y 2))))
248.0ms
(- (* t (log -1)) (+ (* t (log (/ -1 a))) b))

prune1.2s

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 0b

regimes964.0ms

Accuracy

-62.3% (2.8b remaining)

Error of 2.8b against oracle of 0.0b and baseline of 1.7b

bsearch1.9s

end0.0ms

sample17.3s

Algorithm
intervals