Dimitri Pourbaix wrote: > Hi, > >>> 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.
OK, so you are OK with keeping the fix version for this one at 2.2? Phil > > Regards, > Dim. > ---------------------------------------------------------------------------- > > 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 > http://sb9.astro.ulb.ac.be/~pourbaix * mailto:[email protected] > > --------------------------------------------------------------------- > To unsubscribe, e-mail: [email protected] > For additional commands, e-mail: [email protected] > --------------------------------------------------------------------- To unsubscribe, e-mail: [email protected] For additional commands, e-mail: [email protected]
