c++ - 简单的程序来移动一组共面点不起作用

我正在使用简单的叉积计算一组原始点的一组偏移点。后来我比较看看是否得到的向量集是平行的。我看到比较测试在几乎所有情况下都失败了。有人可以告诉我为什么吗？

// Here osg::Vec3 is a container for three floats
std::vector<osg::Vec3 > originalPoints;
originalPoints.push_back(osg::Vec3(10.0,11.0,13.0));
originalPoints.push_back(osg::Vec3(-10.0,12.0,13.0));
originalPoints.push_back(osg::Vec3(5.0,6.0,13.0));
originalPoints.push_back(osg::Vec3(6.75,11.0,13.0));
originalPoints.push_back(osg::Vec3(8.6,-11.0,13.0));

// Here the "^" operator represents cross product
osg::Vec3 tileUpVec = originalPoints[0]^originalPoints[1];
tileUpVec.normalize();

std::vector<osg::Vec3> shift;

osg::Vec3 lineVec1, lineVec2. resultVec;
lineVec1  = originalPoints[0] - originalPoints[1];
resultVec = lineVec1 ^ tileUpVec;
resultVec.normalise();
shift.push_back(resultVec);

for(unsigned int i = 0; i < originalPoints.size(); i++)
{
  lineVec1 = originalPoints[i-1] - originalPoints[i];
  lineVec2 = originalPoints[i] - originalPoints[i+1];
  resultVec = (lineVec + lineVec2) ^ tileUpVec;
  resultVec.normalise();
  shift.push_back(resultVec);
}

lineVec1  = originalPoints[i-1] - originalPoints[i];
resultVec = lineVec1 ^ tileUpVec;
resultVec.normalise();
shift.push_back(resultVec);

ASSERT_TRUE(shift.size() == originalPoints.size());

std::vector<osg::Vec3> shiftedPoints;

for(unsigned int i = 0; i < shift.size(); i++)
{
  shiftedPoints.push_back((shift[i]+originalPoints[i]));
}

float result;

// Here "*" operator represnts dot product
for(unsigned int j = 0; j < (originalPoints.size() - 1); j++)
{
  lineVec1 = shiftedPoints[j] - shiftedPoints[j+1];
  lineVec1.normalize();
  lineVec2 = originalPoints[j] - originalPoints[j+1];
  lineVec2.normalize();
  result = lineVec1 * lineVec2;
  std::cout << result << std::endl;
  EXPECT_TRUE(result > 0.995);
}

EXPECT_TRUE始终失败。这意味着它们不是平行的。什么原因？

最佳答案

您正在创建一个充满所有不同移位的向量。当您将旧集合中的点移动一堆不同的偏移向量时，您不应期望新点和旧点之间的向量相同。 （如果您做过一些数学运算，认为可以保证达到目标，则应在问题中加入该数值，以便我们进行检查。）

如果您确信数学是正确的，并且问题出在代码中，那么您似乎在for循环的第一次迭代中出现了界限错误。查找originalPoints[-1]。

编辑添加：

在第一个for循环中，您定义

linevec1 = original_points[i-1] - original_points[i] linevec2 = original_points[i] - original_points[i+1]

那么您只使用linevec1+linevec2，为什么不只说

linevec = original_points[i-1]-original_points[i+1]？

关于c++ - 简单的程序来移动一组共面点不起作用，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/23982297/