avoid using roots(multiple=True) if possible

- add test for DiracDelta simplification involving several roots