torch.symeig solves a regular symmetric eigenvalue problem but not the generalized. In my case I have essentially Ax = lambda*Bx where A and B are symmetric PD matrices.  I guess if B is non-singular I could solve this by computing inv(B)*A and solving that regular eigenvalue problem but there are no longer guarantees on the symmetry of that matrix and the inverse is generally unstable.