multithread fix

Author: D-McDougall

matexp/__init__.py

@@ -18,6 +18,7 @@
 from .optimizer import Optimize1D, Optimize2D
 from pathlib import Path
 import multiprocessing
+import threadpoolctl
 import numpy as np
 import dill
 import time
@@ -50,9 +51,7 @@ def _initialize_worker_process(derivative_pickle):
     global _derivative
     _derivative = dill.loads(derivative_pickle)
     # Disable automatic multithreading
-    os.environ['OPENBLAS_NUM_THREADS'] = '1'
-    os.environ['MKL_NUM_THREADS'] = '1'
-    os.environ['OMP_NUM_THREADS'] = '1'
+    threadpoolctl.ThreadpoolController().limit(limits=1)
 
 def main(nmodl_filename, inputs, time_step, temperature,
          error, target,

pyproject.toml

@@ -27,5 +27,5 @@
         "Operating System :: OS Independent",
     ]
     requires-python = ">=3.8"
-    dependencies = ["neuron", "numpy", "scipy", "dill"]
+    dependencies = ["neuron", "numpy", "scipy", "dill", "threadpoolctl"]
     scripts = {matexp = "matexp:__main__._placeholder", matexp-manual = "matexp:manual._placeholder"}