矩阵linesP0为3xN。我想从3x1的 vector planeP0中减去它。有没有更聪明,更快捷的方法来做到这一点?

现在,我正在使用for循环。下面的示例代码:

MatrixXf temp(linesP0.rows(), linesP0.cols());
for (int i = 0; i < linesP0.cols(); i++)
{
    temp.col(i) = planeP0 - linesP0.block<3, 1>(0, i);
}

我尝试使用colwise(),但没有用。

最佳答案

您可以使用.colwise()来做到这一点,只需要一点点创意即可。

Vector3d v      = Vector3d(1.0, 2.0, 3.0);
Matrix3d m      = Matrix3d::Random();
Matrix3d result = (-m).colwise() + v;
std::cout << result << std::endl;

样本结果:
v = [1 2 3]' (3x1)
m = [1 1 1; 2 2 2]' (3x2)
result = [0 1 2; -1 0 1]' (3x2)

08-25 06:28