Can you please share your script, to show how do you implement the Numba function which calls the sorted_eigenvalues function provided by Jérôme. And have you impelement the parallel part?