I propose to work on refactoring np.random.RandomState usages in the code base to use the preferred np.random.Generator available starting from NumPy 1.17.

It was discussed in some issue (I cannot find it now) that we could start dropping NumPy 1.16 for the next release of SciPy.

This work would be a follow up of #13863 and #13778