ediff1d is unexpectedly slow. Initially I thought it was due to an unneeded copy operation (I was calling it with to_begin). Then I noticed that a flat iterator is much slower than a ravel array, which was the main reason for slower than expected performance. I plan on submitting a pull request to address this.