Agreed. The problem is deeper. I am out of ideas.
On Sun, Apr 7, 2013 at 5:00 AM, Sean Owen <sro...@gmail.com> wrote: > Oh, I thought you were stil referring to the toy example in this > thread. Yes, k=10 is definitely larger than its rank. This is why I > said Y'*Y is certainly not invertible in this case. I was responding > to you showing that ALS still worked -- indeed it does but it never > inverts / solves a system involving Y'*Y. It's later, in fold-in. > > All that said I don't think inverting is the issue here. Using the SVD > to invert didn't change things, and neither did actually solving the > Ax=b problem instead of inverting A by using Householder reflections. > > > > On Sun, Apr 7, 2013 at 1:45 AM, Koobas <koo...@gmail.com> wrote: > > Okay, you do have a problem. > > Y'*Y is 10x10, but it's rank is 5. > > Has to have something to do with the input data. > > > > > > > > On Sat, Apr 6, 2013 at 7:47 PM, Sean Owen <sro...@gmail.com> wrote: > > > >> For example, here's Y: > >> > >> Y = > >> > >> -0.278098 -0.256438 0.127559 -0.045869 -0.769172 -0.255599 > >> 0.150450 -0.436548 0.209881 -0.526238 > >> 0.613175 -0.600739 -0.291662 -1.142282 0.277204 -0.296846 > >> -0.175122 0.031656 -0.202138 -0.254480 > >> -0.187816 -0.889571 0.052191 -0.304053 0.498097 -0.049822 > >> -0.972282 -0.240532 0.155711 -0.627668 > >> -0.065179 -0.055424 0.977480 0.104342 0.594501 0.033205 > >> -0.896222 -0.345715 -0.371288 -0.489602 > >> -0.434807 -0.403650 0.264583 -0.110285 -1.318951 -0.452470 > >> 0.274445 -0.755704 0.313150 -0.903234 > >> > >> and R from the QR decomposition of Y' * Y: > >> > >> R = > >> > >> 2.56259 -1.35164 -2.43837 1.27844 -0.17692 -0.30514 1.09366 > >> -0.84664 0.58601 1.06875 > >> 0.00000 1.03316 2.61600 -0.46070 -1.46785 -0.10841 0.24828 > >> -2.32186 -2.00163 -0.71470 > >> 0.00000 0.00000 2.11507 1.15523 1.10757 0.36407 -0.31567 > >> 2.77361 0.77367 -0.84055 > >> 0.00000 0.00000 0.00000 0.54242 -0.01545 0.21761 0.26630 > >> 0.13972 0.44089 0.02783 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00000 -0.00000 > >> 0.00000 0.00000 -0.00000 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 > >> -0.00000 0.00000 0.00000 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 > >> 0.00000 0.00000 -0.00000 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 > >> 0.00000 0.00000 -0.00000 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 > >> 0.00000 0.00000 0.00000 > >> 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 > >> 0.00000 0.00000 0.00000 > >> > >> > >> Separately I tried avoiding the inverse altogether here and just using > >> the QR decomposition to solve a system where necessary. Probably a > >> better move anyway. But same result. I think I'm not really > >> quantifying the problem properly, but it's not really a matter of > >> condition number or machine precision. Condition numbers are >1 in > >> these cases but not that large. > >> > >> > >> On Sun, Apr 7, 2013 at 12:19 AM, Koobas <koo...@gmail.com> wrote: > >> > I don't see why the inverse of Y'*Y does not exist. > >> > What Y do you end up with? > >> >