


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
Results
if t < 7.021435509051398e-308Initial program 5.1
rmApplied associate-/l*5.1
rmApplied add-log-exp5.1
if 7.021435509051398e-308 < t < 1.20691627027635e-103Initial program 5.5
rmApplied associate-/l*6.8
rmApplied flip-+10.6
Applied frac-sub10.6
Applied associate-*r/10.6
Applied frac-sub10.7
Simplified4.3
if 1.20691627027635e-103 < t Initial program 2.2
rmApplied associate-/l*0.2
rmApplied add-cbrt-cube0.2
rmApplied pow1/30.2
Final simplification2.3
herbie shell --seed 2019010
(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)))))))))))
Time bar (total: 1.2m)Debug log
| 288.0ms | (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 4.4b
Found 4 expressions with local error:
| 3.3b | (/ (* z (sqrt (+ t a))) t) |
| 3.2b | (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 0.3b | (/ 2.0 (* t 3.0)) |
| 0.2b | (* z (sqrt (+ t a))) |
| 20× | frac-sub |
| 9× | flip-- |
| 9× | flip3-- |
| 8× | associate-*r/ |
| 8× | frac-times |
| 6× | add-log-exp |
| 6× | associate--r+ |
| 6× | add-cube-cbrt |
| 6× | associate-*r* |
| 6× | *-un-lft-identity |
| 6× | add-sqr-sqrt |
| 5× | flip-+ |
| 5× | flip3-+ |
| 4× | add-exp-log |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 3× | distribute-lft-in |
| 3× | sub-neg |
| 3× | distribute-rgt-in |
| 2× | associate-/l/ |
| 2× | div-inv |
| 2× | associate--l+ |
| 2× | associate-*l/ |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | associate-/l* |
| 1× | *-commutative |
| 1× | diff-log |
| 1× | associate-/r* |
| 9.0ms | (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 4.0ms | (/ (* z (sqrt (+ t a))) t) |
| 2.0ms | (* z (sqrt (+ t a))) |
| 0.0ms | (/ 2.0 (* t 3.0)) |
| 169.0ms | (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 55.0ms | (/ (* z (sqrt (+ t a))) t) |
| 42.0ms | (* z (sqrt (+ t a))) |
| 13.0ms | (/ 2.0 (* t 3.0)) |
| 824.0ms | (* t (* (+ b c) (* (- a (/ 5.0 6.0)) (* t 3.0)))) |
| 506.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))))) |
| 425.0ms | (- (* (* z (sqrt (+ t a))) (* (+ 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 (* (- (* b b) (* c c)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3))))) |
| 315.0ms | (- (* (* z (sqrt (+ t a))) (+ (* (+ 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 (* (- b c) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3))))) |
| 313.0ms | (- (* (* z (sqrt (+ t a))) (* (+ (* b b) (+ (* c c) (* b c))) (* (+ (* a a) (- (* (/ 5.0 6.0) (/ 5.0 6.0)) (* a (/ 5.0 6.0)))) (* t 3.0)))) (* t (* (- (pow b 3) (pow c 3)) (- (* (+ (pow a 3) (pow (/ 5.0 6.0) 3)) (* t 3.0)) (* (+ (* a a) (- (* (/ 5.0 6.0) (/ 5.0 6.0)) (* a (/ 5.0 6.0)))) 2.0))))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 0.5b
Found 4 expressions with local error:
| 3.2b | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 2.0b | (/ z (/ t (sqrt (+ t a)))) |
| 0.4b | (/ t (sqrt (+ t a))) |
| 0.3b | (/ 2.0 (* t 3.0)) |
| 20× | frac-sub |
| 13× | associate-/r* |
| 9× | flip-- |
| 9× | flip3-- |
| 8× | frac-times |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | associate--r+ |
| 5× | flip-+ |
| 5× | associate-/r/ |
| 5× | flip3-+ |
| 4× | div-inv |
| 4× | add-exp-log |
| 4× | associate-*r/ |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 3× | distribute-lft-in |
| 3× | frac-2neg |
| 3× | sub-neg |
| 3× | clear-num |
| 3× | distribute-rgt-in |
| 2× | associate--l+ |
| 2× | associate-*l/ |
| 1× | diff-log |
| 16.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 1.0ms | (/ t (sqrt (+ t a))) |
| 1.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 1.0ms | (/ 2.0 (* t 3.0)) |
| 180.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 70.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 57.0ms | (/ t (sqrt (+ t a))) |
| 14.0ms | (/ 2.0 (* t 3.0)) |
| 512.0ms | (- (* z (* (+ (* 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 (sqrt (+ t a))) (* (- (pow b 3) (pow c 3)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3))))) |
| 442.0ms | (* (/ t (sqrt (+ t a))) (* (+ (* a a) (- (* (/ 5.0 6.0) (/ 5.0 6.0)) (* a (/ 5.0 6.0)))) (* t 3.0))) |
| 356.0ms | (* (/ t (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))))))) |
| 338.0ms | (- (* z (* (+ (* b b) (+ (* c c) (* b c))) (+ (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) (* (/ t (sqrt (+ t a))) (* (- (pow b 3) (pow c 3)) (- (* (+ a (/ 5.0 6.0)) (+ a (/ 5.0 6.0))) (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))))))) |
| 328.0ms | (- (* z (* (+ 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 (sqrt (+ t a))) (* (- (* b b) (* c c)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3))))) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0.5b
Found 4 expressions with local error:
| 3.2b | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 2.0b | (/ z (/ t (sqrt (+ t a)))) |
| 1.2b | (log (exp (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))) |
| 0.4b | (/ t (sqrt (+ t a))) |
| 20× | frac-sub |
| 12× | associate-/r* |
| 9× | flip-- |
| 9× | add-cube-cbrt |
| 9× | flip3-- |
| 9× | *-un-lft-identity |
| 9× | add-sqr-sqrt |
| 8× | frac-times |
| 6× | add-log-exp |
| 6× | associate--r+ |
| 5× | flip-+ |
| 5× | associate-/r/ |
| 5× | flip3-+ |
| 5× | log-prod |
| 5× | pow1 |
| 4× | log-pow |
| 4× | add-exp-log |
| 4× | associate-*r/ |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | sub-neg |
| 3× | div-inv |
| 3× | sqrt-prod |
| 3× | distribute-lft-in |
| 3× | exp-prod |
| 3× | distribute-rgt-in |
| 2× | associate--l+ |
| 2× | associate-*l/ |
| 2× | exp-sum |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | rem-log-exp |
| 1× | diff-log |
| 1× | exp-diff |
| 1× | log-div |
| 16.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 14.0ms | (log (exp (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))) |
| 1.0ms | (/ t (sqrt (+ t a))) |
| 1.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 172.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 167.0ms | (log (exp (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))) |
| 59.0ms | (/ t (sqrt (+ t a))) |
| 52.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 642.0ms | (* (/ t (sqrt (+ t a))) (* (+ (* b b) (+ (* c c) (* b c))) (+ (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) |
| 404.0ms | (* (/ t (sqrt (+ t a))) (* (+ (* a a) (- (* (/ 5.0 6.0) (/ 5.0 6.0)) (* a (/ 5.0 6.0)))) (* t 3.0))) |
| 375.0ms | (- (* z (* (+ (* 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 (sqrt (+ t a))) (* (- (pow b 3) (pow c 3)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (/ 2.0 (* t 3.0)) 3))))) |
| 355.0ms | (* (/ t (sqrt (+ t a))) (+ (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))) |
| 353.0ms | (- (* z (* (+ (* b b) (+ (* c c) (* b c))) (+ (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))) (* (/ t (sqrt (+ t a))) (* (- (pow b 3) (pow c 3)) (- (* (+ a (/ 5.0 6.0)) (+ a (/ 5.0 6.0))) (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))))))) |
5 alts after pruning (4 fresh and 1 done)
Merged error: 0.5b
Found 4 expressions with local error:
| 38.4b | (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0)))) |
| 3.2b | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0))))))) |
| 2.0b | (/ z (/ t (sqrt (+ t a)))) |
| 0.4b | (/ t (sqrt (+ t a))) |
| 92× | frac-sub |
| 67× | frac-times |
| 49× | cbrt-div |
| 37× | associate-*r/ |
| 37× | associate-*l/ |
| 23× | flip-+ |
| 23× | flip3-+ |
| 21× | flip-- |
| 21× | flip3-- |
| 12× | associate-/r* |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | associate--r+ |
| 5× | associate-/r/ |
| 4× | add-exp-log |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | div-inv |
| 3× | sqrt-prod |
| 3× | distribute-lft-in |
| 3× | sub-neg |
| 3× | distribute-rgt-in |
| 2× | associate--l+ |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | rem-cbrt-cube |
| 1× | pow3 |
| 1× | diff-log |
| 1× | pow1/3 |
| 1× | cbrt-prod |
| 41.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0))))))) |
| 7.0ms | (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0)))) |
| 1.0ms | (/ t (sqrt (+ t a))) |
| 1.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 266.0ms | (- (/ z (/ t (sqrt (+ t a)))) (* (- b c) (- (+ a (/ 5.0 6.0)) (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0))))))) |
| 74.0ms | (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0)))) |
| 56.0ms | (/ z (/ t (sqrt (+ t a)))) |
| 36.0ms | (/ t (sqrt (+ t a))) |
| 759.0ms | (* (/ t (sqrt (+ t a))) (* (+ b c) (* (- a (/ 5.0 6.0)) (cbrt (* (* t 3.0) (* t 3.0)))))) |
| 502.0ms | (* (/ t (sqrt (+ t a))) (* (+ b c) (* (- a (/ 5.0 6.0)) (cbrt (* (* t 3.0) (* t 3.0)))))) |
| 482.0ms | (- (* z (* (+ (* b b) (+ (* c c) (* b c))) (+ (* (+ a (/ 5.0 6.0)) (+ a (/ 5.0 6.0))) (+ (* (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0)))) (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0))))) (* (+ a (/ 5.0 6.0)) (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0))))))))) (* (/ t (sqrt (+ t a))) (* (- (pow b 3) (pow c 3)) (- (pow (+ a (/ 5.0 6.0)) 3) (pow (cbrt (* (* (/ 2.0 (* t 3.0)) (/ 2.0 (* t 3.0))) (/ 2.0 (* t 3.0)))) 3))))) |
| 441.0ms | (* (/ t (sqrt (+ t a))) (* (+ b c) (* (- a (/ 5.0 6.0)) (cbrt (* t 3.0))))) |
| 439.0ms | (* (/ t (sqrt (+ t a))) (* (+ (* a a) (- (* (/ 5.0 6.0) (/ 5.0 6.0)) (* a (/ 5.0 6.0)))) (cbrt (* (* (* t 3.0) (* t 3.0)) (* t 3.0))))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0b
18.3% (2.2b remaining)
Error of 2.3b against oracle of 0.1b and baseline of 2.7b