Re: Orthogonazlie with Method->"Householder"
- To: mathgroup at smc.vnet.net
- Subject: [mg97781] Re: [mg97744] Orthogonazlie with Method->"Householder"
- From: Maris Ozols <marozols at gmail.com>
- Date: Sat, 21 Mar 2009 05:19:18 -0500 (EST)
- References: <200903200740.CAA17270@smc.vnet.net> <49C40120.2020701@wolfram.com>
Sorry for sending such a large numeric matrix as plain text in my e-mail. Singular values mentioned in both replies indicate that you didn't succeed in recovering it correctly. So I sent my notebook in attachment to Daniel Lichtblau and he confirmed the bad behavior of Orthogonalize[M,Method->"Householder"] and filled a bug report. Sorry for the confusion, ~Maris Ozols~ > Consider the singular values of your matrix: > {1.75663*10^15, 2.21874*10^12, > 5.16322*10^10, 488354., 6518.66, 83.9465} > and you will see that 16 decimal digits may lead to problems. > If you calculate with exact numbers you will not have this problem. > But be aware, simply rationalizing your matrix introduces some error and can > make linear depend vectors independent. > Daniel > [I recopy the matrix below, since my email reader showed something that was > messed up a bit. --dl] > > In[2]:= InputForm[SingularValueList[mat]] > Out[2]//InputForm= > {7.053736007427739*^14, 44.4483160745624, 40.569302085579835, > 25.471811796860234, 17.698020291765765}
- References:
- Orthogonazlie with Method->"Householder"
- From: Maris Ozols <marozols@gmail.com>
- Orthogonazlie with Method->"Householder"