title: Do scipy and numpy svd or eig always return the same singular/eigen vector?
body: <p>Since the SVD decomposition is not unique (pairs of left and right singular vectors can have their sign flipped simultaneously), I was wondering to what extent the U and V matrix returned by <code>scipy.linalg.svd()</code> are 'deterministic' / always the same? </p>
  
  <p>I tried it a few times with a random array on my machine and it seems to always return the same thing (fortunately), but could that vary across machines?</p>