I’m using TMatrixDEigen to find the eigenvalues of a TMatrixD
my code looks like this:
TMatrixD S = TMatrixD(3,3); TMatrixDEigen * Sbis = NULL; TMatrixD VP= TMatrixD(3,3); /* put something in S */ Sbis = new TMatrixDEigen(S); VP = Sbis->GetEigenValues();
basically it works, but i sometimes get an error message :
TMatrixDEigen::MakeSchur(…) seems to be the function that triangularize the original matrix via a schur transformation to get the eigenvalues.
my concern is that this function doesn’t give any feedback on what happen in case of such an error and what may have caused it.
In my particular case, “S” is a generalized sphericity matrix, which eigenvalues are used to compute a few event shape variables (since i haven’t found so far something already implemented)
well, can someone help me to get along with this problem ? (avoid makeschurr error or find a nice eventshape program) thanks a lot in advance
all the best