I have added a ``limit`` default argument to limit the number of roots; added a few tests which were
very slow before adding this limitation.