给出一个Vector
表示v
和Matrix
表示M
,在Julia中计算矩阵二次型的最快方法是什么,即v'Mv
?最优雅的是什么?
注意:我希望返回值为标量.有趣的是,如果v = rand(3)
和M = rand(3, 3)
,则v'*M*v
返回包含一个元素而不是标量的向量.我并没有期待这种行为,虽然已经阅读了足够多的github问题页面,怀疑这种行为有充分的理由让我觉得不够聪明.所以,显然(v'*M*v)[1]
会做这个工作,只是想知道是否有更好的方法......
返回标量的一个选项是dot(v, M*v)
.