avoid using roots(multiple=True) if possible

- add test for DiracDelta simplification involving several roots
  and simplify tests in its computation