If you are using PDL, so your data is already in that form, then you should use PDL methods... and PDL has a plethora of matrix operations, as matrices are one of its primary reasons for existing. For example, I think that eigens_sym or eigens is mostly likely going to give you the information you want.