快速计算向量乘法和矩阵生成的对角线？

Question

我想计算y' * diag(A' * B * A) * y

• y是N * 1向量
• A是一个m * N矩阵，N >> m
• B是m * m矩阵，B是对称正定矩阵，如果需要，Cholesky分解B = L * L'是预先计算的。

是否有可能以O(N*m)复杂度计算上述表达式？