Solving eigenvalue problem with root

Hi all,

I have to solve an eigenvalue problem using ROOT. So basically i need to diagonalize a matrix. I don't know how to do that. I checked the user guide for ROOT 5.24, but I couldn't figure it out.

Should I declare the matrix M that is to be diagonalize as

TMatrixDSym(N) M;

if M is an N X N matrix? What happens is when I execute the code any line after this declaration of M is not executed and the pointer comes back to root promt ( root$ )

Then I did the following:

TVectorD eigen;

eigen = M.GetEigenVaues();

Is that correct?

Please guide me. I'm using ROOT 5.22.

Thanks,
Ajinkya
Graduate Student
University of Virginia. Received on Wed Sep 23 2009 - 00:29:55 CEST