MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Linear algebra and machine precision

  • To: mathgroup at
  • Subject: [mg77182] Re: Linear algebra and machine precision
  • From: Jens-Peer Kuska <kuska at>
  • Date: Tue, 5 Jun 2007 06:46:13 -0400 (EDT)
  • Organization: Uni Leipzig
  • References: <f40hja$6lh$>
  • Reply-to: kuska at


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


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

  • Prev by Date: Re: Two Version 6.0 Package Problems
  • Next by Date: Re: Unable to copy graphics from Mathematica 6 to Publicon
  • Previous by thread: Re: Linear algebra and machine precision
  • Next by thread: Very Strange Behaviour about something Very Simple