我最近的一篇文章中的评论者告诉我,我需要更好地利用c++ 11 move-semantics来处理代码中的瓶颈。以下是需要修复的内容的简化版本。
#include <iostream>
#include <Eigen/Dense>
#include <vector>
void makeCopy(std::vector<Eigen::VectorXd> &oldV){
int n = oldV.size();
std::vector<Eigen::VectorXd> mandatoryCopy;
mandatoryCopy.resize(n);
for(int i = 0; i < n; i++){
mandatoryCopy[i] = oldV[i];
}
// swap the two
oldV = mandatoryCopy;
}
int main(int argc, char **argv)
{
// starting vector
int len(1000);
Eigen::VectorXd placeHolder(50);
std::vector<Eigen::VectorXd> v(len, placeHolder);
// copy it a bunch of times
for(int iter = 0; iter < 1000; ++iter){
std::cout << "iter: " << iter << "\n";
makeCopy(v);
}
return 0;
}
问题:在
makeCopy
的for循环中,oldV[i]
是一个左值,那么我该怎么做类似mandatoryCopy[i]&& = oldV[i]
的事情?这是主要的瓶颈,对不对?我在想类似mandatoryCopy[i]&& = std::move(oldV[i])
的东西,但这显然不会编译。编辑
根据@vsoftco的建议,我已经尝试过
std::vector<Eigen::VectorXd> makeCopy2(std::vector<Eigen::VectorXd> oldV){
int n = oldV.size();
std::vector<Eigen::VectorXd> mandatoryCopy;
mandatoryCopy.resize(n);
for(int i = 0; i < n; i++){
mandatoryCopy[i] = oldV[i];
}
return mandatoryCopy;
}
但我发现它比较慢。 @vsoftco和@ggael都提到,返回修改后的复制参数而不是再次复制会更快,我同意,但是我怀疑这对于我的实际代码是否可行。稍后我可以问这个问题,但这将是一个单独的问题/线程。
最佳答案
您没有看正确的线。如果必须复制一份,那么您将无法摆脱它。尽管如此,最好避免使用for循环和right:
std::vector<Eigen::VectorXd> mandatoryCopy = oldV;
另一方面,您可以通过将
oldV=mandatoryCopy
替换为:来省略第二个副本:std::swap(oldV,mandatoryCopy);
这将执行廉价的指针交换。你得到:
void makeCopy(std::vector<Eigen::VectorXd> &oldV){
std::vector<Eigen::VectorXd> V = oldV;
// do something with V
std::swap(oldV,V);
}
对于功能样式,在第二个示例中,您必须直接使用已经是副本的参数:
std::vector<Eigen::VectorXd> makeCopy2(std::vector<Eigen::VectorXd> V){
// do something with V
return V;
}
并用
v = makeCopy2(v);
调用它。不要忘记使用
-std=c++11
进行编译以获得移动语义副本。最后,将
vector<VectorXd>
打包在MatrixXd
中可能更好。这将大大降低了的内存分配数量:void makeCopy3(MatrixXd &oldV){
int n = oldV.cols();
MatrixXd V = oldV;
for(int i = 0; i < n; i++){
V.col(i) *= 0.99;
}
oldV.swap(V); // or oldV = std::move(V); with c++11 enabled
}