


Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
if t < 8.152146896984202e-243Initial program 33.7
rmApplied associate-/l*30.7
rmApplied add-sqr-sqrt30.9
rmApplied add-cube-cbrt31.0
Applied prod-diff31.0
Applied distribute-rgt-in31.0
Simplified31.0
if 8.152146896984202e-243 < t Initial program 31.6
rmApplied associate-/l*28.9
rmApplied sqrt-prod27.2
Final simplification29.3
herbie shell --seed 2019010 +o rules:numerics
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))
Time bar (total: 1.2m)Debug log
| 131.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 31.7b
Found 4 expressions with local error:
| 18.1b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.3b | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.7b | (/ (* l l) Om) |
| 7.3b | (* n (pow (/ l Om) 2)) |
| 20× | associate-*r/ |
| 11× | associate-*r* |
| 10× | flip-- |
| 10× | flip3-- |
| 9× | add-cube-cbrt |
| 9× | *-un-lft-identity |
| 9× | add-sqr-sqrt |
| 8× | frac-sub |
| 7× | distribute-lft-in |
| 7× | distribute-rgt-in |
| 6× | add-exp-log |
| 6× | prod-diff |
| 6× | sqrt-div |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | associate--l+ |
| 4× | unpow-prod-down |
| 4× | sub-neg |
| 4× | expm1-log1p-u |
| 2× | div-inv |
| 2× | *-commutative |
| 1× | unpow2 |
| 1× | associate-/l* |
| 1× | cbrt-unprod |
| 1× | sqrt-prod |
| 1× | prod-exp |
| 1× | pow-prod-down |
| 1× | associate-*l* |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 37.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 36.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 4.0ms | (/ (* l l) Om) |
| 3.0ms | (* n (pow (/ l Om) 2)) |
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 290.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 28.0ms | (* n (pow (/ l Om) 2)) |
| 12.0ms | (/ (* l l) Om) |
| 614.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 588.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 524.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 358.0ms | (sqrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 353.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 23.3b
Found 4 expressions with local error:
| 18.1b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.3b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.3b | (* n (pow (/ l Om) 2)) |
| 6.2b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 22× | associate-*r/ |
| 14× | associate-*r* |
| 11× | flip-- |
| 11× | flip3-- |
| 10× | add-cube-cbrt |
| 10× | *-un-lft-identity |
| 10× | add-sqr-sqrt |
| 9× | distribute-lft-in |
| 9× | distribute-rgt-in |
| 8× | frac-sub |
| 8× | add-exp-log |
| 8× | add-cbrt-cube |
| 8× | pow1 |
| 6× | prod-diff |
| 6× | sqrt-div |
| 6× | sub-neg |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | associate--l+ |
| 4× | unpow-prod-down |
| 4× | expm1-log1p-u |
| 3× | *-commutative |
| 2× | cbrt-unprod |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | associate-*l* |
| 1× | unpow2 |
| 1× | div-inv |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 40.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 12.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 11.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 4.0ms | (* n (pow (/ l Om) 2)) |
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 248.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 95.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 27.0ms | (* n (pow (/ l Om) 2)) |
| 570.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (* t t) (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l))))) (+ (* U U) (+ (* U* U*) (* U U*)))) (* (+ t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- (pow U 3) (pow U* 3))))))) |
| 552.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 529.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 512.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 459.0ms | (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 21.9b
Found 4 expressions with local error:
| 18.1b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 18.1b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.3b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 10.3b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 40× | associate-*r/ |
| 20× | flip-- |
| 20× | flip3-- |
| 16× | frac-sub |
| 14× | distribute-lft-in |
| 14× | distribute-rgt-in |
| 12× | prod-diff |
| 12× | sqrt-div |
| 10× | add-cube-cbrt |
| 10× | *-un-lft-identity |
| 10× | add-sqr-sqrt |
| 8× | add-exp-log |
| 8× | associate--l+ |
| 8× | add-cbrt-cube |
| 8× | sub-neg |
| 8× | pow1 |
| 6× | associate-*r* |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 2× | cbrt-unprod |
| 2× | sqrt-prod |
| 2× | *-commutative |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | associate-*l* |
| 2× | pow1/2 |
| 22.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 21.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 21.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 20.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 1.1s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 326.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 276.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 976.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 663.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 591.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 575.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 528.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 21.9b
Found 4 expressions with local error:
| 18.1b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 18.1b | (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) |
| 10.5b | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 10.5b | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 10× | associate-*r/ |
| 8× | sqrt-div |
| 8× | add-sqr-sqrt |
| 7× | flip-- |
| 7× | add-cube-cbrt |
| 7× | flip3-- |
| 7× | *-un-lft-identity |
| 6× | cbrt-prod |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | sqrt-prod |
| 4× | frac-sub |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | cbrt-div |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 2× | pow1/3 |
| 2× | pow1/2 |
| 1× | flip-+ |
| 1× | flip3-+ |
| 1× | rem-sqrt-square |
| 16.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 7.0ms | (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) |
| 1.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 1.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 2.6s | (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) |
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 239.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 204.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 829.0ms | (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) |
| 637.0ms | (cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3))) |
| 608.0ms | (cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3))) |
| 535.0ms | (sqrt (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))) |
| 500.0ms | (sqrt (* (cbrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) (cbrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))))) |
11 alts after pruning (9 fresh and 2 done)
Merged error: 21.9b
9.3% (8.0b remaining)
Error of 29.3b against oracle of 21.3b and baseline of 30.1b