Added `parallel` option to gp_fit(). NOTE: execution time does not change!!
Modified call to gp_fit adding parallel=False.