wip add test for freethreading

Author: oscarbenjamin

meson.build

@@ -78,6 +78,12 @@ add_project_arguments(
     language : 'cython'
 )
 
+# Enable free-threading if Cython is new enough:
+#cy = meson.get_compiler('cython')
+#if cy.version().version_compare('>=3.1.0')
+#    add_project_arguments('-Xfreethreading_compatible=true', language : 'cython')
+#endif
+
 if get_option('coverage')
    add_project_arguments('-X', 'linetrace=True', language : 'cython')
    add_project_arguments('-DCYTHON_TRACE=1', language : 'c')

src/flint/test/test_all.py

@@ -4,6 +4,8 @@
 import pickle
 import doctest
 import platform
+import random
+from threading import Thread
 
 from flint.utils.flint_exceptions import DomainError, IncompatibleContextError
 
@@ -1615,6 +1617,35 @@ def test_pickling():
         obj2 = pickle.loads(s)
         assert obj == obj2
 
+def test_python_threads():
+    iterations = 10**5
+    threads = 3 + 1
+    size = 3
+    M = flint.fmpz_mat([[0]*size for _ in range(size)])
+
+    def set_values():
+        for i in range(iterations // 5):
+            i = random.randrange(M.nrows())
+            j = random.randrange(M.ncols())
+            if random.uniform(0, 1) > 0.5:
+                # Bigger than 2**62:
+                M[i,j] = 10**128
+            else:
+                # Smaller than 2**62:
+                M[i,j] = 0
+
+    def get_dets():
+        for i in range(iterations):
+            M.det()
+
+    threads = [Thread(target=set_values) for _ in range(threads-1)]
+    threads.append(Thread(target=get_dets))
+
+    for t in threads:
+        t.start()
+    for t in threads:
+        t.join()
+
 def test_fmpz_mod():
     from flint import fmpz_mod_ctx, fmpz, fmpz_mod
 
@@ -4586,6 +4617,7 @@ def test_all_tests():
     test_arb,
 
     test_pickling,
+    test_python_threads,
 
     test_all_tests,
 ]