Numpy V ^ t T V - 寻找更快的计算方式

Question

对于正交V^t T V和三对角对称V的numpy加速T，任何人都有快速的想法吗？