thank you for responding.
Well, I tried transposing as you said too, but then I lost
confidence in the algorithm implemented there, because
for m=[[10,5,-10],[2,-11,10]].to_gsl_matrix, I then get
a matrix u of size (3,2), (it should be 3,3), and u*u.transposed
[ 8.889e-001 -2.222e-001 -2.222e-001
-2.222e-001 5.556e-001 -4.444e-001
-2.222e-001 -4.444e-001 5.556e-001 ],
so nowhere near a unit matrix, either …
With respect to the sizes, I am expecting to be able to keep
(m1,m2) somewhere in the range of 50 or below each.