From: Nuralanur@...
Date: 2006-07-31T07:01:09+09:00
Subject: Re: Setting precision in Ruby-GSL calculations ?

-------------------------------1154296863
Content-Type: text/plain; charset="US-ASCII"
Content-Transfer-Encoding: 7bit

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 

-------------------------------1154296863--