


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 (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) < 4.113440262758482e-308Initial program 56.7
rmApplied pow1/256.7
rmApplied associate-*l*42.3
if 4.113440262758482e-308 < (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) < 1.4037099479477216e+303Initial program 1.2
rmApplied pow1/21.2
if 1.4037099479477216e+303 < (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) Initial program 60.2
rmApplied sqrt-prod54.8
Final simplification28.7
herbie shell --seed 2019022 +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: 45.7s)Debug log
| 70.0ms | 436× | body | 80 | nan |
| 70.0ms | 291× | body | 80 | valid |
| 136.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: 34.4b
Found 4 expressions with local error:
| 21.6b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 8.7b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 8.2b | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.5b | (* n (pow (/ l Om) 2)) |
| 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 |
| 20.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 15.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 6.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 2.0ms | (* n (pow (/ l Om) 2)) |
| 1.3s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 323.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 87.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 40.0ms | (* n (pow (/ l Om) 2)) |
| 587.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 502.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 497.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 399.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))))))) |
| 338.0ms | (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 26.8b
Found 4 expressions with local error:
| 21.6b | (pow (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) 1/2) |
| 8.7b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 8.2b | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.5b | (* n (pow (/ l Om) 2)) |
| 14× | associate-*r* |
| 12× | associate-*r/ |
| 11× | add-exp-log |
| 11× | pow1 |
| 10× | add-cube-cbrt |
| 10× | *-un-lft-identity |
| 10× | add-sqr-sqrt |
| 9× | distribute-lft-in |
| 9× | distribute-rgt-in |
| 8× | add-cbrt-cube |
| 6× | flip-- |
| 6× | prod-diff |
| 6× | flip3-- |
| 6× | sub-neg |
| 5× | unpow-prod-down |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | frac-sub |
| 4× | associate--l+ |
| 4× | expm1-log1p-u |
| 3× | *-commutative |
| 3× | prod-exp |
| 3× | pow-prod-down |
| 2× | cbrt-unprod |
| 2× | pow-exp |
| 2× | associate-*l* |
| 2× | pow-pow |
| 1× | unpow1/2 |
| 1× | unpow2 |
| 1× | div-inv |
| 1× | pow-to-exp |
| 26.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 24.0ms | (pow (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) 1/2) |
| 6.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 2.0ms | (* n (pow (/ l Om) 2)) |
| 461.0ms | (pow (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) 1/2) |
| 262.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 83.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 28.0ms | (* n (pow (/ l Om) 2)) |
| 373.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 363.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 358.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 289.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 280.0ms | (* (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) 1/2) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 26.8b
Found 4 expressions with local error:
| 21.6b | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) 1/2) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 18× | cbrt-div |
| 15× | flip-- |
| 15× | flip3-- |
| 14× | add-exp-log |
| 14× | pow1 |
| 12× | frac-sub |
| 12× | associate-*r/ |
| 9× | cbrt-prod |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | prod-exp |
| 6× | pow-prod-down |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | pow-exp |
| 4× | add-cbrt-cube |
| 4× | pow-pow |
| 4× | expm1-log1p-u |
| 3× | pow1/3 |
| 1× | unpow1/2 |
| 1× | unpow-prod-down |
| 1× | pow-to-exp |
| 124.0ms | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) 1/2) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 359.0ms | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) 1/2) |
| 212.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 209.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 209.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 281.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 249.0ms | (exp (* 1/2 (- (+ (log -2) (* 2 (log (/ -1 Om)))) (+ (* 2 (log (/ -1 l))) (+ (* 2 (log (/ -1 n))) (* 2 (log (/ -1 U)))))))) |
| 242.0ms | (pow (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) 1/2) |
| 207.0ms | (cbrt (sqrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 207.0ms | (cbrt (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 26.8b
Found 4 expressions with local error:
| 21.6b | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (* (cbrt (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))) 1/2) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 13.9b | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 26× | add-exp-log |
| 26× | pow1 |
| 18× | cbrt-div |
| 15× | flip-- |
| 15× | prod-exp |
| 15× | pow-prod-down |
| 15× | flip3-- |
| 12× | frac-sub |
| 12× | associate-*r/ |
| 9× | cbrt-prod |
| 7× | add-cube-cbrt |
| 7× | pow-exp |
| 7× | *-un-lft-identity |
| 7× | pow-pow |
| 7× | add-sqr-sqrt |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | add-cbrt-cube |
| 4× | expm1-log1p-u |
| 3× | pow1/3 |
| 1× | unpow1/2 |
| 1× | unpow-prod-down |
| 1× | pow-to-exp |
| 317.0ms | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (* (cbrt (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))) 1/2) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 5.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 359.0ms | (pow (* (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (* (cbrt (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) (cbrt (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))) 1/2) |
| 206.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 204.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 204.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 505.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 257.0ms | (exp (* 1/2 (- (+ (log -2) (* 2 (log (/ -1 Om)))) (+ (* 2 (log (/ -1 l))) (+ (* 2 (log (/ -1 n))) (* 2 (log (/ -1 U)))))))) |
| 235.0ms | (pow (* (* (* 2 n) U) (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) 1/2) |
| 204.0ms | (cbrt (sqrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 200.0ms | (cbrt (* (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))) |
12 alts after pruning (9 fresh and 3 done)
Merged error: 26.8b
55.6% (3.6b remaining)
Error of 28.7b against oracle of 25.1b and baseline of 33.2b
| 2.8s | 15298× | body | 80 | nan |
| 1.7s | 9221× | body | 80 | valid |