Add a number_of_cores parameter for OpenMP // in Compute_Omega