commons-dev mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From Dimitri Pourbaix <>
Subject Re: [math] Status of MATH-327?
Date Wed, 17 Mar 2010 11:23:09 GMT

>> 28: there is still one spurious non zero singular value
>> (3.547702387229884E-7 which I am going to try to kick out) instead of
>> two.
>> Dim.
> I pushed it to 2.2 because what remained to resolve did not look to
> me like a showstopper for 2.1.  I may be missing something, though,
> regarding severity.  I am fine getting this into 2.1 if you have a fix.

I just checked the code and notice the problem is already present at the
exit of TriDiagonalTransform.  In order to compute SVD, one computes the
eigen decomposition of A^tA which relies upon the Householder tri-diagonal
transformation.  The tri-diagonal maxtrix still has 3 non-zero main diagonal
elements.  The smallest one is about 1.e-13, i.e. slightly too large to be
considered null against the largest, about 900.0

Right now, I see no option but avoiding the computation of A^tA.  It is
fair to ask for 2.2.

Dimitri Pourbaix                         *
Institut d'Astronomie et d'Astrophysique *      Don't worry, be happy
CP 226, office 2.N4.211, building NO     *         and CARPE DIEM.
Universite Libre de Bruxelles            *
Boulevard du Triomphe                    *      Tel : +32-2-650.35.71
  B-1050 Bruxelles                        *      Fax : +32-2-650.42.26     *

To unsubscribe, e-mail:
For additional commands, e-mail:

View raw message