


Bits error versus x



Bits error versus l



Bits error versus t
if t < -5.00212864258164e+19 or -1.271361517563434e-173 < t < -7.128698999489736e-258Initial program 46.1
Simplified46.1
Taylor expanded around -inf 10.3
if -5.00212864258164e+19 < t < -1.271361517563434e-173Initial program 32.2
Simplified32.1
Taylor expanded around inf 12.0
Simplified6.3
rmApplied add-cube-cbrt6.3
Applied associate-*r*6.3
rmApplied add-cube-cbrt6.3
Applied associate-*r*6.4
rmApplied add-sqr-sqrt6.4
Applied associate-*r*6.4
if -7.128698999489736e-258 < t < 2.4158396745379366e+88Initial program 40.3
Simplified40.3
Taylor expanded around inf 18.4
Simplified15.3
rmApplied add-cube-cbrt15.3
Applied associate-*r*15.3
rmApplied add-cube-cbrt15.3
Applied associate-*r*15.4
rmApplied associate-*l*15.3
if 2.4158396745379366e+88 < t Initial program 49.4
Simplified49.4
Taylor expanded around inf 3.0
Simplified3.0
Final simplification9.8
herbie shell --seed 2019022 +o rules:numerics
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
(/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))))
Time bar (total: 36.1s)Debug log
| 90.0ms | 287× | body | 80 | nan |
| 37.0ms | 18× | body | 1280 | nan |
| 29.0ms | 39× | body | 1280 | valid |
| 28.0ms | 163× | body | 80 | valid |
| 14.0ms | 23× | body | 640 | valid |
| 10.0ms | 21× | body | 320 | valid |
| 9.0ms | 17× | body | 640 | nan |
| 3.0ms | 10× | body | 160 | valid |
| 2.0ms | 5× | body | 160 | nan |
| 1.0ms | 2× | body | 320 | nan |
| 470.0ms | (/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 43.7b
Found 4 expressions with local error:
| 21.5b | (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) |
| 12.2b | (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))) |
| 0.5b | (* t (sqrt 2)) |
| 0.0b | (/ (* t (sqrt 2)) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))) |
| 11× | add-sqr-sqrt |
| 10× | add-cube-cbrt |
| 10× | *-un-lft-identity |
| 9× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 6× | associate-/r* |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 3× | associate-*r* |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | fma-udef |
| 1× | *-commutative |
| 1× | div-exp |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 16.0ms | (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) |
| 8.0ms | (/ (* t (sqrt 2)) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))) |
| 2.0ms | (* t (sqrt 2)) |
| 0.0ms | (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))) |
| 322.0ms | (/ (* t (sqrt 2)) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))) |
| 209.0ms | (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) |
| 89.0ms | (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))) |
| 24.0ms | (* t (sqrt 2)) |
| 632.0ms | (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) |
| 628.0ms | (/ (* (* (* t (sqrt 2)) (* t (sqrt 2))) (* t (sqrt 2))) (* (* (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))))) |
| 555.0ms | (/ t (sqrt (* (cbrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) (cbrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))))) |
| 551.0ms | (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l)))) |
| 550.0ms | (/ (sqrt 2) (sqrt (fma (fma (* 2 t) t (* l l)) (/ (+ 1 x) (+ -1 x)) (* l (- l))))) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 4.0b
Found 4 expressions with local error:
| 28.2b | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 3.4b | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 0.5b | (* t (sqrt 2)) |
| 0.2b | (* (* t 4) (/ t x)) |
| 8× | add-sqr-sqrt |
| 7× | add-cube-cbrt |
| 7× | associate-*r* |
| 7× | *-un-lft-identity |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 3× | sqrt-prod |
| 2× | *-commutative |
| 1× | div-inv |
| 1× | fma-udef |
| 1× | cbrt-unprod |
| 1× | associate-*r/ |
| 1× | prod-exp |
| 1× | pow-prod-down |
| 1× | associate-*l* |
| 1× | pow1/2 |
| 1× | rem-sqrt-square |
| 3.0ms | (* (* t 4) (/ t x)) |
| 1.0ms | (* t (sqrt 2)) |
| 1.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 0.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 261.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 49.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 27.0ms | (* t (sqrt 2)) |
| 21.0ms | (* (* t 4) (/ t x)) |
| 478.0ms | (* (* (* (* t 4) (* t 4)) (* t 4)) (* (* (/ t x) (/ t x)) (/ t x))) |
| 141.0ms | (+ (* 2 (pow t 2)) (+ (* 2 (/ (pow l 2) x)) (* 4 (/ (pow t 2) x)))) |
| 121.0ms | (sqrt (* (cbrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) (cbrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))))) |
| 112.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 108.0ms | (+ (* 2 (pow t 2)) (+ (* 2 (/ (pow l 2) x)) (* 4 (/ (pow t 2) x)))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 3.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 3.4b | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 0.3b | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 0.2b | (* (* t 4) (/ t x)) |
| 10× | associate-*r* |
| 9× | add-sqr-sqrt |
| 8× | add-cube-cbrt |
| 8× | add-exp-log |
| 8× | *-un-lft-identity |
| 8× | pow1 |
| 7× | add-cbrt-cube |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 3× | sqrt-prod |
| 3× | cbrt-prod |
| 2× | cbrt-unprod |
| 2× | *-commutative |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | associate-*l* |
| 1× | div-inv |
| 1× | fma-udef |
| 1× | associate-*r/ |
| 1× | pow1/2 |
| 1× | rem-sqrt-square |
| 13.0ms | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 3.0ms | (* (* t 4) (/ t x)) |
| 1.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 0.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 220.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 65.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 27.0ms | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 14.0ms | (* (* t 4) (/ t x)) |
| 835.0ms | (* (* (* (* t 4) (* t 4)) (* t 4)) (* (* (/ t x) (/ t x)) (/ t x))) |
| 165.0ms | (* (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2))))) (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2))))) (sqrt 2)) |
| 110.0ms | (+ (* 2 (pow t 2)) (+ (* 2 (/ (pow l 2) x)) (* 4 (/ (pow t 2) x)))) |
| 96.0ms | (sqrt (* (cbrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) (cbrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))))) |
| 92.0ms | (sqrt (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 3.9b
Found 4 expressions with local error:
| 28.2b | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 3.4b | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 0.4b | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) |
| 0.2b | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 18× | add-exp-log |
| 18× | pow1 |
| 11× | add-cbrt-cube |
| 10× | associate-*r* |
| 9× | cbrt-unprod |
| 9× | prod-exp |
| 9× | pow-prod-down |
| 9× | add-sqr-sqrt |
| 9× | cbrt-prod |
| 8× | add-cube-cbrt |
| 8× | *-un-lft-identity |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 3× | sqrt-prod |
| 2× | *-commutative |
| 2× | associate-*l* |
| 1× | fma-udef |
| 1× | pow1/2 |
| 1× | rem-sqrt-square |
| 57.0ms | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 28.0ms | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) |
| 1.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 0.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 698.0ms | (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) |
| 278.0ms | (sqrt (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x)))) |
| 49.0ms | (fma 2 (fma (/ l x) l (* t t)) (* (* t 4) (/ t x))) |
| 27.0ms | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 372.0ms | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 343.0ms | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 338.0ms | (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2))))) (cbrt (cbrt (sqrt 2)))) |
| 232.0ms | (* (* (* (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2))))) (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2))))) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 219.0ms | (+ (log (* (* t (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (* (cbrt (cbrt (sqrt 2))) (cbrt (cbrt (sqrt 2)))))) (log (cbrt (cbrt (sqrt 2))))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 3.9b
75.1% (6.0b remaining)
Error of 9.8b against oracle of 3.8b and baseline of 27.8b
| 2.5s | 10278× | body | 80 | nan |
| 1.2s | 5110× | body | 80 | valid |
| 994.0ms | 608× | body | 1280 | nan |
| 932.0ms | 886× | body | 1280 | valid |
| 914.0ms | 1063× | body | 640 | valid |
| 572.0ms | 620× | body | 640 | nan |
| 560.0ms | 311× | body | 160 | valid |
| 461.0ms | 630× | body | 320 | valid |
| 176.0ms | 333× | body | 320 | nan |
| 69.0ms | 199× | body | 160 | nan |