5

当我使用 Eigen 创建矩阵时，如下所示：

  Eigen::MatrixXd M(3,3);
  M<< 1.0, 4.0, 7.0, 2.0, 5.0, 8.0, 3.0, 6.0, 9.0 ;
  std::cout<<M<<std::endl;

生产

1 4 7
2 5 8
3 6 9

我可以用指针遍历数据，打印每个元素：

  double * d = M.data();
  for(int i = 0;i<M.size();i++)
  {
    std::cout<<*d<<" ";
    d++;
  }
  std::cout<<std::endl;

生产

1 2 3 4 5 6 7 8 9

我还可以使用 std::copy 将其复制到堆栈上相同类型的数组中，然后打印该数组的元素：

  double data_copy[9];
  std::copy(M.data(),M.data()+M.size(),data_copy);
  for(int i = 0;i<M.size();i++)
  {
    std::cout<<data_copy[i]<<" ";
  }
  std::cout<<std::endl;

生产

1 2 3 4 5 6 7 8 9

但是，我似乎无法使用 memcpy 进行等效复制。这只设法复制第一个元素：

  double data_memcopy[9];
  memcpy(data_memcopy,M.data(),M.size());
  for(int i = 0;i<M.size();i++)
  {
    std::cout<<data_memcopy[i]<<" ";
  }
  std::cout<<std::endl;

生产

1 0 0 0 0 0 0 0 0

memcpy 没有做我期望的事情的任何充分理由？（还是我使用不当？）

4

1 回答 1

8

您需要乘以双精度数：

memcpy(data_memcopy,M.data(),M.size() * sizeof(double));

否则，您只会复制 M.size() 字节，并且每个双精度数在您的机器上都超过一个字节。

因此，您可能只写入了第一个和第二个双精度（它们在您的系统上可能是 8 个字节，因为您成功复制了第一个，而第二个双精度可能只复制了一个可能为零的字节）。

于 2012-02-13T09:52:26.177 回答