Simplification of square root of sum

Question

Using SymPy, can I automatically simplify something like this:

sqrt(a**2 + 2ab + b**2)

to:

[(a+b), -(a+b)]?

Answer 1

Perhaps something like:

>>> var('a b',positive=True)
(a, b)
>>> solve(x**2-(a**2 + 2*a*b + b**2), x)
[-a - b, a + b]