Create an alternative: W(E8) act on positive roots to get some roots
R = RootSystem(['E',8]).root_lattice()
alpha = R.simple_roots();alpha
[w for w in W if w.action(alpha)==(alpha+alpha)]
Output shows : gap: cannot extend the workspace any more!
*Is their any alternative calculation to simplify the above *