Dear Ed,

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

is

[ 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.

Best regards,

Axel