As ++hippo said, you cannot use the methods from one completely separate object type on another.

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.