Eigenvalues and vectors

[RichardBJ]

Hi there, was trying to construct a pymc model that calculates eigenvectors in the middle (and other matrix transforms). PYMC is now using this pytensor library for tensor operations I believe, so I dont think there is a way to do it with pytensor but thought I should check?
Kindest Regards
Richard

[ricardoV94]

We have something called eigvalsh not sure if that's all you need:

pytensor/pytensor/tensor/slinalg.py

Lines 553 to 554 in 7bb18f3

	def eigvalsh(a, b, lower=True):
	return Eigvalsh(lower)(a, b)

You may have more luck getting a reply if you post on our Discourse forum: https://discourse.pymc.io/

[RichardBJ]

Gosh I'm sorry for the delay in reply. And thanks for the suggestion, in fact, I think my problem has an easier solution which I now have running which leaving my matrix code as is and then treating it as a "a “black box” likelihood function" as exemplared...
here
Thanks!
R.