Use primitive instead of clear_denoms in roots_rational
Add option to return only integer roots from roots_rational