MathGroup Archive 1999

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

Search the Archive

RE: Eigensystem kills MathKernel

  • To: mathgroup at smc.vnet.net
  • Subject: [mg17357] RE: [mg17334] Eigensystem kills MathKernel
  • From: "Jean-Marie THOMAS" <jmthomas at agat.net>
  • Date: Mon, 3 May 1999 01:45:57 -0400
  • Sender: owner-wri-mathgroup at wolfram.com

On my machine (Mathematica 3.0.1 for MS Windows, Winnt4SP4) your computation
crashes the kernel too (bug ? very probably)

I've performed the following tests :
g is just a control on time and memory

\!\(maxMemory = 256\ 10\^6; maxTime = 600; \n
  g := TimeConstrained[\
        MemoryConstrained[#, \ maxMemory,
          Print["\<memory used when overflow : \>", MemoryInUse]],
maxTime]&;
  \ng at Eigensystem[matrix]\n (*this command crashes the kernel*)
  reducedMatrix = Map[Chop, matrix, {2}]\n
  g at Eigensystem[reducedMatrix]\n (*this command crashes the kernel*)
  shortenedMatrix = SetPrecision[matrix, 16]\n
  g at Eigensystem[shortenedMatrix]\) (*this command does not crash the
kernel*)

Depending on the precision you want to reach (386 is very high !) you might
build something like :
s[n_]:=SetPrecision[matrix, n]
and then
Do[
Put[{n,g at Eigensystem[s[n]]},"result"],{n,16,386}]
until the kernel crashes

Hope this helps,


****************************************
Jean-Marie THOMAS
mailto:jmthomas at agat.net
Conseil et Audit en Ingenierie de Calcul
Strasbourg, France
http://www.agat.net
****************************************

-----Original Message-----
From: Joeri Kenes [mailto:kenes at imec.be]
To: mathgroup at smc.vnet.net
Subject: [mg17357] [mg17334] Eigensystem kills MathKernel



When I try to calculate the EigenSystem of a 16x16 complex matrix with
high precision (app.380), the Mathematica kernel under Linux just dies
on me without any error message.

On an Apple, with Mathematica 3.0 under System 7, I get the message

Sorry, a system error occured.
"MathKernel"
error type 11
(which I think is a "miscellanous hardware error")

I have attached a notebook with an example of a matrix which produces
this error.

Does anyone have a clue why this happens?

Joeri Kenes





  • Prev by Date: Re: The behaviour of DisplayTogether
  • Next by Date: Phase Plane Diagrams in Mathematica
  • Previous by thread: RE: Groebner bases
  • Next by thread: Phase Plane Diagrams in Mathematica