Re: [Numpy-discussion] test code for user defined types in numpy

2011-12-21 Thread Charles R Harris
On Wed, Dec 21, 2011 at 1:37 PM, Geoffrey Irving wrote: > On Wed, Dec 21, 2011 at 3:56 AM, Charles R Harris > wrote: > > Hi Geoffrey, > > > > On Tue, Dec 20, 2011 at 7:24 PM, Geoffrey Irving wrote: > >> > >> Hello, > >> > >> As a followup to the prior thread on bugs in user defined types in > >

Re: [Numpy-discussion] trouble building numpy 1.6.1 on Scientific Linux 5

2011-12-21 Thread Russell E. Owen
In article , Ralf Gommers wrote: > On Tue, Dec 20, 2011 at 10:52 PM, Russell E. Owen wrote: > > > In article , > > "Russell E. Owen" wrote: > > > > > In article > > > , > > > Ralf Gommers wrote: > > > > > > > On Fri, Dec 9, 2011 at 8:02 PM, Russell E. Owen wrote: > > > > > > > > > I'm tr

Re: [Numpy-discussion] trouble building numpy 1.6.1 on Scientific Linux 5

2011-12-21 Thread Russell E. Owen
In article , Ralf Gommers wrote: > On Tue, Dec 20, 2011 at 10:52 PM, Russell E. Owen wrote: > > > In article , > > "Russell E. Owen" wrote: > > > > > In article > > > , > > > Ralf Gommers wrote: > > > > > > > On Fri, Dec 9, 2011 at 8:02 PM, Russell E. Owen wrote: > > > > > > > > > I'm tr

Re: [Numpy-discussion] test code for user defined types in numpy

2011-12-21 Thread Geoffrey Irving
On Wed, Dec 21, 2011 at 3:56 AM, Charles R Harris wrote: > Hi Geoffrey, > > On Tue, Dec 20, 2011 at 7:24 PM, Geoffrey Irving wrote: >> >> Hello, >> >> As a followup to the prior thread on bugs in user defined types in >> numpy, I converted my rational number class from C++ to C and switched >> to

Re: [Numpy-discussion] test code for user defined types in numpy

2011-12-21 Thread Mark Wiebe
On Tue, Dec 20, 2011 at 10:48 PM, Christopher Jordan-Squire wrote: > On Tue, Dec 20, 2011 at 9:10 PM, Mark Wiebe wrote: > > On Tue, Dec 20, 2011 at 6:24 PM, Geoffrey Irving wrote: > >> > >> Hello, > >> > >> As a followup to the prior thread on bugs in user defined types in > >> numpy, I convert

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

2011-12-21 Thread Lennart Fricke
I read it like that: (**T is the transpose) Let's call M the mass matrix and N the modal mass matrix. Then X**T*M*X=N. If X (matrix of eigenvectors) is normalized with respect to M, N is I (unity) so it just mean that X**T*M*X=I. That is what octave and matlab give you. For this to be true. x**

Re: [Numpy-discussion] test code for user defined types in numpy

2011-12-21 Thread Charles R Harris
Hi Geoffrey, On Tue, Dec 20, 2011 at 7:24 PM, Geoffrey Irving wrote: > Hello, > > As a followup to the prior thread on bugs in user defined types in > numpy, I converted my rational number class from C++ to C and switched > to 32 bits to remove the need for unportable 128 bit numbers. It > shou

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

2011-12-21 Thread Fahreddın Basegmez
According to this page eigenvectors are normalized with respect to the second matrix. Do you guys have any idea how that's done? http://www.kxcad.net/Altair/HyperWorks/oshelp/frequency_response_analysis.htm "If the eigenvectors are normalized with respect to the mass matrix, the modal mass matr

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

2011-12-21 Thread Andrew Jaffe
Just to be completely clear, there is no such thing as a "non-normalized" eigenvector. An eigenvector is only determined *up to a scalar normalization*, which is obvious from the eigenvalue equation: A v = l v where A is the matrix, l is the eigenvalue, and v is the eigenvector. Obviously v is

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

2011-12-21 Thread Olivier Delalleau
Aaah, thanks a lot Lennart, I knew there had to be some logic to Octave's output, but I couldn't see it... -=- Olivier 2011/12/21 Lennart Fricke > Dear Fahreddın, > I think, the norm of the eigenvectors corresponds to some generic > amplitude. But that is something you cannot extract from the s

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

2011-12-21 Thread Lennart Fricke
Dear Fahreddın, I think, the norm of the eigenvectors corresponds to some generic amplitude. But that is something you cannot extract from the solution of the eigenvalue problem but it depends on the initial deflection or velocities. So I think you should be able to use the normalized values just