Added call to use_solver(useUmfpack = False) in fit_surface(), to avoid deprecation warnings about umfpack