


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 U* < 1.069014172200058e+80Initial program 32.3
rmApplied associate-/l*29.8
rmApplied unpow229.8
Applied associate-*r*29.0
rmApplied pow129.0
Applied pow129.0
Applied pow-prod-down29.0
Simplified28.8
rmApplied sub-neg28.8
Applied distribute-rgt-in28.8
Simplified27.1
if 1.069014172200058e+80 < U* Initial program 35.6
rmApplied associate-/l*33.4
rmApplied unpow233.4
Applied associate-*r*32.3
rmApplied associate-*l*31.7
Final simplification28.0
herbie shell --seed 2019008 +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
1 calls. Slowest were:
| 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: 33.0b
Found 4 expressions with local error:
| 19.4b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.9b | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.9b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 7.4b | (/ (* l l) Om) |
4 calls. Slowest were:
| 21.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*)) |
4 calls. Slowest were:
| 1.4s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 348.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 67.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 13.0ms | (/ (* l l) Om) |
58 calls. Slowest were:
| 617.0ms | (+ (log (* n (pow (/ l Om) 2))) (log (- U U*))) |
| 599.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 566.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 22.9b
Found 4 expressions with local error:
| 19.4b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 10.9b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 7.9b | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 4.5b | (* n (pow (/ l Om) 2)) |
4 calls. Slowest were:
| 21.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*))))) |
| 6.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
4 calls. Slowest were:
| 1.2s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))) |
| 311.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))) |
| 90.0ms | (* (* n (pow (/ l Om) 2)) (- U U*)) |
| 44.0ms | (* n (pow (/ l Om) 2)) |
66 calls. Slowest were:
| 622.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 536.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
| 529.0ms | (* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U)) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 22.1b
Found 4 expressions with local error:
| 19.4b | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 10.9b | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 7.9b | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 4.7b | (* n (/ l Om)) |
4 calls. Slowest were:
| 34.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 16.0ms | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 14.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
4 calls. Slowest were:
| 1.1s | (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))))) |
| 331.0ms | (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*)))) |
| 98.0ms | (* (* (* n (/ l Om)) (/ l Om)) (- U U*)) |
| 12.0ms | (* n (/ l Om)) |
117 calls. Slowest were:
| 886.0ms | (sqrt (* (* (* 2 n) U) (- (* (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)) (* (* Om 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) (- (pow U 3) (pow U* 3))))))) |
| 858.0ms | (+ (- t (* 2 (/ l (/ Om l)))) (* (* (* n (/ l Om)) (/ l Om)) (- U U*))) |
| 580.0ms | (* (+ t (* 2 (/ l (/ Om l)))) (* Om (+ (* U U) (+ (* U* U*) (* U U*))))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 21.3b
Found 4 expressions with local error:
| 19.4b | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 10.9b | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 5.6b | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 5.1b | (* (/ l Om) (- U U*)) |
4 calls. Slowest were:
| 5.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 4.0ms | (* (/ l Om) (- U U*)) |
| 2.0ms | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
4 calls. Slowest were:
| 1.5s | (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1)) |
| 318.0ms | (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) |
| 92.0ms | (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))) |
| 57.0ms | (* (/ l Om) (- U U*)) |
48 calls. Slowest were:
| 556.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)))))) |
| 347.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)))) |
| 311.0ms | (sqrt (sqrt (pow (* (* 2 (* U n)) (- t (fma (* (/ l Om) (- U U*)) (* n (/ l Om)) (* (* 2 l) (/ l Om))))) 1))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 18.9b
3.3% (9.3b remaining)
Error of 28.0b against oracle of 18.7b and baseline of 28.3b