


Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
Results
if t < 8.152146896984202e-243Initial program 33.7
rmApplied associate-/l*30.7
rmApplied add-sqr-sqrt30.9
rmApplied add-cube-cbrt30.9
Applied associate-*r*30.9
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
(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: 2.1m)Debug log
| 129.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-- |
| 8× | frac-sub |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | add-exp-log |
| 6× | sqrt-div |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 4× | add-log-exp |
| 4× | associate--l+ |
| 4× | distribute-lft-in |
| 4× | unpow-prod-down |
| 4× | sub-neg |
| 4× | distribute-rgt-in |
| 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 |
| 52.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 21.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 4.0ms | (/ (* l l) Om) |
| 4.0ms | (* n (pow (/ l Om) 2)) |
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 304.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 46.0ms | (* n (pow (/ l Om) 2)) |
| 17.0ms | (/ (* l l) Om) |
| 391.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ (* l l) Om)) 3)) (+ U U*)) (* (+ (* t t) (+ (* (* 2 (/ (* l l) Om)) (* 2 (/ (* l l) Om))) (* t (* 2 (/ (* l l) Om))))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*))))))) |
| 372.0ms | (* (+ t (* 2 (/ (* l l) Om))) (+ U U*)) |
| 318.0ms | (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))) |
| 318.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ (* l l) Om)) 3)) (+ (* U U) (+ (* U* U*) (* U U*)))) (* (+ (* t t) (+ (* (* 2 (/ (* l l) Om)) (* 2 (/ (* l l) Om))) (* t (* 2 (/ (* l l) Om))))) (* (* n (pow (/ l Om) 2)) (- (pow U 3) (pow U* 3))))))) |
| 295.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (* t t) (* (* 2 (/ (* l l) Om)) (* 2 (/ (* l l) Om)))) (+ U U*)) (* (+ t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*))))))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 22.0b
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-- |
| 8× | add-cube-cbrt |
| 8× | frac-sub |
| 8× | add-exp-log |
| 8× | add-cbrt-cube |
| 8× | *-un-lft-identity |
| 8× | pow1 |
| 8× | add-sqr-sqrt |
| 6× | sqrt-div |
| 6× | distribute-lft-in |
| 6× | sub-neg |
| 6× | distribute-rgt-in |
| 4× | add-log-exp |
| 4× | associate--l+ |
| 4× | unpow-prod-down |
| 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 |
| 41.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 17.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 12.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*))))) |
| 280.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 87.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 35.0ms | (* n (pow (/ l Om) 2)) |
| 382.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ (* U U) (+ (* U* U*) (* U U*)))) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (pow U 3) (pow U* 3))))))) |
| 374.0ms | (* (- U*) (* n (pow (/ l Om) 2))) |
| 366.0ms | (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 359.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ U U*)) |
| 345.0ms | (+ (log (* n (pow (/ l Om) 2))) (log (- U U*))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 20.6b
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 |
| 12× | sqrt-div |
| 8× | add-exp-log |
| 8× | associate--l+ |
| 8× | add-cbrt-cube |
| 8× | distribute-lft-in |
| 8× | sub-neg |
| 8× | pow1 |
| 8× | distribute-rgt-in |
| 6× | add-cube-cbrt |
| 6× | associate-*r* |
| 6× | *-un-lft-identity |
| 6× | add-sqr-sqrt |
| 4× | add-log-exp |
| 2× | cbrt-unprod |
| 2× | sqrt-prod |
| 2× | *-commutative |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | associate-*l* |
| 2× | pow1/2 |
| 22.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 20.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 9.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 9.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*))))) |
| 382.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 286.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 600.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))))))) |
| 382.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ U U*)) |
| 378.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ U U*)) |
| 342.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ U U*)) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*))))))) |
| 328.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ (* U U) (+ (* U* U*) (* U U*)))) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (pow U 3) (pow U* 3))))))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 20.6b
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) (* (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) |
| 12.8b | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 12.8b | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 990× | associate-*r/ |
| 894× | cbrt-div |
| 750× | flip-- |
| 750× | flip3-- |
| 600× | frac-sub |
| 540× | frac-times |
| 348× | sqrt-div |
| 90× | associate-*l/ |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 6× | add-sqr-sqrt |
| 6× | cbrt-prod |
| 4× | add-log-exp |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 2× | sqrt-prod |
| 2× | pow1/3 |
| 2× | pow1/2 |
| 117.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) |
| 17.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.0ms | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 10.0ms | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 1.2s | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- 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*))))) |
| 334.0ms | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 314.0ms | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 665.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (* (- (* t t) (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l))))) (+ U U*)) (* (+ t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*)))))) (cbrt (- (* (- (* 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)))))))) (cbrt (- (* (- (* 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)))))))) |
| 543.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ U U*)) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*)))))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (* (- t (* 2 (/ l (/ Om l)))) (- t (* 2 (/ l (/ Om l))))) (* (* (* n (pow (/ l Om) 2)) (- U U*)) (* (* n (pow (/ l Om) 2)) (- U U*))))))) |
| 507.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ U U*)) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*)))))) (cbrt (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (+ U U*)) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*)))))))) (cbrt (- (* (- (* t t) (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l))))) (+ U U*)) (* (+ t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- (* U U) (* U* U*)))))))) |
| 505.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (* (- t (* 2 (/ l (/ Om l)))) (- t (* 2 (/ l (/ Om l))))) (* (* (* n (pow (/ l Om) 2)) (- U U*)) (* (* n (pow (/ l Om) 2)) (- U U*))))) (cbrt (- (* (- (* 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)))))))) (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) |
| 483.0ms | (sqrt (* (* (* (* 2 n) U) (* (cbrt (- (* (- t (* 2 (/ l (/ Om l)))) (- t (* 2 (/ l (/ Om l))))) (* (* (* n (pow (/ l Om) 2)) (- U U*)) (* (* n (pow (/ l Om) 2)) (- U U*))))) (cbrt (- (* (- (* 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)))))))) (cbrt (- (* (- t (* 2 (/ l (/ Om l)))) (- t (* 2 (/ l (/ Om l))))) (* (* (* n (pow (/ l Om) 2)) (- U U*)) (* (* n (pow (/ l Om) 2)) (- U U*))))))) |
13 alts after pruning (11 fresh and 2 done)
Merged error: 20.6b
8.5% (8.9b remaining)
Error of 29.3b against oracle of 20.4b and baseline of 30.1b