


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: 1.0m)Debug log
| 85.0ms | 436× | body | 80 | nan |
| 77.0ms | 291× | body | 80 | valid |
| 83.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 |
| 29.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 26.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 9.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 3.0ms | (* n (pow (/ l Om) 2)) |
| 974.0ms | (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*)))) |
| 79.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 28.0ms | (* n (pow (/ l Om) 2)) |
| 567.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 517.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 474.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 327.0ms | (* (+ t (* 2 (/ (* l l) Om))) (+ (* U U) (+ (* U* U*) (* U U*)))) |
| 326.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))))))) |
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 |
| 20.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 17.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)) |
| 535.0ms | (pow (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) 1/2) |
| 345.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 94.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 33.0ms | (* n (pow (/ l Om) 2)) |
| 926.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 568.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 508.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 461.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 407.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 |
| 186.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) |
| 8.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 8.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 8.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 584.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) |
| 356.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 354.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 350.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 459.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 396.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) |
| 375.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*)))))) |
| 366.0ms | (exp (* 1/2 (- (+ (log -2) (* 2 (log (/ -1 Om)))) (+ (* 2 (log (/ -1 l))) (+ (* 2 (log (/ -1 n))) (* 2 (log (/ -1 U)))))))) |
| 344.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 |
| 762.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) |
| 10.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 9.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 9.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 584.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) |
| 353.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 343.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 316.0ms | (cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 428.0ms | (exp (* 1/2 (- (+ (* 2 (log (/ 1 Om))) (log -2)) (+ (* 2 (log (/ 1 l))) (+ (* 2 (log (/ 1 n))) (* 2 (log (/ 1 U)))))))) |
| 395.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*)))))) |
| 378.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) |
| 342.0ms | (exp (* 1/2 (- (+ (log -2) (* 2 (log (/ -1 Om)))) (+ (* 2 (log (/ -1 l))) (+ (* 2 (log (/ -1 n))) (* 2 (log (/ -1 U)))))))) |
| 321.0ms | (cbrt (- (* (- (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*)))))) |
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
| 5.0s | 15298× | body | 80 | nan |
| 2.2s | 9221× | body | 80 | valid |