Re: Linear algebra and machine precision
- To: mathgroup at smc.vnet.net
- Subject: [mg77182] Re: Linear algebra and machine precision
- From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
- Date: Tue, 5 Jun 2007 06:46:13 -0400 (EDT)
- Organization: Uni Leipzig
- References: <f40hja$6lh$1@smc.vnet.net>
- Reply-to: kuska at informatik.uni-leipzig.de
Hi, for linear algebra Mathematica comes with the MathKernel librarys from Intel and the librarys may behave different on 32 bint and 64 bit machines. But Mathematica itself know nothing about that and report always the same $MachinePrecision Regards Jens Dmitry Garanin wrote: > Hi, > > I have Mathematica 5.0 on my laptop (Win XP pro 32 bit) and Mathematica 5.2 on my workstation (Win XP pro 64 bit). I have a program operating with eigensystems of huge matrices with nearly degenerate eigenvalues. The results become inaccurate on my laptop when I increase the size of the matrix while I always obtain accurate results on the workstation. > > My initial thought was that on the workstation Mathematica uses longer machine numbers and that is why the results are accurate. However, $MachinePrecision outputs a number below 16 on both machines. Also computing 1-Sqrt(1+AVerySmallNumber) gives the same result on both machines. So that evidently Mathematica cannot really use those 64 bits. Or I understand something wrong? > > Another explanation would be that linear algebra in Mathematica 5.0 is buggy. > > BTW my program solving partial differential equations works well with Mathematica 5.0 but gives unstable results with Mathematica 5.2, so here the situation seems to be reversed. > > What do you think? > > Best regards, > > Dmitry > >