Problem
The current default behaviour of linalg.svd computes the full decomposition by setting full_matrices=True by default. On the other hand, the default behaviour of linalg.qr specifies mode="reduced" by default.

Proposal
Change de default to full_matrices=False in linalg.svd. In most cases the reduced SVD is good enough. Even more, it’s better to do less work by default and just make the more expensive operation opt-in.

cc @rgommers @mruberry