Add multithreads to function bilab.geometry.fitting