|
From: | Mathieu Dubois |
Subject: | Re: Memory exhausted when using \ |
Date: | Tue, 17 Jul 2012 13:07:38 +0200 |
User-agent: | Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.1.15) Gecko/20101027 Fedora/3.0.10-1.fc12 Lightning/1.0b2pre Thunderbird/3.0.10 |
Le 17/07/2012 11:48, Martin Helm a écrit :
Am 16.07.2012 17:02, schrieb Mathieu Dubois:AFAI, if the system is not square A\b returns a minimum norm solution... As I want the same output I need to compute the minimum norm solution.Mathieu, did you simply try? X = lsqrSOL(rows(A), columns(A), A, b, 0, 1e-9, 1e-9, 1e8, 1000, 1); with the lsqrSOL I gave the link to, it is just a .m file which you can put in your current working directory, no fancy installation or so needed. Since I am myself interested in that function and if you don't mind you can save your matrix and vector with octave's save command and upload it somewhere, I would then check how this function performs on your data.
Hi Martin, Sorry for the delay (I have a lot of other problems :).I have installed lsqrSOL and used your snippet (except that I set atol = btol = 1e6). It seems to work (but it's rather slow compared to the matlab implementation which use \).
I will try to make sure that the results are similar. Thanks for your help, Mathieu
[Prev in Thread] | Current Thread | [Next in Thread] |