fix: Python 3.10+ compatibility (collections.abc, math.factorial)

Author: basnijholt

adaptive/learner/learnerND.py

@@ -1,5 +1,6 @@
 # -*- coding: utf-8 -*-
-from collections import OrderedDict, Iterable
+from collections import OrderedDict
+from collections.abc import Iterable
 import heapq
 import itertools
 import random
@@ -24,7 +25,7 @@ def volume(simplex, ys=None):
     dim = len(simplex) - 1
 
     # See https://www.jstor.org/stable/2315353
-    vol = np.abs(np.linalg.det(matrix)) / np.math.factorial(dim)
+    vol = np.abs(np.linalg.det(matrix)) / math.factorial(dim)
     return vol
 
 

adaptive/learner/triangulation.py

@@ -1,4 +1,5 @@
-from collections import Counter, Sized, Iterable
+from collections import Counter
+from collections.abc import Sized, Iterable
 from itertools import combinations, chain
 
 import numpy as np
@@ -575,7 +576,7 @@ def add_point(self, point, simplex=None, transform=None):
             return self.bowyer_watson(pt_index, actual_simplex, transform)
 
     def volume(self, simplex):
-        prefactor = np.math.factorial(self.dim)
+        prefactor = factorial(self.dim)
         vertices = np.array(self.get_vertices(simplex))
         vectors = vertices[1:] - vertices[0]
         return float(abs(np.linalg.det(vectors)) / prefactor)