


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
| 84.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 |
| 19.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 12.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 2.0ms | (/ (* l l) Om) |
| 2.0ms | (* n (pow (/ l Om) 2)) |
| 974.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 274.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 36.0ms | (* n (pow (/ l Om) 2)) |
| 10.0ms | (/ (* l l) Om) |
| 513.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)) |
| 442.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 378.0ms | (* (+ t (* 2 (/ (* l l) Om))) (+ (* U U) (+ (* U* U*) (* U U*)))) |
| 339.0ms | (+ (log (* (* 2 n) U)) (log (- (- 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 |
| 45.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 16.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.1s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 334.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 76.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 43.0ms | (* n (pow (/ l Om) 2)) |
| 939.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 596.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 523.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 440.0ms | (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 399.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ (* U U) (+ (* U* U*) (* 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 |
| 42.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 25.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 24.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 24.0ms | (* (* (* 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*))))) |
| 346.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 336.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 734.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ (* U U) (+ (* U* U*) (* U U*)))) |
| 625.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 589.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 547.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 533.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 |
| 14.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))) |
| 13.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U 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*))))) |
| 241.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 202.0ms | (cbrt (- t (* 2 (/ l (/ Om l))))) |
| 605.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))) |
| 604.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)))) |
| 573.0ms | (+ (/ (* (sqrt 2) (* U (* (sqrt -1) (* n l)))) Om) (+ (* 1/2 (/ (* (sqrt 2) (* n (* U* l))) (* (sqrt -1) Om))) (* 1/2 (/ (* (sqrt 2) (* U* l)) (* U (sqrt -1)))))) |
| 572.0ms | (cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3))) |
| 570.0ms | (cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3))) |
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