[英]What is the best way to compute the trace of a matrix product in numpy?
[英]Compute trace of matrix product using numpy/pytorch broadcasting
設A為(nxm)矩陣,M為(mxm)矩陣。 為矩陣的軌跡編寫tr(),我需要計算tr(AM(A ^ T))。 但是,最終的跟蹤操作將丟棄大部分計算。 我可以使用numpy或pytorch的廣播規則來僅計算AM(A ^ T)的必要對角線嗎?
更新:這是我在PyTorch中計算對角線的解決方案:
torch.sum(torch.sum(At()[:,None,:]*M[:,:,None],0)*At(),0)
您將必須計算至少兩個矩陣乘積之一。 隨后,您可以在這里使用答案之一: 在numpy中計算矩陣乘積的跡線的最佳方法是什么?
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.