MathGroup Archive 2001

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

Search the Archive

out of memory, again.

  • To: mathgroup at smc.vnet.net
  • Subject: [mg27473] out of memory, again.
  • From: "Ivan Stegic" <ivan at stegic.com>
  • Date: Tue, 27 Feb 2001 00:37:25 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

hi again...

you may remember the problem i outlined below. i received some help, and
managed to find the eigenvectors. luckily, i am not interested in all
three, rather only the real eigenvector. it turns out that the real eigen
vector is a text file of about 60kb containing cos and sin of two
parameters... ive tried to simplify each component in turn using
FullSimplify, to which the machine crunches for about 26 hours, and spits
out another out of memory exiting message. are there any wise ideas, or
clever tricks i may try? perhaps someone out there has a machine which
could attack the simplification by brute force? i would ofcourse cite your
help, if you like. the eigenvector was too hefty to attach, so i uploaded
it to the web.

1st comp: http://stegic.com/1stcomponent.dat
2nd comp: http://stegic.com/2ndcomponent.dat
3rd comp: http://stegic.com/3rdcomponent.dat (unity!)

and the single eigen vector as a whole, for completeness:
http://stegic.com/N.dat

any advice, help, comments appreciated.

ivan, the desperate.


--- Original Message ---

>Hello everyone...
>
>I've searched the archives extensively, and found various answers to my
>problem. However, none of them helped me. I am trying to find the
>eigenvectors to a 3x3 matrix. I keep running into Out of Memory Exiting
>problems, whether I run it with the frontend or without. I am running
>Mathematica 3 on a RedHat Linux 5.1, Kernel 2.0.35  Pentium II 400 with
>128MB of RAM. The matrix I have contains 2 variables, and I am trying to
>solve the for the Eigenvectors symbolically in terms of these variables.
>The catch is that the variables are trigonometric functions. Following is
>the command I execute. The result is "Out of Memory. Exiting." What can I
>do? Could someone else try this on their more powerful machines, or will
it
>still run into the same problem? Or is there a way of simplifying this? Or
>the matrix?
>
>Thanks very much,
>Ivan.
>
>
>------------------------
>[ivan@apollo ivan]$ more jobs/i.txt
>Eigenvectors[{{Cos[q]^2 - Cos[2*p]*Sin[q]^2,
>    1.*(Sin[2*p]*Sin[2.*p]*Sin[q] + Cos[p]^2*Cos[2.*p]*Sin[2*q]),
>    1.*Cos[2.*p]*Sin[2*p]*Sin[q] - 1.*Cos[p]^2*Sin[2.*p]*Sin[2*q]},
>   {-(Cos[p]^2*Sin[2*q]), 1.*Cos[q]*Sin[2*p]*Sin[2.*p] -
>     1.*Cos[2.*p]*(-(Cos[2*p]*Cos[q]^2) + Sin[q]^2),
>    1.*(Cos[2.*p]*Cos[q]*Sin[2*p] +
>       Sin[2.*p]*(-(Cos[2*p]*Cos[q]^2) + Sin[q]^2))},
>   {Sin[2*p]*Sin[q], -1.*Cos[2.*p]*Cos[q]*Sin[2*p] +
1.*Cos[2*p]*Sin[2.*p],
>    1.*(Cos[2*p]*Cos[2.*p] + Cos[q]*Sin[2*p]*Sin[2.*p])}}]>>>EVec-Symb;
>[ivan@apollo ivan]$
>------------------------
>
>
>(For those who are interested, this matrix is a rotation matrix that takes
>the orthonormal basis of a cone rotating around a replica cone from its
>initial to its final position. p and q are the cone semi-angle, and the
>parametrized position of the rotated cone respectively.)



  • Prev by Date: Re: Re: converting hex strings
  • Next by Date: RE: Re: converting hex strings
  • Previous by thread: Help: Reading unformatted binary fortran files