Option to restore old behaviour for point on variaty (change breaks pseudo-random seeding)

Author: GDeLaurentis

syngular/__init__.py

@@ -18,6 +18,7 @@
 UNICODEPOWERS = True  # noqa, True or False
 NORMALIZE_POWERS_PATTERNS = ()  # eg. re.compile(r"(mt)(\d+)") will force mt2 to be treated as mt^2
 USE_ELLIPSIS_FOR_PRINT = False  # noqa, toggles ellipsis in for str. Use locally for prints only.
+POINT_ON_VARIETY_RANDOM_SOLUTION = True  # noqa, if False, will use the first root, otherwise a (seeded) random one is picked.
 
 __all__ = [
     "__version__",

syngular/variety.py

@@ -251,8 +251,10 @@ def point_on_variety(self, field, base_point={}, directions=None, valuations=tup
             check_solutions(oSemiNumericalIdeal.groebner_basis, root_dicts, field)  # they may be stricter then wanted for mpc.
 
             try:
-                # root_dict = root_dicts[0]
-                root_dict = random.sample(root_dicts, 1)[0]
+                if not syngular.POINT_ON_VARIETY_RANDOM_SOLUTION:
+                    root_dict = root_dicts[0]
+                else:
+                    root_dict = random.sample(root_dicts, 1)[0]
             except (IndexError, ValueError):
                 if not field.is_algebraically_closed:
                     raise RootNotInFieldError(f"Got root_dicts: {root_dicts}, for lex Groebner basis:\n{oSemiNumericalIdeal.groebner_basis}.")