


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 < -7.158607358585534e+32Initial program 33.7
rmApplied associate-/l*31.6
rmApplied unpow231.6
Applied associate-*r*31.1
rmApplied pow131.1
Applied pow131.1
Applied pow-prod-down31.1
Simplified31.1
rmApplied sub-neg31.1
Applied distribute-rgt-in31.1
Simplified29.3
if -7.158607358585534e+32 < t < -1.4586901274772573e-230 or 4.9744455505630306e-176 < t < 8.69812872684918e+198Initial program 31.4
rmApplied associate-/l*28.7
rmApplied unpow228.7
Applied associate-*r*27.6
rmApplied associate-*l*26.8
if -1.4586901274772573e-230 < t < 4.9744455505630306e-176Initial program 36.3
rmApplied associate-/l*33.5
rmApplied unpow233.5
Applied associate-*r*31.7
rmApplied pow131.7
Applied pow131.7
Applied pow-prod-down31.7
Simplified30.7
rmApplied add-sqr-sqrt30.9
if 8.69812872684918e+198 < t Initial program 39.2
rmApplied associate-/l*36.8
rmApplied unpow236.8
Applied associate-*r*36.3
rmApplied sqrt-prod23.1
Final simplification27.7
herbie shell --seed 2019016 +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.3m)Debug log
| 761× | (pre true 80) |
| 464× | (body nan 80) |
| 297× | (body real 80) |
| 68.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: 32.7b
Found 4 expressions with local error:
| 17.9b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.1b | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 9.8b | (/ (* l l) Om) |
| 7.3b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 22× | associate-*r/ |
| 11× | flip-- |
| 11× | flip3-- |
| 9× | distribute-lft-in |
| 9× | distribute-rgt-in |
| 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× | prod-diff |
| 6× | associate-*r* |
| 6× | sqrt-div |
| 6× | sub-neg |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | associate--l+ |
| 4× | expm1-log1p-u |
| 2× | cbrt-unprod |
| 2× | *-commutative |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | associate-*l* |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 21.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 13.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 | (/ (* l l) Om) |
| 887.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 314.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 76.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 16.0ms | (/ (* l l) Om) |
| 633.0ms | (* (* n (pow (/ l Om) 2)) (- U*)) |
| 536.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 482.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 421.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 309.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*))))))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 19.9b
Found 4 expressions with local error:
| 17.9b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.1b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.3b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 5.2b | (* 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 |
| 49.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 26.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 13.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*))))) |
| 325.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 88.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 48.0ms | (* n (pow (/ l Om) 2)) |
| 556.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 507.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 489.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 403.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (+ (* U U) (+ (* U* U*) (* U U*)))) |
| 394.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 (12 fresh and 0 done)
Merged error: 18.9b
Found 4 expressions with local error:
| 17.9b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 10.1b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 7.3b | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 4.4b | (* n (/ l Om)) |
| 104× | associate-*r/ |
| 45× | frac-times |
| 44× | flip-- |
| 44× | frac-sub |
| 44× | flip3-- |
| 30× | associate-*l/ |
| 24× | sqrt-div |
| 11× | add-exp-log |
| 11× | add-cbrt-cube |
| 11× | pow1 |
| 10× | associate-*r* |
| 9× | add-cube-cbrt |
| 9× | distribute-lft-in |
| 9× | *-un-lft-identity |
| 9× | add-sqr-sqrt |
| 9× | distribute-rgt-in |
| 6× | prod-diff |
| 6× | sub-neg |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | cbrt-unprod |
| 4× | associate--l+ |
| 4× | prod-exp |
| 4× | pow-prod-down |
| 4× | expm1-log1p-u |
| 3× | *-commutative |
| 2× | associate-*l* |
| 1× | div-inv |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 35.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 27.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 16.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 2.0ms | (* n (/ l Om)) |
| 1.3s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 360.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 114.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 16.0ms | (* n (/ l Om)) |
| 828.0ms | (* (* (* n l) (/ l Om)) (- (* U U) (* U* U*))) |
| 805.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) (* (+ (* t t) (+ (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l)))) (* t (* 2 (/ l (/ Om l)))))) (* (* (* n l) (/ l Om)) (- (pow U 3) (pow U* 3))))))) |
| 650.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) |
| 554.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) |
| 499.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 (/ l Om)) (/ l Om)) (- (pow U 3) (pow U* 3))))))) |
18 alts after pruning (18 fresh and 0 done)
Merged error: 16.6b
Found 4 expressions with local error:
| 17.9b | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 10.1b | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 5.9b | (* (/ l Om) (- U U*)) |
| 4.7b | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 9× | add-sqr-sqrt |
| 8× | add-cube-cbrt |
| 8× | add-exp-log |
| 8× | add-cbrt-cube |
| 8× | *-un-lft-identity |
| 8× | pow1 |
| 6× | associate-*r* |
| 5× | associate-*l* |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | sqrt-prod |
| 4× | associate-*r/ |
| 4× | distribute-lft-in |
| 4× | sub-neg |
| 4× | expm1-log1p-u |
| 4× | distribute-rgt-in |
| 3× | flip-- |
| 3× | flip3-- |
| 2× | cbrt-unprod |
| 2× | *-commutative |
| 2× | prod-exp |
| 2× | pow-prod-down |
| 2× | frac-times |
| 1× | div-inv |
| 1× | fma-udef |
| 1× | associate-*l/ |
| 1× | pow1/2 |
| 1× | unpow-prod-down |
| 1× | rem-sqrt-square |
| 5.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 4.0ms | (* (/ l Om) (- U U*)) |
| 3.0ms | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 0.0ms | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 1.6s | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 297.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 103.0ms | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 59.0ms | (* (/ l Om) (- U U*)) |
| 451.0ms | (+ (* t t) (+ (* (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om)))) (* t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om)))))) |
| 381.0ms | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 374.0ms | (* (- (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om)))) (* 2 (* U n))) |
| 343.0ms | (sqrt (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1))) |
| 315.0ms | (sqrt (* (cbrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) (cbrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)))) |
17 alts after pruning (17 fresh and 0 done)
Merged error: 15.4b
5.3% (9.8b remaining)
Error of 27.7b against oracle of 18.0b and baseline of 28.3b
| 23967× | (pre true 80) |
| 14776× | (body nan 80) |
| 9191× | (body real 80) |