allow conditions being evaluated in Solve

Author: BlankShrimp

mathics/builtin/numbers/calculus.py

@@ -2246,6 +2246,7 @@ def eval(self, eqs, vars, evaluation: Evaluation):
             eq_list = eqs.elements
         else:
             eq_list = [eqs]
+        sympy_conditions = []
         sympy_eqs = []
         sympy_denoms = []
         for eq in eq_list:
@@ -2254,8 +2255,7 @@ def eval(self, eqs, vars, evaluation: Evaluation):
             elif eq is SymbolFalse:
                 return ListExpression()
             elif not eq.has_form("Equal", 2):
-                evaluation.message("Solve", "eqf", eqs)
-                return
+                sympy_conditions.append(eq.to_sympy())
             else:
                 left, right = eq.elements
                 left = left.to_sympy()
@@ -2269,6 +2269,10 @@ def eval(self, eqs, vars, evaluation: Evaluation):
                 numer, denom = eq.as_numer_denom()
                 sympy_denoms.append(denom)
 
+        if not sympy_eqs:
+            evaluation.message("Solve", "eqf", eqs)
+            return
+
         vars_sympy = [var.to_sympy() for var in vars]
         if None in vars_sympy:
             return
@@ -2333,6 +2337,10 @@ def transform_solution(sol):
             for sol in result:
                 results.extend(transform_solution(sol))
             result = results
+            # filter with conditions before further translation
+            conditions = sympy.And(*sympy_conditions)
+            result = [sol for sol in result if conditions.subs(sol)]
+
             if any(
                 sol and any(var not in sol for var in all_vars_sympy) for sol in result
             ):