diophantine returns permutations when permute=False

[smichr]

For binary_quadratic types, permutations of the solution are returned. This feature may apply to other types, too:

>>> diophantine(x**2+y**2-15**2)
{(-12, -9), (9, -12), (0, 15), (-15, 0), (15, 0), (-9, -12), (0, -15), (-12, 9),
 (-9, 12), (12, -9), (12, 9), (9, 12)}

It seems like only the base solution should be returned in this case and the permutations when permute=True:

>>> set([tuple(sorted([abs(i) for i in i])) for i in _])
{(0, 15), (9, 12)}

[harshavardhan-99]

I would like to work on this issue