


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) |
| 149.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 |
| 23.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 18.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 11.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 4.0ms | (/ (* l l) Om) |
| 1.5s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 298.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 84.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 10.0ms | (/ (* l l) Om) |
| 589.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 512.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 500.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 424.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*))))))) |
| 382.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: 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 |
| 28.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 15.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 6.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 2.0ms | (* n (pow (/ l Om) 2)) |
| 1.1s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 306.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 81.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 35.0ms | (* n (pow (/ l Om) 2)) |
| 985.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 594.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 543.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 411.0ms | (cbrt (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 402.0ms | (+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
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 |
| 55.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 32.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 13.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 1.0ms | (* n (/ l Om)) |
| 1.3s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 300.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 76.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 13.0ms | (* n (/ l Om)) |
| 541.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) |
| 540.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (* t t) (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l))))) (* Om Om)) (* (+ t (* 2 (/ l (/ Om l)))) (* (* (* n l) l) (- U U*)))))) |
| 514.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) |
| 501.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))))))) |
| 461.0ms | (* (* (* n l) (/ l Om)) (- (* U U) (* U* U*))) |
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 |
| 10.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 7.0ms | (* (/ l Om) (- U U*)) |
| 5.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.5s | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 290.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 85.0ms | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 48.0ms | (* (/ l Om) (- U U*)) |
| 731.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 476.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)))))) |
| 319.0ms | (- (+ (* 2 (/ (pow l 2) Om)) (/ (* U (* n (pow l 2))) (pow Om 2))) (/ (* n (* U* (pow l 2))) (pow Om 2))) |
| 302.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)))) |
| 276.0ms | (* (- (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om)))) (* 2 (* U n))) |
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) |