Simplicial interpolation with sorting network

[ivanpribec]

I just remembered I read about sorting networks here:

• Changing std::sort at Google’s Scale and Beyond – Experimental chill

There was also this article from the DeepMind team (hotly debated on HackerNews),

• Faster sorting algorithms discovered using deep reinforcement learning | Nature

where they managed to get rid of one instruction in a sorting network, and provided a patch for this to libc++. It’s relatively easy to use the C++ sort in Fortran, as shown here (with the downside it is linked and not instantiated in-place, e.g. for constant arrays).