Thanks in advance!

]]>Can you confirm you are getting correct results with the what you are doing?

]]>has shifted 2nd and 3rd columns for V and U compared to numpy SVD. That means with the special reflection case you have to multiply the 2nd column of V with -1, not the 3rd. That took me some time to find out.. ]]>