i have a costly calculation to do for fitting some experimental data. The fitting function is a sum over eigenmodes, each of them containing a specific surface integral. As it i
use multiprocessing pool
import multiprocessing as mp p = mp.Pool(10) res = p.map(your_function, range(Mmin, Mmax))