Re: Calculating Condition Numbers

Hi,
	I think you can calculate condition number in this
manner.............. But surely it is not a fast
method............ It need to do a lot of
computations..


                dense1D<int> pivots((*A).nrows(), 0);
                lu_factorize(*A, pivots);
                lu_inverse(*A,pivots,AInv);
                A_norm = mtl::infinity_norm(tempA);
		//tempA is copy of A
                InvA_norm = mtl::infinity_norm(AInv);
                condA = A_norm*InvA_norm;

(This is not at all a good solution for the data you
have.)
Cheers,
Asharam



________________________________________________________________________
Want to chat instantly with your online friends?  Get the FREE Yahoo!
Messenger http://mail.messenger.yahoo.co.uk
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/