[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Octave/C++ matrix Inv() comparison
From: |
CdeMills |
Subject: |
Re: Octave/C++ matrix Inv() comparison |
Date: |
Tue, 9 Jul 2013 06:45:48 -0700 (PDT) |
Hello,
before starting to invest your time into comparing algorithms and toolboxes
... first look at your data. Perform a singular value decomposition, look at
the ratio of the biggest singular value (top left of the S matrix) to the
smallest one (bottom right). A ratio greater than 1/sqrt(eps), around 1e8,
is calling for rounding issues. If this is the case, try understanding why
this ratio is so low: it is the symptom that some of your columns are
strongly linearly correlated. Common cures: compute pseudo-inverse (pinv),
i.e. neglecting "directions" where noise sensitivity is great; or
minimal-norm solutions.
Regards
Pascal
--
View this message in context:
http://octave.1599824.n4.nabble.com/Octave-C-matrix-Inv-comparison-tp4655291p4655433.html
Sent from the Octave - Maintainers mailing list archive at Nabble.com.
- Octave/C++ matrix Inv() comparison, ionone, 2013/07/07
- Re: Octave/C++ matrix Inv() comparison,
CdeMills <=
- Re: Octave/C++ matrix Inv() comparison, ionone, 2013/07/10
- Re: Octave/C++ matrix Inv() comparison, Nir Krakauer, 2013/07/10
- Re: Octave/C++ matrix Inv() comparison, Pascal Dupuis, 2013/07/10
- Re: Octave/C++ matrix Inv() comparison, ionone, 2013/07/10
- Re: Octave/C++ matrix Inv() comparison, CdeMills, 2013/07/10
- Re: Octave/C++ matrix Inv() comparison, ionone, 2013/07/11
- Re: Octave/C++ matrix Inv() comparison, ionone, 2013/07/11
- Re: Octave/C++ matrix Inv() comparison, CdeMills, 2013/07/11
- Re: Octave/C++ matrix Inv() comparison, ionone, 2013/07/12
- Re: Octave/C++ matrix Inv() comparison, CdeMills, 2013/07/12