将Matlab的bsxfun转换为特征

Question

假设我们有维数MxN的矩阵A和维数Mx1的向量a。在Matlab中，将a与'A‘的所有列相乘，我们可以这样做

bsxfun(@times, a, A)

在本征中是否有一种等价的方法，而不必在矩阵的列上循环？

我想做的是

M = bsxfun(@times, a, A) + bsxfun(@times, a2, A2)

并希望艾根的懒惰评估将使它更有效率。

谢谢!

Answer 1

你可以：

M = A.array().colwise()*a.array();

需要使用.array()来重新定义运算符*的语义，使之成为按系数排序的产品(如果A和a是Array<>对象，则不需要)。

在这种特殊情况下，最好将其编写为缩放操作：

M = a.asDiagonal() * A;

在这两种情况下，您将不会得到任何临时的，因为懒惰的评估。