Hello everyone,

Perhaps others have asked this question, and if it's been answered, my apologies.

Has someone developed some code to extract eigenvectors from the eigenvalues computed by hqr? If so, would they kindly point me towards the source?

A related question. I am interested also in Jordan Normal (Canonical) forms. I can compute them "by hand" but I would prefer (who wouldn't?) to compute them with a machine. I have not been able to find any source code that can do this.

Thank you,
Mark Allyn

Hi Mark,

The 3rd edition of NR provides a recipe. Otherwise you can use a package like LAPACK.

Saul Teukolsky

Good morning, Saul,

Thanks for the tip on the 3rd edition. I didn't know it was there. If I'm correct 3rd edition is C++ code, and I don't feel really comfortable with a lot of the baggage that comes with classes. This is why I preferred the 2nd edition. I reckon I'll have to shed my pettiness!

I knew about LAPACK and have used it. Same with EISPACK. But, I like Numerical Recipes because its easy to see inside the box; plus I love the writing style. Excellent.

Regards,

Mark