Restrict domain of nonlinsolve

[Simone Perriello]

Well, the title says pretty much everything. Is it possible to restrict the domain of nonlinsolve to solve, say, only for real numbers? Restrictions on the domain of symbols doesn't seem to work for nonlinsolve, while works well with the simple solve.

[Chris Smith]

You could wrap your solution in a FiniteSet and attempt an Intersection:

>>> nonlinsolve([x**2 - 2*y**2 -2, x*y - 2], [x, y])

{(-2, -1), (2, 1), (-sqrt(2)*I, sqrt(2)*I), (sqrt(2)*I, -sqrt(2)*I)}

>>> _.intersect(S.Naturals**2)

{(2, 1)}