


Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 5.1297397571292906e+138Initial program 50.2
rmApplied times-frac35.0
Applied associate-/l*35.0
rmApplied flip-+35.0
Applied associate-/r/35.0
Applied add-sqr-sqrt35.0
Applied times-frac35.0
Applied *-un-lft-identity35.0
Applied times-frac35.0
Applied times-frac35.0
Simplified35.0
if 5.1297397571292906e+138 < alpha Initial program 62.5
rmApplied times-frac55.2
Applied associate-/l*55.2
rmApplied flip-+55.2
Applied associate-/r/55.2
Applied add-sqr-sqrt55.2
Applied times-frac55.2
Applied *-un-lft-identity55.2
Applied times-frac55.2
Applied times-frac55.2
Simplified55.2
Taylor expanded around 0 48.6
Final simplification37.5
herbie shell --seed 2019016
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:pre (and (> alpha -1) (> beta -1) (> i 1))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))
Time bar (total: 5.6m)Debug log
| 256× | (pre true 80) |
| 256× | (body real 80) |
| 201.0ms | (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 51.1b
Found 4 expressions with local error:
| 45.3b | (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) |
| 32.3b | (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 0.1b | (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) |
| 0.0b | (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) |
| 48× | flip-+ |
| 48× | flip3-+ |
| 44× | times-frac |
| 42× | associate-/r/ |
| 38× | add-exp-log |
| 38× | add-cbrt-cube |
| 32× | frac-times |
| 28× | associate-*r/ |
| 25× | add-cube-cbrt |
| 25× | *-un-lft-identity |
| 25× | add-sqr-sqrt |
| 16× | associate-*l/ |
| 14× | associate-/l* |
| 13× | div-exp |
| 13× | cbrt-undiv |
| 11× | pow1 |
| 10× | cbrt-unprod |
| 10× | prod-exp |
| 9× | associate-/l/ |
| 6× | div-inv |
| 6× | associate-*r* |
| 4× | add-log-exp |
| 4× | associate-*l* |
| 4× | associate-/r* |
| 2× | *-commutative |
| 2× | pow-prod-down |
| 2× | distribute-lft-in |
| 2× | frac-2neg |
| 2× | clear-num |
| 2× | distribute-rgt-in |
| 1× | flip-- |
| 1× | pow-prod-up |
| 1× | flip3-- |
| 1× | pow-plus |
| 1× | pow2 |
| 37.0ms | (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 32.0ms | (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) |
| 11.0ms | (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) |
| 11.0ms | (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) |
| 741.0ms | (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 156.0ms | (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) |
| 82.0ms | (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) |
| 41.0ms | (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) |
| 998.0ms | (/ (- (+ alpha beta) (* 2 i)) (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 910.0ms | (* (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i))) (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i)))) |
| 858.0ms | (* (+ (pow (+ alpha beta) 3) (pow (* 2 i) 3)) (+ (pow (+ alpha beta) 3) (pow (* 2 i) 3))) |
| 634.0ms | (* (* i (+ (+ alpha beta) i)) (* i (+ (+ alpha beta) i))) |
| 618.0ms | (* (* i (+ (+ alpha beta) i)) (* i (+ (+ alpha beta) i))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 32.2b
Found 4 expressions with local error:
| 31.3b | (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) |
| 31.3b | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 26.2b | (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) |
| 4.9b | (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) |
| 1416× | times-frac |
| 530× | add-cube-cbrt |
| 530× | *-un-lft-identity |
| 530× | add-sqr-sqrt |
| 186× | associate-/r/ |
| 99× | div-inv |
| 89× | flip-+ |
| 89× | flip3-+ |
| 71× | associate-/r* |
| 34× | add-exp-log |
| 34× | add-cbrt-cube |
| 20× | div-exp |
| 20× | cbrt-undiv |
| 16× | associate-/l* |
| 7× | associate-/l/ |
| 4× | add-log-exp |
| 4× | frac-2neg |
| 4× | clear-num |
| 4× | pow1 |
| 2× | associate-*r/ |
| 1× | flip-- |
| 1× | div-sub |
| 1× | flip3-- |
| 38.0ms | (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) |
| 32.0ms | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 28.0ms | (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) |
| 16.0ms | (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) |
| 789.0ms | (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) |
| 540.0ms | (/ (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) |
| 68.0ms | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 68.0ms | (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) |
| 975.0ms | (/ 1 (/ (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (* (cbrt (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (cbrt (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))))) |
| 906.0ms | (/ (/ i 1) (/ (* (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) (/ 1 1))) |
| 901.0ms | (/ (/ (* i (+ (+ alpha beta) i)) (- (* (+ alpha beta) (+ alpha beta)) (* (* 2 i) (* 2 i)))) (/ 1 (/ 1 1))) |
| 892.0ms | (/ (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (sqrt (+ (+ alpha beta) (* 2 i))))) |
| 884.0ms | (/ (/ i (* (cbrt (+ (+ alpha beta) (* 2 i))) (cbrt (+ (+ alpha beta) (* 2 i))))) (/ 1 (/ 1 (* (cbrt (+ (+ alpha beta) (* 2 i))) (cbrt (+ (+ alpha beta) (* 2 i))))))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 32.0b
Found 4 expressions with local error:
| 37.4b | (/ (+ (* (+ alpha beta) i) (+ (* i i) (* alpha beta))) (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i))) |
| 35.2b | (sqrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)) |
| 35.2b | (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 8.2b | (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i)) |
| 16× | add-sqr-sqrt |
| 14× | add-cube-cbrt |
| 14× | *-un-lft-identity |
| 9× | times-frac |
| 7× | flip-+ |
| 7× | flip3-+ |
| 6× | associate-/l/ |
| 6× | add-log-exp |
| 6× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 4× | frac-add |
| 4× | associate-*l/ |
| 4× | sqrt-div |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | associate-/r* |
| 2× | flip-- |
| 2× | associate-/r/ |
| 2× | pow1/2 |
| 2× | flip3-- |
| 2× | associate-+l+ |
| 2× | rem-sqrt-square |
| 1× | div-inv |
| 1× | div-exp |
| 1× | distribute-lft-in |
| 1× | frac-2neg |
| 1× | sum-log |
| 1× | clear-num |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 1× | distribute-rgt-in |
| 34.0ms | (/ (+ (* (+ alpha beta) i) (+ (* i i) (* alpha beta))) (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i))) |
| 26.0ms | (sqrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)) |
| 13.0ms | (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 6.0ms | (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i)) |
| 141.0ms | (/ (+ (* (+ alpha beta) i) (+ (* i i) (* alpha beta))) (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i))) |
| 135.0ms | (sqrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)) |
| 112.0ms | (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 63.0ms | (+ (* (+ alpha beta) (+ alpha beta)) (* (* i -4) i)) |
| 557.0ms | (sqrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)) |
| 536.0ms | (sqrt (* (cbrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)) (cbrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0)))) |
| 466.0ms | (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 453.0ms | (sqrt (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 438.0ms | (sqrt (sqrt (- (* (+ (+ alpha beta) (* i 2)) (+ (+ alpha beta) (* i 2))) 1.0))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 32.0b
Found 4 expressions with local error:
| 31.3b | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 26.5b | (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 4.6b | (* (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 0.0b | (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) |
| 162× | times-frac |
| 95× | add-cube-cbrt |
| 95× | *-un-lft-identity |
| 95× | add-sqr-sqrt |
| 54× | associate-*r* |
| 34× | add-exp-log |
| 34× | add-cbrt-cube |
| 22× | associate-/r/ |
| 18× | associate-/l* |
| 16× | flip-+ |
| 16× | flip3-+ |
| 11× | div-inv |
| 11× | pow1 |
| 11× | frac-times |
| 10× | cbrt-unprod |
| 10× | prod-exp |
| 10× | div-exp |
| 10× | cbrt-undiv |
| 8× | associate-*l/ |
| 7× | associate-*r/ |
| 7× | associate-*l* |
| 6× | associate-/r* |
| 4× | add-log-exp |
| 4× | pow-prod-down |
| 3× | associate-/l/ |
| 2× | flip-- |
| 2× | *-commutative |
| 2× | flip3-- |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | distribute-lft-in |
| 1× | distribute-rgt-in |
| 42.0ms | (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 39.0ms | (* (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 17.0ms | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 4.0ms | (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) |
| 735.0ms | (* (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 198.0ms | (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) |
| 102.0ms | (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) |
| 49.0ms | (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) |
| 857.0ms | (/ (/ 1 (+ (+ alpha beta) (* 2 i))) (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 842.0ms | (/ (/ (sqrt (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) 1) (* (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (cbrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))) |
| 682.0ms | (/ (/ 1 (+ (+ alpha beta) (* 2 i))) (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
| 665.0ms | (* (* (* (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i))) (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i)))) (* (/ i (+ (+ alpha beta) (* i 2))) (+ beta (+ alpha i)))) (/ (* (* (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (* (* (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))) |
| 634.0ms | (/ (/ 1 (+ (+ alpha beta) (* 2 i))) (sqrt (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 32.0b
42.1% (1.7b remaining)
Error of 37.5b against oracle of 35.8b and baseline of 38.7b
| 8000× | (pre true 80) |
| 8000× | (body real 80) |