permute_soln.py

# code from sympy.diophantine

            for solution in solns_pell:
                for X_i in [-solution[0], solution[0]]:
                    for Y_i in [-solution[1], solution[1]]:
                        s = P*Matrix([X_i, Y_i]) + Q
                        try:
                            sol.add(tuple([as_int(_) for _ in s]))
                        except ValueError:

for s in solns_pell:
                sol.add(tuple(s))