From: <Pet...@es...> - 2002-03-27 16:25:22
|
> whereas for 30 equations it goes banana (wrong solution): > > vnl_lsqr.cxx : x = 0 is the exact solution. No iterations were performed. > vnl_lsqr.cxx : The iteration limit ITNLIM was reached. > vnl_lsqr.cxx : residual norm estimate = 3.24731 > vnl_lsqr.cxx : result norm estimate = 14.2702 > vnl_lsqr.cxx : condition no. estimate = 106.047 > vnl_lsqr.cxx : iterations = 30 Whenever you see the message "The iteration limit ITNLIM was reached" you can expect that the result is incorrect or at least imprecise. Could you try changing the # iterations to 120 in this case? P.S. I've changed the default in vnl_lsqr.h to 4*n. Peter. |